Re: Interpolator coefficients for Radau method in solve_ivp
Travis Czechorski <[email protected]> Wed, 29 Apr 2020 17:11:53 -0400
| Newsgroups | gmane.comp.python.scientific.user |
|---|---|
| Message-ID | <CAKCV6PdziZv=fw5o-RyJbdKPe_7mank8O9H4AnV1OmO7H_NqLg@mail.gmail.com> |
Hello, I am trying to figure out how to calculate the following coefficients in the source code for the Radau integration method. # 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, On Wed, Apr 29, 2020 at 5:10 PM Travis Czechorski <[email protected]> wrote: > Hello, > > I am trying to figure out how to calculate the following coefficients in > the source code for the Radau integration method. > > # 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]]) > _______________________________________________ SciPy-User mailing list [email protected] https://mail.python.org/mailman/listinfo/scipy-user