Every Lie algebra $\mathrm{\𝔤}$ admits a decomposition into the semi-direct sum $$$\mathrm{\𝔤}\mathit{ }\mathit{\=}\mathit{ }r \mathit{\+}\mathit{ }s\mathit{ }$,where $r$is the radical of $\mathrm{\𝔤}$and $s$ is a semi-simple subalgebra. Such a decomposition is called a Levi decomposition. Since the radical is an ideal we have $\left[r\, r\right] \subset r$, $\left[r\, s\right] \subset r$ and $s\,s\] \subset s\.$The radical$r$is uniquely defined but the semi-simple subalgebra $s$is not.

LeviDecomposition(Alg) calculates a Levi decomposition of the Lie algebra Alg,. If no argument is given, a Levi decomposition of the current algebra is computed.

With saveTemporaryAlgebras = true, the Lie algebras created by the program as part of the Levi decomposition algorithm are not erased and can be initialized for further analysis. Use this option in conjunction with userlevel[LeviDecomposition] := 2. The default is saveTemporaryAlgebras = false.

A list $\left[R\, S\right]$of two lists of vectors is returned, where $R$is a basis for the radical and $S$is a basis for the semi-simple part. Either list may be empty.

This module implements the algorithm for finding the Levi decomposition described in

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

$\mathrm{with}\left(\mathrm{DifferentialGeometry}\right)\:$$\mathrm{with}\left(\mathrm{LieAlgebras}\right)\:$

$\mathrm{L1}≔\mathrm{\_DG}\left(\left[\left[''LieAlgebra''\,\mathrm{Alg1}\,\left[6\right]\right]\,\left[\left[\left[1\,2\,1\right]\,-1\right]\,\left[\left[1\,2\,2\right]\,2\right]\,\left[\left[1\,2\,4\right]\,1\right]\,\left[\left[1\,3\,1\right]\,-1\right]\,\left[\left[1\,3\,3\right]\,1\right]\,\left[\left[1\,3\,4\right]\,1\right]\,\left[\left[1\,3\,6\right]\,-1\right]\,\left[\left[1\,4\,4\right]\,-2\right]\,\left[\left[1\,6\,1\right]\,-1\right]\,\left[\left[1\,6\,4\right]\,-1\right]\,\left[\left[1\,6\,6\right]\,-1\right]\,\left[\left[2\,3\,1\right]\,-1\right]\,\left[\left[2\,3\,3\right]\,1\right]\,\left[\left[2\,3\,4\right]\,1\right]\,\left[\left[2\,3\,6\right]\,-1\right]\,\left[\left[2\,4\,1\right]\,1\right]\,\left[\left[2\,4\,4\right]\,-1\right]\,\left[\left[2\,6\,1\right]\,-1\right]\,\left[\left[2\,6\,2\right]\,2\right]\,\left[\left[2\,6\,3\right]\,1\right]\,\left[\left[2\,6\,4\right]\,1\right]\,\left[\left[2\,6\,6\right]\,-1\right]\,\left[\left[3\,4\,1\right]\,-1\right]\,\left[\left[3\,4\,4\right]\,1\right]\,\left[\left[3\,4\,6\right]\,-1\right]\,\left[\left[3\,5\,3\right]\,1\right]\,\left[\left[3\,6\,1\right]\,-2\right]\,\left[\left[3\,6\,3\right]\,1\right]\,\left[\left[3\,6\,4\right]\,2\right]\,\left[\left[3\,6\,6\right]\,-2\right]\,\left[\left[4\,6\,4\right]\,-2\right]\,\left[\left[5\,6\,1\right]\,-1\right]\,\left[\left[5\,6\,4\right]\,1\right]\,\left[\left[5\,6\,6\right]\,-1\right]\right]\right]\right)$

