Merge branch 'jed/header-symbol-consistency' into 'release'

Remove interface declarations for functions that don't exist

See merge request petsc/petsc!3374

Remove interface declarations for functions that don't exist

MATDAAD: remove vestigial references to deleted DAAD

MATCUSPARSE: Implement fast assembly from COO data

Merge branch 'adams/feature-mat-cuda' into 'master'

Adams/feature mat cuda

See merge request petsc/petsc!3137

add cuda matrix for meta data method, assembly

Merge branch 'jczhang/mat-vec-kokkos' into 'master'

Add vec Kokkos and mat Kokkos

See merge request petsc/petsc!3219

Add MATAIJKOKKOS

Merge branch 'barry/09-08-2020/fix-debug-mpiabort-on-other-ranks' into 'master'

Barry/09 08 2020/fix debug mpiabort on other ranks

See merge request petsc/petsc!3179

Fix error with MatAXPY_Seq/MPIAIJ that lost the layout and MatStencilInfo information from the input matrix. The code has not been touched
for five years so it is odd the bug was never exposed before

Fix error with MatAXPY_Seq/MPIAIJ that lost the layout and MatStencilInfo information from the input matrix. The code has not been touched
for five years so it is odd the bug was never exposed before.

Add a precheck in MatAXPY_SeqAIJ() for identical nonzero pattern to avoid the constant recomputation of the resulting matrix.

Commit-type: optimization, bug-fix
/spend 3h
Reported-by: Ed Bueler

show more ...

Merge branch 'barry/2020-07-12/snes-lag-jacobian-rosw' into 'master'

Barry/2020 07 12/snes lag jacobian rosw

See merge request petsc/petsc!3049

SeqDense treats ILU as LU allowing PC defaults to work when dense matrix
is not yet known

DMSHELL uses dense format if no format provided.

Fix manual pages for MatFactor

Commit-type: feature, test

SeqDense treats ILU as LU allowing PC defaults to work when dense matrix
is not yet known

DMSHELL uses dense format if no format provided.

Fix manual pages for MatFactor

Commit-type: feature, testing-fix, documentation
/spend 1h

show more ...

Merge branch 'barry/2020-08-13/matsetvalues-nopreallocation-seqaij' into 'master'

Added MatSeqAIJSetTotalPreallocation() allows efficiently building a SeqAIJ...

See merge request petsc/petsc!3065

Added MatSeqAIJSetTotalPreallocation() allows efficiently building a SeqAIJ matrix row by row without
preallocating each row instead just provide an upper bound on the entire memory needed.

Commit-t

Added MatSeqAIJSetTotalPreallocation() allows efficiently building a SeqAIJ matrix row by row without
preallocating each row instead just provide an upper bound on the entire memory needed.

Commit-type: optimization
/spend 4h
Reported-by: Sajid Ali <sajidsyed2021@u.northwestern.edu>

show more ...

Merge branch 'barry/2020-08-16/fix-source-formatting' into 'master'

Turn on checkbadSource test to generate an error when found; fix all source code that causes errors

See merge request petsc/petsc

Merge branch 'barry/2020-08-16/fix-source-formatting' into 'master'

Turn on checkbadSource test to generate an error when found; fix all source code that causes errors

See merge request petsc/petsc!3064

show more ...

Turn on checkbadSource test to generate an error when found; fix all source code that causes errors

Commit-type: portability-fix, testing-fix, style-fix, feature, maintainability
/spend 1.5h

Merge branch 'barry/2020-08-19/better-pcfailedreason-handling' into 'master'

Slightly improved the handling of PCFailedReason on multiple ranks

See merge request petsc/petsc!3084

Slightly improved the handling of PCFailedReason on multiple ranks

pc->failedreasonrank now contains the failed reason from the current rank
pc->failedreason continues the unified common value colle

Slightly improved the handling of PCFailedReason on multiple ranks

pc->failedreasonrank now contains the failed reason from the current rank
pc->failedreason continues the unified common value collected in KSPCheckNorm/Dot
This parallel code can check the common value and be sure to do the same action on each process
Also printed message now provide the accurate value instead of the value on rank 0

Changed KSPReasonView() and KSPReasonViewFromOptions() to KSPCOnvergeReasonView() and KSPConvergedReasonViewFromOptions()

Added the PetscViewerFormat argument to the public KSPReasonViewFromOptions() and removed the private function

Removed some off the PCFailedReaon processing code from KSPSolve_PREONLY since buggily assumed the same reason all ranks

The error message on KSP failure now also prints the PCFailedReason if appropriate

Needed to add a new alt file

Reported-by: Yu, Xiangyu <Xiangyu.Yu@bp.com>

show more ...

Merge branch 'barry/2020-06-23/pcfactor-no-ordering-external' into 'master'

Most external factorization packages do not use the ordering provided by...

See merge request petsc/petsc!2913

Most external factorization packages do not use the ordering provided by PETSc, therefor only provide it when needed.

Adds MatOrderingType external to indicate not to generate an ordering and use wh

Most external factorization packages do not use the ordering provided by PETSc, therefor only provide it when needed.

Adds MatOrderingType external to indicate not to generate an ordering and use what the package needs
Now -pc_view will not print the wrong PETSc ordering when the ordering is done externally

Saves some memory and compute time.

Commit-type: optimization
Time: 1.2 hours
Reported-by: Junchao Zhang <junchao.zhang@gmail.com>
Thanks-to: Stefano Zampini <stefano.zampini@gmail.com>

show more ...

Merge branch 'stefanozampini/feature-cusparse-notransgen-herm' into 'master'

feature cusparse notransgen herm

See merge request petsc/petsc!2915

MATSEQAIJCUSPARSE: implement hermitian transpose operation

Simplify the class, using a single routine for all the MatMult variants
Make transpose generation optional (not default)

Merge branch 'jose/scalapack' into 'master'

MATSCALAPACK: new Mat type that wraps a ScaLAPACK matrix

See merge request petsc/petsc!2858

Fix compilation in opensolaris

"/export/home/glci/builds/zxfiRCGf/0/petsc/petsc/include/petscmat.h", line 165: Warning (Anachronism): Formal argument 4 of type extern "C" int(*)(_p_Mat*,MatFactorTyp

Fix compilation in opensolaris

"/export/home/glci/builds/zxfiRCGf/0/petsc/petsc/include/petscmat.h", line 165: Warning (Anachronism): Formal argument 4 of type extern "C" int(*)(_p_Mat*,MatFactorType,_p_Mat**) in call to MatSolverTypeRegister(const char*, const char*, MatFactorType, extern "C" int(*)(_p_Mat*,MatFactorType,_p_Mat**)) is being passed int(*)(_p_Mat*,MatFactorType,_p_Mat**).
"/export/home/glci/builds/zxfiRCGf/0/petsc/petsc/include/petscmat.h", line 167: Error: Formal argument 6 of type extern "C" int(*)(_p_Mat*,MatFactorType,_p_Mat**)* in call to MatSolverTypeGet(const char*, const char*, MatFactorType, PetscBool*, PetscBool*, extern "C" int(*)(_p_Mat*,MatFactorType,_p_Mat**)*) is being passed int(*)(_p_Mat*,MatFactorType,_p_Mat**)*.

This error appears when including petscmat.h from a C++ example in
opensolaris. The fix is a typedef with PETSC_EXTERN_TYPEDEF for the
function argument.

Thanks-To: Barry Smith

show more ...

MATSCALAPACK: new Mat type that wraps a ScaLAPACK matrix