MatSeqDense: added PtAP support for P dense and A either SeqAIJ or SeqDense
add MAT_INPLACE_MATRIX; replace MAT_REUSE_MATRIX with it for MatConvert()
Added missing MatMissingDiagonal() implementations
MatHeaderReplace() corrupted the -objects_dump arrayReported-by: Torquil Macdonald Sørensen <torquil@gmail.com>MatHeaderReplace() and MatHeaderMerge() destroy the second matrix argument therefor
MatHeaderReplace() corrupted the -objects_dump arrayReported-by: Torquil Macdonald Sørensen <torquil@gmail.com>MatHeaderReplace() and MatHeaderMerge() destroy the second matrix argument therefor pass it asa pointer to it can be zeroed and not mistakenly reused.
show more ...
MatSeqDense: fix for MatDiagonalScale when lda != number of rows
Allow user to create the C matrix for the MatMatMult() function in the case when it is dense and use MAT_REUSE_MATRIX instead of requiring the the initial call to MatMatMult() to create for you.
Merge branch 'stefano_zampini/matmumps-schursupport'
Merge branch 'barry/propagate-pcsetup-failures'
MatDENSE: fix string for blas calls
Infrastructure that allows failures in PCSetUp(), PCApply(), MatMult() etc due to, for example a zero pivot or function evaluation outside its domain, to propagate up and become a KSP_DIVERGED inste
Infrastructure that allows failures in PCSetUp(), PCApply(), MatMult() etc due to, for example a zero pivot or function evaluation outside its domain, to propagate up and become a KSP_DIVERGED instead of generating an error that stops the program. In response to Issue 96. This includes failures in MatCreateSNESMF() applications due to domain errors. The mechanism to propagate some errors is by setting Info or Nan into the output vector and using the norm or inner product reductions in SNES or KSP to propagate the error condition to all processes and then handling them immediately after the norm or inner product.This allows, for example, ODE integrators to try again with a smaller time-step if the PCSetUp failed instead of requiring a complete restart of the run with other options.Currently some error conditions, such as function domain error in a line search may not get progated up using the correct SNESConvergedReasonSee src/snes/examples/tests/ex69.c for the handling of several conditions
MatDENSE: this fix the case when MAT_REUSE_MATRIX is specified with a LDA different from the number of rows of the reused matrix.
updated PETSc directory layout to match standard packaging strategiesinclude/petsc finclude,private,mpiunilib/petsc confbin/petsc*
Optimized MatShift() for a completely empty matrix since this is commonly used with SLEPc to produce diagonal matricesReported-by: Jose Roman <jroman@dsic.upv.es>
Merge branch 'maint'
MatLoad_SeqDense() always allocated new space for the data even if the user already provided itContributed-by: Alejandro Lamas Davina <alejandro.lamas@dsic.upv.es>
Merge branch 'stefano_zampini/seqaijtodense'
Merge branch 'barry/patch-flopcounts'Conflicts: src/mat/impls/dense/seq/dense.c
MatLoad_ and VecLoad_ routines that use options database information from the binary file .info files need to PetscViewerSetUp() before accessing the options database because the "lazy" setup of bina
MatLoad_ and VecLoad_ routines that use options database information from the binary file .info files need to PetscViewerSetUp() before accessing the options database because the "lazy" setup of binary viewers defered the opening of the files until binary reading was first done which was after the calls to determine matrix and vector block sizes.
Prevent integer multiplication overflow for some flop loggingReported-by: Stefano Zampini <stefano.zampini@gmail.com>
MATSEQAIJ: added conversion from to MATSEQDENSE
Mat: Added read-only Vec fixes for MATDENSE
Merge branch 'barry/feature-vec-lock'
Added VecLockPush() and VecLockPop() to make vector read onlyUsed VecLockPush() and VecLockPop() in a couple of key locations; PCApply(), KSPSolve(), MatMult()Fixed many uses of VecGetArray() to Ve
Added VecLockPush() and VecLockPop() to make vector read onlyUsed VecLockPush() and VecLockPop() in a couple of key locations; PCApply(), KSPSolve(), MatMult()Fixed many uses of VecGetArray() to VecGetArrayRead()Added VecGet/RestoreArrayPair() for commonly occuring patternAdded const to many locations driven by use of VecGetArrayRead()
bugfix for MatConvert_MPIDense_Elemental(); add ex173.c for testing Matload()+MatCovnert()
manual merge of conflict
1...<<11121314151617181920>>...38