Function extendr_api::prelude::sparse::linalg::triangular_solve::solve_unit_upper_triangular_in_place
pub fn solve_unit_upper_triangular_in_place<I, E>(
u: SparseColMatRef<'_, I, E>,
conj: Conj,
rhs: MatMut<'_, E>,
parallelism: Parallelism<'_>,
)where
I: Index,
E: ComplexField,
Expand description
Assuming self
is a unit upper triangular matrix, solves the equation Op(self) * X = rhs
, and stores the result in rhs
, where Op
is either the conjugate or the identity
depending on the value of conj
.
ยงNote
The matrix indices need not be sorted, but the diagonal element is assumed to be the last stored element in each column.