Plex: Make integration functions public
- Rename
DMPlexComputeResidual_Internal() --> DMPlexComputeResidualForKey()
DMPlexComputeJacobian_Internal() --> DMPlexComputeJacobianByKey()
DMPlexCompu

Plex: Make integration functions public
- Rename
DMPlexComputeResidual_Internal() --> DMPlexComputeResidualForKey()
DMPlexComputeJacobian_Internal() --> DMPlexComputeJacobianByKey()
DMPlexComputeJacobian_Action_Internal() --> DMPlexComputeJacobianActionByKey()
DMPlexComputeResidual_Hybrid_Internal() --> DMPlexComputeResidualHybridByKey()
DMPlexComputeJacobian_Hybrid_Internal() --> DMPlexComputeJacobianHybridByKey()

show more ...

Merge branch 'barry/2025-04-16/docs-converged-normal-equations/release' into 'main'

Improve the clarity of KSP converging for the normal equations.

See merge request petsc/petsc!8327

Improve the clarity of KSP converging for the normal equations.

Reported-by: Mark Adams

Merge branch 'jolivet/cleanup-symbols' into 'main'

Do not expose unneeded symbols

See merge request petsc/petsc!8335

Merge branch 'knepley/feature-plex-transform-cohesive' into 'main'

Plex: Fix cohesive cell integration

See merge request petsc/petsc!8175

Merge branch 'knepley/fix-plex-hdf5-compress' into 'main'

Plex+HDF5: Fix compression handling in HDF5 output

See merge request petsc/petsc!8317

Deprecate PetscSSEIsEnabled()

Plex: Fix hybrid integration
- Add DMFieldCreateDefaultFaceQuadrature() and PetscFEExpandFaceQuadrature()
- Add neighbor geometry to PetscFEIntegrateHybridResidual_Basic(), PetscFEIntegrateHybridJaco

Plex: Fix hybrid integration
- Add DMFieldCreateDefaultFaceQuadrature() and PetscFEExpandFaceQuadrature()
- Add neighbor geometry to PetscFEIntegrateHybridResidual_Basic(), PetscFEIntegrateHybridJacobian_Basic() and PetscFEEvaluateFieldJets_Hybrid_Internal()
- Add cell geometry to PetscFEIntegrateHybridResidual() and PetscFEIntegrateHybridJacobian()

show more ...

Merge branch 'jolivet/feature-pcmatapplytranspose' into 'main'

PC: add PCMatApplyTranspose()

See merge request petsc/petsc!8322

PC: add PCMatApplyTranspose()

Plex+HDF5: Fix compression handling in HDF5 output
- Add compress flag to HDF5 Viewer
- Add PetscViewerHDF5SetCompress() and PetscViewerHDF5GetCompress()
- Initialize IS compressOutput to PETSC_TRUE

Plex+HDF5: Fix compression handling in HDF5 output
- Add compress flag to HDF5 Viewer
- Add PetscViewerHDF5SetCompress() and PetscViewerHDF5GetCompress()
- Initialize IS compressOutput to PETSC_TRUE and check viewer flag as well
- Turn on viewer compression rather than flip IS flag in Plex output

show more ...

Merge branch 'hsuh/taobrgn-type-viewer-fix' into 'main'

TaoBRGNSetRegularizationType

See merge request petsc/petsc!8297

Plex: Add DMPlexTransformCohesiveExtrudeGetUnsplit()

Remove PC_ApplyMultiple

TAOBRGN structure improvement, add TaoBRGNSet,GetRegularizationType

Merge branch 'rmills/petscregressor-update-changelog' into 'main'

Document addition of PetscRegressor and PETSCREGRESSORLINEAR in dev.md

See merge request petsc/petsc!8307

Document addition of PetscRegressor and PETSCREGRESSORLINEAR in dev.md

Merge branch 'tisaac/feature-lmvm-rework' into 'main'

MatLMVM rework

See merge request petsc/petsc!6401

MatLMVM: Complete rework of internals to use LMBasis and LMProducts

I don't think there is a good way to break this up into smaller commits: once I
use LMBasis and LMProducts in Mat_LMVM, all of the

MatLMVM: Complete rework of internals to use LMBasis and LMProducts

I don't think there is a good way to break this up into smaller commits: once I
use LMBasis and LMProducts in Mat_LMVM, all of the implementations that depend on
it are broken until they are changed to also use these new objects. It is probably
better to review all of the MatLMVM implementations (except for denseqn) as new code
rather than trying to review using a diff tool.

This rework achieves:

- Implementations of three classes of algorithms for limitd memory quasi-Newton updates
- MAT_LMVM_MULT_RECURSIVE: Recursive, level 1 BLAS algorithms that closely match
the typical way the methods are presented mathematically.
- MAT_LMVM_MULT_DENSE: level 2 BLAS algorithms that, when possible,
avoid any need for recomputation when the reference Jacobian J0 changes.
- MAT_LMVM_MULT_COMPACT_DENSE: level 2 BLAS algorithms that construct rank-m (or rank-2m)
update representations of the operators, for the best-performance on MatMult() and MatSolve(),
but at the cost of additional setup, and requiring more recomputation when J0 changes
than MAT_LMVM_MULT_DENSE
- Every quasi-Newton method has a dual method (e.g. BFGS and DFP): we exploit this duality
to use only one code path for a dual pair, reducing the numer of implementaitons to maintain
- Special handling of the case J0 = gamma * I, avoiding products that are unnecessary in this case
- Instead of including a MatLMVMDiagBrdn as the J0 matrix when J0 is
variable (which means that both B and J0 of basis vectors to maintain),
a SymBroydenRescale object is shared by MatLMVMDiagBrdn and other formats that use rescaling
- Improved correctness of when internal products need to be recomputed: if the user externally
changes J0, this will be detected

show more ...

Mat: Improve the usability of constantdiagonal and diagonal matrices

Mat_ConstantDiagonal: Add MatConstantDiagonalGetConstant() and implement conjugate and transpose

Mat_Diagonal: Add missing setra

Mat: Improve the usability of constantdiagonal and diagonal matrices

Mat_ConstantDiagonal: Add MatConstantDiagonalGetConstant() and implement conjugate and transpose

Mat_Diagonal: Add missing setrandom, conjugate, and transpose methods

MatView: Don't print MatInfo for constantdiagonal or diagonal matrices

show more ...

Merge branch 'barry/2025-04-08/mv-changed-dev-rst-md' into 'main'

Move changes/dev.rst to changes/dev.md

See merge request petsc/petsc!8294

Move changes/dev.rst to changes/dev.md