DMPlex: Bug fix for hybrid tet refinement
- The new edges for hybrid faces do not depend on the cap face orientation
- This happens because we are not attaching hybrid faces in an order dependent on

DMPlex: Bug fix for hybrid tet refinement
- The new edges for hybrid faces do not depend on the cap face orientation
- This happens because we are not attaching hybrid faces in an order dependent on the cap face orientation. Maybe we should.

show more ...

DMPlex: Bug fix for hybrid tet refinement
- The orientation that we get for hybrid faces depends on both the edge and face orientation
- This happens because we are not attaching hybrid faces in an o

DMPlex: Bug fix for hybrid tet refinement
- The orientation that we get for hybrid faces depends on both the edge and face orientation
- This happens because we are not attaching hybrid faces in an order dependent on the cap face orientation. Maybe we should.

show more ...

Merge branch 'knepley/feature-plex-refine-hex'

* knepley/feature-plex-refine-hex:
DMPlex: Fixes for ViewFromOptions - Passes regression
DM: Moved ViewFromOptions into DMSetUp_DA()
DMPlex: Inte

Merge branch 'knepley/feature-plex-refine-hex'

* knepley/feature-plex-refine-hex:
DMPlex: Fixes for ViewFromOptions - Passes regression
DM: Moved ViewFromOptions into DMSetUp_DA()
DMPlex: Interface fix in refinement
PetscFE: Fixed leak from bad merge
DMPlex: DMPlexCheckSkeleton() now handles hybrid cells
DMPlex ex4: All refinement tests passing
DMPlex: Added hybrid hex refinement
DMPlex: Bug fix for hybrid tet refinement - Wrong offset
DMPlex: Fixed errors in hybrid tet refinement - Use fMax instead of fStart for offset - Use p instead of fEnd for point value - 2 faces are created from each hybrid face, not 4
DMPlex: Corrected refinement coordinates for hybrid faces
DMPlex ex4: Hybrid tet refinement now working
DMPlex: Propagate fixes from tet refinement to hybrid tet refinement

show more ...

Merge branch 'master' into knepley/feature-dmda-section

* master: (319 commits)
SNESGS: Fix for uninitialized variable warning in gssecant.c
DMPlex: Fix for orientation in Gmsh files - Fixed tes

Merge branch 'master' into knepley/feature-dmda-section

* master: (319 commits)
SNESGS: Fix for uninitialized variable warning in gssecant.c
DMPlex: Fix for orientation in Gmsh files - Fixed test output
do not use C99 struct initialization code in PETSc
Nightly Summary: Added time taken to each of the tests in the dashboard.
ASCII printing of floating point numbers must cast to (double) to support quad precision
configure: check fortranDatatypes only when fortran is enabled [this fixes configure crash with --with-fc=0]
don't use non-portable and also dependent on floating point precision system level math constants
picky compilers
fix a couple compiler warnings/errors (depending on compiler) for Matlab engine code
DMPlex: More fixup from bad rebase, moved to DMGet/SetCoordinateSection()
DM: Added DMGet/SetCoordinateSection()
Fuckup: Fix for merge that is still untangling bad rebase from Oct. 17th 2013 - Fixed const in declarations - Fixed merge from PetscFE - Fixed 2 mallocs
replace all echo in test/tutorial makefiles with printf since echo behaves differently on different systems
DMDA: Missing header - Damn rebase
DMPlex ex3: Added 2D Q_1 and 3D P_2 tests
DMDA: Added DMDAProjectFunction() and DMDAComputeL2Diff()
DMDA: Fixed 2D geometry - Should reuse DMPlex routines
DMDA: Fix bug in 2D closure operation for cells
DMDA: Make closure operations also return the size
DMDA: Changed DMDACreateSection() to emulate DMPlexCreateSection()
...

Conflicts:
config/builder.py
include/petscdmda.h
src/dm/dt/interface/dtfe.c
src/dm/impls/da/dageometry.c
src/dm/impls/da/dalocal.c
src/dm/impls/plex/plexgeometry.c

show more ...

Merged petsc/petsc into master

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

* master: (163 commits)
do not use C99 struct initialization code in PETSc
Nightly Summary: Added time taken to each of the tests in the dash

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

