Actual source code: ex171f.F

petsc-3.7.4 2016-10-02
Report Typos and Errors
  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