Digi Area Group - Math Tools for Professionals
   Maple and Mathematica packages - math tools for professionals

atlas™ - modern differential geometry for Maple™

> Features List & Examples
> Template Worksheets
> Screenshots
> Documentation & Downloads
> License & Pricing
> Buy Online

 
 
 
 
Google

Features List & Examples  |  Introduction  |  Dimension  |  Indexing  |  Forms  |  Metric  |  atlasWizard - Maplet™

atlas[Riemann] - calculation of Riemannian tensor

atlas[Ricci] - calculation of Ricci tensor

atlas[RicciScalar] - calculation of Ricci scalar

Calling Sequence:

     Riemann(Id)
     Ricci(Id)
     RicciScalar(Id)

Parameters:

      Id - variable - corresponding identifier

Description:

  • The Riemann  procedure calculates curvature tensor. The procedure is only available if the curvature 2-forms have been calculated (see atlas[Curvature] ).
  • The Ricci  procedure calculates Ricci tensor. The procedure is only available if the curvature 2-forms (see atlas[Curvature] ) has been calculated.
  • The RicciScalar  procedure calculates Ricci scalar. The procedure is only available if metric tensor is definite (see atlas[Metric] ) and the Ricci tensor has been calculated.

Examples:

3-dimensional sphere
restart:
with(atlas):

Declare forms:
Forms(e[j]=1,xi=1);

{xi, e[j]}

Declare vectors:
Vectors(X,Y,Z,E[j]);

{X, Y, Z, E[j]}

Declare constant lambda :
Constants(lambda);

{Catalan, _Z, Pi, I, lambda, -I}

Declare coframe:
Coframe(e[1]=2*d(x)/(1+lambda*(x^2+y^2+z^2)),
e[2]=2*d(y)/(1+lambda*(x^2+y^2+z^2)),e[3]=2*d(z)/(1+lambda*(x^2+y^2+z^2)));

[e[1] = 2*d(x)/(1+lambda*(x^2+y^2+z^2)), e[2] = 2*d(y)/(1+lambda*(x^2+y^2+z^2)), e[3] = 2*d(z)/(1+lambda*(x^2+y^2+z^2))]

Declare frame:
Frame(E[i]);

