Changes: 3.6 — PETSc 3.17.0 documentation

# Changes: 3.6

General:

- User makefiles must be updated. You must change the lines

  - include ${PETSC_DIR}/conf/variables and

  - include ${PETSC_DIR}/conf/rules

  to

  - include ${PETSC_DIR}/lib/petsc/conf/variables and

  - include ${PETSC_DIR}/lib/petsc/conf/rules

- PetscOptionsBool() no longer sets the value of the boolean based on the value passed in

- Script for running MPIUni jobs is now bin/petsc-mpiexec.uni

- Removed all threadcomm support including –with-pthreadclasses and –with-openmpclasses configure arguments

- Updated MUMPS version to 5.0.0 - and updated metis to 5.1.0, parmetis to 4.0.3-p1 - in sync with MUMPS.

- Updated to hypre-2.10.0b, scotch_6.0.3, Elemental-0.85

Configure/Build:

- Add toplevel makefile targets clean and distclean. Target clean deletes petsc libraries built by makefile target all. Target distclean deletes all build files created by configure and [make all]. Target clean in example directories continues to work as before [deletes object files and executables]

IS:

PetscDraw:

- PetscDrawBoxedString() changed to PetscDrawStringBoxed() added PetscDrawBar object for drawing bar graphs added PetscDrawStringCentered()

PF:

Vec:

VecScatter:

PetscSection:

Mat:

- MatGetVecs() replaced with MatCreateVecs() because it actually does create the Vecs that need to be destroyed

- MatCreateMPIAIJConcatenateSeqAIJ and MatCreateMPIBAIJConcatenateSeqBAIJ replaced with MatCreateMPIMatConcatenateSeqMat

- MatGetRedundantMatrix() replaced with MatCreateRedundantMatrix()

- Added support for MUMPS sequential solver in Schur complement mode

- MatGetSubMatricesMPI() will extract parallel submatrices on subcommunicators of the supplied IS objects. Unsorted, but not repeated indices are okay.

- Removed MatMFFDAddNullSpace() just use MatSetNullSpace()

PC:

- Removed -pc_hypre_type euclid due to bit-rot

- Upgraded to SuiteSparse 4.4.1; Cholmod supports using GPUs (developer repository only feature) –with-cuda –download-suitesparse-gpu –with-64-bit-indices; -mat_cholmod_useGPU 1 or 0; 1 by default when configured for it

- Upgraded to Superlu_DIST 4.0; Superlu_DIST supports using GPUs (developer repository only feature) –with-cuda –with-openmp –download-superlu_dist-gpu (developer repository only feature) –with-64-bit-indices -mat_cholmod_useGPU 1 or 0; 1 by default when configured for it

- Added PCBDDCSetChangeOfBasisMat for user defined change of basis

- PCBDDC: added adaptive selection of constraints: it is available by command line options

- PCJacobiSetUseAbs() now takes a PetscBool argument allowing toggling the option

- PCJacobiSetUseRowMax() and PCJacobiSetUseRowSum() have been merged into PCJacobiSetType()

- PCFactorSetUseInplace() now takes a PetscBool argument. Added PCFactorGetUseInplace()

- PCFactorSetAllowDiagonalFill() now takes a PetscBool argument. Added PCFactorGetAllowDiagonalFill()

- PCEisenstatNoDiagonalScaling() changed to PCEisenstatSetNoDiagonalScaling() and takes a boolean

- Added Interface to AMS solver from HYPRE for Nedelec discretizations (-pc_hypre_type ams)

- Added Interface to ADS solver from HYPRE for Raviart-Thomas discretizations (-pc_hypre_type ads)

- Contributed MKL Pardiso interface providing PCCholesky support. Useable with: -pc_type cholesky -mat_type sbaij -pc_factor_mat_solver_package mkl_pardiso

- PCGASMSetTotalSubdomains() can now create multirank subdomains if their number is less than comm size.

- PCGAMGSetReuseProl() is changed to PCGAMGSetReuseInterpolation()

KSP:

- Removed KSPSetNullSpace() always use MatSetNullSpace() provided to the FIRST matrix argument to KSPSetOperators(), or SNESSetJacobian(), or TSSetIJacobian()

- Added MatSetTransposeNullSpace() to indicate the null space of the transpose operator

- Removed KSPSPECEST since its functionality is replaced by the esteig functionality within KSPChebyshev

- Unified KSPChebyshev routines that work with estimating eigenvalues with the prefix name KSPChebyshevEstEig… For example KSPChebyshevSetEstimateEigenvalues() is now KSPChebyshevEstEigSet() and KSPChebyshevSetEstimateEigenvaluesRandom() is now KSPChebyshevEstEigSetRandom() similarly -ksp_chebyshev_estimate_eigenvalues_random is -ksp_chebyshev_esteig_random. The previous for the (GMRES) KSP used to estimate the eigenvalues is now esteig_ not est_

- KSPGetVecs() replaced with KSPCreateVecs() because it actually does create the Vecs that need to be destroyed

SNES:

- Added KSPMonitorSNES() with command line option -ksp_monitor_snes and -ksp_monitor_snes_lg to monitor SNES residual norm at each linear iteration SNESMonitorVI() becomes SNESVIMonitor()

- Added SNESVIMonitorResidual() and -snes_vi_monitor_residual

SNESLineSearch:

TS:

- The postevent routine set with TSSetEventMonitor() now takes an additional PetscBool argument ‘forwardsolve’. The forwardsolve flag is used to distinguish between a TS forward solve (forwardsolve = 1) and adjoint solve (forwardsolve = 0).

- Added TSAdjointSolve() for adjoint sensitivity analysis

- EquationType is introduced to distinguish between ODEs and DAEs. It needs to be set accordingly before solving the problem through TSSetEquationType(). It currently only affects arkimex methods.

DM/DA:

- DMCreateInjection() now returns a Mat, rather than a VecScatter

- The MatType argument is removed from DMCreateMatrix(), you can use DMSetMatType() to indicate the type you want used with a DM, defaults to MATAIJ

- You can now use DMDASetAOType() to indicate the type of AO you want used with a DMDA (defaults to AOBASIC)

DMPlex:

- Can now redistribute parallel meshes

- Can now increase the overlap of parallel meshes

- Added new mesh formats, Fluent CAS and Gmsh

PetscViewer:

- HDF5 viewer will save Vecs and DMDA Vecs with a dimension of 1 if the bs/dof of the vector is 1. To always include this dimension as a specific dimension in the HDF5 file even if if it is of size 1 use PetscViewerHDF5SetBaseDimension2() or -viewer_hdf5_base_dimension2 true

SYS:

- PetscOptionsHead() now takes a PetscOptions argument.

- PetscObjectAddOptionsHandler() now takes an event handler that includes a PetscOptions argument. PetscObjectProcessOptionsHandlers() requires this as input.

AO:

Fortran:

- Fortran include files are now in include/petsc/finclude instead of include/finclude. Thus replace uses of #include “finclude/xxx.h” with #include “petsc/finclude/xxx.h”. Reason for change: to namespace the finclude directory with PETSc for –prefix installs of PETSc and for packaging systems

previous

Changes: 3.7

next

Changes: 3.5