when MPIU_REAL or MPIU_SCALAR is used then MPIU_SUM/MAX/MIN must be used

these bugs are why plex as crashing in quad precision

fixed use of MPI_SUM with MPIU_REAL which fails for quad precision in plexfem.c

Fixed compiler warnings

Plex: Fix compiler warning

Use %D and cast PetscReal to (double) for PetscPrintf() so it works with quad precision and 64 bit indices

Viewer: Now ViewerRead() takes an extra argument to return the number of items read
- An error is thrown for insufficient data unless the extra argument is non-NULL

Merge branch 'knepley/fix-dmload-parallel'

* knepley/fix-dmload-parallel:
Plex: Fixed input/output of vertices - I was using a local vector, which is completely wrong - pEnd needed to be fixed for

Merge branch 'knepley/fix-dmload-parallel'

* knepley/fix-dmload-parallel:
Plex: Fixed input/output of vertices - I was using a local vector, which is completely wrong - pEnd needed to be fixed for other procs on input
Plex ex5: Added a tutorial for loading and saving meshes
Builder: Correct -numProcs handling
Plex ex3: Run dsymutil on the executable

show more ...

Viewer: Fixes for new Read interface
- Fixed example
- Fixed Fortran 90 binding
- Fixed type

Merge branch 'knepley/fix-dmload-parallel' into knepley/feature-viewer-partial-read

* knepley/fix-dmload-parallel:
Plex: Fixed input/output of vertices - I was using a local vector, which is compl

Merge branch 'knepley/fix-dmload-parallel' into knepley/feature-viewer-partial-read

* knepley/fix-dmload-parallel:
Plex: Fixed input/output of vertices - I was using a local vector, which is completely wrong - pEnd needed to be fixed for other procs on input
Plex ex5: Added a tutorial for loading and saving meshes
Builder: Correct -numProcs handling
Plex ex3: Run dsymutil on the executable

# Please enter a commit message to explain why this merge is necessary,
# especially if it merges an updated upstream into a topic branch.
#
# Lines starting with '#' will be ignored, and an empty message aborts
# the commit.

show more ...

Merge branch 'master' into knepley/feature-viewer-partial-read

* master: (48 commits)
ViennaCL: Fixes memory leak when destructing an aijviennacl matrix.
MatLoad_MPIDense() always allocated arra

Merge branch 'master' into knepley/feature-viewer-partial-read

* master: (48 commits)
ViennaCL: Fixes memory leak when destructing an aijviennacl matrix.
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: Reorganized PlexPoint interface - Added DMPlexGetPointLocalField() and DMPlexGetPointGlobalField() - Now the DMPlexGetLocalOffset*_Private() functions return start and end - Uniformly use private functions
Plex: Put in explicit cast
MatSetValuesMPISBAIJ will now print correct global indices when error generated due to new nonzero location
...

# Please enter a commit message to explain why this merge is necessary,
# especially if it merges an updated upstream into a topic branch.
#
# Lines starting with '#' will be ignored, and an empty message aborts
# the commit.

show more ...

Plex: Fixed input/output of vertices
- I was using a local vector, which is completely wrong
- pEnd needed to be fixed for other procs on input

Merge branch 'knepley/fix-plex-offset-calc'

* knepley/fix-plex-offset-calc:
Plex: Reorganized PlexPoint interface - Added DMPlexGetPointLocalField() and DMPlexGetPointGlobalField() - Now the DMPle

Merge branch 'knepley/fix-plex-offset-calc'

* knepley/fix-plex-offset-calc:
Plex: Reorganized PlexPoint interface - Added DMPlexGetPointLocalField() and DMPlexGetPointGlobalField() - Now the DMPlexGetLocalOffset*_Private() functions return start and end - Uniformly use private functions

show more ...

Merge branch 'knepley/fix-plex-geom'

* knepley/fix-plex-geom:
Plex: Put in explicit cast
Plex: Fix geometry struct handling We are treating the structs as arrays of Scalar for storage purposes,

Merge branch 'knepley/fix-plex-geom'

* knepley/fix-plex-geom:
Plex: Put in explicit cast
Plex: Fix geometry struct handling We are treating the structs as arrays of Scalar for storage purposes, so we need to round up in the division

show more ...

Plex: Reorganized PlexPoint interface
- Added DMPlexGetPointLocalField() and DMPlexGetPointGlobalField()
- Now the DMPlexGetLocalOffset*_Private() functions return start and end
- Uniformly use priva

Plex: Reorganized PlexPoint interface
- Added DMPlexGetPointLocalField() and DMPlexGetPointGlobalField()
- Now the DMPlexGetLocalOffset*_Private() functions return start and end
- Uniformly use private functions

show more ...

Merge branch 'knepley/fix-log2'

* knepley/fix-log2:
Plex ex3: Update to PetscLog2Real()
Math: Added test and stub for log2() - FUCK YOU MICROSOFT

Plex: Put in explicit cast

Viewer: PetscViewerRead() should return the number of items actually read
- This also applies to all its implementations

Plex ex3: For want of a semicolon, the kingdom was lost

Plex: Fix geometry struct handling
We are treating the structs as arrays of Scalar for storage purposes, so we need to round up in the division

Plex ex3: Update to PetscLog2Real()

Plex ex5: Added a tutorial for loading and saving meshes

Plex ex3: Run dsymutil on the executable

Merge branch 'knepley/feature-harmonic-advection'

* knepley/feature-harmonic-advection: (86 commits)
SNES: Some docs for NASM
TS: Added view for MIMEX
Plex: Silence compiler warnings
Builder

Merge branch 'knepley/feature-harmonic-advection'

* knepley/feature-harmonic-advection: (86 commits)
SNES: Some docs for NASM
TS: Added view for MIMEX
Plex: Silence compiler warnings
Builder: Check for missing numProcs was wrong
PetscFE: Turn off Fortran binding for PetscFECompositeGetMapping()
Plex: Fix big bug in DMPlexGetGlobalFieldOffset_Internal() - If offset is negative (dof is on another proc), then we must have end < start
Builder: Allow override of numProcs
Plex: Fix Fortran wrappers for DMLabelView() and DMPlexDistribute()
Plex+Label: Added Fortran wrapper for DMPlexCreateLabel()
Plex: Dammit, how did this not get checked in?
Plex: Refiner check for hybrid cells must use hybrid faces as well - Ghost cells in FV method are hybrid
TS: MIMEX now uses work vectors for correct grid
Plex: Adding support for FV in automatic interpolator construction - Notice that the injector cannot do the right thing until we scrap VecScatter
Plex+FEM: Using PetscFVGetDualSpace()
SNES+Doc: Added Options Database keys
PetscFV: Added PetscFVRefine() - This allows interpolation to finer grids needed for FAS
PetscFV: Added a dual space for calculating projections
DM: Better error checking in DMCoarsen()
DMDT+PetscFE: Added PetscQuadratureExpandComposite() - Replaces PetscFECompositeExpandQuadrature() - Added PetscFECompositeGetMapping()
Plex: Fixed complex type error
...

Conflicts:
src/dm/impls/plex/ftn-custom/makefile
src/ts/interface/ts.c

show more ...

Plex: Silence compiler warnings

Plex: Fix Fortran wrappers for DMLabelView() and DMPlexDistribute()