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[d] - exterior differentiation

Calling Sequence:

     d(expr)

Parameters:

      expr - any expression .

Description:

  • The d procedure calculates exterior derivative of a p-form expression.
  • The exterior derivative is operator  d : Omega[p]  -> Omega[p+1] where Omega[p]is p-form and Omega[p+1]is (p+1)-form. The operator has the following properties.
  • For any 0-form f = f(x[1],x[2],`` .. x[n])we have:d(f) = Sum(Diff(f,x[i])*d(x[i]),i = 1 .. n)
  • For any p-forms omega, sigmaand constants alpha, betawe have:d(alpha*omega+beta*sigma) = alpha*d(omega)+beta*d(sigma)
  • If omega[1]is p-form and omega[2]is q-form then exterior derivative of their exterior product is as follows: d(`&^`(omega[1],omega[2])) = `&^`(d(omega[1]),omega[2])+(-1)^p*`&^`(omega[1],d(omega[2]))
  • For any p-form OmegaPoincare's lemma takes place:d(d(Omega)) = 0

Examples:
restart:
with(atlas):

Declare constants:
Constants(alpha,beta,p,q);

{alpha, beta, -I, q, Pi, _Z, I, Catalan, p}

Declare functions:
Functions(F=F(x[1],x[2],x[3]),f=f(x[k]),G=G(x,z),z=z(phi),x[k]=x[k](y[j]));

{F, x[k], f, G, z}

Declare p-forms:
Forms(omega=n,sigma=m,Omega=P,omega[1]=p,omega[2]=q);

{omega, sigma, Omega, omega[1], omega[2]}

Using d- procedure:

f- declared as a function with indefinite number of variables:
'd'(f)=d(f);

d(f) = Sum(Diff(f,x[l[1]])*Sum(Diff(x[l[1]],y[l[2]])*d(y[l[2]]),l[2] = 1 .. dim),l[1] = 1 .. dim)

F- declared as function on three variables:
'd'(F)=d(F);

d(F) = Diff(F,x[1])*Sum(Diff(x[1],y[l[1]])*d(y[l[1]]),l[1] = 1 .. dim)+Diff(F,x[2])*Sum(Diff(x[2],y[l[1]])*d(y[l[1]]),l[1] = 1 .. dim)+Diff(F,x[3])*Sum(Diff(x[3],y[l[1]])*d(y[l[1]]),l[1] = 1 .. dim)
d(F) = Diff(F,x[1])*Sum(Diff(x[1],y[l[1]])*d(y[l[1]]),l[1] = 1 .. dim)+Diff(F,x[2])*Sum(Diff(x[2],y[l[1]])*d(y[l[1]]),l[1] = 1 .. dim)+Diff(F,x[3])*Sum(Diff(x[3],y[l[1]])*d(y[l[1]]),l[1] = 1 .. dim)

G- declared as function G=G(x,z) and z=z(phi) thus:
'd'(G)=d(G);

d(G) = Diff(G,x)*d(x)+Diff(G,z)*Diff(z,phi)*d(phi)

As alpha  declared as constant thus:
'd'(alpha)=d(alpha);

d(alpha) = 0

Exterior derivative is linear operation:
'd'(alpha*omega+beta*sigma)=d(alpha*omega+beta*sigma);

d(alpha*omega+beta*sigma) = alpha*d(omega)+beta*d(sigma)

There were not any declarations about x and y so they are 0-forms by defaults:
'd'(x*y)=d(x*y);

d(x*y) = d(x)*y+x*d(y)

Verify that double-d gives 0 (Poincare's lemma):
'd'(d(Omega))=d(d(Omega));

d(d(Omega)) = 0

As y is 0-form (by defaults) and Omega  has been declared as P-form then:
'd'(y*Omega)=d(y*Omega);

d(y*Omega) = (-1)^P*`&^`(Omega,d(y))+y*d(Omega)

omega[1]and omega[2]  - declared as p- and q-forms respectively:
'd'(omega[1]&^omega[2])=d(omega[1]&^omega[2]);

d(`&^`(omega[1],omega[2])) = (-1)^((1+p)*q)*`&^`(omega[2],d(omega[1]))+(-1)^p*`&^`(omega[1],d(omega[2]))

See Also:

atlas , atlas[Constants] , atlas[Functions] , atlas[Forms] , atlas[`&^`] .