Hello all:
So after some more thought and discussion, here's another shot at
where to go with the AdaptiveTimeSolver classes:
- a previously planned reformulation of AdaptiveTimeSolver to a virtual
base class and a new child called DoublingAdaptiveTimeSolver to have the
same functionality as the current ATS.
- develop a predictor/corrector version called PCAdaptiveTimeSolver that
is 1) core agnostic and 2) determines the LTE based upon the difference
between the predictor/corrector solutions AND their respective orders,
for example:
- the predictor is a first order implicit Euler solver.
- the corrector is a second order trapezoidal solver, which is
initialized with
the solution from the predictor.
- the error is determined directly from the difference between the two
solutions, local extrapolation. In the case where both
solutions are of
the same order, error determination is non-trivial.
--
Nasser Mohieddin Abukhdeir
Graduate Student (Materials Modeling Research Group)
McGill University - Department of Chemical Engineering
http://webpages.mcgill.ca/students/nabukh/web/
http://mmrg.chemeng.mcgill.ca/
-------------------------------------------------------------------------
This SF.Net email is sponsored by the Moblin Your Move Developer's challenge
Build the coolest Linux based applications with Moblin SDK & win great prizes
Grand prize is a trip for two to an Open Source event anywhere in the world
http://moblin-contest.org/redirect.php?banner_id=100&url=/
_______________________________________________
Libmesh-devel mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-devel