Fortran reset

Makefile

@@ -138,7 +138,7 @@ OBJS := $(OBJS) $(SRCDIR)/PETSc_Helper.o \
 		  $(SRCDIR)/PCPFLAREINV.o	
 
 # Define a variable containing all the tests
-export TEST_TARGETS = ex12f ex6f ex6f_getcoeffs ex6 adv_1d adv_diff_2d ex6_cf_splitting adv_diff_cg_supg matrandom
+export TEST_TARGETS = ex12f ex6f ex6f_getcoeffs ex6 adv_1d adv_diff_2d ex6_cf_splitting adv_diff_cg_supg matrandom matrandom_check_reset
 # Include kokkos examples
 ifeq ($(PETSC_HAVE_KOKKOS),1)
 export TEST_TARGETS := $(TEST_TARGETS) adv_1dk

src/C_PETSc_Interfaces.F90

@@ -250,6 +250,16 @@ end subroutine PCGetSetupCalled_c
 
    end interface
 
+   interface   
+
+      subroutine PCMarkNotSetUp_c(A_array) &
+         bind(c, name="PCMarkNotSetUp_c")
+         use iso_c_binding
+         integer(c_long_long) :: A_array
+      end subroutine PCMarkNotSetUp_c         
+
+   end interface   
+
    interface   
 
       subroutine MatGetDiagonalOnly_c(A_array, diag_only) &

src/C_PETSc_Routines.c

@@ -470,6 +470,12 @@ PETSC_INTERN void PCGetSetupCalled_c(PC *pc, PetscInt *setupcalled)
    *setupcalled = (*pc)->setupcalled;
 }
 
+// Mark that we need to do setup again (e.g. after a reset)
+PETSC_INTERN void PCMarkNotSetUp_c(PC *pc)
+{
+   (*pc)->setupcalled = PETSC_FALSE;
+}
+
 // Gets the number of nonzeros in the local 
 PETSC_INTERN PetscErrorCode MatGetNNZs_local_c(Mat *A, PetscInt *nnzs)
 {