Updated TSFunctionDomainError signature

Instead of getting all stages and the stage index, the function gets
directly the current state (e.g. the one to check).

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

* knepley/feature-viewer-partial-read:
Viewer: Fixed initialization
TS: Updated example to new interface
Plex: Fix compiler warning
Viewer:

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

* knepley/feature-viewer-partial-read:
Viewer: Fixed initialization
TS: Updated example to new interface
Plex: Fix compiler warning
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
Viewer: Fixes for new Read interface - Fixed example - Fixed Fortran 90 binding - Fixed type
Builder: Allow a local dictionary for run parameters
Viewer: PetscViewerRead() should return the number of items actually read - This also applies to all its implementations

show more ...

Added domain error checks at each stage

This patch adds two functions:
1 - TSSetFunctionDomainError
2 - TSFunctionDomainError

The first one to set a user function, the second one to call it.
It als

Added domain error checks at each stage

This patch adds two functions:
1 - TSSetFunctionDomainError
2 - TSFunctionDomainError

The first one to set a user function, the second one to call it.
It also adds the calls to TSFunctionDomainError in all implementations calling TSPostStage.

show more ...

Fixed issues in the review: fixed TABs; C89 syntax; TSClone comm; error message; clarified implicit equations in docs

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

ARKIMEX implementation and TSClone notes

added some documetation/fixed some arkimex typos

Cleaned up some issues in TSClone

fixed typo in comment

Replace TSDuplicate with TSClone; all is working but there is still some concern with SNES when the clone is destroyed and a GetSNES and SetSNES is needed

TSComputeRHSJacobian: If TSComputeRHSFunctionLinear is being used, then re-computing the RHS Jacobian (when called from TSComputeRHSFunctionLinear) is necessary, because the saved Jacobian is actuall

TSComputeRHSJacobian: If TSComputeRHSFunctionLinear is being used, then re-computing the RHS Jacobian (when called from TSComputeRHSFunctionLinear) is necessary, because the saved Jacobian is actually scaled-shifted [shift*I - Jac].

show more ...

finalized TSDuplicate implementation and testing

TSDuplicate: Added an initial implementation of TSDuplicate()

added a TS statement to allow saving the last solution to a file

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

TSReset: Call TSAdaptReset()

TSErrorWeightedNorm[2|Infinity]: API changes

* Accept two arguments U,Y instead of always using ts->vec_sol
* Accept argument wnormtype intead of use ts->adapt->wnormtype

Merge branch 'dalcinl/feature-draw'

Merge branch 'barry/feature-bargraph'

missing CHKERRQ; and ugly double declaration of ierr

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

remove redundant brackets

enable users to evaluate cost integral in the forward run

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

TSMonitorLGTimeStep: Don't call PetscDrawLGSetUseMarkers

Otherwise command line option -lg_use_markers false has no effect