Merge branch 'barry/fix-nonew-notcollective/maint' into jed/mat-assembly-perf

Jed this fucking 6+ month running thing in next that is not in master is a royal fucking pain in the ass and
an abuse of

Merge branch 'barry/fix-nonew-notcollective/maint' into jed/mat-assembly-perf

Jed this fucking 6+ month running thing in next that is not in master is a royal fucking pain in the ass and
an abuse of git. I think it is not justified to ever have anything in next for more than a few weeks at most.
Either take it out of next if it is broken or put it in master if it is not broken. Hanging around in next but
not master for ever is not a good policy.

show more ...

Merge branch 'master' into knepley/feature-snes-deflation

* master: (239 commits)
Bib: Moved entry
revert sources to dev - reset PETSC_VERSION_RELEASE - reset changes/dev.html - enable cuda in c

Merge branch 'master' into knepley/feature-snes-deflation

* master: (239 commits)
Bib: Moved entry
revert sources to dev - reset PETSC_VERSION_RELEASE - reset changes/dev.html - enable cuda in configure
release: set v3.6 strings
doc: add PCGAMGSetReuseInterpolation() to changes file
fix: workarround static initialization issue (perhaps compiler bug) on opensolaris c/complex build
fix: ex40,ex41 errors in src/ts/examples/tutorials/ on bsd with c++/complex Perhaps this code nees a separate section when 'enorm == 0) - as pow(0,-0.33) = inf.
manual: fix format
more changes..
PCGASM/MatGetSubMatricesMPI() changes in dev.html.
changes: more updates
Add size checks for static arrays of assumed max size [-Warray-bounds]
Doc: Update to Changes
configure: default to useThreads= 0as the useThreads=1 code currently prints output in jumbled order [producing hard to read configure.log]
Delete PETSc viewer attribute if none are left on the linked list instead of putting an empty attribute
MatGetSubMatricesMPI_MPIAIJ(): guard for an empty B-matrix when delogging memory in MatSetSeqMats().
PCBDDC: fix typos in the specification of MPI operation
MPIU_SUM was used with MPIU_INT when it should only be used with MPIU_REAL or MPIU_SCALAR
fixed bad PetscPrintf() formatting in two examples
Plex: In pointwise callbacks, use embedding dimension
updates for generating Apple installations due to changes in locations of include files
...

show more ...

Merge branch 'master' into mesgarnejad/feature-plex-globaltonatural

* master: (219 commits)
revert sources to dev - reset PETSC_VERSION_RELEASE - reset changes/dev.html - enable cuda in configure

Merge branch 'master' into mesgarnejad/feature-plex-globaltonatural

* master: (219 commits)
revert sources to dev - reset PETSC_VERSION_RELEASE - reset changes/dev.html - enable cuda in configure
release: set v3.6 strings
doc: add PCGAMGSetReuseInterpolation() to changes file
fix: workarround static initialization issue (perhaps compiler bug) on opensolaris c/complex build
fix: ex40,ex41 errors in src/ts/examples/tutorials/ on bsd with c++/complex Perhaps this code nees a separate section when 'enorm == 0) - as pow(0,-0.33) = inf.
manual: fix format
more changes..
PCGASM/MatGetSubMatricesMPI() changes in dev.html.
changes: more updates
Add size checks for static arrays of assumed max size [-Warray-bounds]
Doc: Update to Changes
configure: default to useThreads= 0as the useThreads=1 code currently prints output in jumbled order [producing hard to read configure.log]
Delete PETSc viewer attribute if none are left on the linked list instead of putting an empty attribute
MatGetSubMatricesMPI_MPIAIJ(): guard for an empty B-matrix when delogging memory in MatSetSeqMats().
PCBDDC: fix typos in the specification of MPI operation
MPIU_SUM was used with MPIU_INT when it should only be used with MPIU_REAL or MPIU_SCALAR
fixed bad PetscPrintf() formatting in two examples
Plex: In pointwise callbacks, use embedding dimension
updates for generating Apple installations due to changes in locations of include files
configure: do not fail when checking MPICH version [similar to 2cadc8fc42fbad3c229d43b955581b8d29e02ad6 ]
...

show more ...

Merge branch 'master' of bitbucket.org:petsc/petsc

Merge remote-tracking branch 'petsc/master' into matpartitioning-hierach

Merge remote-tracking branch 'petsc/master' into mat-increaseoverlap

Merge branch 'stefano_zampini/pcbddc-adaptive-selection'

Conflicts:
src/ksp/pc/impls/bddc/bddcprivate.c
src/ksp/pc/impls/bddc/bddcscalingbasic.c

Merge branch 'master' into knepley/solkx

* master: (6933 commits)
Bib: Added reference
Includegraph: Updated for compatibility with new folder layout.
Add SNESLineSearchReason to fortran inclu

Merge branch 'master' into knepley/solkx

* master: (6933 commits)
Bib: Added reference
Includegraph: Updated for compatibility with new folder layout.
Add SNESLineSearchReason to fortran includes Also truncate names to fit 32 char fortran limit
fix compile warnings /usr/home/balay/petsc.clone-2/src/vec/vec/interface/vector.c:1944: warning: division by zero in '1.0e+0 / 0.' /usr/home/balay/petsc.clone-2/src/ksp/ksp/interface/itfunc.c:508: warning: 'vec_rhs' may be used uninitialized in this function
updated output for new trust region initial size
initial size of the trust region is set as a percent of the norm of the initial guess, it should not be related to the initial norm of the function (which is kind of nuts).
The default bounds for SNESVISetVariableBounds() in the manual page were reversed from correct values
Bib: Added refs
Bib: Update MPICH webpage
Plex: Forgot to close file
removed nonexistent ex46f from makefile test rule
Revert "fixed bad merge into master"
fixed bad merge into master
fix memory leak in mkl_pardiso fix some formatting in mkl_pardiso code
mv runex111 from TESTEXAMPLES_C to TESTEXAMPLES_DATAFILESPATH
Mat+Doc: More Fortran docs
Plex ex3: Run tests with the correct number of field components
Bib: Added TetGen cite
parmetis: make sure 'ldd libparmetis.so' points to the correct libmetis.so that its linked with.
small fix as reported by the nightly test
...

