Merge branch 'barry/fix-monitorlg'

previously the KSP and SNES LG norm monitoring routines crashes on -nox
Fixes this by keeping the PetscDraw object available and destroying it rather than obtaining it from the PetscDrawLG object
whi

previously the KSP and SNES LG norm monitoring routines crashes on -nox
Fixes this by keeping the PetscDraw object available and destroying it rather than obtaining it from the PetscDrawLG object
which did not work for -nox

show more ...

Changing -snes_converged_reason and -ksp_converged_reason to take a general viewer argument

Remove XXXRegisterAll() from public interface

This means we are relying on packages to be well behaved, to set
themselves up properly. Passes testexamples in my builds (haven't tried
complex).

Remove XXXRegisterAllCalled from public interface.

All XXXRegisterAll() methods are now idempotent, so a boolean indicating
whether they can be called or not is no longer necessary.

Merge branch 'master' into mlange/refactor-plex-distribute

Conflicts:
src/dm/impls/plex/examples/tests/output/ex12_1.out
src/dm/impls/plex/examples/tests/output/ex12_2.out

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

Merge branch 'barry/fix-setters'

Merge remote-tracking branch 'origin/tisaac/plex-flexible-projection'

* origin/tisaac/plex-flexible-projection:
PetscDualSpace: Fixed function signature
Plex: update libraries for ex3
Plex: mo

Merge remote-tracking branch 'origin/tisaac/plex-flexible-projection'

* origin/tisaac/plex-flexible-projection:
PetscDualSpace: Fixed function signature
Plex: update libraries for ex3
Plex: move DMPlexProjectionFunction() back to libdm
Plex: use -dm_plex_max_projection_height in non-conforming examples
Plex: fixed a typo in DMPlexProjectFunctionLocal()
Plex: added DMPlexComputeProjection3Dto1D_Internal()
Plex: fixes to recent merge
PetscFE: added dim, dimEmbed to PetscFECellGeom
Plex: skip empty loops for empty dual spaces
DualSpace: 0-order Lagrange spaces do not have dual spaces for height > 0.
Plex: protect against field projection with height > 0.
Plex: start implementing projection with height > 0.
PetscDualSpace: added PetscDualSpaceGetHeightSubspace()
Plex: added DMPlexGet/SetMaxProjectionHeight()

Conflicts:
src/dm/impls/plex/plexcreate.c

show more ...

added a variety of getters and fixed calling sequence of a few setters

Merge branch 'barry/fix-ksp-fcg-trunction-type'

shortened spelling of KSP_FCG_TRUNCTION_TYPE to fit within 72 characters in Fortran include file

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

Plex: move DMPlexProjectionFunction() back to libdm

When I first implemented hanging constraints, I moved
DMPlexProjectFunction() to libksp because I used a global solve to
determine the values of d

Plex: move DMPlexProjectionFunction() back to libdm

When I first implemented hanging constraints, I moved
DMPlexProjectFunction() to libksp because I used a global solve to
determine the values of dofs on an anchor point that is not in the
closure of any cell. This was a hack, as those values do not require a
global solve, but can be determined with functionals evaluated on the
anchor point: -dm_plex_max_projection_height > 0 allows functionals to
be evaluated on these anchor points, so the global solve is no longer
necessary.

To see the effect, you can modify the scripts to run ex3 in
src/dm/impls/plex/examples/tests: set -dm_plex_max_projection_height 0
in the runex3_nonconforming_{simplex,tensor}_{2,3} examples, and
projecting a function fails.

I did not move DMPlexProjectField() back into libdm because it evaluates
field jets, not just field values, and doing this on non-cells isn't
straightforward.

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

Merge branch 'master' of bitbucket.org:petsc/petsc into mark/ex49-bs

Merge branch 'barry/docs-matblocksize'

MatGetVecs() and KSPGetVecs() replaced with MatCreateVecs() and KSPCreateVecs()
since they are creators not getters

Merge branch 'psanan/ksp-fcg'

changed some syntax to closer match PETSc style in KSPFCG

Add KSPFCG

This adds the Flexible Conjugate Gradient Method as described and analysed in Notay's 2000 paper. Most of the available options concern storage and use of previous search directions for o

Add KSPFCG

This adds the Flexible Conjugate Gradient Method as described and analysed in Notay's 2000 paper. Most of the available options concern storage and use of previous search directions for orthogonalization. Notay's recommended truncation strategy as well as a standard truncation strategy are included.

Remaining to investigate/implement are more options for the orthogonalization procedure (currently standard Gram-Schmidt is implemented), and efficient implementation of residual monitoring (residuals norms are currently computed explicitly).

show more ...

Move DMPlexProject*() functions to src/ksp to remove ksp dependency in dm

removed extra call to KSPSetFromOptions() on the coarsest level of multigrid
KSP and PC skip displaying other options when preonly and none are selected
fix Chebyshev when using GMRES to estimate eig

removed extra call to KSPSetFromOptions() on the coarsest level of multigrid
KSP and PC skip displaying other options when preonly and none are selected
fix Chebyshev when using GMRES to estimate eigenvalues to skip setting options for PC that is shared with both Cheby and gmres

show more ...

Merge branch 'master' into prbrune/snes-fascontinuation

* master: (855 commits)
fix for moab example had wrong blocksize of mapping
docs:remove src/docs/doxygen/manual [per Barry's instructions]

Merge branch 'master' into prbrune/snes-fascontinuation

* master: (855 commits)
fix for moab example had wrong blocksize of mapping
docs:remove src/docs/doxygen/manual [per Barry's instructions]
gfortran: remove 132 char limit from f90 free-form
SNES: Set up matrix-free *after* copying the user Jacobian to snes->pc
DMPlex-Gmsh: Use trueNumCells when building the coordinate section
Sys: Got guard logic reversed
Sys: The #define guard for PetscMalloc*() must match the one for PetscFree*()
PetscFV: I hate fucking complex
removed ltog scatter from DMDA since functionality can be obtained with gtol scatter and less memory
fix test of DMLocalToGlobalBegin/End and add actual test cases
improve object memory logging for DA, ISLocalToGlobal and ISBLOCK
lower high water mark of constructing parallel scatter by allocating only as much room in svalues as needed
owner in VecScatter construction is the rank, therefor can use PetscMPIInt to hold it do not PetscInt
another reuse of allocated array in 2d/3d DMDA create
improved high water mark of DMDASetups by implementing and using PETSC_USE_POINTER for ISCreateBlock()
PetscFV: Another complex fix
DMPlex+TS: Fixing types
ISIdentity() was wrong for block IS! Indices are NOT scaled by bs
PetscFV: Straighten out complex types
TS ex11: Fix casts and remove unused variable
...

show more ...