[Users] Error Estimation for RK45CK

Erik Schnetter schnetter at gmail.com
Wed Jul 26 15:10:10 CDT 2023


David

This happens in the file "RK45.c" in line 252. The loop there accesses
the "gamma_array" with the error estimation coefficients and then
calls "MoL_LinearCombination" to estimate the error.

-erik

On Wed, Jul 26, 2023 at 3:02 PM Garrison, David <Garrison at uhcl.edu> wrote:
>
> Hello,
>
> Where in the MoL is the error estimates for RK45 calculated?  I can’t seem to find that in the source code.  Somehow they are written to the MOL::ERRORESTIMATE grid function but I can’t see where that is done.
>
> -DG
>
> --
> David Garrison, Ph.D.
> Associate Dean for the College of Science and Engineering,
> Professor and Founding Chair of Physics
> University of Houston-Clear Lake
> Bayou 3611
> Houston, TX 77058
>
> Phone: 281-283-3796
> https://www.uhcl.edu/science-engineering/faculty/garrison-david
> http://www.uhcl.edu/physics
>
> "If we knew what it was we were doing, it would not be called research, would it?" — Albert Einstein.
>
>
>
>
> _______________________________________________
> Users mailing list
> Users at einsteintoolkit.org
> http://lists.einsteintoolkit.org/mailman/listinfo/users



-- 
Erik Schnetter <schnetter at gmail.com>
http://www.perimeterinstitute.ca/personal/eschnetter/


More information about the Users mailing list