DMPlex: Removed debugging prints

DMPlex: Added DMPlexOrient() which correctly orients a set of adjacent cells
- You need this when extracting an embedded submesh

Merge branch 'knepley/plex'

Improved support for discretizing PDEs. Cell geometry and FVM. Little to
no topology.

* knepley/plex:
SNES ex12: Fix leak
TS ex11: error if dim != DIM because memory

Merge branch 'knepley/plex'

Improved support for discretizing PDEs. Cell geometry and FVM. Little to
no topology.

* knepley/plex:
SNES ex12: Fix leak
TS ex11: error if dim != DIM because memory corruption occurs otherwise
DMPlex: All quadrilateral geometry tests working
DMPlex: remove unused variable
DMPlex: Hexahedron geometry tests now pass
DMPlex: More geometry work
DMPlex: use PETSC_HAVE_TGAMMA to get the build working on windows
configure: add test for tgamma() [this is missing on windows]
DMPlex: Added random transform tests for tets - Fixed memory leak
DMPlex: Fixed 3D geometry
DMPlex: FVM geometry tests now work for triangles
DMPlex: one more win32/complex fix
DMPlex: Enhanced geometry tests - Added random transforms of the reference triangle in 2D
DMPlex: Fixed complex problems in geometry calculation
DMPlex: Fixed Fortran test output
DMDT: M_PI --> PETSC_PI
DMPlex: Reorganized geometry testing
TS ex11: Put in 3D inflow and outflow ids from the command line
DMPlex: Added DMPlexRestoreCone/Support() to F90 interface definition header - Also added to test ex1f90
TS ex11: Need to make these ids input from the command line
DMPlex: Add Fortran bindings
DMPlex: Reorder tet output since I order tets differently
Builder: Link Fortran examples with Fortran linker
DMPlex: Fixed leak in example
DMPlex: Tutorial now shows how to output VTK
PetscViewer: Added VTK type for Fortran
DMPlex: Remove parameter names from prototype
DMPlex: Geometry fixes for C++ complex
DMAKKT: Remove from build
DMPlex: Now DMPlexCreateSection() only has an F90 interface
SNES: add Fortran bindings to DMSNESSet{Function,Jacobian}
SNES: remove superfluous casting in DM local Fortran bindings
SNES: compile zdmlocalsnesf.c
DMDT: Fix Fortran wrapper
DMPlex: Added Fortran interface for DMPlexMatSetClosure()
PetscSection: Added docs, enabled Fortran wrappers
DM: Added Fortran bindings for DMSNESSetFunction/JacobianLocal()
SNES ex12: Looks like 3D Neumann conditions are working
DMPlex: Fix volume determination for 2D in 3D
DMPlex: Better error reporting
SNES ex12: Boundary integration seems to work in 2D
DMPlex+FEM: Fix for boundary integration
TS ex11: Better error reporting
DMPlex ex8: Added FVM geometry tests
DMPlex: Fixed FVM geometry for 2D in 3D
SNES ex62: Now using PetscDT quadrature
DMDT: Translated FIAT's Gauss-Jacobi quadrature
DMPlex: Fix geometry tests - Fix tests for 2D projection
DMPlex: Geometry fixes
DMPlex: Small fixes for geometry
DMPlex: Fixed damn prototype
PetscSection: Added VecSetValuesSectionF90() - Fixed bug in Vec Fortran header
DMPlex: Added Fortran defines for additional insert modes
DMPlex: Silence warning in DMPlexInvertCells_Internal()
DMPlex: Functions passed to evaluation routines now return void and pass results in arguments
DMPlex: Now we give an explicit embedding dimension to the FVM geometry methods
DMPlex: Fixed complication with complex - Started to add 3D geometry stuff
TS ex11: Replace edge geometry with call to DMPlexComputeCellGeometryFVM()
DMPlex: Added normal argument to DMPlexComputeCellGeometryFVM() - This is only calculated for faces - Added calculation for faces in 2D
DMPlex: Damn makefiles
DMPlex: Added declaration for DMPlexCreateCGNS()
TS ex11: Reorganized computation of geometry in preparation for 3D
Builder: Now individual tests can have requirements - The key is 'requires', and it takes a list of package names - Also fixed up showSingleRun()
TS ex11: Use new DMPlexComputeCellGeometryFVM()
DMPlex: Added DMPlexComputeCellGeometryFVM() - Added internal volume methods
DMPlex: fix const and int vs. PetscInt
DMPlex: Better error reporting
DMPlex+ExodusII: Fixed reading of quads - Was broken by 2e1b13c25062c3c40593ce7412c5cd227259ade7
DMPlex: Hex cell geometry was broken - I don't see how this test passed before
DMPlex: Regression cleanup
DMPlex: Regression cleanup
SNES ex12: Fixed sign for Neumann BC
DMPlex+FEM: The dimension does not change for boundary elements
DMPlex: Geometry now works for boundary elements
DMPlex+FEM: Fixed quadrature coordinate handling for boundary integrals
SNES ex12: Reorganized tests, still working on Neumann conditions - Fixed boundary face label
DMPlex: Turn off CGNS by default
DMPlex: Added 1D cell geometry
DMPlex: DMPlexCreateCGNS()
DMPlex: Updated test ex1 output
DMPlex: Fixed compiler problems with complex
SNES ex12: Now does boundary integration - Still not verifying exact solution for Neumann conditions
DMPlex: Fix DMPlexComputeCellGeometry() to handle lower dimensional cells
DMPlex: Fixed error in projection from 3D to 2D - Need to check for case where normal is already z
DMPlex: Adding boundary integration
DMPlex: Fixed cell inversion for TetGen to commute with mesh interpolation
SNES ex12: Cleaning up ex12 testing
DMPlex: Added integration over boundary to DMPlexComputeResidualFEM() - Added quadBd, f0/1BdFuncs to PetscFEM struct - Needs to be tested
FEM Generation: We now allow a *_bd.h header that holds boundary discreization
Builder: Fixed regression requirements
Builder: Fixed specification of multiple tests to run
DMPlex ex7: Orientations now working properly
DMPlex: Changed orientation convention and fixed many bugs
SNEX ex12: Ignore generated headers
SNEX ex12 and ex52: Remove CTetGen warnings from test output - Fixed this in CTetGen repository
DMPlex: Fix default |J| for DMPlexComputeCellGeometry()
DMPlex: Fix for cell geometry (unbelievable screwup)
DMPlex ex7: Check that |J| is nonzero for each cell in interpolated mesh
DMPlex: Allow DMPlexComputeCellGeometry() to work with interpolated meshes
DMPlex: Added prototype for DMPlexCopyCoordinates(), and docs for that and DMPlexInterpolate()
DMPlex: Added F90 interface for DMPlexComputeCellGeometry()

