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@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@einsteintoolkit.org http://lists.einsteintoolkit.org/mailman/listinfo/users