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
21

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,
) {
    // Older faer releases only accept CSC input, so compute y^T = x^T X.
    faer::sparse::linalg::matmul::dense_sparse_matmul(
        output.as_mat_mut().transpose_mut(),
        Accum::Replace,
        input.as_mat().transpose(),
        matrix.as_ref(),
        F::one(),
        parallelism,
    );
}