1: !
  2: !   This program tests MatNullSpaceCreate()
  3: !
  4: program main
  5: #include <petsc/finclude/petscmat.h>
  6:   use petscmat
  7:   implicit none

  9:   PetscErrorCode ierr
 10:   MatNullSpace nsp
 11:   Vec v(1)
 12:   PetscInt nloc, on
 13:   PetscScalar one
 14:   PetscReal norm
 15:   Vec, pointer :: vnsp(:)

 17:   PetscCallA(PetscInitialize(ierr))

 19:   nloc = 12
 20:   on = 1
 21:   PetscCallA(VecCreateFromOptions(PETSC_COMM_WORLD, PETSC_NULL_CHARACTER, on, nloc, PETSC_DETERMINE, v(1), ierr))
 22:   one = 1.0
 23:   PetscCallA(VecSet(v(1), one, ierr))
 24:   PetscCallA(VecNormalize(v(1), norm, ierr))
 25:   PetscCallA(MatNullSpaceCreate(PETSC_COMM_WORLD, PETSC_FALSE, on, [v], nsp, ierr))
 26:   PetscCallA(MatNullSpaceGetVecs(nsp, PETSC_NULL_BOOL, PETSC_NULL_INTEGER, vnsp, ierr))
 27:   PetscCallA(MatNullSpaceRestoreVecs(nsp, PETSC_NULL_BOOL, PETSC_NULL_INTEGER, vnsp, ierr))
 28:   PetscCallA(MatNullSpaceDestroy(nsp, ierr))
 29:   PetscCallA(VecDestroy(v(1), ierr))
 30:   PetscCallA(PetscFinalize(ierr))
 31: end

 33: !/*TEST
 34: !
 35: !   test:
 36: !      output_file: output/empty.out
 37: !
 38: !TEST*/