Skip to main content

SparseSolver

Struct SparseSolver 

Source
pub struct SparseSolver { /* private fields */ }
Expand description

The one linear solver Carapace ships — a concrete struct, not a trait or enum, since there is exactly one implementation to commit to.

faer’s sparse LU (with built-in COLAMD/AMD fill-reducing ordering), not a dense nalgebra solve. A dense O(n^3) LU doesn’t scale to real models. It builds and solves correctly under wasm32.

Implementations§

Source§

impl SparseSolver

Source

pub fn new() -> Self

Source

pub fn factor( &self, k: &SparseColMat<usize, f64>, ) -> Result<SparseFactorization, AnalysisError>

Factor k, without solving — the half of a solve that TangentStrategy::ReuseAtStepStart/Initial (docs/algorithms.md §3.1) let a caller skip repeating across iterations. faer‘s Lu is a fully owned value (no lifetime tied back to k), so the result outlives the matrix it was formed from — see SparseFactorization. Symbolic analysis is reused while the full CSC pattern matches; values (including zeros) never participate in cache invalidation. Numeric factorization is always recomputed here, independently of callers’ numeric tangent-reuse caches.

Source

pub fn solve( &self, k: &SparseColMat<usize, f64>, r: &DVector<f64>, ) -> Result<DVector<f64>, AnalysisError>

Factor-then-solve in one call — for Algorithm::Linear and anywhere else a single ad-hoc solve is enough. Equivalent to self.factor(k)? .solve(r), kept as its own method so a one-off caller doesn’t have to name SparseFactorization just to discard it immediately.

Trait Implementations§

Source§

impl Default for SparseSolver

Source§

fn default() -> Self

Returns the “default value” for a type. Read more

Auto Trait Implementations§

Blanket Implementations§

Source§

impl<T> Any for T
where T: 'static + ?Sized,

Source§

fn type_id(&self) -> TypeId

Gets the TypeId of self. Read more
Source§

impl<T> Borrow<T> for T
where T: ?Sized,

Source§

fn borrow(&self) -> &T

Immutably borrows from an owned value. Read more
Source§

impl<T> BorrowMut<T> for T
where T: ?Sized,

Source§

fn borrow_mut(&mut self) -> &mut T

Mutably borrows from an owned value. Read more
§

impl<T> ByRef<T> for T

§

fn by_ref(&self) -> &T

Source§

impl<T> From<T> for T

Source§

fn from(t: T) -> T

Returns the argument unchanged.

§

impl<T, U> Imply<T> for U
where T: ?Sized, U: ?Sized,

Source§

impl<T, U> Into<U> for T
where U: From<T>,

Source§

fn into(self) -> U

Calls U::from(self).

That is, this conversion is whatever the implementation of From<T> for U chooses to do.

Source§

impl<T> Same for T

Source§

type Output = T

Should always be Self
§

impl<SS, SP> SupersetOf<SS> for SP
where SS: SubsetOf<SP>,

§

fn to_subset(&self) -> Option<SS>

The inverse inclusion map: attempts to construct self from the equivalent element of its superset. Read more
§

fn is_in_subset(&self) -> bool

Checks if self is actually part of its subset T (and can be converted to it).
§

fn to_subset_unchecked(&self) -> SS

Use with care! Same as self.to_subset but without any property checks. Always succeeds.
§

fn from_subset(element: &SS) -> SP

The inclusion map: converts self to the equivalent element of its superset.
Source§

impl<T, U> TryFrom<U> for T
where U: Into<T>,

Source§

type Error = Infallible

The type returned in the event of a conversion error.
Source§

fn try_from(value: U) -> Result<T, <T as TryFrom<U>>::Error>

Performs the conversion.
Source§

impl<T, U> TryInto<U> for T
where U: TryFrom<T>,

Source§

type Error = <U as TryFrom<T>>::Error

The type returned in the event of a conversion error.
Source§

fn try_into(self) -> Result<U, <U as TryFrom<T>>::Error>

Performs the conversion.