Lines Matching refs:i
96 PetscInt i; in spbas_allocate_data() local
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()
123 for (i = 0; i < nrows; i++) { in spbas_allocate_data()
124 r_nnz = result->row_nnz[i]; in spbas_allocate_data()
125 PetscCall(PetscMalloc1(r_nnz, &result->icols[i])); in spbas_allocate_data()
128 for (i = 0; i < nrows; i++) { in spbas_allocate_data()
129 r_nnz = result->row_nnz[i]; in spbas_allocate_data()
130 PetscCall(PetscMalloc1(r_nnz, &result->values[i])); in spbas_allocate_data()
182 PetscInt i, i1, i2; /* Loop counters for (partly) sorted arrays */ in spbas_mergesort_icols() local
216 for (i = istart; i < i2end; i++) { in spbas_mergesort_icols()
218 ihlp1[i] = ihlp2[i1]; in spbas_mergesort_icols()
221 ihlp1[i] = ihlp2[i2]; in spbas_mergesort_icols()
224 ihlp1[i] = ihlp2[i1]; in spbas_mergesort_icols()
238 for (i = 0; i < nrows; i++) isort[i] = ihlp2[i]; in spbas_mergesort_icols()
261 PetscInt i, j; in spbas_compress_pattern() local
276 for (i = 0; i < nrows; i++) B->icol0[i] = icol_in[irow_in[i]]; in spbas_compress_pattern()
284 for (i = 0; i < nrows; i++) { in spbas_compress_pattern()
285 B->row_nnz[i] = irow_in[i + 1] - irow_in[i]; in spbas_compress_pattern()
286 isort[i] = i; in spbas_compress_pattern()
287 ipoint[i] = i; in spbas_compress_pattern()
295 for (i = 1; i < nrows; i++) { in spbas_compress_pattern()
296 …if (spbas_row_order_icol(isort[i - 1], isort[i], irow_in, icol_in, col_idx_type) == 0) ipoint[isor… in spbas_compress_pattern()
300 for (i = 0; i < nrows; i++) used[ipoint[i]] = PETSC_TRUE; in spbas_compress_pattern()
304 for (i = 0; i < nrows; i++) { in spbas_compress_pattern()
305 if (used[i]) B->n_alloc_icol += B->row_nnz[i]; in spbas_compress_pattern()
311 for (i = 0; i < B->nrows; i++) { in spbas_compress_pattern()
312 if (used[i]) { in spbas_compress_pattern()
313 B->icols[i] = &B->alloc_icol[ptr]; in spbas_compress_pattern()
314 icols = &icol_in[irow_in[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()
328 for (i = 0; i < nrows; i++) B->icols[i] = B->icols[ipoint[i]]; in spbas_compress_pattern()
352 PetscInt i; in spbas_delete() local
359 for (i = 0; i < matrix.nrows; i++) PetscCall(PetscFree(matrix.icols[i])); in spbas_delete()
362 for (i = 0; i < matrix.nrows; i++) PetscCall(PetscFree(matrix.values[i])); in spbas_delete()
381 PetscInt i, j, r_nnz, i0; in spbas_matrix_to_crs() local
403 for (i = 0; i < nrows; i++) { in spbas_matrix_to_crs()
404 r_nnz = matrix_A.row_nnz[i]; in spbas_matrix_to_crs()
405 i0 = irow[i]; in spbas_matrix_to_crs()
406 irow[i + 1] = i0 + r_nnz; in spbas_matrix_to_crs()
407 icol_A = matrix_A.icols[i]; in spbas_matrix_to_crs()
410 val_A = matrix_A.values[i]; in spbas_matrix_to_crs()
420 for (j = 0; j < r_nnz; j++) icol[i0 + j] += i; in spbas_matrix_to_crs()
422 i0 = matrix_A.icol0[i]; in spbas_matrix_to_crs()
439 PetscInt i, j, k; in spbas_transpose() local
457 for (i = 0; i < nrows; i++) result->row_nnz[i] = 0; in spbas_transpose()
459 for (i = 0; i < ncols; i++) { in spbas_transpose()
460 r_nnz = in_matrix.row_nnz[i]; in spbas_transpose()
461 irow = in_matrix.icols[i]; in spbas_transpose()
465 for (j = 0; j < r_nnz; j++) result->row_nnz[i + irow[j]]++; in spbas_transpose()
467 icol0 = in_matrix.icol0[i]; in spbas_transpose()
476 for (i = 0; i < nrows; i++) result->row_nnz[i] = 0; in spbas_transpose()
480 for (i = 0; i < ncols; i++) { in spbas_transpose()
481 r_nnz = in_matrix.row_nnz[i]; in spbas_transpose()
482 irow = in_matrix.icols[i]; in spbas_transpose()
483 val = in_matrix.values[i]; in spbas_transpose()
486 else if (col_idx_type == SPBAS_DIAGONAL_OFFSETS) icol0 = i; in spbas_transpose()
487 else if (col_idx_type == SPBAS_OFFSET_ARRAY) icol0 = in_matrix.icol0[i]; in spbas_transpose()
490 result->icols[k][result->row_nnz[k]] = i; in spbas_transpose()
496 for (i = 0; i < ncols; i++) { in spbas_transpose()
497 r_nnz = in_matrix.row_nnz[i]; in spbas_transpose()
498 irow = in_matrix.icols[i]; in spbas_transpose()
501 else if (col_idx_type == SPBAS_DIAGONAL_OFFSETS) icol0 = i; in spbas_transpose()
502 else if (col_idx_type == SPBAS_OFFSET_ARRAY) icol0 = in_matrix.icol0[i]; in spbas_transpose()
506 result->icols[k][result->row_nnz[k]] = i; in spbas_transpose()
529 PetscInt i, i1, i2; /* Loop counters for (partly) sorted arrays */ in spbas_mergesort() local
573 for (i = istart; i < i2end; i++) { in spbas_mergesort()
575 ihlp1[i] = ihlp2[i1]; in spbas_mergesort()
576 vhlp1[i] = vhlp2[i1]; in spbas_mergesort()
579 ihlp1[i] = ihlp2[i2]; in spbas_mergesort()
580 vhlp1[i] = vhlp2[i2]; in spbas_mergesort()
583 ihlp1[i] = ihlp2[i1]; in spbas_mergesort()
584 vhlp1[i] = vhlp2[i1]; in spbas_mergesort()
589 for (i = istart; i < i2end; i++) { in spbas_mergesort()
591 ihlp1[i] = ihlp2[i1]; in spbas_mergesort()
594 ihlp1[i] = ihlp2[i2]; in spbas_mergesort()
597 ihlp1[i] = ihlp2[i1]; in spbas_mergesort()
615 for (i = 0; i < nnz; i++) icol[i] = ihlp2[i]; in spbas_mergesort()
617 for (i = 0; i < nnz; i++) val[i] = vhlp2[i]; in spbas_mergesort()
632 PetscInt i, j, ip; in spbas_apply_reordering_rows() local
646 for (i = 0; i < nrows; i++) { in spbas_apply_reordering_rows()
647 ip = permutation[i]; in spbas_apply_reordering_rows()
648 if (do_values) vals[i] = matrix_A->values[ip]; in spbas_apply_reordering_rows()
649 icols[i] = matrix_A->icols[ip]; 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()
670 PetscInt i, j; in spbas_apply_reordering_cols() local
680 for (i = 0; i < nrows; i++) { in spbas_apply_reordering_cols()
681 icols = matrix_A->icols[i]; in spbas_apply_reordering_cols()
682 row_nnz = matrix_A->row_nnz[i]; in spbas_apply_reordering_cols()
683 if (do_values) vals = matrix_A->values[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()
706 PetscInt i, j, i0, r_nnz; in spbas_pattern_only() local
719 for (i = 0; i < nrows; i++) retval.row_nnz[i] = ai[i + 1] - ai[i]; in spbas_pattern_only()
722 for (i = 0; i < retval.nrows; i++) { in spbas_pattern_only()
723 i0 = ai[i]; in spbas_pattern_only()
724 r_nnz = ai[i + 1] - i0; in spbas_pattern_only()
726 for (j = 0; j < r_nnz; j++) retval.icols[i][j] = aj[i0 + j] - i; in spbas_pattern_only()
744 PetscInt i, j, nnz; in spbas_mark_row_power() local
751 for (i = 0; i < nnz; i++) { in spbas_mark_row_power()
752 j = row + in_matrix->icols[row][i]; in spbas_mark_row_power()
754 …PetscCall(spbas_mark_row_power(iwork, row + in_matrix->icols[row][i], in_matrix, marker / 2, minmr… in spbas_mark_row_power()
760 for (i = 0; i < nnz; i++) { in spbas_mark_row_power()
761 j = row + in_matrix->icols[row][i]; in spbas_mark_row_power()
780 PetscInt i, j, kend; in spbas_power() local
807 for (i = 1; i < power; i++) marker *= 2; in spbas_power()
809 for (i = 0; i < nrows; i++) { in spbas_power()
812 nnz = in_matrix.row_nnz[i]; in spbas_power()
813 kend = i + in_matrix.icols[i][nnz - 1]; in spbas_power()
815 PetscCall(spbas_mark_row_power(iwork, i, &in_matrix, marker, i, maxmrk)); in spbas_power()
819 for (j = i; j < maxmrk; j++) nnz += (iwork[j] != 0); in spbas_power()
822 retval.row_nnz[i] = nnz; in spbas_power()
823 PetscCall(PetscMalloc1(nnz, &retval.icols[i])); in spbas_power()
827 for (j = i; j < maxmrk; j++) { in spbas_power()
829 retval.icols[i][inz] = j - i; in spbas_power()
847 PetscInt i, j; in spbas_keep_upper() local
852 for (i = 0; i < inout_matrix->nrows; i++) { in spbas_keep_upper()
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 …inout_matrix->icols[i] = (PetscInt *)realloc((void *)inout_matrix->icols[i], inout_matrix->row_nnz… in spbas_keep_upper()
865 …rix->values) inout_matrix->values[i] = (PetscScalar *)realloc((void *)inout_matrix->values[i], ino… in spbas_keep_upper()