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.