solvr 0.2.0-beta.2

Advanced computing library for real-world problem solving - optimization, differential equations, interpolation, statistics, and more
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63

//! CUDA implementation of finite difference PDE solvers.

use numr::runtime::cuda::{CudaClient, CudaRuntime};
use numr::tensor::Tensor;

use crate::pde::error::PdeResult;
use crate::pde::impl_generic::{heat_2d_impl, heat_3d_impl, poisson_impl, wave_impl};
use crate::pde::traits::FiniteDifferenceAlgorithms;
use crate::pde::types::{
    BoundarySpec, FdmOptions, FdmResult, Grid2D, Grid3D, TimeDependentOptions, TimeResult,
};

impl FiniteDifferenceAlgorithms<CudaRuntime> for CudaClient {
    fn fdm_poisson(
        &self,
        f: &Tensor<CudaRuntime>,
        grid: &Grid2D,
        boundary: &[BoundarySpec<CudaRuntime>],
        options: &FdmOptions,
    ) -> PdeResult<FdmResult<CudaRuntime>> {
        poisson_impl(self, f, grid, boundary, options)
    }

    fn fdm_heat_2d(
        &self,
        u0: &Tensor<CudaRuntime>,
        alpha: f64,
        source: Option<&Tensor<CudaRuntime>>,
        grid: &Grid2D,
        boundary: &[BoundarySpec<CudaRuntime>],
        time_opts: &TimeDependentOptions,
        options: &FdmOptions,
    ) -> PdeResult<TimeResult<CudaRuntime>> {
        heat_2d_impl(self, u0, alpha, source, grid, boundary, time_opts, options)
    }

    fn fdm_heat_3d(
        &self,
        u0: &Tensor<CudaRuntime>,
        alpha: f64,
        source: Option<&Tensor<CudaRuntime>>,
        grid: &Grid3D,
        boundary: &[BoundarySpec<CudaRuntime>],
        time_opts: &TimeDependentOptions,
        options: &FdmOptions,
    ) -> PdeResult<TimeResult<CudaRuntime>> {
        heat_3d_impl(self, u0, alpha, source, grid, boundary, time_opts, options)
    }

    fn fdm_wave(
        &self,
        u0: &Tensor<CudaRuntime>,
        v0: &Tensor<CudaRuntime>,
        c: f64,
        source: Option<&Tensor<CudaRuntime>>,
        grid: &Grid2D,
        boundary: &[BoundarySpec<CudaRuntime>],
        time_opts: &TimeDependentOptions,
        options: &FdmOptions,
    ) -> PdeResult<TimeResult<CudaRuntime>> {
        wave_impl(self, u0, v0, c, source, grid, boundary, time_opts, options)
    }
}