Merge branch 'dalcinl/bt-byte' into 'main'

Implement PetscBT with new PetscByte datatype as unsigned char

See merge request petsc/petsc!9046

Merge remote-tracking branch 'origin/release'

Fix MPI_Comm_rank/size API mixups in helpers (#1865,#1866)

Switch some defined(PETSC_ to PetscDefined(

Implement PetscBT with new PetscByte datatype as unsigned char

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

Merge branch 'release'

Minor cleanup: use MPI_IN_PLACE when possible and avoid extra allocation

Remove legacy '+1' allocations in XXXDiagonalScaleLocalSetUp() routines for AIJ and BAIJ parallel matrices

This commit removes '+1's in some memory allocations that are believed
to be leftover from

Remove legacy '+1' allocations in XXXDiagonalScaleLocalSetUp() routines for AIJ and BAIJ parallel matrices

This commit removes '+1's in some memory allocations that are believed
to be leftover from long ago to handle corner cases with malloc() calls
of size 0, which some implementations would not handle correctly. This
ceased to be an issue long ago.

Thanks-to: Pierre Jolivet <pierre@joliv.et>, Barry Smith <bsmith@petsc.dev>

show more ...

Fix long-standing bug in MatMPIAIJDiagonalScaleLocalSetUp()

This commit fixes a problem I discovered in October 2017 and then forgot about
after some unsuccessful initial attempts to find the proble

Fix long-standing bug in MatMPIAIJDiagonalScaleLocalSetUp()

This commit fixes a problem I discovered in October 2017 and then forgot about
after some unsuccessful initial attempts to find the problem. The MPIAIJ
implementation for MatDiagonalScaleLocalSetUp() was broken when the
local-to-global mappings in PETSc were refactored to always use
block indices; the MPIAIJ implementation of
MatDiagonalScaleLocalSetUp(), however, used element-wise, not block
indexing, and this was not refactored at the time, as it was overlooked
due to having no test case.
My commit finally fixes this; it also fixes several memory leaks in the
AIJ and BAIJ implementations of MatDiagonalScaleLocal() which were
overlooked, again, due to having no test cases.

show more ...

Set MAT_NO_OFF_PROC_ENTRIES when appropriate to avoid MPI_Allreduce()

MatView(): do not print block size twice

Rework MAT_SYMMETRIC and MAT_HERMITIAN options

VecSetValues() instead

MatNorm: fix some non-scalable implementations for NORM_1

Reported-by: Vincent Robert <vincent.robert@cea.fr>

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 ...

Merge remote-tracking branch 'origin/release'