Changes between Version 2 and Version 3 of Ticket #28707, comment 1


Ignore:
Timestamp:
11/09/19 09:31:37 (23 months ago)
Author:
egourgoulhon
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • Ticket #28707, comment 1

    v2 v3  
    33[[Image(test_geod_bad.png)]]
    44
    5 The geodesics 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).
     5The 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).
    66This issue is due to the default values of 1.49012e-8 for the parameters `rtol` and `atol` of
    77[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