using PETSC_NULL_OBJECT instead of dummy where possible

code/spectral_damage.f90

@@ -81,7 +81,6 @@ subroutine spectral_damage_init()
  DM :: damage_grid
  Vec :: uBound, lBound
  PetscErrorCode :: ierr
- PetscObject    :: dummy
  character(len=100) :: snes_type
 
  external :: &
@@ -124,7 +123,8 @@ subroutine spectral_damage_init()
  CHKERRQ(ierr)
  call SNESSetDM(damage_snes,damage_grid,ierr); CHKERRQ(ierr)                                        !< connect snes to da
  call DMCreateGlobalVector(damage_grid,solution,ierr); CHKERRQ(ierr)                                !< global solution vector (grid x 1, i.e. every def grad tensor)
- call DMDASNESSetFunctionLocal(damage_grid,INSERT_VALUES,spectral_damage_formResidual,dummy,ierr)   !< residual vector of same shape as solution vector
+ call DMDASNESSetFunctionLocal(damage_grid,INSERT_VALUES,spectral_damage_formResidual,&
+                                                                            PETSC_NULL_OBJECT,ierr) !< residual vector of same shape as solution vector
  CHKERRQ(ierr) 
  call SNESSetFromOptions(damage_snes,ierr); CHKERRQ(ierr)                                           !< pull it all together with additional cli arguments
  call SNESGetType(damage_snes,snes_type,ierr); CHKERRQ(ierr)

code/spectral_mech_AL.f90

@@ -116,7 +116,6 @@ subroutine AL_init
    temp33_Real = 0.0_pReal
 
  PetscErrorCode :: ierr
- PetscObject :: dummy
  PetscScalar, pointer, dimension(:,:,:,:) :: xx_psc, F, F_lambda
  integer(pInt), dimension(:), allocatable :: localK  
  integer(pInt) :: proc
@@ -165,9 +164,9 @@ subroutine AL_init
  CHKERRQ(ierr)
  call SNESSetDM(snes,da,ierr); CHKERRQ(ierr)
  call DMCreateGlobalVector(da,solution_vec,ierr); CHKERRQ(ierr)
- call DMDASNESSetFunctionLocal(da,INSERT_VALUES,AL_formResidual,dummy,ierr)
+ call DMDASNESSetFunctionLocal(da,INSERT_VALUES,AL_formResidual,PETSC_NULL_OBJECT,ierr)
  CHKERRQ(ierr)
- call SNESSetConvergenceTest(snes,AL_converged,dummy,PETSC_NULL_FUNCTION,ierr)
+ call SNESSetConvergenceTest(snes,AL_converged,PETSC_NULL_OBJECT,PETSC_NULL_FUNCTION,ierr)
  CHKERRQ(ierr)
  call SNESSetFromOptions(snes,ierr); CHKERRQ(ierr)
 

code/spectral_mech_Basic.f90

@@ -105,7 +105,6 @@ subroutine basicPETSc_init
    temp33_Real = 0.0_pReal
 
  PetscErrorCode :: ierr
- PetscObject :: dummy 
  PetscScalar, pointer, dimension(:,:,:,:)   ::  F
 
  integer(pInt), dimension(:), allocatable :: localK  
@@ -153,10 +152,10 @@ subroutine basicPETSc_init
  CHKERRQ(ierr)
  call SNESSetDM(snes,da,ierr); CHKERRQ(ierr)
  call DMCreateGlobalVector(da,solution_vec,ierr); CHKERRQ(ierr)                                     ! global solution vector (grid x 9, i.e. every def grad tensor)
- call DMDASNESSetFunctionLocal(da,INSERT_VALUES,BasicPETSC_formResidual,dummy,ierr)                 ! residual vector of same shape as solution vector
+ call DMDASNESSetFunctionLocal(da,INSERT_VALUES,BasicPETSC_formResidual,PETSC_NULL_OBJECT,ierr)     ! residual vector of same shape as solution vector
  CHKERRQ(ierr) 
  call SNESSetDM(snes,da,ierr); CHKERRQ(ierr)                                                        ! connect snes to da
- call SNESSetConvergenceTest(snes,BasicPETSC_converged,dummy,PETSC_NULL_FUNCTION,ierr)              ! specify custom convergence check function "_converged"
+ call SNESSetConvergenceTest(snes,BasicPETSC_converged,PETSC_NULL_OBJECT,PETSC_NULL_FUNCTION,ierr)  ! specify custom convergence check function "_converged"
  CHKERRQ(ierr)
  call SNESSetFromOptions(snes,ierr); CHKERRQ(ierr)                                                  ! pull it all together with additional cli arguments
 

code/spectral_mech_Polarisation.f90

@@ -116,7 +116,6 @@ subroutine Polarisation_init
    temp33_Real = 0.0_pReal
 
  PetscErrorCode :: ierr
- PetscObject :: dummy
  PetscScalar, pointer, dimension(:,:,:,:) :: xx_psc, F, F_tau
  integer(pInt), dimension(:), allocatable :: localK  
  integer(pInt) :: proc
@@ -165,9 +164,9 @@ subroutine Polarisation_init
  CHKERRQ(ierr)
  call SNESSetDM(snes,da,ierr); CHKERRQ(ierr)
  call DMCreateGlobalVector(da,solution_vec,ierr); CHKERRQ(ierr)
- call DMDASNESSetFunctionLocal(da,INSERT_VALUES,Polarisation_formResidual,dummy,ierr)
+ call DMDASNESSetFunctionLocal(da,INSERT_VALUES,Polarisation_formResidual,PETSC_NULL_OBJECT,ierr)
  CHKERRQ(ierr)
- call SNESSetConvergenceTest(snes,Polarisation_converged,dummy,PETSC_NULL_FUNCTION,ierr)
+ call SNESSetConvergenceTest(snes,Polarisation_converged,PETSC_NULL_OBJECT,PETSC_NULL_FUNCTION,ierr)
  CHKERRQ(ierr)
  call SNESSetFromOptions(snes,ierr); CHKERRQ(ierr)
 

code/spectral_thermal.f90

@@ -86,7 +86,6 @@ subroutine spectral_thermal_init
  DM :: thermal_grid
  PetscScalar,  dimension(:,:,:), pointer     :: x_scal
  PetscErrorCode :: ierr
- PetscObject    :: dummy
 
  external :: &
    SNESCreate, &
@@ -123,7 +122,8 @@ subroutine spectral_thermal_init
  CHKERRQ(ierr)
  call SNESSetDM(thermal_snes,thermal_grid,ierr); CHKERRQ(ierr)                                      ! connect snes to da
  call DMCreateGlobalVector(thermal_grid,solution        ,ierr); CHKERRQ(ierr)                       ! global solution vector (grid x 1, i.e. every def grad tensor)
- call DMDASNESSetFunctionLocal(thermal_grid,INSERT_VALUES,spectral_thermal_formResidual,dummy,ierr) ! residual vector of same shape as solution vector
+ call DMDASNESSetFunctionLocal(thermal_grid,INSERT_VALUES,spectral_thermal_formResidual,&
+                                                                            PETSC_NULL_OBJECT,ierr) ! residual vector of same shape as solution vector
  CHKERRQ(ierr) 
  call SNESSetFromOptions(thermal_snes,ierr); CHKERRQ(ierr)                                          ! pull it all together with additional cli arguments