Changes: 3.20 — PETSc 3.25.1 documentation

Skip to main content
`Ctrl`+`K`

# Changes: 3.20

General:

- Add `PetscOptionsBegin()`, `PetscOptionsEnd()`, `PetscOptionsInt()`, `PetscOptionsBool()`, `PetscOptionsIntArray()`, `PetscOptionsReal()`, `PetscOptionsRealArray()`, and `PetscOptionsScalar()` for Fortran

- Add `PetscAssertPointer()` as a replacement for `PetscValidPointer()`, `PetscValidCharPointer()`, `PetscValidIntPointer()`, `PetscValidInt64Pointer()`, `PetscValidCountPointer()`, `PetscValidBoolPointer()`, `PetscvalidScalarPointer()`, and `PetscValidRealPointer()`

- Convert standard header-guards in PETSc header-files to `#pragma once`. Users relying on specific preprocessor definitions to determine whether particular PETSc headers (e.g. `petscdm.h` used to be guarded by `PETSCDM_H`) have been included will find that this no longer works. While header-guards are technically “public” (insofar that they appear in public header-files), their name and value are considered implementation details

Configure/Build:

- Add support for external-packages to prefer tarball download for regular use - as currently all packages prefer git clones. MPICH is a package using this feature. Here MPICH tarball will be downloaded for regular use. However on providing `--download-mpich-commit=main` option - configure download and builds from MPICH git repository

- Add support for external package SLATE, optional dependency for STRUMPACK

- Add support for external package ZFP, optional dependency for STRUMPACK

- Add support for external package ButterflyPACK, optional dependency for STRUMPACK

Sys:

- Add `PetscDeviceContextGetStreamHandle()` to return a handle to the stream the current device context is using

- Add utility functions `PetscGetCurrentCUDAStream()` and `PetscGetCurrentHIPStream()` to return the GPU stream the current device context is using

- Add `PetscStrcmpAny()` to compare against multiple non-empty strings

- Change arguments 4 and 5 of `PetscViewerBinaryReadAll()` and `PetscViewerBinaryWriteAll()` to `PetscInt64`

- Add `PetscIsCloseAtTolScalar()`

- Remove `PetscTellMyCell()` and related runtime option

- Remove `PetscTextBelt()` and related runtime option

- Remove deprecated `-malloc [no]` startup option

- Remove deprecated `-malloc_info` startup option

- Remove deprecated `-log_summary` option

- Remove `PetscURLShorten()`, it has not worked since 2019

- Move `PetscIntStackCreate()`, `PetscIntStackDestroy()`, `PetscIntStackPush()`, `PetscIntStackPop()`, and `PetscIntStackEmpty()` declarations to public API in `petsclog.h`

- Add `-on_error_malloc_dump` option

Event Logging:

- Add `PetscLogState` interface for describing profiling events and stages

- Add `PetscLogHandler` interface that interprets the profiling events of a `PetscLogState`

- Add `PETSCLOGHANDLERDEFAULT` implementation of `PetscLogHandler`

- Add `PETSCLOGHANDLERNESTED` implementation of `PetscLogHandler`

- Add `PETSCLOGHANDLERTRACE` implementation of `PetscLogHandler`

- Add `PetscLogHandlerCreateTrace()`

- Add `PETSCLOGHANDLERMPE` implementation of `PetscLogHandler`

- Add `PETSCLOGHANDLERPERFSTUBS` implementation of `PetscLogHandler`

- Add `PETSCLOGHANDLERLEGACY` implementation of `PetscLogHandler`

- Add `PetscLogHandlerCreateLegacy()`

- Add `PetscLogPerfstubsBegin()` for log handling with PerfStubs/TAU

- Add `PetscLogLegacyCallbacksBegin()` to transition old log handler callbacks to the `PetscLogHandler`-based approach

- Add `PetscLogHandlerStart()` to connect a `PetscLogHandler` to PETSc’s global logging events

- Add `PetscLogStageGetName()`

- Add `PetscLogEventGetName()`

- Add `PetscLogEventsPause()` and `PetscLogEventsResume()` to isolate a logging event from unrelated in-process events

- Add `PetscLogClassGetClassId()`

- Add `PetscLogClassIdGetName()`

- Deprecate `PetscLogClassPerfInfoClear()` , `PetscLogClassPerfLogCreate()`, and `PetscLogClassPerfLogDestroy()`

- Deprecate `PetscClassPerfLogEnsureSize()`, `PetscClassRegInfoDestroy()`, `PetscClassRegLogCreate()`, `PetscClassRegLogDestroy()`, `PetscClassRegLogGetClass()`, and `PetscClassRegLogRegister()`

