Actual source code: ex171f.F

petsc-3.12.2 2019-11-22
Report Typos and Errors
  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*/