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


Ignore:
Timestamp:
11/09/19 09:29:30 (2 years ago)
Author:
egourgoulhon
Comment:

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):
     1In 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):
    22
    33[[Image(test_geod_bad.png)]]