TSEvent: Add option -ts_event_monitor to monitor the event detection algorithm.

TSEvent: Fix TSEvent when an event occurs at the last time-step.

TSEvent used to not work correctly when the event was detected at the
final TS time step. This was because TS set TS_CONVERGED_TIME f

TSEvent: Fix TSEvent when an event occurs at the last time-step.

TSEvent used to not work correctly when the event was detected at the
final TS time step. This was because TS set TS_CONVERGED_TIME flag
causing the time-loop to exit inspite of the event being detected.
A reproducible example for this scenario (pointed by Hong) was
./ex40 -ts_final_time 4.1 -ts_exact_final_time MATCHSTEP. An event occured
at the last step that was not getting correctly processed by TSEventMonitor.

Fixed this issue by setting ts->reason to TS_CONVERGED_ITERATING whenever
a step is rollbacked in TSEventMonitor.

show more ...

TSEvent: Modified TSEvent event incidence condition.

Previously, for an event to occur two conditions had to be satisfied:
(a) a change in function sign between consecutive time-steps, AND
(b) the e

TSEvent: Modified TSEvent event incidence condition.

Previously, for an event to occur two conditions had to be satisfied:
(a) a change in function sign between consecutive time-steps, AND
(b) the event function value within the specified tolerance.

The new implementation only uses condition (b), i.e., an event
is detected if its function value is within the specified tolerance.
Doing so, eliminates cases where TS used to take ridiculously small
time steps (~ 1e-15) just because the function value had not changed
sign.

An additional check had to be added to the event location code
so that the same event is not detected again.

show more ...

Merge branch 'master' into knepley/feature-composable-callbacks

* master: (151 commits)
PetscMalloc: Fix handling of zero-sized allocations
add note to DMGetLocal/GlobalVector() about intended s

Merge branch 'master' into knepley/feature-composable-callbacks

* master: (151 commits)
PetscMalloc: Fix handling of zero-sized allocations
add note to DMGetLocal/GlobalVector() about intended short-time usage of the vectors.
sf: use patched openmpi tarball that has a potential fix for sf example hang
SNES problems that set VI bounds must use a VI solver since the other solvers ignore the bounds
pounders viewer did not cast PetscReal to double in printf nor use %D instead of %d
fix ts ex20adj
fix ts theta
fixed another variable declaration in the middle of a code block
example main() had PetscInt return type and argument that should be int
DMDAVTKWriteAll_VTS: Fix handling of higher-dimension coordinates
fixed example that declared variables in the middle of a function
C programs do not return PetscInt nor do they take an argument that is PetscInt, fixed example that did
You cannot declare variables in the middle of functions in PETSc.
add ierr checking to PetscOptionsBegin/End() in example
configure: add error check if user specifies PETSC_ARCH=-dbg [i.e starts with a '-']
scale norm of vector difference by norm of vector to insure it is small enough to consider the two vectors equal
missing PetscViewerDestroy() added
MatNestGetISs: use custom fortranstub to handle PETSC_NULL_OBJECT
update location of open-mx website on FAQ page
If the user sets bounds and calls a solver that does not support them, throw an error.
...

show more ...

Merge branch 'pr192/Fande-Kong/scalable-matincreaseoverlap/master' of https://bitbucket.org/petsc/petsc

Merged petsc/petsc into master

Merge branch 'master' into stefano_zampini/pcbddc-adaptive-selection

Conflicts:
src/ksp/pc/impls/bddc/bddcprivate.c
src/ksp/pc/impls/bddc/bddcscalingbasic.c
src/mat/impls/aij/mpi/mumps/mumps.c
s

Merge branch 'master' into stefano_zampini/pcbddc-adaptive-selection

Conflicts:
src/ksp/pc/impls/bddc/bddcprivate.c
src/ksp/pc/impls/bddc/bddcscalingbasic.c
src/mat/impls/aij/mpi/mumps/mumps.c
src/mat/impls/is/matis.h

show more ...

Merged master into semihozmen/ksp-fcg-eigen-comp-support

Merged petsc/petsc into master

Merge branch 'master' into pr192/Fande-Kong/scalable-matincreaseoverlap/master

Needed to merge master to update to new calling sequence PetscCommBuildTwoSided

fixed missing CHKERRQ() detected by clang's --analyze

Merged petsc/petsc into master

Merge branch 'barry/fix-directory-layout'

updated PETSc directory layout to match standard packaging strategies

include/petsc finclude,private,mpiuni
lib/petsc conf
bin/petsc*

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

Merge branch 'master' into tisaac/dmforest

If we're going to start developing on this again, the recent changes to the
location of configuration files will make it annoying to switch to and from
thi

Merge branch 'master' into tisaac/dmforest

If we're going to start developing on this again, the recent changes to the
location of configuration files will make it annoying to switch to and from
this branch. Matt's and my fixes in DMPlex will also be relevant.

* master: (1557 commits)
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 ex3: For want of a semicolon, the kingdom was lost
MatSetValuesMPISBAIJ will now print correct global indices when error generated due to new nonzero location
completed basic use of MatSetValuesBlocked_SeqBAIJ_Inlined() by MatSetValuesBlocked_MPIBAIJ() to allow proper global indices being printed when error generated by new nonzero locations.
Added inline version of MatSetValuesBlocked_SeqBAIJ() to be used by MatSetValuesBlocked_MPIBAIJ() to print error message for nonzero block with global indices
...

show more ...

Merge branch 'shri/adjoints-and-events'

Merge branch 'master' into shri/ts-is-for-differential-variables

Conflicts:
include/petsc-private/tsimpl.h
src/ts/impls/eimex/eimex.c
src/ts/interface/ts.c

ADJOINTS+EVENTS: Should not use PetscAbsScalar on real numbers.

Merge branch 'shri/adjoints-and-events'

ADJOINTS+EVENTS: Initialize TSEvent at the beginning of each TSSolve.

ADJOINTS+EVENTS: The cost function and the sensitivity matches for the example with events (ex3adj_events) compared with the original adjoint example (ex3adj)

i) Added TSAdjointEventMonitor that si

ADJOINTS+EVENTS: The cost function and the sensitivity matches for the example with events (ex3adj_events) compared with the original adjoint example (ex3adj)

i) Added TSAdjointEventMonitor that simply checks if an event had occured in the forward solve at the current time. If true, calls the user postevent routine with the forwardsolve flag set to FALSE.
ii) Added stepnum field in event recorder
iii) Fixed the time interval of the disturbance in ex3adj.
iv) PostStepFunction in ex3adj_events (also added to ex3adj) displays value of delta after every time-step.
v) PostEventFunction in ex3adj_events now handles both the forward and the adjoint solve.
vi) Have commented out TSPostStep() from TSAdjointSolve as we don't need it in the adjoint solve (I think).

show more ...

Merge branch 'master' into jed/chem
Note merging a branch with good stuff in it for a year and a half is insane!

Conflicts:
src/ts/examples/tutorials/makefile
src/ts/interface/ts.c

ADJOINTS+EVENTS: Preliminary implementation of event recording.

ADJOINTS+EVENTS: Bug fixes in TSEventMonitor.

i) Avoid taking very small time-step after PostEvent.
ii) TS step counters need to be decremented on rolling back step.