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

//! Unified client trait for stiff ODE solvers (BDF, Radau).
//!
//! This module provides a single `StiffSolverClient` trait that conditionally
//! includes `IterativeSolvers` and `SparseOps` when the sparse feature is enabled,
//! avoiding trait duplication between solvers.
use crate::DType;

use numr::ops::{LinalgOps, ScalarOps, TensorOps, UtilityOps};
use numr::runtime::{Runtime, RuntimeClient};

#[cfg(feature = "sparse")]
use numr::algorithm::iterative::IterativeSolvers;
#[cfg(feature = "sparse")]
use numr::ops::IndexingOps;
#[cfg(feature = "sparse")]
use numr::sparse::SparseOps;

/// Client trait for stiff ODE solvers requiring linear algebra operations.
///
/// When the `sparse` feature is enabled, this trait also requires `IterativeSolvers`
/// for GMRES-based sparse linear system solving, `SparseOps` for dense-to-CSR
/// conversion, and `IndexingOps` for `gather_2d` used by the direct sparse solver.
#[cfg(feature = "sparse")]
pub trait StiffSolverClient<R: Runtime<DType = DType>>:
    TensorOps<R>
    + ScalarOps<R>
    + LinalgOps<R>
    + UtilityOps<R>
    + RuntimeClient<R>
    + IterativeSolvers<R>
    + SparseOps<R>
    + IndexingOps<R>
{
}

#[cfg(feature = "sparse")]
impl<R, T> StiffSolverClient<R> for T
where
    R: Runtime<DType = DType>,
    T: TensorOps<R>
        + ScalarOps<R>
        + LinalgOps<R>
        + UtilityOps<R>
        + RuntimeClient<R>
        + IterativeSolvers<R>
        + SparseOps<R>
        + IndexingOps<R>,
{
}

/// Client trait for stiff ODE solvers requiring linear algebra operations.
///
/// Without the `sparse` feature, only dense linear algebra is available.
#[cfg(not(feature = "sparse"))]
pub trait StiffSolverClient<R: Runtime<DType = DType>>:
    TensorOps<R> + ScalarOps<R> + LinalgOps<R> + UtilityOps<R> + RuntimeClient<R>
{
}

#[cfg(not(feature = "sparse"))]
impl<R, T> StiffSolverClient<R> for T
where
    R: Runtime<DType = DType>,
    T: TensorOps<R> + ScalarOps<R> + LinalgOps<R> + UtilityOps<R> + RuntimeClient<R>,
{
}