atlas 2 for Maple

atlas package Overview

Declaration operators

Calculation operators

Standard DG operators

Utility operators

Features List & Examples

atlasWizard - Maplet™

Template Worksheets

Licenses & Pricing

atlas[div] - divergence operator 

Calling Sequence: 

    div(expr) 

Parameters: 

       expr - any vector expression. 

Description: 

  • The div - procedure allows one to calculate the divergence operator on a vector expression.
 

Examples: 

> restart:
with(atlas):
 

Declare constants:  

> Constants(Lambda);
 

{`+`(`-`(I)), I, Pi, _Z, Catalan, Lambda}(2.1)
 

Declare functions:  

> Functions(f=f(x,y),h=h(f));
 

{f, h}(2.2)
 

Declare p-forms:  

> Forms(e[k]=1,xi=1);
 

{xi, e[k]}(2.3)
 

Declare vectors:  

> Vectors(E[j],X,Y,Z);
 

{X, Y, Z, E[j]}(2.4)
 

By definition:
'div(X+Y)'=div(X+Y);
'div(F*X)'=div(F*X);
 

 

div(`+`(X, Y)) = `+`(div(X), div(Y))
div(`*`(F, `*`(X))) = `+`(`*`(F, `*`(div(X))), iota[X](d(F)))(2.5)
 

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

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

Declare frame:
Frame(E[k]); 

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

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))))(2.8)
 

Calculate connection:
Connection(omega); 

omega[i, j](2.9)
 

Divergence of frame vectors:
'div(E[1])'=div(E[1]); 

div(E[1]) = `+`(`-`(`/`(`*`(4, `*`(Lambda, `*`(x))), `*`(`+`(1, `*`(Lambda, `*`(`^`(x, 2))), `*`(Lambda, `*`(`^`(y, 2))))))))(2.10)
 

Divergence of vector X:
'div(X)'=div(ToBasis(X)); 

div(X) = `+`(`-`(`/`(`*`(`+`(`*`(4, `*`(iota[X](e[1]), `*`(Lambda, `*`(x)))), `-`(iota[E[1]](d(iota[X](e[1])))), `-`(`*`(iota[E[1]](d(iota[X](e[1]))), `*`(Lambda, `*`(`^`(x, 2))))), `-`(`*`(iota[E[1]]...
div(X) = `+`(`-`(`/`(`*`(`+`(`*`(4, `*`(iota[X](e[1]), `*`(Lambda, `*`(x)))), `-`(iota[E[1]](d(iota[X](e[1])))), `-`(`*`(iota[E[1]](d(iota[X](e[1]))), `*`(Lambda, `*`(`^`(x, 2))))), `-`(`*`(iota[E[1]]...
div(X) = `+`(`-`(`/`(`*`(`+`(`*`(4, `*`(iota[X](e[1]), `*`(Lambda, `*`(x)))), `-`(iota[E[1]](d(iota[X](e[1])))), `-`(`*`(iota[E[1]](d(iota[X](e[1]))), `*`(Lambda, `*`(`^`(x, 2))))), `-`(`*`(iota[E[1]]...
div(X) = `+`(`-`(`/`(`*`(`+`(`*`(4, `*`(iota[X](e[1]), `*`(Lambda, `*`(x)))), `-`(iota[E[1]](d(iota[X](e[1])))), `-`(`*`(iota[E[1]](d(iota[X](e[1]))), `*`(Lambda, `*`(`^`(x, 2))))), `-`(`*`(iota[E[1]]...
(2.11)
 

Divergence of "rotation" vector:
'div(y*E[1]-x*E[2])'=div(y*E[1]-x*E[2]); 

div(`+`(`*`(y, `*`(E[1])), `-`(`*`(x, `*`(E[2]))))) = 0(2.12)
 

>
 

See Also:  

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