| atlas[div] - divergence operator Calling Sequence: div(expr) Parameters: expr - any vector expression. Description: - The div procedure calculates the divergence operator on a vector expression.
Examples: restart: with(atlas): Declare constants: Constants(Lambda); 
Declare functions: Functions(f=f(x,y),h=h(f)); 
Declare p-forms: Forms(e[k]=1,xi=1); ![{e[k], xi}](prod/atlas/help/images/div3.gif)
Declare vectors: Vectors(E[j],X,Y,Z); ![{E[j], X, Y, Z}](prod/atlas/help/images/div4.gif)
By definition: 'div(X+Y)'=div(X+Y); 'div(F*X)'=div(F*X); 
)](prod/atlas/help/images/div6.gif)
Declare coframe: Coframe(e[1]=d(x),e[2]=d(y)); ![[e[1] = d(x), e[2] = d(y)]](prod/atlas/help/images/div7.gif)
Declare frame: Frame(E[k]); ![[E[1] = Diff(``,x), E[2] = Diff(``,y)]](prod/atlas/help/images/div8.gif)
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](prod/atlas/help/images/div9.gif)
Calculate connection: Connection(omega); ![omega[i,j]](prod/atlas/help/images/div10.gif)
Divergence of frame vectors: 'div(E[1])'=div(E[1]); ![div(E[1]) = -4*Lambda*x/(1+Lambda*x^2+Lambda*y^2)](prod/atlas/help/images/div11.gif)
Divergence of vector X: 'div(X)'=div(ToBasis(X)); *Lambda*x/(1+Lambda*x^2+Lambda*y^2)+iota[E[1]](d(iota[X](e[1])))-4*iota[X](e[2])/(1+Lambda*x^2+Lambda*y^2)*Lambda*y+iota[E[2]](d(iota[X](e[2])))](prod/atlas/help/images/div12.gif)
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](prod/atlas/help/images/div13.gif)
See Also: atlas , atlas[d] , atlas[cov] , atlas[`&.`] , atlas[`&^`] , atlas[iota] . |