Problem with integrated_geodesic on open subsets on Manifolds. U.integrated_geodesic(g, (s, 0, 1500), v0, chart=S)

This isn't working. As far as I have got it should produce the same as the Sagemanifolds Schwarzschild example. The only difference is I have introduced open subsets for the manifold otherwise it can't find a basis. I want different charts for future development. It shouldn't produce the error as far as I can see.

I'd really appreciate some help.

I have just found that if I take out the restriction for the phi range then the geodesic computes and is the same as the example. As it was it only ever computes one orbit. I can't see why this should be so I think it is a bug.

edit retag close merge delete

The link https://cocalc.com/projects/0a54691c-2a01-4d4d-8f94-d0785ae70fa0/files/KM_simple_geod_Kerr.ipynb is not public. Can you fix this? (click on Files, check the file, click on Share and then Make item public).

( 2018-03-29 13:43:55 +0200 )edit

... and then replace the link here with the share url that CoCalc will provide.

( 2018-04-06 09:13:11 +0200 )edit

You could also post the code here in your question!

( 2018-04-06 09:13:42 +0200 )edit

To display inline code, like z = x*y, use backticks.

To display blocks of code or error messages, skip a line above and below, and do one of the following (all give the same result):

• indent all code lines with 4 spaces
• select all code lines and click the "code" button (the icon with '101 010')
• select all code lines and hit ctrl-K

For instance, typing

If we define f by

def f(x, y):
return (x, y)

then f(2, 3) returns (2, 3) but f(2) gives:

TypeError: f() takes exactly 2 arguments (1 given)


produces:

If we define f by

def f(x, y):
return (x, y)


then f(2, 3) returns (2, 3) but f(2) gives:

TypeError: f() takes exactly 2 arguments (1 given)

( 2018-04-06 09:14:09 +0200 )edit