PCBDDC: optimized eigenvalue problem

explicit inversions of sum(S_E) and sum(St_E) blocks have been removed and the eigenvalue problem has been recasted
remove unneded transpositions for the symmetr

PCBDDC: optimized eigenvalue problem

explicit inversions of sum(S_E) and sum(St_E) blocks have been removed and the eigenvalue problem has been recasted
remove unneded transpositions for the symmetric case

show more ...

PCBDDC: restore lines erroneusly removed in previous commit

PCBDDC: get rid of invert_stildas branch

PCBDDC: reuse interior solver computed by MUMPS when computing the Schur complement

PCBDDC: add option to select faster application of deluxe scaling

QR factorization of constraints resulting in a unitary transformation is mandatory if a change of basis is requested with faster del

PCBDDC: add option to select faster application of deluxe scaling

QR factorization of constraints resulting in a unitary transformation is mandatory if a change of basis is requested with faster deluxe

show more ...

PCBDDC: handle the case for 1 subset only

eigenvalue problem still need to be investigated.

Merge branch 'master' into tisaac/dmforest

If we're going to start developing on this again, the recent changes to the
location of configuration files will make it annoying to switch to and from
thi

Merge branch 'master' into tisaac/dmforest

If we're going to start developing on this again, the recent changes to the
location of configuration files will make it annoying to switch to and from
this branch. Matt's and my fixes in DMPlex will also be relevant.

* master: (1557 commits)
MatLoad_MPIDense() always allocated array space even if it was previously provided
Clarified lumping options for PCFieldSplitSchurPreType SELFP.
fixed up PetscInfo in GAMG
TSErrorWeightedNorm[2|Infinity]: API changes
MatLoad_SeqDense() always allocated new space for the data even if the user already provided it
Fixed DMCompositeGetAccessArray() to support readonly vectors
Fix obviously-wrong aspects of documentation on SNESVINEWTONRSLS and SNESVINEWTONSSLS.
fixed example to compile - demonstrates a bug
var name 'len' was F90 key word
For user provided Jacobian functions they should always assembly into the second matrix, the one from which the preconditioner is built because the first may be used as matrix-free with -snes_mf_operator
Suggest --with-64-bit-indices if users mistakenly use --with-64bit-indices
Fix configure logging of integer size.
missing underscore in #define of vecgetarrayread translation of fortran stub for VECGETARRAYREAD
updated test output files to reflect differences introduced previously in PCView output
-ksp_type preonly does not support -ksp_monitor_short so remove from test line to prevent extranious warning output
missing CHKERRQ; and ugly double declaration of ierr
Plex ex3: For want of a semicolon, the kingdom was lost
MatSetValuesMPISBAIJ will now print correct global indices when error generated due to new nonzero location
completed basic use of MatSetValuesBlocked_SeqBAIJ_Inlined() by MatSetValuesBlocked_MPIBAIJ() to allow proper global indices being printed when error generated by new nonzero locations.
Added inline version of MatSetValuesBlocked_SeqBAIJ() to be used by MatSetValuesBlocked_MPIBAIJ() to print error message for nonzero block with global indices
...

show more ...

Merge branch 'master' into shri/ts-is-for-differential-variables

Conflicts:
include/petsc-private/tsimpl.h
src/ts/impls/eimex/eimex.c
src/ts/interface/ts.c

PCBDDC: get rid of parallel application of deluxe scaling

code has been greatly simplified
introduced subcommunicators in PCBDDCSubSchursSetUp to improve performances in multilevel cases
adaptivity

PCBDDC: get rid of parallel application of deluxe scaling

code has been greatly simplified
introduced subcommunicators in PCBDDCSubSchursSetUp to improve performances in multilevel cases
adaptivity is supported only in combination with deluxe scaling

show more ...

PCBDDC: multilevel guard on S_all

Merge branch 'master' into jed/chem
Note merging a branch with good stuff in it for a year and a half is insane!

Conflicts:
src/ts/examples/tutorials/makefile
src/ts/interface/ts.c

PCBDDC: optimized computation of Stildas matrices

Better to compute once the inverse of the local Schur complement and then extract principal minors
In the previous version, "n_subs" problems, each

PCBDDC: optimized computation of Stildas matrices

Better to compute once the inverse of the local Schur complement and then extract principal minors
In the previous version, "n_subs" problems, each of size "local_size-subset_size" where factorized and applied to "subset_size" right hand sides

show more ...

PCBDDC: added code for loop on rhs instead of computing explicitly the operator

PCBDDC: fixes for complex builds

PCBDDC: when bs>1 better to convert to SeqAIJ

added guard in PCBDDCComputeExplicitSchur when ksp has 0 size

PCBDDC: remove MatReuse argument from MatGetSubMatrixUnsorted

PETSCBLASLAPACK: added bindings for POTRI

updated sub_schurs to use potri
(bug fixes from previous commit for cum variable)

PCBDDC: reuse workspace for GetSubmatrix

lapack workspace renamed to Bwork to avoid conflicts

PCBDDC: manage symmetry for implicit schur complements

PCBDDC: optimize PCBDDCComputeExplicitSchur

PCBDDC: no need to store all S_Ej_tildas

PCBDDC: added slots for hermitianess and positive definitess in SubSchurs

local matrix used in SubSchurs should be converted to SeqAIJ

Merge branch 'master' into barry/saws-options-merge-master

Totally confusing manual merge, could be wrong

Conflicts:
include/petscoptions.h
src/dm/impls/da/dacreate.c
src/sys/objects/aoptions.c

Merge branch 'master' into barry/saws-options-merge-master

Totally confusing manual merge, could be wrong

Conflicts:
include/petscoptions.h
src/dm/impls/da/dacreate.c
src/sys/objects/aoptions.c
src/sys/threadcomm/interface/threadcomm.c

show more ...

Merge branch 'master' into pr252/dmay/PetscViewerBinarySkipHeaderVecViewLoad
Manual merge of makefile

Conflicts:
src/sys/classes/viewer/examples/tutorials/makefile

PCBDDC: multilevel fixes to SubSchurs