Conflicts:
include/petscdmplex.h

show more ...

Merge branch 'knepley/pylith'

* knepley/pylith:
DMPlex: Add missing declaration - This should merge alright with the other branch which fixes it
DMPlex: Submesh fixes - Coordinates must use a fi

Merge branch 'knepley/pylith'

* knepley/pylith:
DMPlex: Add missing declaration - This should merge alright with the other branch which fixes it
DMPlex: Submesh fixes - Coordinates must use a field - Allow same subpointMap to be specified again
DMPlex: Fix submesh marking for uninterpolated meshes - This now works for embedded submeshes
DMPlex: Interpolation corner cases - Do not interpolate a 1D mesh - No need to copy coordinates of the same DM
DMPlex: Interpolation of a submesh now starts at the faces - Do not know a better way to do this right now
DMPlex: Fix error in work space allocation
ISieve: Coordinate conversion should add a field

show more ...

DMPlex: Fix error in work space allocation

Merge branch 'knepley/pylith'

* knepley/pylith:
DMPlex: Fixed DMPlexCreateSohesiveSubmesh() - We should take the fault face directly from the cohesive cell rather than using the face insertion rou

Merge branch 'knepley/pylith'

* knepley/pylith:
DMPlex: Fixed DMPlexCreateSohesiveSubmesh() - We should take the fault face directly from the cohesive cell rather than using the face insertion routine
DMPlex: Add argument to DMPlexCreateSubmesh() that allows user to choose stratum - Added Fortran binding - Updated test
DMPlex: Added DMPlexCreateCohesiveSubmesh() - Needed to have submesh faces in the same order as cohesive cells
DMPlex: Fix improper preallocation - The size for a closure is the sum
Sieve: Missing VecDestroy() in conversion to DM
DMPlex: Submesh extraction now works without coordinates
ISieve: Fix conversion of orientations to DMPlex

show more ...

Merge branch 'jed/matmatmult-sort'

New version provides in the neighborhood of 20% speedup for MatMatMult
as used in PCGAMG.

* jed/matmatmult-sort:
MatMatMultSymbolic_SeqAIJ_SeqAIJ: switch to "so

Merge branch 'jed/matmatmult-sort'

New version provides in the neighborhood of 20% speedup for MatMatMult
as used in PCGAMG.

* jed/matmatmult-sort:
MatMatMultSymbolic_SeqAIJ_SeqAIJ: switch to "sort" version as default
segbuffer: g++ does not like casting int*restrict* to void*, void** okay
MatMatMult_SeqAIJ_SeqAIJ_Sort: new implementation using char flags and sorting
segbuffer: add PetscSegBufferGetInts() to encourage use of PETSC_RESTRICT
segbuffer: create wrapper structure
segbuffer: add PetscSegBufferGetSize() and PetscSegBufferUnuse()
segbuffer: fix bug in which seg->used was not reset

