Merge branch 'jczhang/feature-kokkos-spgemm' into 'main'

Add initial MatProduct support for MATAIJKOKKOS

See merge request petsc/petsc!4303

MAT: use MatSeqAIJ{Get,Restore}Array{Read,Write} to hide explicit use of offloadmask

Merge branch 'rmills/feature-dm-bind-below' into 'main'

Rmills/feature dm bind below

See merge request petsc/petsc!3455

MatBindToCPU_MPIAIJ() now also binds the 'lvec' and 'diag' vectors used in matrix-vector multiplication.

Merge remote-tracking branch 'origin/release'

Merge branch 'jczhang/fix-PCTelescope-GPU' into 'release'

Fix PCTelescope GPU matrices

See merge request petsc/petsc!4471

Merge remote-tracking branch 'origin/release'

Merge branch 'jczhang/fix-mat-gpu-bugs' into 'release'

Fix several matrix GPU related bugs

See merge request petsc/petsc!4302

MatProduct: need to check symmetry between symbolic and numeric

Suppose we have code:
MatTransposeMatMult(A,B,MAT_INITIAL_MATRIX,fill,C);
MatTransposeMatMult(E,B,MAT_REUSE_MATRIX,fill,C);

If A

MatProduct: need to check symmetry between symbolic and numeric

Suppose we have code:
MatTransposeMatMult(A,B,MAT_INITIAL_MATRIX,fill,C);
MatTransposeMatMult(E,B,MAT_REUSE_MATRIX,fill,C);

If A is symmetric and E is not, but E has the same nonzero pattern as A, then above code is legitimate.
If we change the MatProductType in MAT_INITIAL_MATRIX from MATPRODUCT_AtB to MATPRODUCT_AB, then
in MAT_REUSE_MATRIX, we need to redo symbolic otherwise C=E^B would be wronlgy computed as C=EB.

show more ...

MPIAIJ: use type of the input mat to determine type of the output mat in MatCreateMPIMatConcatenateSeqMat_MPIAIJ

Merge branch 'jczhang/vecscatter-view-as-mat' into 'main'

Add an option to view the Mvctx of a sparse matrix in matlab, draw or binary

See merge request petsc/petsc!4427

Mat: add code to view Mvctx

Merge branch 'psanan-release-patch-82973' into 'main'

Developer docs, style: remove rule about using "if (!rank)" and replace all "!rank" usage with "rank == 0"

See merge request petsc/petsc!4146

Style: replace "!rank" with "rank == 0"

In SF tutorial ex1 (src/vec/is/sf/tutorials/ex1.c), add parentheses around "rank == 0". This is relevant because ! binds more tightly than +, but == binds les

Style: replace "!rank" with "rank == 0"

In SF tutorial ex1 (src/vec/is/sf/tutorials/ex1.c), add parentheses around "rank == 0". This is relevant because ! binds more tightly than +, but == binds less
tightly, so the result of the computation would otherwise change.

show more ...

Merge remote-tracking branch 'origin/release'

Merge branch 'jczhang/fix-MatCreateMPIAIJWithSplitArrays' into 'release'

Revise code and doc of MatCreateMPIAIJWithSplitArrays

See merge request petsc/petsc!4279

Merge remote-tracking branch 'origin/jose/fix-manpages'

Fix manpages to avoid errors reported by doctext

See merge request petsc/petsc!4309

Merge branch 'rmills/feature-matgetcolumnreductions' into 'main'

Add MatGetColumnReductions()

See merge request petsc/petsc!4272

Add MatGetColumnSums() and MatGetColumnMeans() (and real and imaginary part variants)

MatGetColumnReductions() now takes a PetscInt so ReductionType or NormType members can be passed.
ReductionType

Add MatGetColumnSums() and MatGetColumnMeans() (and real and imaginary part variants)

MatGetColumnReductions() now takes a PetscInt so ReductionType or NormType members can be passed.
ReductionType no longer duplicates entries of NormType.
Replace MatGetColumnNorms_SeqBAIJ() with MatGetColumnReductions_SeqBAIJ()

show more ...

Fix manpages: mismatching number of parameters

Fix manpages: Input/Output Parameters --> Parameter

Fix manpages: Input/Output Parameter --> Parameters

Added MatGetColumnReductions(). MatGetColumnNorms() is now a wrapper that calls MatGetColumnReductions().

Made changes requested by Richard.

Merge branch 'jose/spelling-errors' into 'main'

Fix many spelling errors in manpages and comments

See merge request petsc/petsc!4283

Fix spelling errors in manpages and comments