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.


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


The link 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)


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