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/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 'stefanozampini/submat-zerorows' into 'main'

implement ZeroRowsLocal and ZeroRowsColumnsLocal for local references

See merge request petsc/petsc!8705

Mat: add ops for zerorowscolumnslocal

Merge branch 'jolivet/clang-format-21' into 'main'

clang-format-21

See merge request petsc/petsc!8202

Remove unnecessary braces around one-liners

git grep -lE "[ ]*(if|for|while) \(.*\) {[^;]*;[^;]*}$" -- '*.c' '*.cxx' '*.cu' '*.h' '*.hpp' '*.cpp' | xargs sed -i '' -E 's#([ ]*)(if|for|while) \((.*)\

Remove unnecessary braces around one-liners

git grep -lE "[ ]*(if|for|while) \(.*\) {[^;]*;[^;]*}$" -- '*.c' '*.cxx' '*.cu' '*.h' '*.hpp' '*.cpp' | xargs sed -i '' -E 's#([ ]*)(if|for|while) \((.*)\) {([^;]*);([^;]*)}$#\1\2 \(\3\)\4;\5#'

show more ...

One-liners from petsc/petsc!5344 and petsc/petsc!5557

Slightly reworked regular expression

git ls-files -z -- '*.c' '*.cxx' '*.cu' '*.h' '*.hpp' '*.cpp' | while IFS= read -r -d '' file; do
cat

One-liners from petsc/petsc!5344 and petsc/petsc!5557

Slightly reworked regular expression

git ls-files -z -- '*.c' '*.cxx' '*.cu' '*.h' '*.hpp' '*.cpp' | while IFS= read -r -d '' file; do
cat $file | tr '\n' '\r' | sed -E 's/\r([ ]*)(for|if|while|else) ([^\r]*)\{\r[ ]*Petsc([a-zA-Z]*)\(([^\r]*)\);\r[ ]*\}\r/\r\1\2 \3Petsc\4(\5);\r/g' | tr '\r' '\n' > ${file}.joe; mv ${file}.joe ${file}
done

show more ...

Merge branch 'MarDiehl/fix-typo' into 'main'

missing "s" for isascii and issundials

See merge request petsc/petsc!8559

missing "s" for isascii and issundials

Merge branch 'lindad/hypre-control-through-pmat' into 'main'

Change hypre PC option defaults based on hypre matrix memory location

See merge request petsc/petsc!8505

Make MatGetCurrentMemType first class API

Merge branch 'jolivet/cleanup' into 'main'

Housekeeping

See merge request petsc/petsc!8392

MatOps cleanup

Merge branch 'stefanozampini/matoption-no-error' into 'main'

MatSetOption: ignore unused options instead of erroring out

See merge request petsc/petsc!8179

MatSetOption: ignore unused options instead of erroring out

Merge branch 'lindad/new-aij-from-hash' into 'main'

Add MatCopyHashToXAIJ() to copy hash table entries into an XAIJ matrix type

See merge request petsc/petsc!7939

Add MatCopyHashToXAIJ() to copy hash table entries into an XAIJ matrix type

Merge branch 'stefanozampini/useless-cast' into 'main'

Remove useless cast

See merge request petsc/petsc!7894

Brain dead fixes for useless casts