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[cov] - covariant derivative

Calling Sequence:

     cov(V,expr)

Parameters:

         expr - any expression.
        V - vector field.

Description:

  • The cov  procedure calculates the covariant derivative on an expression along given vector field. The derivative has the following properties.
  • For any vector field X and 0-form f  we have:    cov(X,f) = iota[X](d(f))
  • For vector fields X, Y and Z we have: cov(X,Y+Z) = cov(X,Y)+cov(X,Z)  and cov(f*X+h*Y,Z) = f*cov(X,Z)+h*cos(Y,Z)  for any functions f  and h .
  • For any vector field X and tensor fields Omega  and T  the Leibniz rule for the derivative takes place: cov(X,`&.`(Omega,T)) = `&.`(cov(X,Omega),T)+`&.`(Omega,cov(X,T))

Examples:
restart:
with(atlas):

Declare constants:
Constants(Lambda);

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

Declare functions:
Functions(f=f(x,y),h=h(z));

{h, f}

Declare p-forms:
Forms(e[k]=1,xi=1);

{xi, e[k]}

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

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

Declare tensors:
Tensors(T=[n,k],Omega=[l,m]);

{T, Omega}

Using cov - procedure:

As h is 0-form by defaults then:
'cov(X,h)'=cov(X,h);

cov(X,h) = Diff(h,z)*iota[X](d(z))

f- declared as function f=f(x,y) thus:
'cov(X,f)'=cov(X,f);

cov(X,f) = Diff(f,x)*iota[X](d(x))+Diff(f,y)*iota[X](d(y))

Obviously that:
'cov(X,x^2*sin(x))'=cov(X,x^2*sin(x));

cov(X,x^2*sin(x)) = 2*x*sin(x)*iota[X](d(x))+x^2*cos(x)*iota[X](d(x))

As Lambda  declared as constant thus:
'cov(X,Lambda)'=cov(X,Lambda);

cov(X,Lambda) = 0

Lie derivative is linear with respect to any argument:
'cov(X,Y+Z)'=cov(X,Y+Z);

cov(X,Y+Z) = cov(X,Y)+cov(X,Z)

And:
'cov(X+Y,Z)'=cov(X+Y,Z);

cov(X+Y,Z) = cov(X,Z)+cov(Y,Z)

More complex examples:
'cov(f*X+h*Y,Z)'=cov(f*X+h*Y,Z);

cov(f*X+h*Y,Z) = f*cov(X,Z)+h*cov(Y,Z)

And:
'cov(Z,f*X+h*Y)'=cov(Z,f*X+h*Y);

cov(Z,f*X+h*Y) = (Diff(f,x)*iota[Z](d(x))+Diff(f,y)*iota[Z](d(y)))*X+f*cov(Z,X)+Diff(h,z)*iota[Z](d(z))*Y+h*cov(Z,Y)
cov(Z,f*X+h*Y) = (Diff(f,x)*iota[Z](d(x))+Diff(f,y)*iota[Z](d(y)))*X+f*cov(Z,X)+Diff(h,z)*iota[Z](d(z))*Y+h*cov(Z,Y)

Verify that Leibniz rule takes place for tensor product:
'cov(X,T&.Omega)'=cov(X,T&.Omega);

cov(X,`&.`(T,Omega)) = `&.`(cov(X,T),Omega)+`&.`(T,cov(X,Omega))

Examples of calculations with definite metric

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

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

Declare frame:
Frame(E[k]);

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

Declare metric:
Metric(g=4*(d(x)&.d(x)+d(y)&.d(y))/(1+Lambda*(x^2+y^2))^2);

g = 4*(`&.`(e[1],e[1])+`&.`(e[2],e[2]))/(1+Lambda*(x^2+y^2))^2

Calculate connection:
Connection(theta);

theta[i,j]

For frame vectors E[j]:
'cov(E[j],sin(x)*cos(y))'=cov(E[j],sin(x)*cos(y));

cov(E[j],sin(x)*cos(y)) = cos(x)*cos(y)*delta[1,j]-sin(x)*sin(y)*delta[2,j]

As g is a metric then for any vector field X:
'cov(X,g)'=cov(X,g);

cov(X,g) = 0

By definition:
'cov(E[i],E[1])'=cov(E[i],E[1]);

cov(E[i],E[1]) = (-2*Lambda*x/(1+Lambda*x^2+Lambda*y^2)*delta[1,i]-2/(1+Lambda*x^2+Lambda*y^2)*Lambda*y*delta[2,i])*E[1]+(2/(1+Lambda*x^2+Lambda*y^2)*Lambda*y*delta[1,i]-2*Lambda*x/(1+Lambda*x^2+Lambda...
cov(E[i],E[1]) = (-2*Lambda*x/(1+Lambda*x^2+Lambda*y^2)*delta[1,i]-2/(1+Lambda*x^2+Lambda*y^2)*Lambda*y*delta[2,i])*E[1]+(2/(1+Lambda*x^2+Lambda*y^2)*Lambda*y*delta[1,i]-2*Lambda*x/(1+Lambda*x^2+Lambda...

For coframe 1-forms:
'cov(E[k],e[2])'=cov(E[k],ToBasis(e[2]));

cov(E[k],e[2]) = -(2/(1+Lambda*x^2+Lambda*y^2)*Lambda*y*delta[1,k]-2*Lambda*x/(1+Lambda*x^2+Lambda*y^2)*delta[2,k])*e[1]-(-2*Lambda*x/(1+Lambda*x^2+Lambda*y^2)*delta[1,k]-2/(1+Lambda*x^2+Lambda*y^2)*La...
cov(E[k],e[2]) = -(2/(1+Lambda*x^2+Lambda*y^2)*Lambda*y*delta[1,k]-2*Lambda*x/(1+Lambda*x^2+Lambda*y^2)*delta[2,k])*e[1]-(-2*Lambda*x/(1+Lambda*x^2+Lambda*y^2)*delta[1,k]-2/(1+Lambda*x^2+Lambda*y^2)*La...

For exterior product:
'cov(E[k],e[1]&^e[2])'=cov(E[k],e[1]&^e[2]);

cov(E[k],`&^`(e[1],e[2])) = 4*Lambda*(x*delta[1,k]+y*delta[2,k])/(1+Lambda*x^2+Lambda*y^2)*`&^`(e[1],e[2])

For tensor product:
'cov(E[k],e[2]&.E[2])'=cov(E[k],e[2]&.E[2]);

cov(E[k],`&.`(e[2],E[2])) = -2*Lambda*(y*delta[1,k]-x*delta[2,k])/(1+Lambda*x^2+Lambda*y^2)*`&.`(e[1],E[2])-2*Lambda*(y*delta[1,k]-x*delta[2,k])/(1+Lambda*x^2+Lambda*y^2)*`&.`(e[2],E[1])

See Also:

atlas , atlas[d] , atlas[cov] , atlas[`&.`] , atlas[`&^`] , atlas[iota] .