Merge branch 'jacobf/2021-12-09/attribute-format-mat' into 'main'

Feature: Attribute Format - Mat

See merge request petsc/petsc!4634

applied attribute format fixes to mat headers

Merge branch 'jed/thrust-nosync' into 'main'

thrust: upgrade to 1.16 prerelease for thrust::cuda::par_nosync in vector operations

See merge request petsc/petsc!4676

thrust: upgrade to 1.16 prerelease for async on streams

Headers in this version are less promiscuous so we need more includes.

Merge branch 'stefanozampini/petsc4py' into 'main'

Some work on petsc4py

See merge request petsc/petsc!4644

Merge remote-tracking branch 'origin/release'

Merge branch 'jczhang/revise-a-comment' into 'release'

Revise a comment

See merge request petsc/petsc!4675

Revise a comment

Mat: add missing general fallbacks for matrices flagged as symmetric

Fix MatMatSolve_Basic using MatDense API to prevent from unneeded device2host copies

Merge branch 'adams/landau-cuda-opt' into 'main'

Landau optimizations + GPU assembly bug fix

See merge request petsc/petsc!4358

add hook for CPU cusparse solves

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

Add initial MatProduct support for MATAIJKOKKOS

See merge request petsc/petsc!4303

MATSEQAIJCUSPARSE: implement MatSeqAIJ{Get,Restore}Array{Read,Write}

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

Rmills/feature dm bind below

See merge request petsc/petsc!3455

Prevent unsafe setting of inode.use to true in MatBindToCPU() for cuSPARSE and ViennaCL matrices.

This fixes an issue caused by Stefano's commit c58ef05e, "MatInode: set it to
false when using devic

Prevent unsafe setting of inode.use to true in MatBindToCPU() for cuSPARSE and ViennaCL matrices.

This fixes an issue caused by Stefano's commit c58ef05e, "MatInode: set it to
false when using device matrices", which introduced flipping the inode.use flag
inside MatBindToCPU_SeqAIJCUSPARSE() and MatBindToCPU_SeqAIJViennaCL(). I have
found cases in which this would cause "Corrupted Petsc object, Missing Inode
Structure" errors. Only flipping inode.use to true if inode.size is nonzero
eliminates these errors.

show more ...

MatGetFactor_seqaijcusparse_cusparse() now sets ops to AIJ versions if matrix is bound to CPU.

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

MAT: access spptr after MatSeqAIJCUSPARSECopyToGPU, but not before

MAT: form explicit transpose whenever MatSeqAIJCUSPARSEFormExplicitTransposeForMult is called

Merge branch 'jacobf/2021-09-04/petscdevice-cupminitialization' into 'main'

Replace current CUPM initialization with unified PetscDevice version

See merge request petsc/petsc!4310

- Replace the current CUPM initialization infrastructure with PetscDevice.
- Remove -petsc_use_default_null_stream and -[cuda|hip]_synchronize options
- Remove PetscCUDASynchronize and PetscHIPSynchr

- Replace the current CUPM initialization infrastructure with PetscDevice.
- Remove -petsc_use_default_null_stream and -[cuda|hip]_synchronize options
- Remove PetscCUDASynchronize and PetscHIPSynchronize
- Remove PetscCUDAInitialize(), PetscCUDAInitializeCheck(), PetscHIPInitialize(), and
PetscHIPInitializeCheck()
- Remove PetscCUBLASInitializeHandle(), PetscCUSOLVERDnInitializeHandle(),
PetscHIPBLASInitializeHandle(), and PetscHIPSOLVERInitializeHandle().
- Remove petsc_gputimer_begin and petsc_gputimer_begin
- Add -device_enable, -device_select and -device_view options
- Replace -[cuda|hip]_device with split options -device_enable_[cuda|hip] and -device_select_[cuda|hip]
- Replace -[cuda|hip]_view with -device_view_[cuda|hip]
- Add PetscDeviceInitType
- Add PetscDeviceInitialize() and PetscDeviceInitialized()
- Add PetscDeviceView()

show more ...

Merge branch 'stefanozampini/hypre-gpu' into 'main'

Interface HYPRE GPU solvers

See merge request petsc/petsc!4248

MATSEQAIJCUSPARSE: add missing man pages