|
tensor[geodesic_eqns] - generate the Euler-Lagrange equations for the geodesic curves
|
|
Calling Sequence
|
|
geodesic_eqns(coord, param, Cf2)
|
|
|
Parameters
|
|
|
coord
|
-
|
list of coordinate names
|
|
param
|
-
|
name of the variable to parametrize the curves with
|
|
Cf2
|
-
|
Christoffel symbols of the second kind
|
|
|
|
|
|
Description
|
|
|
•
|
The function geodesic_eqns(coord, Tau, Cf2) generates (but does not solve) the Euler-Lagrange equations of the geodesics for a metric with Christoffel symbols of the second kind Cf2 and coordinate variables coord. The equations are written in terms of the coordinate variable names as functions of the given parameter Tau. They are returned in the format of a list of equations.
|
|
•
|
Cf2 should be indexed using the cf2 indexing function provided by the tensor package. It can be computed using the Christoffel2 routine.
|
|
|
|
Examples
|
|
|
>
|
|
Determine the geodesic equations for the Poincare half-plane. The coordinates are:
|
>
|
![coord := [u, v]](/support/helpjp/helpview.aspx?si=5520/file04532/math90.png)
|
|
![coord := [u, v]](/support/helpjp/helpview.aspx?si=5520/file04532/math93.png)
| (1) |
The metric is:
|
>
|
![g_compts := array(symmetric, sparse, 1 .. 2, 1 .. 2, [(1, 1) = 1/v^2, (2, 2) = 1/v^2])](/support/helpjp/helpview.aspx?si=5520/file04532/math99.png)
|
|
>
|
![g := create([-1, -1], eval(g_compts))](/support/helpjp/helpview.aspx?si=5520/file04532/math103.png)
|
|
![g := TABLE([index_char = [-1, -1], compts = matrix([[1/v^2, 0], [0, 1/v^2]])])](/support/helpjp/helpview.aspx?si=5520/file04532/math106.png)
| (2) |
|
>
|
|
|
>
|
 
|
|
>
|
|
|
>
|
|
|
>
|
|
|
| (3) |
Now generate the geodesic equations:
|
>
|
|
|
| (4) |
How about Euclidean 3-space in Cartesian coordinates?
|
>
|
![coord := [x, y, z]](/support/helpjp/helpview.aspx?si=5520/file04532/math160.png)
|
|
![coord := [x, y, z]](/support/helpjp/helpview.aspx?si=5520/file04532/math163.png)
| (5) |
|
>
|
![g_compts := array(symmetric, sparse, 1 .. 3, 1 .. 3, [(1, 1) = 1, (2, 2) = 1, (3, 3) = 1])](/support/helpjp/helpview.aspx?si=5520/file04532/math167.png)
|
|
>
|
![g := create([-1, -1], eval(g_compts))](/support/helpjp/helpview.aspx?si=5520/file04532/math171.png)
|
|
![g := TABLE([index_char = [-1, -1], compts = matrix([[1, 0, 0], [0, 1, 0], [0, 0, 1]])])](/support/helpjp/helpview.aspx?si=5520/file04532/math174.png)
| (6) |
|
>
|
|
|
>
|
 
|
|
>
|
|
|
>
|
|
|
>
|
|
|
| (7) |
|
>
|
|
|
| (8) |
|
>
|
|
|
| (9) |
and in spherical-polar coordinates?
|
>
|
![coord := [r, theta, phi]](/support/helpjp/helpview.aspx?si=5520/file04532/math227.png)
|
|
![coord := [r, theta, phi]](/support/helpjp/helpview.aspx?si=5520/file04532/math230.png)
| (10) |
The metric is:
|
>
|
![g_compts := array(symmetric, sparse, 1 .. 3, 1 .. 3, [(1, 1) = 1, (2, 2) = r^2, (3, 3) = r^2*sin(theta)^2])](/support/helpjp/helpview.aspx?si=5520/file04532/math236.png)
|
|
>
|
![g := create([-1, -1], eval(g_compts))](/support/helpjp/helpview.aspx?si=5520/file04532/math240.png)
|
|
![g := TABLE([index_char = [-1, -1], compts = matrix([[1, 0, 0], [0, r^2, 0], [0, 0, r^2*sin(theta)^2]])])](/support/helpjp/helpview.aspx?si=5520/file04532/math243.png)
| (11) |
|
>
|
|
|
>
|
 
|
|
>
|
|
|
>
|
|
|
>
|
|
|
| (12) |
Now generate the geodesic equations:
|
>
|
|
|
| (13) |
|
|
