DMDA: Fixed FE projection functions

Merged upstream/master to feature-dmmoab

Merge remote-tracking branch 'upstream/master'

Merge branch 'jed/dmda-memusage'

remove DMDAGetGlobalIndices() since handled by DMGetLocalToGlobalMapping() followed by ISLocalToGlobalMappingGetIndices()
fixed cases where bs of ISLocalToGlobalMapping in Mat did not match matrix bl

remove DMDAGetGlobalIndices() since handled by DMGetLocalToGlobalMapping() followed by ISLocalToGlobalMappingGetIndices()
fixed cases where bs of ISLocalToGlobalMapping in Mat did not match matrix block size

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

Merge branch 'master' into barry/saws-options

Merge branch 'jed/dmda-memusage'

removed ltog scatter from DMDA since functionality can be obtained with gtol scatter and less memory

Merge branch 'master' into stefano_zampini/pcbddc-multilevel-improvement

Conflicts:
src/docs/website/documentation/changes/dev.html
src/ksp/ksp/examples/tutorials/ex59.c
src/ksp/pc/impls/bddc/bdd

Merge branch 'master' into stefano_zampini/pcbddc-multilevel-improvement

Conflicts:
src/docs/website/documentation/changes/dev.html
src/ksp/ksp/examples/tutorials/ex59.c
src/ksp/pc/impls/bddc/bddcprivate.c

show more ...

Fixing merge conflicts in dmregall.

Merged petsc/petsc into master

Merge branch 'master' into shri/ts-events

Conflicts:
src/ts/interface/ts.c

Merge branch 'master' into shri/feature-dm-circuit

Conflicts:
src/dm/impls/makefile

Merge branch 'knepley/fix-fem-multifield'

* knepley/fix-fem-multifield:
DMDT: Fixed Fortran PetscQuadrature definition
DMPlex ex3: Added Q2 tests - Made DA tests only 2x2 - Changed output format

Merge branch 'knepley/fix-fem-multifield'

* knepley/fix-fem-multifield:
DMDT: Fixed Fortran PetscQuadrature definition
DMPlex ex3: Added Q2 tests - Made DA tests only 2x2 - Changed output format - Quad tests now actually work
DMDA: Added DMDAComputeL2GradientDiff()
PetscSpace: Fixed polynomial derivative evaluation for quads
DMDA: Fix closure operations in 2D
DMPlex: Interpolation cleanup - Remove unused PetscSpace - Fix uninitialized variable
PetscFE: Fix PetscBool handling
SNES ex62: Updated output
DMPlex ex3: Use PetscFECreateDefault() and correct tensor tests
PetscFE: Fix PetscFECreateDefault() - The space should default to tensor on a non-simplex - The quadrature should be tensor Gauss on a non-simplex
DT: Added PetscDTGaussTensorQuadrature()
DT: Added citation for Gloub+Welsch method
PetscSection: Now global section can also be viewed
DMPlex: Allow NULL function for ProjectFunction()
PetscDualSpace: Added doc for PetscDualSpaceApply() - Cherry-picked from df2d8bb
PetscFE: Fixed PetscFECreateDefault() - Setup FE - Minimum quad order - Destroy quad after setting
PetscFE: PetscFECreateDefault() now takes a number of components and a reference cell type - Cherry-picked from a255314
PetscFE: Added PetscFECreateDefault() - Cherry-picked from 8919f5a
PetscFE: Allow lowest order cell Lagrange dual basis functions
DMPlex: Fix preallocation with cell variables - I had an incorrect loop optimization

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

show more ...

Merge branch 'master' into prbrune/matcoloring-greedy

DMDA: Added DMDAComputeL2GradientDiff()

Merge branch 'master' into knepley/fix-fem-multifield

* master:
tests: move ex151 out of complex tests to avoid compile errors on windows
DM: Fixed pointer comparison
Microsoft C++ compiler ne

Merge branch 'master' into knepley/fix-fem-multifield

* master:
tests: move ex151 out of complex tests to avoid compile errors on windows
DM: Fixed pointer comparison
Microsoft C++ compiler needs extra care with complex numbers and casting
DMPlex: C++ has fucking namespaces, use them assholes
Field has 3 components therefor need to initialize three components, Microsoft C compiler demands it
DMPlex ex3: Added interpolation test - Only works for 2D simplices - Shows the restriction only gets constants right - Other small cleanup
DMPlex ex3: Change to interpolated meshes by default
DMPlex: Fix interpolation of multicomponent fields
DM: Fix reference counting in DMSetField()
DMPlex: Fix interpolation scaling
DMPlex: Fix allocation in DMPlexMatSetClosureRefined()
DMDA: Fixed signature of DMDACreateSection()
When configuring with MPIUni remove any MPI includes in the $PETSC_DIR/PETSC_ARCH/include directory that may have been put there by earlier ./configure runs.

show more ...

Merge branch 'knepley/feature-interpolation-testing'

* knepley/feature-interpolation-testing:
DM: Fixed pointer comparison
DMPlex ex3: Added interpolation test - Only works for 2D simplices - Sh

Merge branch 'knepley/feature-interpolation-testing'

