id,summary,reporter,owner,description,type,status,priority,milestone,component,resolution,keywords,cc,merged,author,reviewer,upstream,work_issues,branch,commit,dependencies,stopgaps
25936,Improvements to integrated curves and geodesics on manifolds,gh-FlorentinJ,,"This ticket implements some improvements regarding integrated curves on manifolds, first appearing in #22951.
It adds `scipy.integrate.ode` as an integration solver. The main advantage is that it supports basic event handling, making it possible to detect frontiers of charts. This allows for integration on multiple charts, as shown in [http://nbviewer.jupyter.org/github/FlorentinJ/Sage_examples/blob/master/solve_across_charts%20example.ipynb this] notebook, which can be particularly useful to avoid singularities, or in manifold which cannot be represented by a single chart, like in [http://nbviewer.jupyter.org/github/FlorentinJ/Sage_examples/blob/master/solve_across_charts%20example%202.ipynb this one] (not commented).
A side effect of this new method is that the time needed for the integration can be reduced (a lot!) by using objects created with `fast_callable` instead of symbolic expressions. This is used in [http://nbviewer.jupyter.org/github/FlorentinJ/Sage_examples/blob/master/kerr_geodesic_minimal.ipynb this] notebook to compute the shape of a geodesic in a Kerr spacetime. Speed test shows a x700 speedup compared to the default integration method ""rk4_maxima"".
",enhancement,closed,major,sage-8.5,geometry,fixed,"manifold, geodesic, integrated curve, integration, chart",egourgoulhon karimvanaelst,,Florentin Jaffredo,"Travis Scrimshaw, Eric Gourgoulhon",N/A,,61061c152b58324f182dccc2d89a35d8fe045709,61061c152b58324f182dccc2d89a35d8fe045709,,