Expose MatSeqBAIJ[Get|Restore]Array

Replace PetscInt with PetscMPIInt when returning a rank

Improves MatConvert([MPI|Seq][S]BAIJ)

MatMatMult_SeqSBAIJ_SeqDense

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

Double semicolons

Merge branch 'wence/fix/freeze-petsc-layout' into 'master'

freeze PetscLayout after PetscLayoutSetUp

Closes #356

See merge request petsc/petsc!2028

fixes for c89

PetscLayout: forbid calling setup twice with different sizes

* PetscLayout objects should be considered immutable, especially because they can be shared.
* Make use of PetscLayoutCreateFromSizes().

PetscLayout: forbid calling setup twice with different sizes

* PetscLayout objects should be considered immutable, especially because they can be shared.
* Make use of PetscLayoutCreateFromSizes().
* setupcalled is now the single indicator that PetscLayout is ready.
* Insure the rest of variables (e.g. map->rstart) are calculated even if n, N and range are set.

show more ...

Same logic between MPIAIJ and MPI[SBAIJ,BAIJ].

MatMPISBAIJSetPreallocationCSR: Fix for roworiented input

add fast path insertion for bs==1 for the other block classes

Merge branch 'barry/2019-08-01/fix-matmpisbaijsetpreallocationcsr' [PR #1941]

* barry/2019-08-01/fix-matmpisbaijsetpreallocationcsr:
Fix bug in MatMPISBAIJSetPreallocationCSR()
Add test case that de

Merge branch 'barry/2019-08-01/fix-matmpisbaijsetpreallocationcsr' [PR #1941]

* barry/2019-08-01/fix-matmpisbaijsetpreallocationcsr:
Fix bug in MatMPISBAIJSetPreallocationCSR()
Add test case that detected bug
Fix over-preallocation that occured in MatMPISBAIJSetPreallocationCSR() and MatSeqSBAIJSetPreallocationCSR()
Added documentation MATMPI[B]AIJ and MATSEQ[B]AIJ noting their support of two features

show more ...

Merge branch 'jolivet/matconcatenate-matsbaij-for-real' [PR #1934]

* jolivet/matconcatenate-matsbaij-for-real:
Fixes MatCreateMPIMatConcatenateSeqMat_MPISBAIJ for good.

Missing needed <immintrin.h> file

Commit-type: bug-fix, portability-fix
Reported-by: Lisandro Dalcin <dalcinl@gmail.com>

Comments from Barry.

Fix bug in MatMPISBAIJSetPreallocationCSR()

Add test case that detected bug

Fix over-preallocation that occured in MatMPISBAIJSetPreallocationCSR() and MatSeqSBAIJSetPreallocationCSR()

Added docum

Fix bug in MatMPISBAIJSetPreallocationCSR()

Add test case that detected bug

Fix over-preallocation that occured in MatMPISBAIJSetPreallocationCSR() and MatSeqSBAIJSetPreallocationCSR()

Added documentation MATMPI[B]AIJ and MATSEQ[B]AIJ noting their support of two features

Commit-type: optimization, bug-fix, testing, documentation, example
Reported-by: Pierre Jolivet <pierre.jolivet@enseeiht.fr>

show more ...

Fix manual pages for MatXXXYBAIJSetPreallocationCSR() routines

Commit-type: documentation

Add PetscIntCast() to check for integer overflow

Used in MatSeqXAIJSetPreallocation() to determine if the number of nonzeros in the matrix
(and hence the ai[] column indices and lengths of a

Add PetscIntCast() to check for integer overflow

Used in MatSeqXAIJSetPreallocation() to determine if the number of nonzeros in the matrix
(and hence the ai[] column indices and lengths of aj[] columns indices and a[] numerical values) fit
in PetscInt.

Commit-type: error-checking
Reported-by: Karl Lin <karl.linkui@gmail.com>

show more ...

Merge branch 'barry/2019-07-19/PetscKernel-macros-to-functions' [PR #1894]

* barry/2019-07-19/PetscKernel-macros-to-functions:
Replace macros with inline static functions for BAIJ operations

Remove stray charactor > from source file

Commit-type: bug-fix
Reported-by: nightly tests

Replace macros with inline static functions for BAIJ operations

This removes warning about unused error checking, warning appeared in commit 580bdb3

Commit-type: portability-fix
Reported-by: Jed Br

Replace macros with inline static functions for BAIJ operations

This removes warning about unused error checking, warning appeared in commit 580bdb3

Commit-type: portability-fix
Reported-by: Jed Brown <jed@jedbrown.org>

show more ...

Replaced use of private PETSc kernels with MatInvertBlockDiagonal().

Mat: expose PetscLINPACKgefa and gedi --with-visibility and separate libraries

Merge branch 'dalcinl/fix-headers' [PR #1897]

* dalcinl/fix-headers:
Minor fixes to headers.

MATSEQBAIJ: Add missing CHKERRQ() to silence -Wunused-value warning

/home/devel/petsc/dev/src/mat/impls/baij/seq/baijfact81.c: In function ‘MatLUFactorNumeric_SeqBAIJ_9_NaturalOrdering’:
/home/devel

MATSEQBAIJ: Add missing CHKERRQ() to silence -Wunused-value warning

/home/devel/petsc/dev/src/mat/impls/baij/seq/baijfact81.c: In function ‘MatLUFactorNumeric_SeqBAIJ_9_NaturalOrdering’:
/home/devel/petsc/dev/include/petscsys.h:1771:77: warning: value computed is not used [-Wunused-value]
1771 | #define PetscArraycpy(str1,str2,cnt) ((sizeof(*(str1)) != sizeof(*(str2))) || PetscMemcpy(str1,str2,(cnt)*sizeof(*(str1))))

show more ...