Actual source code: ex171f.F
petsc-3.7.2 2016-06-05
1: !
2: ! This program tests MatNullSpaceCreate()
3: !
4: program main
5: implicit none
7: #include <petsc/finclude/petscsys.h>
8: #include <petsc/finclude/petscis.h>
9: #include <petsc/finclude/petscvec.h>
10: #include <petsc/finclude/petscmat.h>
11: #include <petsc/finclude/petscviewer.h>
13: PetscErrorCode ierr
14: MatNullSpace nsp
15: Vec v(1)
16: PetscInt nloc
17: PetscScalar one
18: PetscReal norm
20: call PetscInitialize(PETSC_NULL_CHARACTER,ierr)
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: call MatNullSpaceCreate(PETSC_COMM_WORLD,PETSC_FALSE,1,v(1), &
28: & nsp,ierr)
29: call MatNullSpaceDestroy(nsp,ierr)
30: call VecDestroy(v(1),ierr)
31: call PetscFinalize(ierr)
32: end