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[`&/`] - restriction of a [0, k] tensor field under a mapping

Calling Sequence:

     Expr &/ MapId

Parameters:

      MapId - variable  - the mapping identifier  i.e. MapId : dom1 ---> dom2
      Expr - expression - a tensor expression which has to be restricted

Description:

  • The &/ procedure calculates restriction of a tensor field under a mapping. The restriction is linear  operation defined on [0,k] tensors only. The definition is as follows.
  • Let M and N be manifolds of dimensions m = dim(M), n = dim(N) . Let F be mapping between the manifolds: F: proc (M) options operator, arrow; N end proc  defined by functions:

    PIECEWISE([``, y[1] = y[1](x[1],x[2],`` .. x[m])],[``, y[2] = y[2](x[1],x[2],`` .. x[m])],[``, `.........`],[``, y[n] = y[n](x[1],x[2],`` .. x[m])])

    where
    {x[1], x[2], `` .. x[m]}  are local coordinates on M and {y[1], y[2], `` .. y[n]}  are local coordinates on N (in some domains).
  • For any [0,1] tensor field T on N the restriction of  T  under F is tensor field Omega = `&/`(T,F)  on M with components Omega[j] = Sum(T[i](y[k](x[1],x[2],`` .. x[m]))*Diff(y[i],x[j]),i = 1 .. n)  in local coordinates.
  • For tensor product of any [0, k] tensor fields T[1], T[2]   on N the following formula takes place: `&/`(`&.`(T[1],T[2]),F) = `&.`(`&/`(T[1],F),`&/`(T[2],F))
  • The formulas considered above defined linear restriction operator &/ completely.
  • According to the definition it is necessary to calculate the restrictions on the domain M. Use atlas[Domain]  procedure to jump on M manifold if needed.

Examples:

The following example shows how the restriction operator can be used.
Let M be 2-dimentional sphere
S^2  and N be 3-dimensional Euclidean space R^3 . Let F:proc (M) options operator, arrow; N end proc   be standard embedding of sphere S^2  into R^3 .

restart:
with(atlas):

This procedure is presented just for appropriate simplification (see atlas[simp] ).
`atlas/simp`:=proc(a) normal(a);subs({cos(theta)^2=1-sin(theta)^2,cos(phi)^2=1-sin(phi)^2},%);normal(%); factor(%) end;

`atlas/simp` := proc (a) normal(a); subs({cos(theta)^2 = 1-sin(theta)^2, cos(phi)^2 = 1-sin(phi)^2},%); normal(%); factor(%) end proc
`atlas/simp` := proc (a) normal(a); subs({cos(theta)^2 = 1-sin(theta)^2, cos(phi)^2 = 1-sin(phi)^2},%); normal(%); factor(%) end proc
`atlas/simp` := proc (a) normal(a); subs({cos(theta)^2 = 1-sin(theta)^2, cos(phi)^2 = 1-sin(phi)^2},%); normal(%); factor(%) end proc

Declare 1-forms e[j] and u[k] for corresponding coframes:
Forms(e[j]=1,u[k]=1);

{e[j], u[k]}

Declare vectors for corresponding frames:
Vectors(E[j],U[k]);

{E[j], U[k]}

Declare Euclidean space - R^3 :
Domain(R^3);

R^3

Declare coframe on R^3 :
Coframe(e[1]=d(x),e[2]=d(y),e[3]=d(z));

[e[1] = d(x), e[2] = d(y), e[3] = d(z)]

Declare frame on R^3 :
Frame(E[j]);

[E[1] = Diff(``,x), E[2] = Diff(``,y), E[3] = Diff(``,z)]

Declare metric on R^3  (standard flat metric):
Metric(g=d(x)&.d(x)+d(y)&.d(y)+d(z)&.d(z));

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

Declare sphere - S^2 :
Domain(S^2);

S^2

Declare coframe on S^2 :
Coframe(u[1]=d(theta),u[2]=d(phi));

[u[1] = d(theta), u[2] = d(phi)]

Declare frame on S^2 :
Frame(U[j]);

[U[1] = Diff(``,theta), U[2] = Diff(``,phi)]

Declare definite mapping F:proc (S) options operator, arrow; R^3 end proc :
Mapping(F,S^2,R^3,
        x=sin(theta)*cos(phi),
        y=sin(theta)*sin(phi),
        z=cos(theta));

F

S^2.`--->`.(R^3)

Who(F);


F: mapping