* master: (163 commits)
do not use C99 struct initialization code in PETSc
Nightly Summary: Added time taken to each of the tests in the dashboard.
ASCII printing of floating point numbers must cast to (double) to support quad precision
configure: check fortranDatatypes only when fortran is enabled [this fixes configure crash with --with-fc=0]
don't use non-portable and also dependent on floating point precision system level math constants
picky compilers
fix a couple compiler warnings/errors (depending on compiler) for Matlab engine code
DMPlex: More fixup from bad rebase, moved to DMGet/SetCoordinateSection()
DM: Added DMGet/SetCoordinateSection()
Fuckup: Fix for merge that is still untangling bad rebase from Oct. 17th 2013 - Fixed const in declarations - Fixed merge from PetscFE - Fixed 2 mallocs
replace all echo in test/tutorial makefiles with printf since echo behaves differently on different systems
DMDA: Missing header - Damn rebase
DMPlex ex3: Added 2D Q_1 and 3D P_2 tests
DMDA: Added DMDAProjectFunction() and DMDAComputeL2Diff()
DMDA: Fixed 2D geometry - Should reuse DMPlex routines
DMDA: Fix bug in 2D closure operation for cells
DMDA: Make closure operations also return the size
DMDA: Changed DMDACreateSection() to emulate DMPlexCreateSection()
DMDA: Added functions which emulate DMPlex functionality
DMDA: Now DMDAGetNumCells() returns the cells in each direction as well
...

Conflicts:
config/builder.py

show more ...

Merge branch 'madams/gamg-destroy' into prbrune/pcgamg-classicalinterpolationstrategies

Merge branch 'master' into madams/gamg-destroy

Conflicts:
src/ksp/pc/impls/gamg/gamg.c

DMPlex: Interface fix in refinement

DMPlex: More fixup from bad rebase, moved to DMGet/SetCoordinateSection()

Merge branch 'master' into barry/reduce-dmsetup-da-memoryusage

Conflicts:
src/dm/examples/tests/ex15.c
src/dm/examples/tutorials/ex3.c
src/dm/impls/da/da2.c
src/dm/impls/da/da3.c

Merged master into barry/update-xxxviewfromoptions

Merge branch 'master' into barry/xcode

DMPlex: Added hybrid hex refinement

DMPlex: Bug fix for hybrid tet refinement
- Wrong offset

DMPlex: Fixed errors in hybrid tet refinement
- Use fMax instead of fStart for offset
- Use p instead of fEnd for point value
- 2 faces are created from each hybrid face, not 4

DMPlex: Corrected refinement coordinates for hybrid faces

Merge branch 'barry/saws-push-header-body' into barry/saws-options

DMPlex: Propagate fixes from tet refinement to hybrid tet refinement

Merge branch 'knepley/feature-plex-hybrid-3d'

* knepley/feature-plex-hybrid-3d: (45 commits)
DMPlex ex5: Check faces on all meshes
DMPlex ex5: Fix consistency of original meshes - Fix orientatio

Merge branch 'knepley/feature-plex-hybrid-3d'

* knepley/feature-plex-hybrid-3d: (45 commits)
DMPlex ex5: Check faces on all meshes
DMPlex ex5: Fix consistency of original meshes - Fix orientations of tet meshes - Redo hex test to automatically interpolate - Replaced output
DMPlex ex5: Check original meshes for consistency
DMPlex: make non-namespaced CellRefiner functions static
DMPlex: Fix hybrid bounds in DMPlexCheckFaces()
DMPlex ex5: Now all test working - Two faults working correctly
DMPlex: Fix for hybrid meshes with multiple surfaces - Move the original hybrid cells to the end of the numbering, so that split+hybrid points are inserted as a block - Use DMPlexShiftPoint_Internal() uniformly
SNES ex4: Fix for DMCheckSkeleton()
DMPlex ex4: All non-hybrid refinement test pass (Fuck Yeah)
DMPlex: In tet refinement, fixed edge orientation in internal face
DMPlex ex4: Correct output after cross-tet edge orientation fix
DMPlex: In tet refinement, correct internal tet face orientations
DMPlex: In tet refinement, correct internal face cones
DMPlex: In tet refinement, correct check for cross-tet edge
DMPlex: Now DMPlexCheckFaces() also checks all subfaces
DMPlex: Added 1D stuff to DMPlexGetRawFaces_Internal()
DMPlex ex4: hex refinement now fixed
DMPlex: Fixed split face supports for tet refinement
DMPlex: Multiple refinement for hex now working - Fixed face orientations for hexes - Fixed split face orientations - Fixed split face supports - Fixed internal face orientations - Fixed internal edge supports
DMPlex: Added GetQuadEdgeInverse_Static() and GetQuadSubfaceInverse_Static() for refinement
...

