Merge branch 'barry/2026-01-13/ftn-generation-for-context-functions' into 'main'

Fortran generation of interface definitions for functions that take and return context variables (void *ctx and PetC

Merge branch 'barry/2026-01-13/ftn-generation-for-context-functions' into 'main'

Fortran generation of interface definitions for functions that take and return context variables (void *ctx and PetCtx ctx)

See merge request petsc/petsc!8956

show more ...

Add automatic generated fortran support for setting and returning contexts

Major cleanup of the code for setting and returning application (previously called user) contexts
- PetscCtx ctx is now t

Add automatic generated fortran support for setting and returning contexts

Major cleanup of the code for setting and returning application (previously called user) contexts
- PetscCtx ctx is now the standard for passing in contexts
- PetscCtxRt ctx is the standard for passing out contexts, including context destroy routines based on PetscCtxDestroyFn

Both of these are typedef to void*

Also removed most use of user and userctx for context arguments now consistently named ctx

Improved documentation on PetscCtx and PetscCtxRt and on how to use contexts (including getting them) from Fortran

show more ...

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

[Mat|Vec]Conjugate(): no need to check for USE_COMPLEX in implementations

See merge request petsc/petsc!8937

MatConjugate(): no need to check for USE_COMPLEX in implementations

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

Minor cleanup: use MPI_IN_PLACE when possible and avoid extra allocation

See merge request petsc/petsc!8925

Minor cleanup: use MPI_IN_PLACE when possible and avoid extra allocation

Merge branch 'jolivet/fix-unnecessary-reduce' into 'main'

Mat: set appropriate NO_OFF_PROC_ENTRIES to bypass MPI_Allreduce()

See merge request petsc/petsc!8872

MatAssemblyEnd(): fix missing space in -info

Merge branch 'jolivet/safeguard-pchpddm-schur-ci' into 'main'

Housekeeping + more safeguards in PCHPDDM

See merge request petsc/petsc!8859

Housekeeping

Merge branch 'barry/2025-10-18/optimize-aij-ops' into 'main'

Refactor handling of diagonal marking in *AIJ and SELL matrices

See merge request petsc/petsc!8797

Optimizations on Seq_[SB]AIJ operations that use MatGetDiagonalMarkers.

- Replace MatMissingDiagonal() and MatMarkDiagonals_SeqXXX() with MatGetDiagonalMarkers_SeqXXX()

- Mat_SeqXXX->diag is not automatically created during MatAssemblyEnd() saving memory and time

- Replace MatMissingDiagonal() and MatMarkDiagonals_SeqXXX() with MatGetDiagonalMarkers_SeqXXX()

- Mat_SeqXXX->diag is not automatically created during MatAssemblyEnd() saving memory and time

- Accessing Mat_SeqXXX->diag now requires the use of MatGetDiagonalMarkers_SeqXXX() except when the current values are known to be correct; for example during numerical factorizations and solves

- Mat_SeqXXX->diag is now never shared among matrices; hence the free_diag flag is gone. That was always a risky proposition since any of the owning matrices could chang the values thus making them incorrect for other owners.

show more ...

This MR optimize some MAT[SB]AIJ operations related to the diagonal

The longer term goal is to refactor some Mat_SeqXXX non-numeric code to make it easier to also implement it on GPUs. For example,

This MR optimize some MAT[SB]AIJ operations related to the diagonal

The longer term goal is to refactor some Mat_SeqXXX non-numeric code to make it easier to also implement it on GPUs. For example, MatGetDiagonal() and friends

This commit
- MatMarkDiagonal_SeqAIJ/ELL()
- now use change in mat->nonzerostate to determine if it needs to recheck the diagonal locations
- sets diagDense flag for a complete diagonal
- MatInvertDiagonal_SeqAIJ/ELL()
- now uses change in mat->state to determine if diagonal entries and inverses need to be recomputed
- name changed to MatInvertDiagonalForSOR_SeqAIJ/ELL() for code clarity
- Added MatInvertDiagonal_SeqAIJ_Inode() for use by MatSOR_SeqAIJ_Inode()
- now uses change in mat->state to determine if block diagonal entries and inverses need to be recomputed

show more ...

Merge branch 'barry/2025-09-30/add-matgetdiagonal-cuda' into 'main'

Implements MatGetDiagonal for CUDA to take place on device.

See merge request petsc/petsc!8756

Implements MatGetDiagonal for CUDA to take place on device.

Reported-by: Pierre.LEDAC@cea.fr
Thanks-to: Junchao Zhang

Merge remote-tracking branch 'origin/release'

Merge branch 'jolivet/fix-arrays' into 'release'

Fix error in MatTransposeMatMultSymbolic() in case of empty rows

See merge request petsc/petsc!8773

Merge branch 'jczhang/2025-09-19/precision-mumps' into 'main'

Support different precisions with MUMPS

See merge request petsc/petsc!8715

Enable ScaLAPACK and MUMPS for other precisions

Merge remote-tracking branch 'origin/release'

Avoid extra allocation

Merge branch 'barry/2025-10-07/fix-matsetpreallocationcoo-mpiaij-bad-error-check/release' into 'release'

Fix illegal access of i1[n1 - 1] when n1 is zero.

See merge request petsc/petsc!8768

Fix illegal access of i1[n1 - 1] when n1 is zero.

Co-authored-by: Junchao Zhang <jczhang@mcs.anl.gov>

Merge branch 'stefanozampini/submat-zerorows' into 'main'

implement ZeroRowsLocal and ZeroRowsColumnsLocal for local references

See merge request petsc/petsc!8705