Actual source code: ex171f.F
petsc-3.12.1 2019-10-22
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
16: call PetscInitialize(PETSC_NULL_CHARACTER,ierr)
17: if (ierr .ne. 0) then
18: print*,'Unable to initialize PETSc'
19: stop
20: endif
22: nloc = 12
23: call VecCreateMPI(PETSC_COMM_WORLD,nloc,PETSC_DETERMINE,v(1),ierr)
24: one = 1.0
25: call VecSet(v(1),one,ierr)
26: call VecNormalize(v(1),norm,ierr)
27: on = 1
28: call MatNullSpaceCreate(PETSC_COMM_WORLD,PETSC_FALSE,on,v(1), &
29: & nsp,ierr)
30: call MatNullSpaceDestroy(nsp,ierr)
31: call VecDestroy(v(1),ierr)
32: call PetscFinalize(ierr)
33: end
35: !/*TEST
36: !
37: ! test:
38: !
39: !TEST*/