Actual source code: ex171f.F90
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: PetscCallA(PetscInitialize(ierr))
18: nloc = 12
19: PetscCallA(VecCreateMPI(PETSC_COMM_WORLD,nloc,PETSC_DETERMINE,v(1),ierr))
20: one = 1.0
21: PetscCallA(VecSet(v(1),one,ierr))
22: PetscCallA(VecNormalize(v(1),norm,ierr))
23: on = 1
24: PetscCallA(MatNullSpaceCreate(PETSC_COMM_WORLD,PETSC_FALSE,on,v(1),nsp,ierr))
25: PetscCallA(MatNullSpaceDestroy(nsp,ierr))
26: PetscCallA(VecDestroy(v(1),ierr))
27: PetscCallA(PetscFinalize(ierr))
28: end
30: !/*TEST
31: !
32: ! test:
33: !
34: !TEST*/