Changes: Development

General:

• Change behavior of -options_left when set to true: it no longer triggers a call to PetscOptionsView()

• Change \(PETSC_DIR/lib/petsc/bin/petscfreesharedmemory to \)PETSC_DIR/lib/petsc/bin/petscfreesharedmemory.sh

Configure/Build:

• Make SYCL a devicePackage, i.e., builds --with-sycl now have PETSC_HAVE_DEVICE defined

• Add the option --with-devicelanguage to compile PetscDevice code using either a C or C++ compiler

• Add Caliper, an instrumentation and performance profiling library that can be used to profile Hypre.

• Add typing stubs for the Python extension module petsc4py.PETSc.

Sys:

• Add PetscCallHYPRE() to check HYPRE error codes and print error messages on failure

• Add PetscBTCountSet() to count set bits in PetscBT

Event Logging:

• Add two approaches for GPU energy monitoring: -log_view_gpu_energy and -log_view_gpu_energy_meter

• Add API PetscLogGpuEnergy(), PetscLogGpuEnergyMeter(), PetscLogGpuEnergyMeterBegin() and PetscLogGpuEnergyMeterEnd() for GPU energy monitoring

PetscViewer:

• Change the final argument of PetscViewerGLVisSetFields() to PetscCtxDestroyFn *. This means the destroy function must dereference the argument before operating on it

PetscDraw:

AO:

IS:

VecScatter / PetscSF:

PF:

Vec:

PetscSection:

• Add PetscSectionMigrateData(), akin to DMPlexDistributeData()

PetscPartitioner:

Mat:

• Change the destroy() function argument of MatShellSetMatProductOperation() to type PetscCtxDestroyFn *. This means the destroy function must dereference the argument before operating on it

• Remove MatMissingDiagonal(). Developers should use MatGetDiagonalMarkers_SeqXXX() when the functionality is needed

• Change MatSetOption(A, MAT_HERMITIAN, PETSC_TRUE) for MatSBAIJ to no longer automatically set the option MAT_SYMMETRIC to PETSC_FALSE. It is now the duty of the user to call MatSetOption(A, MAT_SYMMETRIC, PETSC_FALSE) if a MatSBAIJ is Hermitian but not symmetric

• Deprecate -matmatmult_Bbn in favor of -matproduct_batch_size

MatCoarsen:

PC:

• Add multi-precision support for MUMPS. One could use -pc_precision <single, double> to set the precision to be used by MUMPS, which can be different from PetscScalar’s precision

• Add support for MUMPS out-of-core facility with the option -mat_mumps_ooc_tmpdir <dir> and new functions MatMumpsSetOocTmpDir(), MatMumpsGetOocTmpDir()

KSP:

• Remove KSPHPDDMPrecision in favor of PetscPrecision

SNES:

• Change the destroy() function argument of SNESSetConvergenceTest() to type PetscCtxDestroyFn *. This means the destroy function must dereference the argument before operating on it

• Add SNESSetObjectiveDomainError()

• Change SNES_DIVERGED_FNORM_NAN to SNES_DIVERGED_FUNCTION_NANORINF

• Add SNES_DIVERGED_OBJECTIVE_NANORINF

• Add SNES_DIVERGED_OBJECTIVE_DOMAIN

• Add developer functions SNESCheckFunctionDomainError(), SNESLineSearchCheckFunctionDomainError(), SNESCheckObjectiveDomainError(), SNESLineSearchCheckObjectiveDomainError(), SNESCheckJacobianDomainError(), and SNESLineSearchCheckJacobianDomainError()

SNESLineSearch:

TS:

• Add TSPseudoComputeFunction() to get nonlinear residual while avoiding recalculation if possible

• Remove unused TSPseudoVerifyTimeStepDefault()

• Remove TSPseudoComputeTimeStep() and TSPseudoVerifyTimeStep()

• Change the destroy() function argument of TSTrajectorySetTransform() to type PetscCtxDestroyFn *. This means the destroy function must dereference the argument before operating on it

• Correct option -ts_max_reject to -ts_max_step_rejections

• Correct option -ts_dt to -ts_time_step

TAO:

PetscRegressor:

DM/DA:

• Change the final argument of DMShellSetDestroyContext() to PetscCtxDestroyFn *. This means the destroy function must dereference the argument before operating on it

• Add DMLabelGetValueISGlobal() to get globally consistent IS of values in a DMLabel

DMSwarm:

DMPlex:

• Add DMPlexVecGetClosureAtDepth()

• Add an extra communicator argument to DMPlexFilter() to allow extracting local meshes

• Add DMPlexGetLETKFLocalizationMatrix to compute localization weight matrix for LETKF

• Change verticesAdjSaved parameter in DMPlexCreateFromCell*Parallel*() functions to be allocated by function rather than by user

FE/FV:

DMNetwork:

DMStag:

DT:

Fortran:

• Replace ./configure option --with-mpi-f90module-visibility with --with-mpi-ftn-module=<mpi or mpi_f08>

• Add PETSC_INT_KIND and PETSC_MPIINT_KIND

• Fortran code should now use MPIU_Comm instead of MPI_Comm, and similarly for other MPI types, see section “Fortran and MPI” in the users guide

• Fortran interface definitions are now automatically generated for all functions that take context variable arguments, represented in the C source code with a type of PetscCtx, allowing the use of any Fortran derived type (or PETSc object) as the context

• For all PETSc functions XXXGetYYY() that return a context variable as an argument, represented in the C source code with an argument type of PetscCtxRt, a macro is generated used with Interface_XXXGetYYY(AppCtx) which tells the Fortran compiler that a pointer to that derived type type(AppCtx) is returned from the Fortran version of XXXGetYYY(). See src/snes/tutorials/ex5f90.F90`