[E[1] = (1/2+1/2*lambda*x^2+1/2*lambda*y^2+1/2*lambda*z^2)*Diff(``,x), E[2] = (1/2+1/2*lambda*x^2+1/2*lambda*y^2+1/2*lambda*z^2)*Diff(``,y), E[3] = (1/2+1/2*lambda*x^2+1/2*lambda*y^2+1/2*lambda*z^2)*Di...
[E[1] = (1/2+1/2*lambda*x^2+1/2*lambda*y^2+1/2*lambda*z^2)*Diff(``,x), E[2] = (1/2+1/2*lambda*x^2+1/2*lambda*y^2+1/2*lambda*z^2)*Diff(``,y), E[3] = (1/2+1/2*lambda*x^2+1/2*lambda*y^2+1/2*lambda*z^2)*Di...

 

d(x);

1/2*e[1]+1/2*e[1]*lambda*x^2+1/2*e[1]*lambda*y^2+1/2*e[1]*lambda*z^2

Declare metric on   S^3  (see atlas[Metric] ):
Metric(g=e[1]&.e[1]+e[2]&.e[2]+e[3]&.e[3]);

g = `&.`(e[1],e[1])+`&.`(e[2],e[2])+`&.`(e[3],e[3])

Connection calculation:
Connection(omega);

omega[i,j]

Curvature calculation:
Curvature(Omega);

Omega[i,j]

Riemannian tensor calculation:
Riemann(R);

R = lambda*`&.`(`&^`(e[1],e[2]),`&^`(e[1],e[2]))+lambda*`&.`(`&^`(e[1],e[3]),`&^`(e[1],e[3]))+lambda*`&.`(`&^`(e[2],e[3]),`&^`(e[2],e[3]))
R = lambda*`&.`(`&^`(e[1],e[2]),`&^`(e[1],e[2]))+lambda*`&.`(`&^`(e[1],e[3]),`&^`(e[1],e[3]))+lambda*`&.`(`&^`(e[2],e[3]),`&^`(e[2],e[3]))

Ricci tensor calculation:
Ricci(r);

r = 2*lambda*`&.`(e[1],e[1])+2*lambda*`&.`(e[2],e[2])+2*lambda*`&.`(e[3],e[3])

Ricci scalar calculation:
RicciScalar(s);

s = 6*lambda

Example 2
restart:
with(atlas):

Declare forms:
Forms(e[j]=1,xi=1);

{e[j], xi}

Declare vectors:
Vectors(X,Y,Z,E[j]);

{X, Y, Z, E[j]}

Declare coframe:
Coframe(e[1]=x*d(x)+y*d(y),e[2]=x*d(y)-y*d(x));

[e[1] = x*d(x)+y*d(y), e[2] = x*d(y)-y*d(x)]

Declare frame:
Frame(E[i]);

[E[1] = 1/(y^2+x^2)*x*Diff(``,x)+1/(y^2+x^2)*y*Diff(``,y), E[2] = -1/(y^2+x^2)*y*Diff(``,x)+1/(y^2+x^2)*x*Diff(``,y)]

Connection definition:
omega[1,1]:=x*e[1];

omega[1,1] := x*e[1]

omega[2,2]:=y*e[2];

omega[2,2] := e[2]*y

omega[1,2]:=y*e[1];

omega[1,2] := y*e[1]

omega[2,1]:=-x*e[2];

omega[2,1] := -x*e[2]

Connection declaration:
Connection(omega);

omega[i,j]

Curvature calculation:
Curvature(Omega);

Omega[i,j]

Riemann(R);

R = 1/2*y*(-1+x^3+x*y^2)/(y^2+x^2)*`&.`(E[1],e[1],`&^`(e[2],e[1]))-1/2*x*(-3+x^3+x*y^2)/(y^2+x^2)*`&.`(E[2],e[1],`&^`(e[2],e[1]))-1/2*(-x+y^4+y^2*x^2)/(y^2+x^2)*`&.`(E[1],e[2],`&^`(e[2],e[1]))-1/2*y*(3...
R = 1/2*y*(-1+x^3+x*y^2)/(y^2+x^2)*`&.`(E[1],e[1],`&^`(e[2],e[1]))-1/2*x*(-3+x^3+x*y^2)/(y^2+x^2)*`&.`(E[2],e[1],`&^`(e[2],e[1]))-1/2*(-x+y^4+y^2*x^2)/(y^2+x^2)*`&.`(E[1],e[2],`&^`(e[2],e[1]))-1/2*y*(3...

Ricci calculation:
Ricci(r);

r = -y*(-1+x^3+x*y^2)/(y^2+x^2)*`&.`(e[1],e[2])-x*(-3+x^3+x*y^2)/(y^2+x^2)*`&.`(e[1],e[1])+(-x+y^4+y^2*x^2)/(y^2+x^2)*`&.`(e[2],e[2])-y*(3+x^3+x*y^2)/(y^2+x^2)*`&.`(e[2],e[1])
r = -y*(-1+x^3+x*y^2)/(y^2+x^2)*`&.`(e[1],e[2])-x*(-3+x^3+x*y^2)/(y^2+x^2)*`&.`(e[1],e[1])+(-x+y^4+y^2*x^2)/(y^2+x^2)*`&.`(e[2],e[2])-y*(3+x^3+x*y^2)/(y^2+x^2)*`&.`(e[2],e[1])

 

RicciScalar(s);


Warning, There is no actual metric tensor



See Also:

atlas , atlas[Connection] , atlas[Curvature] , atlas[Metric] .