36
36
void Paso_SystemMatrix_nullifyRows(Paso_SystemMatrix* A, double* mask_row, double main_diagonal_value) {
40
40
if (A ->col_block_size==1 && A ->row_block_size ==1) {
41
41
if (A->type & MATRIX_FORMAT_CSC) {
44
44
} else if (A->type & MATRIX_FORMAT_TRILINOS_CRS) {
49
49
Paso_SystemMatrix_startRowCollect(A,mask_row);
50
50
Paso_SparseMatrix_nullifyRows_CSR_BLK1(A->mainBlock,mask_row,main_diagonal_value);
51
51
Paso_SparseMatrix_nullifyRows_CSR_BLK1(A->col_coupleBlock,mask_row,0.);