- Deprecate `PetscEventPerfInfoClear()` and `PetscEventPerfInfoCopy()`

- Deprecate `PetscEventPerfLogActivate()`, `PetscEventPerfLogActivateClass()`, `PetscEventPerfLogCreate()`, `PetscEventPerfLogDeactivate()`, `PetscEventPerfLogDeactivateClass()`, `PetscEventPerfLogDeactivatePop()`, `PetscEventPerfLogDeactivatePush()`, `PetscEventPerfLogDestroy()`, `PetscEventPerfLogEnsureSize()`, `PetscEventPerfLogGetVisible()`, and `PetscEventPerfLogSetVisible()`

- Deprecate `PetscEventRegLogCreate()`, `PetscEventRegLogDestroy()`, `PetscEventRegLogGetEvent()`, and `PetscEventRegLogRegister()`

- Deprecate `PetscStageLogCreate()`, `PetscStageLogDestroy()`, `PetscStageLogGetActive()`, `PetscStageLogGetClassPerfLog()`, `PetscStageLogGetClassRegLog()`, `PetscStageLogGetEventRegLog()`, `PetscStageLogGetStage()`, `PetscStageLogGetVisible()`, `PetscStageLogPop()`, `PetscStageLogPush()`, `PetscStageLogRegister()`, `PetscStageLogSetActive()`, `PetscStageLogSetVisible()`, `PetscStageLogGetCurrent()`, and `PetscStageLogGetEventPerfLog()`

- Deprecate `PetscLogGetStageLog()`

- Deprecate `PetscLogPushCurrentEvent_Internal()`

- Deprecate `PetscLogPopCurrentEvent_Internal()`

- Deprecate `PetscLogAllBegin()`

- Deprecate `PetscLogSet()`

- Add `PetscLogHandlerGetEventPerfInfo()`

- Add `PetscLogHandlerSetLogActions()`

- Add `PetscLogHandlerSetLogObjects()`

- Add `PetscLogHandlerLogObjectState()`

- Add `PetscLogHandlerGetNumObjects()`

- Add `PetscLogHandlerEventDeactivatePush()`

- Add `PetscLogHandlerEventDeactivatePop()`

- Add `PetscLogHandlerEventsPause()`

- Add `PetscLogHandlerEventsResume()`

- Add `PetscLogHandlerDump()`

- Add `PetscLogHandlerStageSetVisible()`

- Add `PetscLogHandlerStageGetVisible()`

- Add `PetscLogHandlerGetStagePerfInfo()`

- Add `PetscLogStageGetPerfInfo()`

- Add `PETSCLOGHANDLERNVTX` implementation of `PetscLogHandler`

PetscViewer:

- Add `PetscViewerASCIIOpenWithFileUnit()` and `PetscViewerASCIISetFileUnit()`

PetscDraw:

AO:

IS:

VecScatter / PetscSF:

PF:

Vec:

- Add `VecErrorWeightedNorms()` to unify weighted local truncation error norms used in `TS`

- Add CUDA/HIP implementations for `VecAbs()`, `VecSqrtAbs()`, `VecExp()`, `VecLog()`, `VecPointwiseMax()`, `VecPointwiseMaxAbs()`, and `VecPointwiseMin()`

- Add `VecMAXPBY()`

- Deprecate `VecChop()` in favor of `VecFilter()`

- Add `VecCreateFromOptions()`

PetscSection:

- Add `PetscSectionResetClosurePermutation()`

PetscPartitioner:

Mat:

- Add `MatCreateDenseFromVecType()`

- Add support for calling `MatDuplicate()` on a matrix preallocated via `MatSetPreallocationCOO()`, and then `MatSetValuesCOO()` on the new matrix

- Remove `MATSOLVERSPARSEELEMENTAL` since it is no longer functional

- Add `MATSELLCUDA`. It supports fast `MatMult()`, `MatMultTranspose()`, and `MatMultAdd()` on GPUs

- Add support for `MAT_FACTOR_LU` and `MAT_FACTOR_CHOLESKY` with `MATSOLVERMUMPS` for `MATNEST`

- `MatGetFactor()` can now return `NULL` for some combinations of matrices and solvers types. This is to support those combinations that can only be inspected at runtime (i.e. MatNest with AIJ blocks vs MatNest with SHELL blocks)

- Remove `MatSetValuesDevice()`, `MatCUSPARSEGetDeviceMatWrite()`, `MatKokkosGetDeviceMatWrite`

- Add `MatDenseCUDASetPreallocation()` and `MatDenseHIPSetPreallocation()`

- Add support for KOKKOS in `MATH2OPUS`