Conflicts:
config/builder.py

show more ...

Merged petsc/petsc into master

Merge branch 'master' into tisaac/dmforest

directory layout changed again: merging to switching branches back and forth
doesn't cause reconfigure to have a hard time

* master: (357 commits)
remov

Merge branch 'master' into tisaac/dmforest

directory layout changed again: merging to switching branches back and forth
doesn't cause reconfigure to have a hard time

* master: (357 commits)
removed nonexistent ex46f from makefile test rule
Revert "fixed bad merge into master"
fixed bad merge into master
fix memory leak in mkl_pardiso fix some formatting in mkl_pardiso code
mv runex111 from TESTEXAMPLES_C to TESTEXAMPLES_DATAFILESPATH
Mat+Doc: More Fortran docs
Bib: Added TetGen cite
parmetis: make sure 'ldd libparmetis.so' points to the correct libmetis.so that its linked with.
small fix as reported by the nightly test
add a test for MatPtAP()
small cleanup; add tests
cleanup MatPtAPSymbolic_MPIAIJ_MPIAIJ()
Fix nonzerostate tracking in all MATMPI types.
update location of test matrices in faq.html file
unused variable still declared
Fix for MatLoad_MPIDense() from sparsely stored matrix also fix for MatLoad_MPIDense_DenseInFile() when matrix has not had its local sizes preset by user
Preserve the B nonzerostate during disassembly.
TSEvent: Fix arch-mswin nightly errors due to declaring variables in the middle of a function.
Fix error in parallel loading of dense format matrix
removed KSPSetNullSpace() since mat already carries that information and having two locations for it is confusing, unnecessary and prone to error added MatSetTransposeNullSpace() for providing the null space of the transpose of the matrix and added support in KSPSolve() to remove this null space from the right handside of the linear system thus making the system consistent.
...

show more ...

Merge branch 'barry/propagate-pcsetup-failures'

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 SNESConvergedReason

See src/snes/examples/tests/ex69.c for the handling of several conditions

show more ...

Merge branch 'master' into knepley/feature-composable-callbacks

* master: (151 commits)
PetscMalloc: Fix handling of zero-sized allocations
add note to DMGetLocal/GlobalVector() about intended s

Merge branch 'master' into knepley/feature-composable-callbacks

* master: (151 commits)
PetscMalloc: Fix handling of zero-sized allocations
add note to DMGetLocal/GlobalVector() about intended short-time usage of the vectors.
sf: use patched openmpi tarball that has a potential fix for sf example hang
SNES problems that set VI bounds must use a VI solver since the other solvers ignore the bounds
pounders viewer did not cast PetscReal to double in printf nor use %D instead of %d
fix ts ex20adj
fix ts theta
fixed another variable declaration in the middle of a code block
example main() had PetscInt return type and argument that should be int
DMDAVTKWriteAll_VTS: Fix handling of higher-dimension coordinates
fixed example that declared variables in the middle of a function
C programs do not return PetscInt nor do they take an argument that is PetscInt, fixed example that did
You cannot declare variables in the middle of functions in PETSc.
add ierr checking to PetscOptionsBegin/End() in example
configure: add error check if user specifies PETSC_ARCH=-dbg [i.e starts with a '-']
scale norm of vector difference by norm of vector to insure it is small enough to consider the two vectors equal
missing PetscViewerDestroy() added
MatNestGetISs: use custom fortranstub to handle PETSC_NULL_OBJECT
update location of open-mx website on FAQ page
If the user sets bounds and calls a solver that does not support them, throw an error.
...

show more ...

PCBDDC: reuse as much as possible the MUMPS Schur solver

optimizations still need to be applied

PCBDDC: Implement caching strategy for dirichlet dofs in graph

Merge branch 'pr192/Fande-Kong/scalable-matincreaseoverlap/master' of https://bitbucket.org/petsc/petsc

Merged petsc/petsc into master

Merge branch 'master' into stefano_zampini/pcbddc-adaptive-selection

Conflicts:
src/ksp/pc/impls/bddc/bddcprivate.c
src/ksp/pc/impls/bddc/bddcscalingbasic.c
src/mat/impls/aij/mpi/mumps/mumps.c
s

Merge branch 'master' into stefano_zampini/pcbddc-adaptive-selection

Conflicts:
src/ksp/pc/impls/bddc/bddcprivate.c
src/ksp/pc/impls/bddc/bddcscalingbasic.c
src/mat/impls/aij/mpi/mumps/mumps.c
src/mat/impls/is/matis.h

show more ...

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

Merged master into semihozmen/ksp-fcg-eigen-comp-support

Merged petsc/petsc into master

Merge branch 'master' into pr192/Fande-Kong/scalable-matincreaseoverlap/master

Needed to merge master to update to new calling sequence PetscCommBuildTwoSided

fixed missing CHKERRQ() detected by clang's --analyze

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