#include : Over-estimate of the number of nonzeros needed toierr=PetscFree(lvec);CHKERRQ(ierr); ierr=PetscFree(val);CHKERRQ(ierr);
ierr = spbas_cholesky_garbage_collect( &retval, nrows, &n_row_alloc_ok, &n_alloc_used, max_row_nnz);CHKERRQ(ierr); ierr=PetscFree(max_row_nnz);CHKERRQ(ierr);
/* Place the inverse of the diagonals in the matrix */ for (i=0; i<nrows; i++) { r_nnz = retval.row_nnz[i]; retval.values[i][r_nnz-1] = 1.0 / diag[i]; for (j=0; j<r_nnz-1; j++) { retval.values[i][j] *= -1; } } ierr=PetscFree(diag);CHKERRQ(ierr); *matrix_L = retval; PetscFunctionReturn(0); }
Level:none
Location:src/mat/impls/aij/seq/bas/spbas_cholesky.h
Index of all Mat routines
Table of Contents for all manual pages
Index of all manual pages