[SciPy-User] Interpolator coefficients for Radau method in solve_ivp

classic Classic list List threaded Threaded
1 message Options
Reply | Threaded
Open this post in threaded view
|

[SciPy-User] Interpolator coefficients for Radau method in solve_ivp

Travis Czechorski
This post was updated on .
Hello,

I am trying to figure out how to calculate the following coefficients in the source code for the Radau integration method. The source code is here https://github.com/scipy/scipy/blob/v1.4.1/scipy/integrate/_ivp/radau.py#L180-L536

# Interpolator coefficients.
P = np.array([ [13/3 + 7*S6/3, -23/3 - 22*S6/3, 10/3 + 5 * S6],
                     [13/3 - 7*S6/3, -23/3 + 22*S6/3, 10/3 - 5 * S6],
                     [1/3, -8/3, 10/3]])  

I read that it uses a cubic polynomial to interpolate, can someone provide me with the process to calculate these coefficients?

Thanks,
___________
SciPy-User mailing list
SciPy-User@python.org
https://mail.python.org/mailman/listinfo/scipy-user