mod_PETSc.F90

Module petscmanagement

      type userctx
#include <finclude/petscsysdef.h>
#include <finclude/petscvecdef.h>
#include <finclude/petscdmdef.h>
        petscint xs,xe,xm,gxs,gxe,gxm,mx
        petscmpiint rank
        petscmpiint num_procs
      end type userctx

End Module
 



Module global_x
Implicit None
#include <finclude/petscsysdef.h>
#include <finclude/petscsnesdef.h>
#include <finclude/petscvecdef.h>    

snes     :: snes                 !nonlinear solver context
vec      :: x                    !array of unknowns  
vec      :: r                    !array of residuals
petscint :: ngrid                !number of grid points
petscint :: ngp                  !number of ghost points
 
DOUBLE PRECISION :: timer        ! timer
DOUBLE PRECISION :: timer_old    ! timer
DOUBLE PRECISION :: total_time   ! timer


End Module global_x



Module f90moduleinterfaces
Use petscmanagement
#include <finclude/petscsnesdef.h>

      Interface snessetapplicationcontext
        Subroutine snessetapplicationcontext(snes,ctx,ierr)
        use petscmanagement
          snes snes
          type(userctx) ctx
          petscerrorcode ierr
        End Subroutine
      End Interface snessetapplicationcontext

      Interface snesgetapplicationcontext
        Subroutine snesgetapplicationcontext(snes,ctx,ierr)
        use petscmanagement
          snes snes
          type(userctx), pointer :: ctx
          petscerrorcode ierr
        End Subroutine
      End Interface snesgetapplicationcontext

End Module f90moduleinterfaces