TABLE([coframe = {e[1] = cos(theta)*u[1]*cos(phi)-sin(theta)*sin(phi)*u[2], e[2] = cos(theta)*u[1]*sin(phi)+sin(theta)*cos(phi)*u[2], e[3] = -sin(theta)*u[1]}, manifolds = [S^2, R^3], equations = [x = ...
TABLE([coframe = {e[1] = cos(theta)*u[1]*cos(phi)-sin(theta)*sin(phi)*u[2], e[2] = cos(theta)*u[1]*sin(phi)+sin(theta)*cos(phi)*u[2], e[3] = -sin(theta)*u[1]}, manifolds = [S^2, R^3], equations = [x = ...
TABLE([coframe = {e[1] = cos(theta)*u[1]*cos(phi)-sin(theta)*sin(phi)*u[2], e[2] = cos(theta)*u[1]*sin(phi)+sin(theta)*cos(phi)*u[2], e[3] = -sin(theta)*u[1]}, manifolds = [S^2, R^3], equations = [x = ...
TABLE([coframe = {e[1] = cos(theta)*u[1]*cos(phi)-sin(theta)*sin(phi)*u[2], e[2] = cos(theta)*u[1]*sin(phi)+sin(theta)*cos(phi)*u[2], e[3] = -sin(theta)*u[1]}, manifolds = [S^2, R^3], equations = [x = ...
TABLE([coframe = {e[1] = cos(theta)*u[1]*cos(phi)-sin(theta)*sin(phi)*u[2], e[2] = cos(theta)*u[1]*sin(phi)+sin(theta)*cos(phi)*u[2], e[3] = -sin(theta)*u[1]}, manifolds = [S^2, R^3], equations = [x = ...
TABLE([coframe = {e[1] = cos(theta)*u[1]*cos(phi)-sin(theta)*sin(phi)*u[2], e[2] = cos(theta)*u[1]*sin(phi)+sin(theta)*cos(phi)*u[2], e[3] = -sin(theta)*u[1]}, manifolds = [S^2, R^3], equations = [x = ...
TABLE([coframe = {e[1] = cos(theta)*u[1]*cos(phi)-sin(theta)*sin(phi)*u[2], e[2] = cos(theta)*u[1]*sin(phi)+sin(theta)*cos(phi)*u[2], e[3] = -sin(theta)*u[1]}, manifolds = [S^2, R^3], equations = [x = ...
TABLE([coframe = {e[1] = cos(theta)*u[1]*cos(phi)-sin(theta)*sin(phi)*u[2], e[2] = cos(theta)*u[1]*sin(phi)+sin(theta)*cos(phi)*u[2], e[3] = -sin(theta)*u[1]}, manifolds = [S^2, R^3], equations = [x = ...

Verify that we are on the sphere:
Domain();

S^2

Calculate metric induced on the sphere using &/ operator:
Metric(G = g &/ F);

G = `&.`(u[1],u[1])+sin(theta)^2*`&.`(u[2],u[2])

One can calculate restriction of any [0,k] tensor field on R^3  under the mapping:

Restriction of coframe 1-forms:
'e[1]&/F'=e[1]&/F;
'e[2]&/F'=e[2]&/F;

`&/`(e[1],F) = cos(theta)*u[1]*cos(phi)-sin(theta)*sin(phi)*u[2]

`&/`(e[2],F) = cos(theta)*u[1]*sin(phi)+sin(theta)*cos(phi)*u[2]

Restriction of 0-forms (scalars):
'(x^2+y^2)&/F'=(x^2+y^2)&/F;
'(y/x)&/F'=(y/x)&/F;

`&/`(x^2+y^2,F) = sin(theta)^2

`&/`(y/x,F) = sin(phi)/cos(phi)

Restriction of "rotation" 1-form:
'(x*d(y)-y*d(x))&/F'=(x*d(y)-y*d(x))&/F;

`&/`(x*d(y)-y*d(x),F) = sin(theta)^2*u[2]

Restriction of tensor product d(x)&.d(z):
'(d(x)&.d(z))&/F'=(d(x)&.d(z))&/F;

`&/`(`&.`(d(x),d(z)),F) = -sin(theta)*cos(theta)*cos(phi)*`&.`(u[1],u[1])+sin(theta)^2*sin(phi)*`&.`(u[2],u[1])

Restriction of exterior product d(x)&^d(y):
'(d(x)&^d(y))&/F'=(d(x)&^d(y))&/F;

`&/`(`&^`(d(x),d(y)),F) = cos(theta)*sin(theta)*`&^`(u[1],u[2])

Some more examples

Declare abstract mapping between S^2  and R^3 :
Mapping(Phi,S^2,R^3);

Phi

S^2.`--->`.(R^3)

Restriction of exterior product d(x)&^d(y) under abstract mapping Phi :
'(d(x)&^d(y))&/Phi'=(d(x)&^d(y))&/Phi;

`&/`(`&^`(d(x),d(y)),Phi) = (Diff(x,theta)*Diff(y,phi)-Diff(x,phi)*Diff(y,theta))*`&^`(u[1],u[2])

Restriction of coframe 1-forms
'e[1]&/Phi'=e[1]&/Phi;
'e[2]&/Phi'=e[2]&/Phi;

`&/`(e[1],Phi) = Diff(x,theta)*u[1]+Diff(x,phi)*u[2]

`&/`(e[2],Phi) = Diff(y,theta)*u[1]+Diff(y,phi)*u[2]

Who is who?
Who();

PIECEWISE([{R^3, S^2}, Domains],[{F, Phi}, Mappings],[{e[j], u[k], g, E[j], U[k], G}, Tensors],[{e[j], u[k]}, Forms],[{_Z, Pi, Catalan, I, -I}, Constants],[{}, Functions])

See Also:

atlas , atlas[Mapping] , atlas[Domain] , atlas[Metric] .