- Add `-pc_precision single` option for use with `MATSOLVERSUPERLU_DIST`

- Add `MATDIAGONAL` which can be created with `MatCreateDiagonal()`

- Add `MatDiagonalGetDiagonal()`, `MatDiagonalRestoreDiagonal()`, `MatDiagonalGetInverseDiagonal()`, and `MatDiagonalRestoreInverseDiagonal()`

- Add support for `MatLoad()` and `MatView()` to load and store `MATMPIAIJ` matrices that have more than `PETSC_INT_MAX` nonzeros, so long as each rank has fewer than `PETSC_INT_MAX`

- Add `MatLRCSetMats()` and register creation routine for `MatLRC`

- Add CUDA/HIP implementation for `MatGetDiagonal()`

- Add a Boolean parameter to `MatEliminateZeros()` to force the removal of zero diagonal coefficients

- Expose `MatComputeVariableBlockEnvelope()` in public headers

- Add `MatEliminateZeros()` implementations for `MatBAIJ` and `MatSBAIJ`

- Deprecate `MatChop()` in favor of `MatFilter()`, with two additional Boolean parameters to compress the underlying storage and keep or discard near-zero diagonal coefficients

- Add `MatCreateFromOptions()`

- Add `MatSTRUMPACKGetReordering()`, `MatSTRUMPACKSetGeometricNxyz()`, `MatSTRUMPACKSetGeometricComponents()`, `MatSTRUMPACKSetGeometricWidth()`, `MatSTRUMPACKGetColPerm()`, `MatSTRUMPACKSetGPU()`, `MatSTRUMPACKGetGPU()`, `MatSTRUMPACKSetCompression()`, `MatSTRUMPACKGetCompression()`, `MatSTRUMPACKSetCompRelTol()`, `MatSTRUMPACKGetCompRelTol()`, `MatSTRUMPACKSetCompAbsTol()`, `MatSTRUMPACKGetCompAbsTol()`, `MatSTRUMPACKSetCompMinSepSize()`, `MatSTRUMPACKGetCompMinSepSize()`, `MatSTRUMPACKSetCompLeafSize()`, `MatSTRUMPACKGetCompLeafSize()`, `MatSTRUMPACKSetCompLossyPrecision()`, `MatSTRUMPACKGetCompLossyPrecision()`, `MatSTRUMPACKSetCompButterflyLevels()`, `MatSTRUMPACKGetCompButterflyLevels()`

- Add `MAT_STRUMPACK_GEOMETRIC`, `MAT_STRUMPACK_AMD`, `MAT_STRUMPACK_MMD`, `MAT_STRUMPACK_AND`, `MAT_STRUMPACK_MLF`, `MAT_STRUMPACK_SPECTRAL` to `MatSTRUMPACKReordering`

- Add `MatSTRUMPACKCompressionType`

- Remove `MatSTRUMPACKSetHSSLeafSize()`, `MatSTRUMPACKSetHSSMaxRank()`, `MatSTRUMPACKSetHSSMinSize()`, `MatSTRUMPACKSetHSSMinSepSize()`, `MatSTRUMPACKSetHSSAbsTol()`, `MatSTRUMPACKSetHSSRelCompTol()`, `MatSTRUMPACKSetHSSRelTol()`

MatCoarsen:

PC:

- Add `PCMatGetApplyOperation()` and `PCMatSetApplyOperation()`

- Add `PCReduceFailedReason()`

- Add `PCSetKSPNestLevel()` and `PCSetKSPNestLevel()`

- Refactor `PCMPI` to be a private system used automatically when `-mpi_linear_solver_server` is used. The `KSP` and `PC` objects that solve the system now inherit any prefix provided initially with `KSPSetPrefix()` and do not require the previously required `mpi_` prefix

- Add option `-fieldsplit_1_pc_hpddm_schur_precondition` to use `PCHPDDM` on the Schur complements from `PCFIELDSPLIT`

- Add `PCGAMGSetAggressiveSquareGraph()` to use square graph method for aggressive coarsening that was the previous default coarsening approach before release 3.19

- Add `PCGAMGSetAggressiveMISk()` to set the number of levels (k) of aggressive MIS-k coarsening (2 is very common)

- Add `PCGAMGSetMinDegreeOrderingMISk()` to use a minimum degree ordering for the (greedy) MIS-k algorithm

- Change `PCGAMGSetUseParallelCoarseGridSolve()` to `PCGAMGSetParallelCoarseGridSolve()`

- Add `PCGAMGSetRecomputeEstEig()` to set flag to have Chebyshev recompute its eigen estimates (default set to true)

KSP:

- Add `KSPSetMinimumIterations()` and `KSPGetMinimumIterations()`

