lazymatrix 0.4.0

Lazy normalized design matrices: apply column centering/scaling without materializing (X - 1cᵀ)S⁻¹, preserving sparsity.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20

use super::{faer, faer_traits};
use crate::Scalar;
use faer::sparse::SparseColMat;
use faer::{Accum, Col, Par};

pub(super) fn multiply<F: Scalar + faer_traits::ComplexField>(
    matrix: &SparseColMat<usize, F>,
    input: &Col<F>,
    output: &mut Col<F>,
    parallelism: Par,
) {
    faer::sparse::linalg::matmul::sparse_dense_matmul(
        output.as_mat_mut(),
        Accum::Replace,
        matrix.as_ref().transpose(),
        input.as_mat(),
        F::one(),
        parallelism,
    );
}