Ask Your Question
1

christoffel symbols of the second kind: how exactly to calculate?

asked 2023-05-25 16:29:18 +0200

pfeifhns gravatar image

Hello, I have a function f(u,v) = ( cos(v)sin(u), sin(u)sin(v), cos(u) ) which describes a sphere in R3.

How can I calculate the Christoffel symbols of the second kind? I didn't understand the examples I saw so far. I have no idea on how to proceed. Thank you for your efforts, Best wishes Hans

edit retag flag offensive close merge delete

Comments

1 Answer

Sort by ยป oldest newest most voted
1

answered 2023-05-27 11:37:30 +0200

achrzesz gravatar image

updated 2023-05-27 12:29:18 +0200

From: sage/src/sage/calculus/test.py

def christoffel(i,j,k,vars,g):
     s = 0
     ginv = g^(-1)
     for l in range(g.nrows()):
         s = s + (1/2)*ginv[k,l]*(g[j,l].diff(vars[i])+g[i,l].diff(vars[j])-\
                                  g[i,j].diff(vars[l]))
     return s

var('th ph')
g = matrix(SR, [[1,0],[0,sin(th)^2]])
for i in [0,1]:
    for j in [0,1]:
        for k in[0,1]:
            print((i,j,k),christoffel(i,j,k, [th,ph], g))

(0, 0, 0) 0
(0, 0, 1) 0
(0, 1, 0) 0
(0, 1, 1) cos(th)/sin(th)
(1, 0, 0) 0
(1, 0, 1) cos(th)/sin(th)
(1, 1, 0) -cos(th)*sin(th)
(1, 1, 1) 0

SageManifolds version:

S2=manifolds.Sphere(2)     
Phi=S2.embedding()                      
sph.<th,ph>=S2.spherical_coordinates()  
E=S2.ambient()                          
g=E.metric() 
g1 = S2.metric('g1')
g1.set( Phi.pullback(g))
nab = g1.connection()
nab.display(coordinate_labels=False,only_nonzero=False)
edit flag offensive delete link more

Comments

Explanations in the comment above

achrzesz gravatar imageachrzesz ( 2023-05-27 11:41:06 +0200 )edit

Thank you very very much for your efforts1 Best wishes Hans

pfeifhns gravatar imagepfeifhns ( 2023-05-29 13:13:05 +0200 )edit

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools

Stats

Asked: 2023-05-25 16:29:18 +0200

Seen: 458 times

Last updated: May 29 '23