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

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

Needed to merge master to update to new calling sequence PetscCommBuildTwoSided

Merged petsc/petsc into master

Merge branch 'barry/fix-directory-layout'

updated PETSc directory layout to match standard packaging strategies

include/petsc finclude,private,mpiuni
lib/petsc conf
bin/petsc*

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

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

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

Merge branch 'barry/remove-unneeded-paren-petscmalloc'

Conflicts:
src/mat/impls/baij/seq/baij.c

convert some PetscMalloc() to PetscMalloc1()
remove unneeded () in first argument to PetscMalloc/Calloc1()
fixes for MATLAB examples

Merge branch 'master' into knepley/fix-configure-petsc-dep

* master: (210 commits)
error out on special cases where DMLocalToGlobal_DA won't work
remove unneed uses of local DM vectors and DMLoc

Merge branch 'master' into knepley/fix-configure-petsc-dep

* master: (210 commits)
error out on special cases where DMLocalToGlobal_DA won't work
remove unneed uses of local DM vectors and DMLocalToGlobalBegin/End in examples removed to prevent cut and paste repeats of their use
openmpi: update to version 1.8.3
added more movies and pubs
more movies
MatGetDiagonalBlock() can be called before matrix is assembled
added a couple more movies/images
PetscLogEventGetId() was missing public prototype
simply output from make streams benchmark test
c2html: fixup so that C2HTML make variable is set [similar to cmake.py]
chaco: current blas test does'nt work with Cray CC [which automatically links with blas] and linux - where -lm needs to be listed explicitly. Switch the test to check if we are using PETSc modifled chaco
builddist: always install c2html/sowing
argonne logo for users manual cover
Remove PETSC_HAVE_CUSP_SMOOTHED_AGGREGATION, use PETSC_HAVE_CUSP instead.
mpich: use older tarball with cygwin/gnu compilers
updated information about packages PETSc can use in users manual
Updated users manual to latest ANL Report format requirements
mpich: update to mpich-3.1.3.tar.gz - remove -download-mpich-mpe [mpich no longer packages mpe] - remove -download-mpich-shared. It was added as a split from --with-mpi-shared to --known-mpi-shared and --download-mpich-shared in 50237a55163f90ac7341b723c6422dfa1dbce792 removing this gives back the control of shared=0/1 to package.py - remove download_solaris as the new tarball works with both gcov and on solaris - remove MPICHInstall() [and use GNUPackage.Install()] as the new tarball fixes 'make clean' requiring sudo
added more images and PETSc pubs
Doc: More refs
...

show more ...

Merge branch 'tisaac/feature-fe-geom-dim' into tisaac/plex-flexible-projection

Use the new PetscFECellGeom

Conflicts:
include/petscdmplex.h
include/petscfe.h
src/dm/impls/plex/plexfem.c

Merge branch 'tisaac/thplex'

This merge brings in related work in the following areas:
- Added a local constraint matrix to DM (see DMSetDefaultConstraints())
- Added point-to-anchor-point constrain

Merge branch 'tisaac/thplex'

This merge brings in related work in the following areas:
- Added a local constraint matrix to DM (see DMSetDefaultConstraints())
- Added point-to-anchor-point constraints in DMPlex (see
DMPlexSetAnchors()), which build on the DM constraints
- Added hierarchically-constrained points to DMPlex (see
DMPlexSetReferenceTree(), DMPlexSetTree()), which builds on the anchor
constraints
- Tests for these constraints in src/dm/impls/plex/examples/ex3.c (try
runex3_* examples in that directory)
- Added quadrature order (see PetscQuadratureGet/SetOrder()) (work by
Matt Knepley)
- Added default return values for DMGetCoordinateDim() (defaults to
DMGetDimension())

Conflicts:
include/petscdmplex.h
src/dm/dt/interface/dt.c
src/dm/impls/plex/examples/tests/ex3.c
src/dm/impls/plex/makefile
src/dm/impls/plex/plexfem.c

show more ...

DM+Plex: change the collective status of default constraints, anchors, and reference tree

The default constraints, anchors, and reference tree are changed to have
PETSC_COMM_SELF communicators.

The

DM+Plex: change the collective status of default constraints, anchors, and reference tree

The default constraints, anchors, and reference tree are changed to have
PETSC_COMM_SELF communicators.

The reasoning is that we would like tree constraints to cause only
minimal changes to the solution workflow demonstrated in
snes/examples/tutorials/ex12.c. We would like to insert only the
following code:

/* ... (construct/choose a reference tree) */
DMPlexSetReferenceTree(dm,refTree);
/* ... (construct the tree structure) */
DMPlexSetTree(dm,...);

This means that the presence of some type of constraints may not be
detected until DMPlexMatSetClosure(), within a cell loop. If that
happens, the following call graph occurs:

DMPlexMatSetClosure(dm,...) calls
DMGetDefaultConstraints(dm,...) calls
*DMCreateDefaultConstraints_Plex(dm,...) calls
DMPlexGetAnchors(dm,...) calls
*DMPlexCreateAnchors_Tree(dm,...)
DMPlexComputeAnchorMatrix_Tree_FromReference(dm,...) calls
DMGetDefaultConstraints(refTree,...) calls
*DMCreateDefaultConstraints_Plex(refTree,...) calls
DMPlexGetAnchors(refTree,...) calls
*DMPlexCreateAnchors_Tree(refTree,...)

The objects created in this cascade have to be local.

show more ...

Plex: complete rename of constraints to anchors, make matrix construction automatic.

Now all the user has to do is set the reference tree, set the tree, and
the constraint matrix is automatically co

Plex: complete rename of constraints to anchors, make matrix construction automatic.

Now all the user has to do is set the reference tree, set the tree, and
the constraint matrix is automatically constructed in the call to
DMGetDefaultConstraints()

show more ...

Plex: added createanchors method.

Now DMPlexGetAnchors() is like DMGetConstraints(): if a method exists to
create the constraints, it is called.

Plex: rename Constraints as Anchors for clarity.

Now that there are DMGet/SetConstraints() functions, some clarity is
needed:

DMGet/SetConstraints() are any type of linear constraints that are
cons

Plex: rename Constraints as Anchors for clarity.

Now that there are DMGet/SetConstraints() functions, some clarity is
needed:

DMGet/SetConstraints() are any type of linear constraints that are
considered part of the global-to-local map.

DMPlexGet/SetAnchors() are linear constraints that make a point's dofs
linear dependent on the dofs of a set of anchor points.

DMPlexGet/SetTree() are anchor constraints that derive from a
hierarchically refined mesh.

show more ...

Plex: Fix compiler warning

PlexTree: Use correct MPI type

Plex: Complex fixes

Plex: Complex fix

Plex: Fix calls to DMPlexGet/SetDimension()

Plex: Complex fixes