Lines Matching refs:nnz

45   else memreq += matrix.nnz * sizeof(PetscInt);  in spbas_memory_requirement()
51 else memreq += matrix.nnz * sizeof(PetscScalar); in spbas_memory_requirement()
97 PetscInt nnz = result->nnz; in spbas_allocate_data() local
106 result->n_alloc_icol = nnz; in spbas_allocate_data()
108 PetscCall(PetscMalloc1(nnz, &result->alloc_icol)); in spbas_allocate_data()
115 result->n_alloc_val = nnz; in spbas_allocate_data()
117 PetscCall(PetscMalloc1(nnz, &result->alloc_val)); in spbas_allocate_data()
252 PetscInt nnz = irow_in[nrows]; in spbas_compress_pattern() local
253 size_t mem_orig = (nrows + nnz) * sizeof(PetscInt); in spbas_compress_pattern()
268 B->nnz = nnz; in spbas_compress_pattern()
330 …PetscCall(PetscInfo(NULL, " (%g nonzeros per row)\n", (double)((PetscReal)nnz / (PetscReal… in spbas_compress_pattern()
380 PetscInt nnz = matrix_A.nnz; in spbas_matrix_to_crs() local
392 PetscCall(PetscMalloc1(nnz, &icol)); in spbas_matrix_to_crs()
396 PetscCall(PetscMalloc1(nnz, &val)); in spbas_matrix_to_crs()
436 PetscInt nnz = in_matrix.nnz; in spbas_transpose() local
449 result->nnz = nnz; in spbas_transpose()
526 static PetscErrorCode spbas_mergesort(PetscInt nnz, PetscInt *icol, PetscScalar *val) in spbas_mergesort() argument
541 PetscCall(PetscMalloc1(nnz, &ialloc)); in spbas_mergesort()
546 PetscCall(PetscMalloc1(nnz, &valloc)); in spbas_mergesort()
552 for (istep = 1; istep < nnz; istep *= 2) { in spbas_mergesort()
562 for (istart = 0; istart < nnz; istart += 2 * istep) { in spbas_mergesort()
566 if (i1end > nnz) i1end = nnz; in spbas_mergesort()
569 if (i2end > nnz) i2end = nnz; 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()
712 retval.nnz = ai[nrows]; in spbas_pattern_only()
744 PetscInt i, j, nnz; in spbas_mark_row_power() local
747 nnz = in_matrix->row_nnz[row]; in spbas_mark_row_power()
751 for (i = 0; i < nnz; i++) { in spbas_mark_row_power()
760 for (i = 0; i < nnz; i++) { in spbas_mark_row_power()
781 PetscInt nnz, inz; in spbas_power() local
795 retval.nnz = 0; 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()
818 nnz = 0; 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()
834 retval.nnz += nnz; in spbas_power()
866 inout_matrix->nnz -= jstart; in spbas_keep_upper()