Changes between Version 2 and Version 3 of Ticket #28707, comment 1
- Timestamp:
- 11/09/19 09:31:37 (3 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
Ticket #28707, comment 1
v2 v3 3 3 [[Image(test_geod_bad.png)]] 4 4 5 The geodesic sarises from the upper right (''(x,y)'' approx. (10, 5)), circles twice around the black hole (the grey area) and then departs away to the lower left (`rkf45` "exact" solution) or to the lower right (`ode_int` bad solution).5 The geodesic arises from the upper right (''(x,y)'' approx. (10, 5)), circles twice around the black hole (the grey area) and then departs away to the lower left (`rkf45` "exact" solution) or to the lower right (`ode_int` bad solution). 6 6 This issue is due to the default values of 1.49012e-8 for the parameters `rtol` and `atol` of 7 7 [https://docs.scipy.org/doc/scipy/reference/generated/scipy.integrate.odeint.html scipy.integrate.odeint] being not sufficient in the present case. Now, with the Sage 8.9 implementation of `IntegratedCurve.solve()`, there is no way to change the values of `rtol` and `atol`. Thanks to the new argument `**control_param` introduced here this becomes possible. Morevoer, the default values of