grid solver adjustments for PETSc 3.22.x

and a few smaller adjustments

Author: MarDiehl

CMakeLists.txt

@@ -11,7 +11,7 @@ endif()
 project(Prerequisites LANGUAGES)
 set(ENV{PKG_CONFIG_PATH} "$ENV{PETSC_DIR}/$ENV{PETSC_ARCH}/lib/pkgconfig:$ENV{PKG_CONFIG_PATH}")
 pkg_check_modules(PETSC_MIN REQUIRED PETSc>=3.12.0 QUIET) #CMake does not support version range
-pkg_check_modules(PETSC REQUIRED PETSc<3.22.0)
+pkg_check_modules(PETSC REQUIRED PETSc<3.23.0)
 
 pkg_get_variable(CMAKE_Fortran_COMPILER PETSc fcompiler)
 pkg_get_variable(CMAKE_C_COMPILER PETSc ccompiler)

src/CLI.f90

@@ -6,7 +6,7 @@
 !> @brief Parse command line interface for PETSc-based solvers
 !--------------------------------------------------------------------------------------------------
 #define PETSC_MINOR_MIN 12
-#define PETSC_MINOR_MAX 21
+#define PETSC_MINOR_MAX 22
 
 module CLI
   use, intrinsic :: ISO_fortran_env

src/grid/grid_mech_spectral_Galerkin.f90

@@ -135,6 +135,7 @@ subroutine MatShellSetOperation(mat,op_num,op_callback,ierr)
     end subroutine MatShellSetOperation
   end interface MatShellSetOperation
 
+#if (PETSC_VERSION_MAJOR==3 && PETSC_VERSION_MINOR<22)
   interface SNESSetJacobian
     subroutine SNESSetJacobian(snes_mech,A,P,jac_callback,ctx,ierr)
       use petscsnes
@@ -145,6 +146,7 @@ subroutine SNESSetJacobian(snes_mech,A,P,jac_callback,ctx,ierr)
       PetscErrorCode :: ierr
     end subroutine SNESSetJacobian
   end interface SNESSetJacobian
+#endif
 
   interface SNESSetUpdate
     subroutine SNESSetUpdate(snes_mech,upd_callback,ierr)
@@ -678,14 +680,17 @@ end subroutine GK_op
 subroutine set_F_aim(snes, step, ierr)
   SNES      :: snes
   PetscInt  :: step ! curr completed petsc iter
-  PetscErrorCode :: ierr
+  PetscErrorCode, intent(out) :: ierr
+
 
   real(pREAL), dimension(3,3) :: &
     deltaF_aim
 
   deltaF_aim = math_mul3333xx33(S, P_av - P_aim)
   F_aim = F_aim - deltaF_aim
 
+  ierr = 0_pPETSCERRORCODE
+
 end subroutine set_F_aim
 
 end module grid_mechanical_spectral_Galerkin

src/prec.f90

@@ -23,10 +23,10 @@ module prec
   integer,     parameter :: pI32       = selected_int_kind(9)                                       !< number with at least up to +-1e9 (typically 32 bit)
   integer,     parameter :: pI64       = selected_int_kind(18)                                      !< number with at least up to +-1e18 (typically 64 bit)
 #ifdef PETSC
-  PetscInt,    private   :: dummy_int
-  integer,     parameter :: pPETSCINT  = kind(dummy_int)
-  PetscScalar, private   :: dummy_scalar
-  real(pREAL), parameter, private :: pPETSCSCALAR = kind(dummy_scalar)
+  PetscInt,        private   :: dummy_int
+  PetscErrorCode,  private   :: dummy_error_code
+  integer,         parameter :: pPETSCINT = kind(dummy_int)
+  integer,         parameter :: pPETSCERRORCODE = kind(dummy_error_code)
 #endif
   integer,     parameter :: pSTRLEN = 256                                                           !< default string length
   integer,     parameter :: pPATHLEN   = 4096                                                       !< maximum length of a path name on linux
@@ -251,10 +251,11 @@ subroutine prec_selfTest()
   real(pREAL),   dimension(1) :: f
   integer(pI64), dimension(1) :: i
   real(pREAL),   dimension(2) :: r
+#ifdef PETSC
+  PetscScalar :: dummy_scalar
 
 
-#ifdef PETSC
-  if (pREAL /= pPETSCSCALAR)                error stop 'PETSc and DAMASK scalar datatypes do not match'
+  if (pREAL /= kind(dummy_scalar))          error stop 'PETSc and DAMASK scalar datatypes do not match'
 #endif
   realloc_lhs_test = [1,2]
   if (any(realloc_lhs_test/=[1,2]))         error stop 'LHS allocation'