show more ...

Merge branch 'jed/matmatmult-sort' into jed/fix-optionsinsertfile-large

Branch contains segbuffer enhancements in addition to
MatMatMult_SeqAIJ_SeqAIJ_Sort. I need the segbuffer enhancements, one of

Merge branch 'jed/matmatmult-sort' into jed/fix-optionsinsertfile-large

Branch contains segbuffer enhancements in addition to
MatMatMult_SeqAIJ_SeqAIJ_Sort. I need the segbuffer enhancements, one of
which is a portability fix that was applied after the
MatMatMult_SeqAIJ_SeqAIJ_Sort commit.

* jed/matmatmult-sort:
segbuffer: g++ does not like casting int*restrict* to void*, void** okay
MatMatMult_SeqAIJ_SeqAIJ_Sort: new implementation using char flags and sorting
segbuffer: add PetscSegBufferGetInts() to encourage use of PETSC_RESTRICT
segbuffer: create wrapper structure
segbuffer: add PetscSegBufferGetSize() and PetscSegBufferUnuse()
segbuffer: fix bug in which seg->used was not reset

show more ...

DMPlex: Fix improper preallocation
- The size for a closure is the sum

DMPlex: Added Fortran interface for DMPlexMatSetClosure()

DMPlex: Silence warning in DMPlexInvertCells_Internal()

Merged petsc/petsc into rmills/petsc master

Merge branch 'jed/fix-plex-const' into knepley/plex

* jed/fix-plex-const:
DMPlex: fix const and int vs. PetscInt

DMPlex: fix const and int vs. PetscInt

Merge branch 'jed/sftypes'

* jed/sftypes:
PetscSF: split petscsftypes.h out of petscsf.h

PetscSF: split petscsftypes.h out of petscsf.h

Callers of PetscSF functions must now explicitly include petscsf.h.

Merge "DMPlex: Made DMPlexGetNumFaceVertices() externally visible"

Merge branch 'barry/rm-xxxregisterdynamic'

* barry/rm-xxxregisterdynamic:
Registration: remove stale 'XXRegisterDynamic)' entries in man pages
TS examples: fix use of PetscFunctionList and add t

Merge branch 'barry/rm-xxxregisterdynamic'

* barry/rm-xxxregisterdynamic:
Registration: remove stale 'XXRegisterDynamic)' entries in man pages
TS examples: fix use of PetscFunctionList and add to nightlies
Changes: Note updates to XRegisterDynamic/PetscObjectComposeFunctionDynamic
PetscObjectComposeFunctionDynamic: remove stale docs and usage
developers.tex: remove complications from function composition with dlls
removed string version of function name for XXXRegister(), PetscFunctionListAdd() and PetscObjectComposeFunction()
changes: document PetscFunctionListAdd() API change
developers.tex: update documentation of PetscObjectComposeFunction
removed path and MPI_Comm arguments from PetscFunctionListFind/Add()
removed path argument to XXXInitializePackage() and XXXRegister()
removed XXXRegisterDynamic() but kept the APIs for everything else underneath the same phase I of the update to handling registering function pointers

show more ...

DMPlex: Fixed compiler problems with complex

DMPlex: Fixed cell inversion for TetGen to commute with mesh interpolation

Merged petsc/petsc into rmills/petsc master

Merge branch 'barry/rm-xxxregisterdynamic' into jed/ts-eimex

PetscObjectComposeFunctionDynamic() and TSRegisterDynamic() were
replaced by PetscObjectComposeFunction() and TSRegister(), both of which

Merge branch 'barry/rm-xxxregisterdynamic' into jed/ts-eimex

PetscObjectComposeFunctionDynamic() and TSRegisterDynamic() were
replaced by PetscObjectComposeFunction() and TSRegister(), both of which
drop the string name argument.

* barry/rm-xxxregisterdynamic: (82 commits)
...

Conflicts:
src/ts/interface/tsregall.c

show more ...

removed path argument to XXXInitializePackage() and XXXRegister()

segbuffer: add PetscSegBufferGetInts() to encourage use of PETSC_RESTRICT

segbuffer: create wrapper structure

The old structure exposed the linked list property of the implementation.
Now we have a simple wrapper structure so that implementation detail is
not leaked. This

segbuffer: create wrapper structure

The old structure exposed the linked list property of the implementation.
Now we have a simple wrapper structure so that implementation detail is
not leaked. This will also make it possible to implement some thread-safe
operations.

show more ...