Update-TSEvent: Use higher event location tolerance for single precision.

Update-TSEvent: Use PetscRealPart() for supporting complex builds.

Update-TSEvent: Removed the limit on events recorded.

Removed the limit on events recorded. The initial array sizes in the event recorder
is set to 12 (changed by option -ts_event_recorder_initial_s

Update-TSEvent: Removed the limit on events recorded.

Removed the limit on events recorded. The initial array sizes in the event recorder
is set to 12 (changed by option -ts_event_recorder_initial_size) Each time limit is
reached, the array sizes are doubled (see TSEventRecorderResize)

show more ...

Update-TSEvent: Major overhaul of the event location logic.

Modified the code to use the Anderson-Bjorck variant of the regular falsi method for event location. Preliminary tests show this variant t

Update-TSEvent: Major overhaul of the event location logic.

Modified the code to use the Anderson-Bjorck variant of the regular falsi method for event location. Preliminary tests show this variant takes less iterations than pure regular falsi (previously implemented)

show more ...

Update-TSEvent: Added counter to monitor number of iterations.

Added iteration counter in TSEvent.
Updated the event monitoring messages.

TSEvent: Remove redundant initialization

TSEvent: Refactor TSPostEvent() and make it private

TSEvent: Fix time step resetting to play well with adaptivity

* After finishing with event handling, the time step has to be reset to
some previous value. Instead of using the initial time step, u

TSEvent: Fix time step resetting to play well with adaptivity

* After finishing with event handling, the time step has to be reset to
some previous value. Instead of using the initial time step, use the
time step that was current before entering the event handling business.

show more ...

TSAdapt: Minor additional cleanups (mostly coding conventions)

TSEvent: Rename TSSetEventMonitor() -> TSSetEventHandler()

TSEvent: Various fixes and cleanups

* Rename most of TSEventMonitorXXX() to TSEvent[Handler]XXX
* Rename the private TSEvent structure as well as some slots within it

Introduce MPIU_Allreduce() to detect cases where all processes that share an MPI comm do not call the same sequence of MPI_Allreduce()
This is a bug that can be difficult to track down.

Suggested-by

Introduce MPIU_Allreduce() to detect cases where all processes that share an MPI comm do not call the same sequence of MPI_Allreduce()
This is a bug that can be difficult to track down.

Suggested-by: Eric Chamberland <Eric.Chamberland@giref.ulaval.ca>
Suggested-by: Patrick Lacasse <placasse@giref.ulaval.ca>

show more ...

fix violations of PETSc style guide: Usage of SETERRQ and NULL

TSEvent: PETSc blame digest (next) 2015-05-19 fix => PetscOptionsBegin needs CHKERRQ(ierr)

TSEvent: PETSc blame digest (next) 2015-05-19 fix => Cast to TSEventStatus

TSEvent: Major fix to the TSEvent code to handle multiple events located within an interval.

The previous TSEvent code did not properly handle cases where multiple events were located
within an inte

TSEvent: Major fix to the TSEvent code to handle multiple events located within an interval.

The previous TSEvent code did not properly handle cases where multiple events were located
within an interval. For instance, if two events were located in an interval with the first one
expected to cross zero at an earlier time compared to the second one, TSEvent would incorrectly
detect the second one if the second event function value was within the tolerance earlier. To avoid
this issue, the logic for TSEvent has been reorganized now to locate other events even if a single (or multiple)
events have been detected.

show more ...

TSEvent: Added TSSetEventTolerances() - Support for setting tolerances for each individual event.

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