- Add `KSPSetNestLevel()` and `KSPGetNestLevel()`

- Support `KSPSetInitialGuessNonzero()` with `KSPPREONLY` and `PCDISTRIBUTE` when it is called on both the outer and inner `KSP`

SNES:

- Add a convenient, developer-level `SNESConverged()` function that runs the convergence test and updates the internal converged reason

- Swap the order of monitor and convergence test. Now monitors are always called after a convergence test

- Deprecate option `-snes_ms_norms` in favor of `-snes_norm_schedule always`

SNESLineSearch:

TS:

- Remove `TSErrorWeightedNormInfinity()`, `TSErrorWeightedNorm2()`, `TSErrorWeightedENormInfinity()`, `TSErrorWeightedENorm2()` since the same functionality can be obtained with `VecErrorWeightedNorms()`

- Add support for time-dependent solvers with varying solution size using `TSSetResize()`

- Add support for Diagonally Implicit Runge-Kutta methods with the new `TSDIRK` type

TAO:

- Add `TaoADMMGetRegularizerCoefficient()`

- Add `TAOBNCG`, `TaoBNCGGetType()` and `TaoBNCGSetType()`

DM/DA:

- Add support for `DMDAGetElements()` for Fortran

- Add support for clearing named vectors with `DMClearNamedGlobalVectors()` and `DMClearNamedLocalVectors()`

- Add `DMPrintCellIndices()` for debugging

- Add `DMUseTensorOrder()`

- Add `DMCeedCreate()`, `DMCeedComputeGeometry()`, and `DMCeedDestroy()`

DMSwarm:

- Add scatter mode to `DMSwarmProjectFields()` and no longer create vectors

- Promote `DMSwarmDataFieldGetEntries()`, `DMSwarmDataFieldRestoreEntries()`, `DMSwarmDataBucketGetDMSwarmDataFieldByName()`, `DMSwarmDataBucketGetDMSwarmDataFieldIdByName()`, and `DMSwarmDataBucketQueryDMSwarmDataFieldByName()` to public interface

DMPlex:

- Add `DMPlexTransformExtrudeGetPeriodic()` and `DMPlexTransformExtrudeSetPeriodic()`

- Replace `DMPlexGetGhostCellStratum()` with `DMPlexGetCellTypeStratum()`

- Add `DMPlexSetUseMatClosurePermutation()` and `DMPlexGetUseMatClosurePermutation()`

- Add `DMPlexGetUseCeed()`, `DMPlexSetUseCeed()`, and `DMPlexSNESComputeResidualCEED()`

- Add `DMPlexSNESComputeResidualDS()`

FE/FV:

- Add `PLEXFE_QFUNCTION` macro to create CEED QFunctions from Plex pointwise functions

DMNetwork:

- Add `DMNetworkViewSetShowRanks()`, `DMNetworkViewSetViewRanks()`, `DMNetworkViewSetShowGlobal()`, `DMNetworkViewSetShowVertices()`, `DMNetworkViewSetShowNumbering()`

- Add `-dmnetwork_view_all_ranks` `-dmnetwork_view_rank_range` `-dmnetwork_view_no_vertices` `-dmnetwork_view_no_numbering` for viewing DMNetworks with the Matplotlib viewer

- Add `-dmnetwork_view_zoomin_vertices` `-dmnetwork_view_zoomin_vertices_padding` for viewing DMNetworks with the Matplotlib viewer

DMStag:

DT:

- Add `PETSCFEVECTOR` for making a finite element that is a vector of copies of another finite element

- Add `PetscFECreateVector()`

- Add `PETSCDUALSPACESUM` for a finite element dual basis that is the sum of other dual bases

- Add `PetscDualSpaceCreateSum()`

- Add `PetscDualSpaceSumSetNumSubspaces()`

- Add `PetscDualSpaceSumGetNumSubspaces()`

- Add `PetscDualSpaceSumSetSubspace()`

- Add `PetscDualSpaceSumGetSubspace()`

- Add `PetscDualSpaceSumSetConcatenate()`

- Add `PetscDualSpaceSumGetConcatenate()`

- Add `PetscDualSpaceSumSetInterleave()`

- Add `PetscDualSpaceSumGetInterleave()`

- Add `PetscSpaceSumSetInterleave()`

- Add `PetscSpaceSumGetInterleave()`

- Add `PetscDualSpaceGetInteriorSection()`

Fortran:

- Add `PetscCheck()` and `PetscCheckA()` for Fortran

- Change `PETSC_HAVE_FORTRAN` to `PETSC_USE_FORTRAN_BINDINGS` to indicate if PETSc is built with Fortran bindings