bugfix for MatGetRowMaxAbs_MPIBAIJ(); add more tests for BAIJ and SBAIJ

Merge branch 'dalcinl/fix-null' into 'master'

Fix 0 -> NULL for pointers

See merge request petsc/petsc!3086

Mat: Replace 0 -> NULL for pointers

Merge remote-tracking branch 'origin/maint'

Merge branch 'barry/2020-07-19/bugs-flopcount-view/maint' into 'maint'

Limit chance of integer overflow in calls to PetscLogFlops()

See merge request petsc/petsc!2986

Limit chance of integer overflow in calls to PetscLogFlops()

Don't view inner KSP when it has not yet been created.

Commit-type: bug-fix
Funded-by: ECP
/spend 20m
Reported-by: Karl Yang <y.juntao@

Limit chance of integer overflow in calls to PetscLogFlops()

Don't view inner KSP when it has not yet been created.

Commit-type: bug-fix
Funded-by: ECP
/spend 20m
Reported-by: Karl Yang <y.juntao@hotmail.com>

show more ...

Merge branch 'jolivet/feature-kspmatsolve-deprecate-ksphpddmmatsolve' into 'master'

Change how solves with multiple RHS work

See merge request petsc/petsc!2811

Change how solves with multiple RHS work

* Deprecate KSPHPDDMMatSolve
* Add KSPMatSolve, previously known as KSPHPDDMMatSolve
* New option -ksp_matsolve_block_size
* New routines MatDense[Get|Restor

Change how solves with multiple RHS work

* Deprecate KSPHPDDMMatSolve
* Add KSPMatSolve, previously known as KSPHPDDMMatSolve
* New option -ksp_matsolve_block_size
* New routines MatDense[Get|Restore]SubMatrix

show more ...

Merge branch 'stefanozampini/feature-add-hpackages-rebased' into 'master'

Fixes to MatProduct + many bug fixes + new hierarchical matrix support

See merge request petsc/petsc!2800

MatProduct: proper cleanup of additional data structure

Do not mess with view/destroy routine, add an extra slot to MatProduct struct to stash the data, destroy/view routines
Deprecate MatFreeInterm

MatProduct: proper cleanup of additional data structure

Do not mess with view/destroy routine, add an extra slot to MatProduct struct to stash the data, destroy/view routines
Deprecate MatFreeIntermediateStructures in favor of MatProductClear
Clean up base headers for matrix classes
Update tests (remove some todos from ex70 because now the tests work properly)
MatProductSetType should not set the setfromoptions operation!
MatProductClear should reset operations
Remove no longer needed function composition
Add few more PetscInfo messages
Single dispatch routine
MatProductSetFromOptions routines are now allowed to fail
MatProductSetFromOptions_SeqXBAIJ_SeqDense: support for MATPRODUCT_AtB
MatTransposeMatMultNumeric_SeqAIJ_SeqDense: fix bug with leading dimensions
Make internal functions PETSC_INTERN
Fix a couple of places where the destroy routine of the matrix was changed
Products should all use internal data (or composed data attached to the internal data)
Fix few bugs in MPIAIJ_MPIDENSE product with stages

show more ...

MATSEQSBAIJ: fixes for missing row

there could be other places to fix

MatMatMultNumeric_SeqSBAIJ_SeqDense: fix bug

Merge remote-tracking branch 'origin/maint'

Merge branch 'stefanozampini/fix-matmatmult-prj' into 'maint'

Fix MatMatMult old API when reusing dense output matrix

See merge request petsc/petsc!2761

MatMatMult support for SBAIJ HERMITIAN case

MATSEQSBAIJ: fix bug with empty rows

Merge remote-tracking branch 'origin/hzhng/matproduct-api'

This branch reorganizes the API of mat-mat operations: A*B, At*B, A*Bt, PtAP, RARt, A*B*C
- See matproduct.c for detailed description of new API
- removed API function XXXSymbolic() and XXXNumeric()

This branch reorganizes the API of mat-mat operations: A*B, At*B, A*Bt, PtAP, RARt, A*B*C
- See matproduct.c for detailed description of new API
- removed API function XXXSymbolic() and XXXNumeric()
- removed subroutines MatMatMult_xxx, MatPtAP_xxx, ...
- modified MatxxxSymbolic_XXX() prototype (product is an input now)
- selection of algorithm is in MatProductSetFromOptions_xxx()

show more ...

Merge branch 'stefanozampini/update-suitesparse' into 'master'

Update SuiteSparse to use github releases + fixes to SBAIJ hermitian complex case

See merge request petsc/petsc!2305

MATSBAIJ: Fix Hermitian MatMult for Seq and MPI code paths

we can reuse the same MatMult operation with no impact for non-complex builds
Remove improper check for the hermitian complex case from Mat

MATSBAIJ: Fix Hermitian MatMult for Seq and MPI code paths

we can reuse the same MatMult operation with no impact for non-complex builds
Remove improper check for the hermitian complex case from MatGetFactor and move it where it actually belongs
Fix conversion routines from SBAIJ when the hermitian flag is set

show more ...

Merge branch 'jolivet/fix-upstream-hpddm' into 'master'

Multiple fixes to PCHPDDM + MatMatMult_SeqSBAIJ_SeqDense

See merge request petsc/petsc!2123

MatMatMult_SeqSBAIJ_SeqDense

Merge remote-tracking branch 'origin/maint'

Merge branch 'barry/2019-10-15/bug-gamg-complexity/maint' into 'maint'

Bug fix: computing MG complexity had integer overflow. Would only occur for very large problems

See merge request petsc/petsc!

Merge branch 'barry/2019-10-15/bug-gamg-complexity/maint' into 'maint'

Bug fix: computing MG complexity had integer overflow. Would only occur for very large problems

See merge request petsc/petsc!2179

show more ...

Bug fix: computing MG complexity had integer overflow. Would only occur for very large problems

Also fixed code that computes parallel matrix storage information to communicate with PetscLogDouble i

Bug fix: computing MG complexity had integer overflow. Would only occur for very large problems

Also fixed code that computes parallel matrix storage information to communicate with PetscLogDouble instead of PetscReal.

Commit-type: bug-fix
Reported-by: Mark Lohry <mlohry@gmail.com>

show more ...