In the LieAlgebras package, the command DGsetup is used to initialize a Lie algebra, that is, to define the basis elements for the Lie algebra and its dual and to store the structure constants for the Lie algebra in memory.  The first argument for DGsetup is a Lie algebra data structure which contains the structure constants in a standard format used by the LieAlgebras package.

An important class of Lie algebras are provided by the infinitesimal generators (vector fields) for any r-parameter transformation group. These vector fields form a Lie algebra with respect the LieBracket operation. The program LieAlgebraData(Gamma, AlgName) will create the Lie algebra data structure for the Lie algebra of vector fields Gamma.

This program uses the GetComponents command with the option method = "real" to calculate the Lie brackets of the vector fields in Gamma as constant, real, linear combinations of the vectors in Gamma. When the component functions of the vector fields in Gamma are transcendental functions, it is best to specify a small list of points at which these component functions and their derivatives are easily evaluated.

The command LieAlgebraData is part of the DifferentialGeometry:-LieAlgebras package.  It can be used in the form LieAlgebraData(...) only after executing the commands with(DifferentialGeometry) and with(LieAlgebras), but can always be used by executing DifferentialGeometry:-LieAlgebras:-LieAlgebraData(...).

with(DifferentialGeometry): with(LieAlgebras):

DGsetup([x, y], M1);

$\mathrm{frame\; name:\; M1}$

Gamma1 := evalDG([D_x, D_y, x*D_x, y*D_x]);

$\mathrm{\Γ1}:=\left[\mathrm{D\_x}\,\mathrm{D\_y}\,x\mathrm{D\_x}\,y\mathrm{D\_x}\right]$

LieAlgebraData(Gamma1, Alg1);

$\left[\left[\mathrm{e1}\,\mathrm{e3}\right]\=\mathrm{e1}\,\left[\mathrm{e2}\,\mathrm{e4}\right]\=\mathrm{e1}\,\left[\mathrm{e3}\,\mathrm{e4}\right]\=-\mathrm{e4}\right]$

DGsetup([x, y, z], M2);

$\mathrm{frame\; name:\; M2}$

Gamma2 := evalDG([D_y, cos(y)*D_x-cos(x)/sin(x)*sin(y)*D_y+sin(y)/sin(x)*D_z, -sin(y)*D_x-cos(x)/sin(x)*cos(y)*D_y+cos(y)/sin(x)*D_z]);

$\mathrm{\Γ2}:=\left[\mathrm{D\_y}\,\cos \left(y\right)\mathrm{D\_x}-\frac{\cos \left(x\right)\sin \left(y\right)\mathrm{D\_y}}{\sin \left(x\right)}\+\frac{\sin \left(y\right)\mathrm{D\_z}}{\sin \left(x\right)}\,-\sin \left(y\right)\mathrm{D\_x}-\frac{\cos \left(x\right)\cos \left(y\right)\mathrm{D\_y}}{\sin \left(x\right)}\+\frac{\cos \left(y\right)\mathrm{D\_z}}{\sin \left(x\right)}\right]$

$\mathrm{\Γ2}:=\left[\mathrm{D\_y}\,\cos \left(y\right)\mathrm{D\_x}-\frac{\cos \left(x\right)\sin \left(y\right)\mathrm{D\_y}}{\sin \left(x\right)}\+\frac{\sin \left(y\right)\mathrm{D\_z}}{\sin \left(x\right)}\,-\sin \left(y\right)\mathrm{D\_x}-\frac{\cos \left(x\right)\cos \left(y\right)\mathrm{D\_y}}{\sin \left(x\right)}\+\frac{\cos \left(y\right)\mathrm{D\_z}}{\sin \left(x\right)}\right]$

LieAlgebraData(Gamma2, Alg2, initialpointlist = [[x = Pi/2, y = 0, z = 0]]);

$\left[\left[\mathrm{e1}\,\mathrm{e2}\right]\=\mathrm{e3}\,\left[\mathrm{e1}\,\mathrm{e3}\right]\=-\mathrm{e2}\,\left[\mathrm{e2}\,\mathrm{e3}\right]\=\mathrm{e1}\right]$