Lines Matching refs:row_nnz
67 PetscCall(PetscMalloc1(nrows, &result->row_nnz)); in spbas_allocate_pattern()
111 for (i = 1; i < nrows; i++) result->icols[i] = result->icols[i - 1] + result->row_nnz[i - 1]; in spbas_allocate_data()
120 … for (i = 1; i < nrows; i++) result->values[i] = result->values[i - 1] + result->row_nnz[i - 1]; in spbas_allocate_data()
124 r_nnz = result->row_nnz[i]; in spbas_allocate_data()
129 r_nnz = result->row_nnz[i]; in spbas_allocate_data()
257 PetscInt row_nnz; in spbas_compress_pattern() local
285 B->row_nnz[i] = irow_in[i + 1] - irow_in[i]; in spbas_compress_pattern()
305 if (used[i]) B->n_alloc_icol += B->row_nnz[i]; in spbas_compress_pattern()
315 row_nnz = B->row_nnz[i]; in spbas_compress_pattern()
317 for (j = 0; j < row_nnz; j++) B->icols[i][j] = icols[j]; in spbas_compress_pattern()
319 for (j = 0; j < row_nnz; j++) B->icols[i][j] = icols[j] - i; in spbas_compress_pattern()
321 for (j = 0; j < row_nnz; j++) B->icols[i][j] = icols[j] - icols[0]; in spbas_compress_pattern()
323 ptr += B->row_nnz[i]; in spbas_compress_pattern()
366 PetscCall(PetscFree(matrix.row_nnz)); in spbas_delete()
404 r_nnz = matrix_A.row_nnz[i]; in spbas_matrix_to_crs()
457 for (i = 0; i < nrows; i++) result->row_nnz[i] = 0; in spbas_transpose()
460 r_nnz = in_matrix.row_nnz[i]; in spbas_transpose()
463 for (j = 0; j < r_nnz; j++) result->row_nnz[irow[j]]++; in spbas_transpose()
465 for (j = 0; j < r_nnz; j++) result->row_nnz[i + irow[j]]++; in spbas_transpose()
468 for (j = 0; j < r_nnz; j++) result->row_nnz[icol0 + irow[j]]++; in spbas_transpose()
476 for (i = 0; i < nrows; i++) result->row_nnz[i] = 0; in spbas_transpose()
481 r_nnz = in_matrix.row_nnz[i]; in spbas_transpose()
490 result->icols[k][result->row_nnz[k]] = i; in spbas_transpose()
491 result->values[k][result->row_nnz[k]] = val[j]; in spbas_transpose()
492 result->row_nnz[k]++; in spbas_transpose()
497 r_nnz = in_matrix.row_nnz[i]; in spbas_transpose()
506 result->icols[k][result->row_nnz[k]] = i; in spbas_transpose()
507 result->row_nnz[k]++; in spbas_transpose()
634 PetscInt *row_nnz; in spbas_apply_reordering_rows() local
643 PetscCall(PetscMalloc1(nrows, &row_nnz)); in spbas_apply_reordering_rows()
650 row_nnz[i] = matrix_A->row_nnz[ip]; in spbas_apply_reordering_rows()
651 for (j = 0; j < row_nnz[i]; j++) icols[i][j] += ip - i; in spbas_apply_reordering_rows()
656 PetscCall(PetscFree(matrix_A->row_nnz)); in spbas_apply_reordering_rows()
660 matrix_A->row_nnz = row_nnz; in spbas_apply_reordering_rows()
672 PetscInt row_nnz; in spbas_apply_reordering_cols() local
682 row_nnz = matrix_A->row_nnz[i]; in spbas_apply_reordering_cols()
685 for (j = 0; j < row_nnz; j++) icols[j] = permutation[i + icols[j]] - i; in spbas_apply_reordering_cols()
686 PetscCall(spbas_mergesort(row_nnz, icols, vals)); in spbas_apply_reordering_cols()
719 for (i = 0; i < nrows; i++) retval.row_nnz[i] = ai[i + 1] - ai[i]; in spbas_pattern_only()
747 nnz = in_matrix->row_nnz[row]; in spbas_mark_row_power()
812 nnz = in_matrix.row_nnz[i]; in spbas_power()
822 retval.row_nnz[i] = nnz; in spbas_power()
853 …for (jstart = 0; (jstart < inout_matrix->row_nnz[i]) && (inout_matrix->icols[i][jstart] < 0); jsta… in spbas_keep_upper()
855 …for (j = 0; j < inout_matrix->row_nnz[i] - jstart; j++) inout_matrix->icols[i][j] = inout_matrix->… in spbas_keep_upper()
858 …for (j = 0; j < inout_matrix->row_nnz[i] - jstart; j++) inout_matrix->values[i][j] = inout_matrix-… in spbas_keep_upper()
861 inout_matrix->row_nnz[i] -= jstart; in spbas_keep_upper()
863 …s[i] = (PetscInt *)realloc((void *)inout_matrix->icols[i], inout_matrix->row_nnz[i] * sizeof(Petsc… in spbas_keep_upper()
865 … = (PetscScalar *)realloc((void *)inout_matrix->values[i], inout_matrix->row_nnz[i] * sizeof(Petsc… in spbas_keep_upper()