Evaluating the tangent vector of a previously saved and loaded integrated geodesic

edit

Hello Bubbz, thank you for your question. It seems indeed that the contents of the attributes _solutions and _interpolations of geod are lost in the process save/load, although all the other attributes are well preserved (e.g. _initial_tangent_vector or _parameters).

In the case of the attribute _interpolations, this seems to be expected because pickling does not work on Spline objects: calling e.g. save(geod._interpolations, "geod_interps.sobj") yields the error message TypeError: cannot pickle 'sage.calculus.interpolation.Spline' object.

In the case of the attribute _solutions, the reason why it is lost is probably tricky, as this is easily worked around with the idea above: you could save and reload the attribute _solutions independently from geod itself, by executing e.g. save(geod._solutions, "geod_sols.sobj"), so that you could execute geod._solutions = load("geod_sols.sobj") after having reloaded geod itself at any later time.