* knepley/feature-interpolation-testing:
DM: Fixed pointer comparison
DMPlex ex3: Added interpolation test - Only works for 2D simplices - Shows the restriction only gets constants right - Other small cleanup
DMPlex ex3: Change to interpolated meshes by default
DMPlex: Fix interpolation of multicomponent fields
DM: Fix reference counting in DMSetField()
DMPlex: Fix interpolation scaling
DMPlex: Fix allocation in DMPlexMatSetClosureRefined()
DMDA: Fixed signature of DMDACreateSection()

show more ...

DMDA: Fixed signature of DMDACreateSection()

Merge branch 'knepley/feature-plex-fas'

* knepley/feature-plex-fas: (56 commits)
DMPlex: Removed unused variable
DMPlex: Fix for C++ build, small fix for PetscFE
DMPlex: Removed leak from Tria

Merge branch 'knepley/feature-plex-fas'

* knepley/feature-plex-fas: (56 commits)
DMPlex: Removed unused variable
DMPlex: Fix for C++ build, small fix for PetscFE
DMPlex: Removed leak from Triangle
SNES ex12: Two working unstructured FAS solves for the unregularized p-Laplacian
DMPlex: Cleaned up some memory issues
DMPlex: Implement injection - I don't think this should be required
DMPlex: DM interpolation requires the scaling vector for FAS
DMPlex: Can now create a map from coarse to fine vertices - This is currently necessary for FAS, which I think is wrong
DMPlex: Turn off Fortran stub
DMPlex: Fix DMPlexMatSetClosureRefined() - Reject duplicate fine points - Screen out point with dof = 0 - Fix the allocation for fine points
DMPlex: My brain did not work for several days - I have fixed the interpolation code - I think I can reuse that other code for something else
DMPlex: Fix setting Plex options on automatically refined meshed
SNES ex12/31/52/62/72: Use new SNESComputeJacobian() interface
DMPlex: Remove MatStructure from DMPlexComputeJacobianFEM() - Add missing declarations
Sys Logging: Fixes for log detail - Must preserve reduction count - Output objects and memory - Output stage summary
SNES ex12: Fixing things for multilevel - Need to set PetscFE on all coarse grids - Need application context for interpolator construction
SNES ex12: Use PetscFESetUp(), and fix quadrature ref handling
SNES ex52/62: Use PetscFESetUp(), and fix quadrature ref handling
DMPlex ex3: Uses automatic section creation, PetscFESetUp(), and fix quadrature ref handling
DMPlex: Fixing interpolator construction - Fixed definition of rectangular sizes - Use PetscFERefine() - Use special MatSet
...

show more ...

Merge branch 'master' into barry/google

Merge branch 'master' into knepley/feature-plex-fas

* master: (863 commits)
Revert "mpich: use latest mpich tarball with gcov but use '--enable-coverage' mpich configure option."
MatCreateTransp

Merge branch 'master' into knepley/feature-plex-fas

* master: (863 commits)
Revert "mpich: use latest mpich tarball with gcov but use '--enable-coverage' mpich configure option."
MatCreateTranspose: fix for lazy block size
mpich: use latest mpich tarball with gcov but use '--enable-coverage' mpich configure option.
configure: add FC_DEFINE_FLAG to petscvariables
added a citation for pounders, Stefan finally came through
update to latest netcdf fixes bug that did not allow building on recent Mac OS X
DMPlex: Fixes for CGNS - Fixed loop bounds since CGNS uses closed intervals - Bug with mixed cell types - Typo with coordinates
final updates from nightly tests
fix for countcfunction
removed outdated docs and comments on MatStructure usage removed SAME_PRECONDITIONER since no longer used removed unneed use of MatStructure in Schur complement matrix
removed MatStructure arguments from Tao
fixes for Sundials and fun3d cannot use MPIU_SUM for integer arrays
inplace ILU factorization must update state of inplace matrix to prevent infinite recursion of updating PC
fixes for BDDC by the way not having BDDC always in the testing is going to lead to trouble
the compute Jacobian functions for SNES and TS now do not take a MatStructure flag
KSPSetOperators() no longer has the MatStructure argument. The Mat objects now track that information themselves. Use KPS/PCSetReusePreconditioner() to prevent the recomputation of the preconditioner if the operator changed in the way that SAME_PRECONDITIONER did with KSPSetOperators()
tracks when Mat's have a change in nonzero structure this information can be used by PC to determine how preconditioners must be rebuilt currently only does an error check if the user said same nonzero pattern but actually changed the pattern
completed cleanup of removal of matrix pointers to compute jacobian routines
makefile: Added 'make includegraph' to create an inclusion dependency graph.
Contrib: Updated include-graph.sh script to ignore *.hh files
...

Conflicts:
src/dm/dt/interface/dt.c

show more ...

Merging master with feature-dmmoab. Only conflict was with SOURCECXX specification under impls/moab/makefile.

removed unneeded includes in include files
added findvisibleprivatefunctions for finding symbols in libraries that might not belong
removed a couple of global symbols
named spaced a couple of Tao fun

removed unneeded includes in include files
added findvisibleprivatefunctions for finding symbols in libraries that might not belong
removed a couple of global symbols
named spaced a couple of Tao functions

show more ...