solvr 0.2.0

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
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92

//! CUDA implementation of matrix equation solvers.

#[cfg(feature = "cuda")]
use crate::linalg::impl_generic::{
    continuous_lyapunov_impl, discrete_lyapunov_impl, solve_care_impl, solve_care_iterative_impl,
    solve_dare_impl, solve_dare_iterative_impl, solve_discrete_lyapunov_iterative_impl,
    sylvester_impl,
};
#[cfg(feature = "cuda")]
use crate::linalg::traits::matrix_equations::MatrixEquationAlgorithms;
#[cfg(feature = "cuda")]
use numr::error::Result;
#[cfg(feature = "cuda")]
use numr::runtime::cuda::{CudaClient, CudaRuntime};
#[cfg(feature = "cuda")]
use numr::tensor::Tensor;

#[cfg(feature = "cuda")]
impl MatrixEquationAlgorithms<CudaRuntime> for CudaClient {
    fn solve_sylvester(
        &self,
        a: &Tensor<CudaRuntime>,
        b: &Tensor<CudaRuntime>,
        c: &Tensor<CudaRuntime>,
    ) -> Result<Tensor<CudaRuntime>> {
        sylvester_impl(self, a, b, c)
    }

    fn solve_continuous_lyapunov(
        &self,
        a: &Tensor<CudaRuntime>,
        q: &Tensor<CudaRuntime>,
    ) -> Result<Tensor<CudaRuntime>> {
        continuous_lyapunov_impl(self, a, q)
    }

    fn solve_discrete_lyapunov(
        &self,
        a: &Tensor<CudaRuntime>,
        q: &Tensor<CudaRuntime>,
    ) -> Result<Tensor<CudaRuntime>> {
        discrete_lyapunov_impl(self, a, q)
    }

    fn solve_care(
        &self,
        a: &Tensor<CudaRuntime>,
        b: &Tensor<CudaRuntime>,
        q: &Tensor<CudaRuntime>,
        r: &Tensor<CudaRuntime>,
    ) -> Result<Tensor<CudaRuntime>> {
        solve_care_impl(self, a, b, q, r)
    }

    fn solve_dare(
        &self,
        a: &Tensor<CudaRuntime>,
        b: &Tensor<CudaRuntime>,
        q: &Tensor<CudaRuntime>,
        r: &Tensor<CudaRuntime>,
    ) -> Result<Tensor<CudaRuntime>> {
        solve_dare_impl(self, a, b, q, r)
    }

    fn solve_care_iterative(
        &self,
        a: &Tensor<CudaRuntime>,
        b: &Tensor<CudaRuntime>,
        q: &Tensor<CudaRuntime>,
        r: &Tensor<CudaRuntime>,
    ) -> Result<Tensor<CudaRuntime>> {
        solve_care_iterative_impl(self, a, b, q, r)
    }

    fn solve_dare_iterative(
        &self,
        a: &Tensor<CudaRuntime>,
        b: &Tensor<CudaRuntime>,
        q: &Tensor<CudaRuntime>,
        r: &Tensor<CudaRuntime>,
    ) -> Result<Tensor<CudaRuntime>> {
        solve_dare_iterative_impl(self, a, b, q, r)
    }

    fn solve_discrete_lyapunov_iterative(
        &self,
        a: &Tensor<CudaRuntime>,
        q: &Tensor<CudaRuntime>,
    ) -> Result<Tensor<CudaRuntime>> {
        solve_discrete_lyapunov_iterative_impl(self, a, q)
    }
}