#include "petscmat.h" PetscErrorCode MatFDColoringGetPerturbedColumns(MatFDColoring coloring,PetscInt *n,const PetscInt *cols[])Not Collective
coloring | - coloring context created with MatFDColoringCreate() |
n | - the number of local columns being perturbed | |
cols | - the column indices, in global numbering |
Note: IF the matrix type is BAIJ, then the block column indices are returned
#include <petsc/finclude/petscmat.h>
use petscmat
PetscInt, pointer :: array(:)
PetscErrorCode ierr
MatFDColoring i
call MatFDColoringGetPerturbedColumnsF90(i,array,ierr)
use the entries of array ...
call MatFDColoringRestorePerturbedColumnsF90(i,array,ierr)