Add subtype concept to SeqAIJ matrices

Introduced -mat_seqaij_type <seqaijperm, ...> MatSeqAIJSetType() to allow running code written for AIJ to work
with any subtype. See MatSeqAIJRegister() and Ma

Add subtype concept to SeqAIJ matrices

Introduced -mat_seqaij_type <seqaijperm, ...> MatSeqAIJSetType() to allow running code written for AIJ to work
with any subtype. See MatSeqAIJRegister() and MatSeqAIJRegisterAll() for registering new subtypes

Use

export PETSC_OPTIONS="-mat_seqaij_type seqaijperm -malloc_test"
make alltests

to run the entire test suite using an alternative subclass, in this case seqaijperm. This will produce output regarding the different name for
the matrix class but otherwise should produce no errors.

This was implemented to 1) allow easy testing of an MKL subtype and 2) to allow users to easily use the MKL subtype
without changing their code.

Also fixes several bugs in seqaijperm, if it isn't tested its broken.

Commit-type: testing-fix,style-fix, feature
Funded-by: ECP
Project: MKL
Time: 32 hours
Reported-by: Irina Sokolova <irina.sokolova@intel.com>

show more ...

configure: add option --with-fortran-bindings

Now one can use --with-fc=gfortran --with-fortran-bindings=0 to detect fortran compiler properly
[for externalpacakge use - for eg: scalapack] but avoid

configure: add option --with-fortran-bindings

Now one can use --with-fc=gfortran --with-fortran-bindings=0 to detect fortran compiler properly
[for externalpacakge use - for eg: scalapack] but avoid building PETSc fortran bindings in library and modules.

Also disable sowing when with-fortran-bindings is disabled

Fix makefiles (corresponding to fortran bindings) to consistantly use PETSC_HAVE_FORTRAN

Also remove unused PETSC_USE_FORTRAN_DATATYPES flag
Reported-by: Patrick Farrell <patrick.farrell@maths.ox.ac.uk>

show more ...

MatCopy implementations need to update the state since values have changed

Commit-type: bug-fix
Development Tools: Vim, Emacs, Eclipse

Merge remote-tracking branch 'origin/barry/update-petscfileretrieve'

Merge remote-tracking branch 'origin/barry/consistent-formating-ascii-solver-viewers'

update help menu for MatRARt() -- MatRARt_mpiaij_mpiaij is supported now, but not efficient.

update help menu for MatRARt as suggested by Franck Houssen

update help menu for MatRARt()

add default RARt=R*A*R^T as suggested by Jed

Function does actually exist for MPIAIJ

Commit-type: documentation
Reported-by: Franck Houssen <franck.houssen@inria.fr>

Mat: add MAT_SPD case to MatGetOption

MAT: add SeqDense to the documentation of MatMatTransposeMult

Merge branch 'stefano_zampini/feature-improvematfactorschur'

Implement Barry's comments

I didn't realize there was a hacky special case for dense B matrix

Commit-type: documentation
Reported-by: Pierre Jolivet <Pierre.Jolivet@enseeiht.fr>

Fix docs for MatMatMult() that C matrix can only be created by user for sequential case.

Commit-type: documentation
Time: .2 hours
Reported-by: Pierre Jolivet <pierre.jolivet@enseeiht.fr>

add MAT_STRUCTURE_ONLY to MatOption for create matrix with structure only

fix block sizes of MatCreateMPIAIJWithSeqAIJ()

Added another test case for extchem

Switched extchem to use dense matrix since MatSetValues() is much much faster
fix outdate MatOptions[] definitions

Commit-type: bug-fix, feature, example
Funded-

Added another test case for extchem

Switched extchem to use dense matrix since MatSetValues() is much much faster
fix outdate MatOptions[] definitions

Commit-type: bug-fix, feature, example
Funded-by: IMEX
Project: chemistry
Time: 3 hours
Thanks-to: "Aithal, Shashikant M." <aithal@cels.anl.gov>

show more ...

Mat: silent overzelous error

we should be able to duplicate matrices if MAT_COPY_VALUES has not been requested

Mat: allow inplace factorization of the Schur complement

Mat: Hermitian and Symmetric are synonyms with reals

Mat: Added schur and schur_status to the private implementation

This permits to simplify a lot the code for handling Schur complements computations
from factorization packages.
Added factorschur.c t

Mat: Added schur and schur_status to the private implementation

This permits to simplify a lot the code for handling Schur complements computations
from factorization packages.
Added factorschur.c to store intermediate routines
Removed MatFactorSetSchurComplementSolver that was horrible

show more ...

Mat: added interface function for MatMatSolveTranspose

Mat: added MatFactorSchurStatus to enable the user to progress factorization/inversion of the Schur complement