Conflicts:
src/dm/impls/plex/examples/tests/ex4.c
src/dm/impls/plex/plexsubmesh.c

show more ...

Merge branch 'knepley/fix-plex-submesh-parallel'

* knepley/fix-plex-submesh-parallel:
DMPlex: Fixed leak in SF creation for submesh - Changed name of subpointIS to avoid confusion
DMPlex: Fixed

Merge branch 'knepley/fix-plex-submesh-parallel'

* knepley/fix-plex-submesh-parallel:
DMPlex: Fixed leak in SF creation for submesh - Changed name of subpointIS to avoid confusion
DMPlex: Fixed SF creation in hybrid tet refinement, added better check
DMPlex: Guard for empty submesh in parallel during subSF creation
DMPlex: Make DMPlexCreateCohesiveSubmesh_Interpolated() also create PetscSF
PetscSection: Fix for PetscSectionCreateGlobalSectionCensored() - Offset was wrong
DM: Added check in DMCreateDefaultSF() - Check for overrun of largest offset
PetscSection: Fixed PetscSectionGetOffsetRange() for empty section
DMPlex: Small doc fix for DMPlexConstructCohesiveCells()
DMPlex: Protect another instance of empty IS
DMPlex: Protect submesh selection when some procs have empty labels
DMPlex: Turn off check for cells bordering a cohesive face - In parallel, the negative side cell for a cohesive face can be on another process - We might want to consider communicating these

show more ...

DMPlex: make non-namespaced CellRefiner functions static

Merge branch 'jed/malloc-array'

Type arguments dropped from PetscMalloc[2-7], PetscNew, and PetscNewLog,
added PetscMalloc1 for allocating typed arrays, add PetscCalloc[1-7] for
allocating cleared (

Merge branch 'jed/malloc-array'

Type arguments dropped from PetscMalloc[2-7], PetscNew, and PetscNewLog,
added PetscMalloc1 for allocating typed arrays, add PetscCalloc[1-7] for
allocating cleared (zeroed) memory.

* jed/malloc-array:
Sys: Add pointer casts from (void **) in calls to PetscMalloc1()
Sys: drop explicit type arguments from PetscNew() and PetscNewLog()
Sys: add PetscCalloc[1-7]
Sys: add PetscMalloc1 macro, array allocation without redundant types
PetscMalloc[2-7]: remove type arguments, infer from pointer type
SNESComputeJacobianDefaultColor: fix uninitialized variable

show more ...

Merge commit '58723a9' into knepley/feature-plex-hybrid-3d

* commit '58723a9':
DMPlex: Fix DMPlexCheckSkeleton() for submeshes
DMPlex ex5: Added second fault for quad test
DMPlex: Fix setting

Merge commit '58723a9' into knepley/feature-plex-hybrid-3d

* commit '58723a9':
DMPlex: Fix DMPlexCheckSkeleton() for submeshes
DMPlex ex5: Added second fault for quad test
DMPlex: Fix setting of hybrid bounds for multiple faults
DMPlex ex5: Turn on second fault in quad test
DMPlex: Removed compiler warnings
DMPlex ex5: Cleaned up 2D tensor tests to rely on interpolation
DMPlex ex5: Cleaned up 2D simplex tests to rely on interpolation

Conflicts:
include/petscdmplex.h
src/dm/impls/plex/examples/tests/ex4.c
src/dm/impls/plex/plexrefine.c

show more ...

Sys: add PetscMalloc1 macro, array allocation without redundant types

The type is inferred from the pointer return type. This patch is
automated via the following script:

git grep -l 'PetscMalloc(

Sys: add PetscMalloc1 macro, array allocation without redundant types

The type is inferred from the pointer return type. This patch is
automated via the following script:

git grep -l 'PetscMalloc(.*sizeof' src | xargs perl -pi -e 's@PetscMalloc\(([^,;]*[^,; ]) *\* *sizeof\([^,;()]+\),@PetscMalloc1($1,@'

This commit contains an additional bug-fix in csrperm.c, fixing pointer
arity. The code was introduced in 2006, but the allocation could not
have been correct at any time. This probably means that
MatDuplicate_SeqAIJPERM has never been tested.

a54129beb540034ba105796c682d589e7e1111f2
Richard Tran Mills <rmills@ornl.gov>

Added MATSEQCSRPERM support for MatDuplicate() and conversion to/from
MATSEQAIJ. Note that these changes are not quite debugged.

show more ...