Method run

src/modules/linear_systems/sparse_linear_systems.cpp:208–272

* @brief Perform the sparse linear system transition step */

 * @brief Perform the sparse linear system transition step
 */
AnyType
sparse_direct_linear_system_transition::run(AnyType &args) {

    SparseDirectLinearSystemTransitionState<MutableArrayHandle<double> > 
                                                              state = args[0];
    int32_t row_id = args[1].getAs<int32_t>();
    int32_t col_id = args[2].getAs<int32_t>();
    double value = args[3].getAs<double>();
    double _b = args[4].getAs<double>();


    
    // When the segment recieves the first non-zero in the sparse matrix
    // we initialize the state
    if (state.nnz_processed == 0) {
        
        int32_t num_equations = args[5].getAs<int32_t>();
        int32_t num_vars = args[6].getAs<int32_t>();
        int32_t total_nnz = args[7].getAs<int32_t>();
        int algorithm = args[8].getAs<int>();
        state.initialize(*this, 
                        num_vars,
                        num_equations,
                        total_nnz
                        );
        state.algorithm = algorithm;
        state.NNZA = total_nnz;
        state.numVars = num_vars;
        state.numEquations = num_equations;
        state.b_stored.setZero();
        state.b.setZero();

    }

	// Now do the transition step
	// First we create a block of zeros in memory
	// and then add the vector in the appropriate position
    ColumnVector bvec(static_cast<uint32_t>(state.numEquations));
    ColumnVector rvec(static_cast<uint32_t>(state.NNZA));
    ColumnVector cvec(static_cast<uint32_t>(state.NNZA));
    ColumnVector vvec(static_cast<uint32_t>(state.NNZA));
    
    bvec.setZero();
    rvec.setZero();
    cvec.setZero();
    vvec.setZero();

    rvec(state.nnz_processed) = row_id;
    cvec(state.nnz_processed) = col_id;
    vvec(state.nnz_processed) = value;

    if (state.b_stored(row_id) == 0) {
        bvec(row_id) = _b;
        state.b_stored(row_id) = 1; 
        state.b += bvec;
    }

    // Build the vector & matrices based on row_id