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

asked 2018-03-28 09:56:51 -0500

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.

link:[text]://cocalc.com/projects/0a54691c-2a01-4d4d-8f94-d0785ae70fa0/files/KM_simple_geod_Kerr.ipynb?session=

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 flag offensive close merge delete

Comments

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).

eric_g gravatar imageeric_g ( 2018-03-29 06:43:55 -0500 )edit

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

slelievre gravatar imageslelievre ( 2018-04-06 02:13:11 -0500 )edit

You could also post the code here in your question!

slelievre gravatar imageslelievre ( 2018-04-06 02:13:42 -0500 )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)
slelievre gravatar imageslelievre ( 2018-04-06 02:14:09 -0500 )edit