Changes between Version 1 and Version 2 of Ticket #28707, comment 1
 Timestamp:
 11/09/19 09:29:30 (3 years ago)
Legend:
 Unmodified
 Added
 Removed
 Modified

Ticket #28707, comment 1
v1 v2 1 In Sage 8.9, when integrating a null geodesic with an impact parameter close to critical (all details are in [https://nbviewer.jupyter.org/github/egourgoulhon/SageMathTest/blob/master/Notebooks/test_geodesics.ipynb this notebook]), the solution obtained with `solve(method='ode_int')` (green in the figure below) differs significantly from those obtained by means of other solvers, for instance the `rkf45` solver (red in the figure below):1 In Sage 8.9, when integrating a null geodesic in Schwarzschild spacetime with an impact parameter close to critical (all details are in [https://nbviewer.jupyter.org/github/egourgoulhon/SageMathTest/blob/master/Notebooks/test_geodesics.ipynb this notebook]), the solution obtained with `solve(method='ode_int')` (green in the figure below) differs significantly from those obtained by means of other solvers, for instance the `rkf45` solver (red in the figure below): 2 2 3 3 [[Image(test_geod_bad.png)]]