$\mathrm{L1}≔\left[\left[\mathrm{e1}\,\mathrm{e2}\right]\=-\mathrm{e1}\+2\mathrm{e2}\+\mathrm{e4}\,\left[\mathrm{e1}\,\mathrm{e3}\right]\=-\mathrm{e1}\+\mathrm{e3}\+\mathrm{e4}-\mathrm{e6}\,\left[\mathrm{e1}\,\mathrm{e4}\right]\=-2\mathrm{e4}\,\left[\mathrm{e1}\,\mathrm{e6}\right]\=-\mathrm{e1}-\mathrm{e4}-\mathrm{e6}\,\left[\mathrm{e2}\,\mathrm{e3}\right]\=-\mathrm{e1}\+\mathrm{e3}\+\mathrm{e4}-\mathrm{e6}\,\left[\mathrm{e2}\,\mathrm{e4}\right]\=\mathrm{e1}-\mathrm{e4}\,\left[\mathrm{e2}\,\mathrm{e6}\right]\=-\mathrm{e1}\+2\mathrm{e2}\+\mathrm{e3}\+\mathrm{e4}-\mathrm{e6}\,\left[\mathrm{e3}\,\mathrm{e4}\right]\=-\mathrm{e1}\+\mathrm{e4}-\mathrm{e6}\,\left[\mathrm{e3}\,\mathrm{e5}\right]\=\mathrm{e3}\,\left[\mathrm{e3}\,\mathrm{e6}\right]\=-2\mathrm{e1}\+\mathrm{e3}\+2\mathrm{e4}-2\mathrm{e6}\,\left[\mathrm{e4}\,\mathrm{e6}\right]\=-2\mathrm{e4}\,\left[\mathrm{e5}\,\mathrm{e6}\right]\=-\mathrm{e1}\+\mathrm{e4}-\mathrm{e6}\right]$

$\mathrm{DGsetup}\left(\mathrm{L1}\right)\:$

$\mathrm{LD}≔\mathrm{LeviDecomposition}\left(\right)$

$\mathrm{LD}≔\left[\left[\mathrm{e1}-\mathrm{e4}\+\mathrm{e6}\,\mathrm{e3}\,\mathrm{e5}\right]\,\left[\mathrm{e1}\,\mathrm{e2}\,\mathrm{e4}\right]\right]$

$R≔\mathrm{LD}\left[1\right]\;$$S≔\mathrm{LD}\left[2\right]$

$R≔\left[\mathrm{e1}-\mathrm{e4}\+\mathrm{e6}\,\mathrm{e3}\,\mathrm{e5}\right]$

$S≔\left[\mathrm{e1}\,\mathrm{e2}\,\mathrm{e4}\right]$

$\mathrm{Radical}\left(\right)$

$\left[\mathrm{e3}\,\mathrm{e1}-\mathrm{e4}\+\mathrm{e6}\,\mathrm{e5}\right]$

$\mathrm{Query}\left(R\,''Ideal''\right)$

$\mathrm{true}$

$\mathrm{Query}\left(R\,''Solvable''\right)$

$\mathrm{true}$

$\mathrm{Query}\left(S\,''Subalgebra''\right)$

$\mathrm{true}$

$\mathrm{Query}\left(S\,''Semisimple''\right)$

$\mathrm{true}$

$B≔\mathrm{map}\left(\mathrm{op}\,\mathrm{LD}\right)$

$B≔\left[\mathrm{e1}-\mathrm{e4}\+\mathrm{e6}\,\mathrm{e3}\,\mathrm{e5}\,\mathrm{e1}\,\mathrm{e2}\,\mathrm{e4}\right]$

$\mathrm{L3}≔\mathrm{LieAlgebraData}\left(B\,\mathrm{Alg3}\right)$

$\mathrm{L3}≔\left[\left[\mathrm{e1}\,\mathrm{e3}\right]\=\mathrm{e1}\,\left[\mathrm{e1}\,\mathrm{e4}\right]\=\mathrm{e1}\,\left[\mathrm{e1}\,\mathrm{e5}\right]\=\mathrm{e1}-\mathrm{e2}\,\left[\mathrm{e2}\,\mathrm{e3}\right]\=\mathrm{e2}\,\left[\mathrm{e2}\,\mathrm{e4}\right]\=\mathrm{e1}-\mathrm{e2}\,\left[\mathrm{e2}\,\mathrm{e5}\right]\=\mathrm{e1}-\mathrm{e2}\,\left[\mathrm{e2}\,\mathrm{e6}\right]\=-\mathrm{e1}\,\left[\mathrm{e4}\,\mathrm{e5}\right]\=-\mathrm{e4}\+2\mathrm{e5}\+\mathrm{e6}\,\left[\mathrm{e4}\,\mathrm{e6}\right]\=-2\mathrm{e6}\,\left[\mathrm{e5}\,\mathrm{e6}\right]\=\mathrm{e4}-\mathrm{e6}\right]$

$\mathrm{DGsetup}\left(\mathrm{L3}\right)$

$\mathrm{Lie\; algebra:\; Alg3}$

$\mathrm{MultiplicationTable}\left(''LieBracket''\right)$

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