graphblas_sparse_linear_algebra 0.44.0

mod doc_test;

// #[cfg(test)]
mod tests {
    use graphblas_sparse_linear_algebra::collections::sparse_matrix::operations::{
        GetSparseMatrixElement, SetSparseMatrixElement,
    };
    use rayon::prelude::*;
    use std::sync::Mutex;

    use graphblas_sparse_linear_algebra::collections::sparse_matrix::{
        Coordinate, GetMatrixElementValue, MatrixElement, Size, SparseMatrix,
    };
    use graphblas_sparse_linear_algebra::context::Context;
    use graphblas_sparse_linear_algebra::operators::binary_operator::Plus;
    use graphblas_sparse_linear_algebra::operators::element_wise_multiplication::{
        ApplyElementWiseMatrixMultiplicationBinaryOperator,
        ElementWiseMatrixMultiplicationBinaryOperator,
    };
    use graphblas_sparse_linear_algebra::operators::mask::SelectEntireMatrix;
    use graphblas_sparse_linear_algebra::operators::options::OptionsForOperatorWithMatrixArguments;

    #[test]
    fn parallel_calls_to_graphblas() {
        let context = Context::init_default().unwrap();

        let number_of_matrices = 100;

        let matrix_size = Size::new(10, 5);
        let mut matrices: Vec<SparseMatrix<i32>> = (0..number_of_matrices)
            .into_par_iter()
            .map(|_| SparseMatrix::<i32>::new(&context, &matrix_size).unwrap())
            .collect();

        matrices.par_iter_mut().for_each(|matrix| {
            matrix
                .set_matrix_element(&MatrixElement::from_triple(1, 2, 3))
                .unwrap()
        });

        let add_operator = Plus::<i32>::new();
        let options = OptionsForOperatorWithMatrixArguments::new_default();
        let result_matrix = Mutex::new(SparseMatrix::<i32>::new(&context, &matrix_size).unwrap());

        let element_wise_matrix_add_operator = ElementWiseMatrixMultiplicationBinaryOperator::new();

        matrices.par_iter().for_each(|matrix| {
            element_wise_matrix_add_operator
                .apply(
                    &*matrix,
                    &add_operator,
                    &*matrix,
                    &add_operator,
                    &mut *result_matrix.lock().unwrap(),
                    &SelectEntireMatrix::new(&context),
                    &options,
                )
                .unwrap();
        });

        let result_matrix = result_matrix.into_inner().unwrap();

        assert_eq!(
            600,
            result_matrix
                .element(Coordinate::new(1, 2))
                .unwrap()
                .unwrap()
                .value()
        );
    }
}