Merge remote-tracking branch 'origin/release'

Switch some defined(PETSC_ to PetscDefined(

Merge branch 'jose/fix-hipsparse-csrmv-trans' into 'release'

hipSPARSE: workaround for issue with transpose matmult in ROCm 7.2.0

See merge request petsc/petsc!9056

HIPSPARSE: workaround for issue with transpose matmult in ROCm 7.2.0

When built with ROCm 7.2.0 MatMultTranspose() fails with:

[0]PETSC ERROR: --------------------- Error Message ------------------

HIPSPARSE: workaround for issue with transpose matmult in ROCm 7.2.0

When built with ROCm 7.2.0 MatMultTranspose() fails with:

[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
[0]PETSC ERROR: GPU error
[0]PETSC ERROR: hipSPARSE errorcode 3 (HIPSPARSE_STATUS_INVALID_VALUE)
[0]PETSC ERROR: See https://petsc.org/release/faq/ for trouble shooting.
[0]PETSC ERROR: PETSc Development Git Revision: v3.24.4-488-g2c105492db5 Git Date: 2026-02-19 16:00:40 +0000
[0]PETSC ERROR: ../test1 with 1 MPI process(es) and PETSC_ARCH arch-gnu-c-opt-hip Fri Feb 20 16:47:04 2026
[0]PETSC ERROR: Configure options: --with-hip --with-debugging=0 COPTFLAGS=-O FOPTFLAGS=-O CXXOPTFLAGS=-O HIPOPTFLAGS=-O
[0]PETSC ERROR: #1 MatMultAddKernel_SeqAIJHIPSPARSE() at petsc/src/mat/impls/aij/seq/seqhipsparse/aijhipsparse.hip.cxx:3116
[0]PETSC ERROR: #2 MatMultTranspose_SeqAIJHIPSPARSE() at petsc/src/mat/impls/aij/seq/seqhipsparse/aijhipsparse.hip.cxx:3010
[0]PETSC ERROR: #3 MatMultTranspose() at petsc/src/mat/interface/matrix.c:2741
[0]PETSC ERROR: #4 MatMult_Transpose() at petsc/src/mat/impls/transpose/transm.c:9

The issue can be attributed to ROCm 7.2.0 using four minimal examples
for CSR mat-vec: level2 and generic interfaces, both transpose and
non-transpose. The generic-transpose version does not work, it fails
with code: 10, reason: function is not implemented. The workaround in
PETSc is to fall back to the level2 interface if ROCm version is 7.2.0.

The second change is to avoid a compile warning.

show more ...

Fix compile warning with recent version of ROCm

src/mat/impls/aij/seq/seqhipsparse/aijhipsparse.hip.cpp:4199:27: warning: 'identity<int>' is deprecated: use ::internal::identity instead [-Wdeprecate

Fix compile warning with recent version of ROCm

src/mat/impls/aij/seq/seqhipsparse/aijhipsparse.hip.cpp:4199:27: warning: 'identity<int>' is deprecated: use ::internal::identity instead [-Wdeprecated-declarations]
4199 | auto pred = thrust::identity<int>();
| ^
/opt/rocm-7.2.0/include/thrust/functional.h:1182:8: note: 'identity<int>' has been explicitly marked deprecated here
1182 | struct THRUST_DEPRECATED_BECAUSE("use ::internal::identity instead") identity
| ^
/opt/rocm-7.2.0/include/thrust/detail/config/deprecated.h:41:44: note: expanded from macro 'THRUST_DEPRECATED_BECAUSE'
41 | # define THRUST_DEPRECATED_BECAUSE(MSG) [[deprecated(MSG)]]
| ^

show more ...

Change *.cpp to *.cxx

Remove the last of the use of horizontal lines --- as seperators. With modern code development systems they are not needed and just an eye-sore

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

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

Minor cleanup: use MPI_IN_PLACE when possible and avoid extra allocation

Kokkos: update to v5.0.0

CI: remove kokkos from linux-cuda-pkgs (in favor of linux-kokkos-cuda) as v5 breaks with many external pkgs,
CI: continue to use kokkos v4.7.01 for hip builds - as they curr

Kokkos: update to v5.0.0

CI: remove kokkos from linux-cuda-pkgs (in favor of linux-kokkos-cuda) as v5 breaks with many external pkgs,
CI: continue to use kokkos v4.7.01 for hip builds - as they currently use rocm-5.4.3

Co-authored-by: Satish Balay <balay@mcs.anl.gov>

show more ...

Merge remote-tracking branch 'origin/release'

Fix builds against CUDA 13.1

- Fix usage of deprecated thrust:: APIs in CCCL 3.1
- Explicitly include "thrust/copy.h" to fix builds w/ CUDA 13.1
- configure: use SuperLU_DIST 9.2.1 for CUDA 13.1
- C

Fix builds against CUDA 13.1

- Fix usage of deprecated thrust:: APIs in CCCL 3.1
- Explicitly include "thrust/copy.h" to fix builds w/ CUDA 13.1
- configure: use SuperLU_DIST 9.2.1 for CUDA 13.1
- CI: update linux-cuda-latest to test with CUDA 13.1

show more ...

Rework MAT_SYMMETRIC and MAT_HERMITIAN options

Merge branch 'release'

MatAssemblyEnd(): fix missing space in -info

Fix typos

Merge branch 'jose/warning-hip' into 'main'

Fix HIP compile warning

See merge request petsc/petsc!8861

Fix HIP compile warning

src/mat/impls/aij/seq/seqhipsparse/aijhipsparse.hip.cpp:1401:15: error: unused variable 'i' [-Werror,-Wunused-variable]
1401 | PetscInt i;
| ^
src/m

Fix HIP compile warning

src/mat/impls/aij/seq/seqhipsparse/aijhipsparse.hip.cpp:1401:15: error: unused variable 'i' [-Werror,-Wunused-variable]
1401 | PetscInt i;
| ^
src/mat/impls/aij/seq/seqhipsparse/aijhipsparse.hip.cpp:1647:15: error: unused variable 'i' [-Werror,-Wunused-variable]
1647 | PetscInt i;
| ^

show more ...

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

MatSeqBAIJ: add MatProduct_AtB implementation (B of type MatSeqDense)

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