MVR5
I&k_scl%boe:\mathware\vessiot\programs\vessiot_source_code\LieAlgebras\Chevalle
y\G6",
I*Chevalley=F$%&falseGE\[l?%;classify_4dim_solvable_430GR6#%$algG66%"kG%$kvbG%#
daG%*alg_classG%*new_basisG%)alg_infoG%$H_1G%$X_1G%$Y_1G%$Z_1G%$cbmG%0eigen_inf
o_listG%&eigenG%*eval_listG%+evect_listG%"iG%"jG%(xy_valsG%)P_matrixG%$nilGF$F$
C4@$/9#"""-%6change_Lie_algebra_toG6#&9"6#FF>8%7#-%$seqG6$-%%vectG6#-%$catG6$%"
zG8$/FZ;FF-%3frameBaseDimensionGF$>8&-%5bracket_of_subspacesG6$FNFN>8*-%#opG6#-
%4complementary_basisG6$FjnFN>8/-&%'linalgG6#%-eigenvectorsG6#-%.restricted_adG
6$F_oFjn@$/-%%nopsG6#7#FgoFF>FgoFfp>817">82Fjp?&80Fgo%%trueGC$>Fip7$-Fao6#Fip&F
^qFL>F\q7$-Fao6#F\q-Fao6#&F^q6#""$@$/-FdpFiqFF>8'Q'L(4,6)F$@$/Far""#>FcrQ'L(4,4
)F$@$/FarF^r@%/-%(convertG6$-%$mapG6$%#ImGFip%$setG<#""!>FcrQ'L(4,2)F$>FcrQ'L(4
,5)F$@$/FcrFisC+?(83FFFF-FdpFdqF_q@$/-%$absG6#&Fip6#F`t-%$maxG6#-Fao6#-Fbs6$Fet
Fip%&breakG>F_o-%&&multG6$,$*$Fgt!""FguF_o>8+-%&v_zipG6%-F_s6$&F\qFht%%listGFjn
%%plusG>8,-F[v6%-F_s6$&-%'removeG6%%$hasGF\qF_vFLF`vFjnFav>8--F[v6%-F_s6$&Fiv6#
FgrF`vFjnFav>85-Fjv6%%#isG-Fbs6$Fao-Fbs6%%-linear_comboG-Fbs6%%,Lie_bracketG7$F
cvF^wF_oFbxFgs@%-Fiw6#2&FfwFdw&FfwFL>8(7&FiuFcvF^wF_oC$>Fjx7&FiuF^wFcvF_o>Ffw7$
FgxFhx@$/Ffw7$FguFguC$>Fjx7&FcvF^wFiu-Fcu6$FguF_o>Ffw7$FFFgu>8)7%%+winternitzG7
$""%FgrFfw@$/FcrFirC/-&Fjo6#%'jordanG6$F^p86>877#-FQ6$-F[v6%7#-FQ6$&Fiz6$F`t84/
F`t;FFF^rFjnFav/Ff[lFh[l>FiuFjp?(F`tFFFF-Fdp6#F[[lF_q@%/-Fdp6#-Fjv6%F\w-F^x6$-F
ax6$F_o&F[[lFhtF[[lFgsFF>Fiu7$-Fao6#FiuFh\l>FcvFh\l>F[[l7%&FiuFLFcv&FiuFdw@$/&-
F^x6$-Fax6$F_oFcvF[[lFLFgs>F[[l7%Fa]lFcvF`]l>Fiu&F[[lFL>F^w&F[[lF]r>F_o-Fcu6$*$
&-F^x6$-Fax6$F_oFiuF[[lFLFguF_o>Fcv-Fcu6$*$Fd]lFguFcv>Fjx7&FiuFcvF^wFgy>Ffw-F^x
6$-Fax6$&FjxF]r&Fjx6#F`z7#Fc_l>F\z7%F^z7$F`zF`zFfw@$/FcrF[tC*?(F`tFFFFFatF_q@$/
-FdsFftFgsF`u>FiuFju?(Ff[lFFFFFatF_q@$-Fiw6#2Fgs-Fds6#&Fip6#Ff[lF`u>Fcv-F[v6%-F
_s6$-Fbs6$%#ReG&F\qFj`lF`vFjnFav>F^w-F[v6%-F_s6$-Fbs6$FdsFcalF`vFjnFav>F_o-Fcu6
$,$*$Fg`lFguFguF_o>Ffw7$*&FgtFFFg`lFgu*&-FbalFh`lFFFg`lFgu@%-Fiw6#2FgsFhxC$>Fjx
F[y>F\z7%F^z7$F`z""&FfwC$>Fjx7&Fiu-Fcu6$FguF^w-Fcu6$FguFcvFgy>F\z7%F^zF]cl7$,$F
hxFgu,$FgxFgu@$/FcrFdrC*Fdz>F[[lF\[l>FiuF\^l>Fcv-Fcu6$&FipFL&F[[lFdw>F^w-Fcu6$*
$FcdlFgrF^^l>F_o-Fcu6$,$*$FcdlFguFguF_o>FjxF[y>F\z7$F^z7$F`z""'>8.-&Fjo6#%*tran
sposeG6#-%'matrixG6#-Fbs6%F^xFjxFN-%'RETURNG6#7$-&Fjo6#%(inverseG6#FdelF\zF$F$F
$%5classify_3dim_simpleGRF+65%%kdimGF/%-basis_vectorG%"HG%#CHG%+eigen_infoGF:F<
F;%(counterGF5F6%&_betaG%#HpG%#XpG%#YpGF2F=F8F$F$C2@$FD-%.change_alg_toGFI>FZFg
n>FN7#-FQ6$-FT6#(FYFfw/Ffw;FFFZ?&FjnFNF_qC'>FcrFjn>Fjx-Fdo6$7#FcrFN>F\z7#-Fio6#
-F_p6$Fcr7$Fcr-Fao6#Fjx>Fcv-Fjv6%F\w7#-FQ6$&F\z6$FfwFF/Ffw;FF-Fdp6#F\zFgs@$0Fcv
FjpF`u>FiuFjp?&F_oF\zF_q@$0&F_oFLFgs>Fiu7$F[]l&F_o6$F^rFF@$/-Fbs6$FdsFcv7$FgsFg
sC+>Fdel-F[v6%7#-FQ6$&F`]l6#F>/F>Fbhl7%Fcr&FjxFL&FjxFdwFav>Fgo-F[v6%7#-FQ6$&Fa]
l6#F./F.FbhlFi[mFav>Fip-Fcu6$*$&FcvFLFguFcr>F^q-Fao6#-F^x6$-Fax6$FdelFgo7#Fip>F
\q-Fcu6$*$F^qFguFdel>F`tFgo>Ff[l7%F`tFip-Fcu6$FgrF\q>Fiz-Ffel6#-F[fl6#-Fbs6%F^x
Ff[lFN-F`fl6#7$-%&evalmG6#-Fdfl6#Fiz7$F^z7$F^rF^cl?(F^wFFFF-Fdp6#FcvF_q@$-Fiw6$
-Fds6#&Fcv6#F^w%)positiveGF`u>Fip-Fcu6$*$Fb_mFguFcr>Fdel-F[v6%-F_s6$-Fbs6$Fbal&
FiuFe_mF`vF`ilFav>Fgo-F[v6%-F_s6$-Fbs6$FdsFb`mF`vF`ilFav>F^q*$-%%sqrtG6#-Fet6#F
[]mFgu>Ff[l7%-Fcu6$F^qFdel-Fcu6$,$F^qFguFgoFip@$/-F^x6$-Fax6$&Ff[lFL&Ff[lFdw7#&
Ff[lF]r7#FguC$>Fiz-Ff^m6#-%#&*G6$F\^m-F[fl6#7%7%FgsFFFgs7%FFFgsFF7%FFFgsFguFb^m
@$/Fjam7#FFC$>FizF\^m-F`fl6#7$Fe^m7$F^z7$F^rFbel-%&ERRORG6#QinIf~you~get~this~m
essage,~there~is~a~big~problem~with~the~code.F$F$F$F$%+decom_alg3GR6&%)alg_data
G%"AG%#f1G%#f2G60%#r1G%#r2G%"GGF]gl%#N1G%#N2G%"MG%"rG%#d1G%#d2G%#s1G%#s2G%,stru
c_constG%(new_algGF$F$C/-%&gcdexG6'9&9'%(_lambdaGFZFN>F_o-Ff^m6#-%%subsG6$/Fiem
9%-%'expandG6#*&FgemFFFZFF>Fjx-%3canonicalNullspaceG6#F_o>F\z-Fhfm6#-Ff^m6#,&FF
FFF_oFgu>Fjn-Ffel6#-F[fl6#7$Fail-FaoF^jl>Fcr-Fdfl6#Fjn>Fcv-FdpFbil>F^w,&&&&9$FL
F]rFLFFFcvFgu>Fip-%9change_Lie_algebra_basisG6%Fcr%'factorGFahm>FdelFjp>FgoFjp?
&F^q&FipFdwF_q@%1&FeqFLFcv>Fdel7$-FaoFgflF^q>Fgo7$-Fao6#Fgo7$7%,&F]imFFFcvFgu,&
&FeqFdwFFFcvFgu,&&FeqF]rFFFcvFgu&F^qFdw7%7$7%%(Lie_algG-FW6$&F`hmFdwFF7#FcvFdel
7$7%F`jm-FW6$FcjmFgr7#F^wFgo-Ff^m6#FcrF$F$F$%<classify_4dim_solvable_4310GRF+69
F/F0F1F2F3F4F5F6F7F8F9F:F;F<%#UVG%&x_valGF=F>%&elem1G%&elem2GFA%'_alphaGFaglF$F
$C8@$FDFG>FZ7#-FQ6$-FT6#-FW6$FYFf[l/Ff[lFfn>FN-F\o6$FZFZ>88-%6L31_to_canonical_
formG6#FN>F\z-Fao6#-Fdo6$Fc\nFZ>89&-F^x6$-Fax6$&Fc\nF]rF\zFc\nFL>8:&-F^x6$-Fax6
$&Fc\nFdwF\zFc\nFL>F\z-%'&minusG6$-%&&plusG6$F\z-Fcu6$F]]nF[^n-Fcu6$Fe]nFc]n>Fd
el-Fio6#-F_p6$F\zFc\n@$/-Fdp6#7#FdelFF>FdelF`_n>F^qFjp>FipFjp?&FgoFdelF_q@$0&Fg
oFLFgsC$>F^q7$-Fao6#F^qFg_n>Fip7$Fcq-Fao6#&FgoF]r@$/FatF^r@%/-F_s6$-Fbs6$FdsF^q
FesFfs>FjnQ'L(4,9)F$>FjnQ(L(4,12)F$@$/FatFgrC$@%0Ff`nFfs>FjnQ(L(4,11)F$>FjnQ'L(
4,8)F$@$/-Fdp6#-%'centerGF$Fgs>FjnQ(L(4,10)F$@$/FjnFfanC.>Fdel-Fio6#-F_p6$F\z7$
F[^nFc]n>F^qFjp>FipFjp?&FgoFdelF_qC$>F^qFj_n>FipF^`n>Fiu-F[v6%-F_s6$FcdlF`vFgbn
Fav>Fcv-F[v6%-F_s6$FjhmF`vFgbnFav>F_o&Fc\nFL>Fc\n-Fe\n6#7%F_oFiuFcv>F]]n-Fao6#-
F^x6$Fi]n7#F[^n>F\z-Fcu6$*$F]]nFguF\z>Fcr7$-Fao6#Fc\nF\z>Fjx7$F^z7$F`z"")@$/Fjn
F[anC->F_o-F[v6%FacnFc\nFav>Fiu-F[v6%FfcnFc\nFav>Fcv-F[v6%-F_s6$&FipF]rF`vFc\nF
av>Fc\nF[dn>F]]n-Fbs6$Fao7$Fadn-F^x6$Fa]n7#Fc]n@$-Fiw6#2-Fet6#&F]]nFL-Fet6#&F]]
nFdw>Fc\n7%Ficn-Fcu6$FguFc]nF[^n>F]]nF_dn>F\zFedn>F`t-Fao6#Fdfn>FcrFidn>Fjx7%F^
z7$F`z""*7#F`t@$/FjnF^bnC.>FipFcbn>Fiu-F[v6%-F_s6$-Fao6#F^fnF`vFgbnFav>Fcv-Fao6
#-Fdo6$7#FiuFgbn>F_o-Fax6$FiuFcv>Fc\nF[dn>F]]nF_dn>F\zFedn>Fe]n&-F^x6$Fa]nFgbnF
L>Fc\n7%FicnF[^n-Fcu6$*$Fe]nFguFc]n>Ficn-Fax6$F[^nFc]n>FcrFidn>Fjx7$F^z7$F`z"#5
@$/FjnFdanC.>FdelFcbn>F^qFjp>FipFjp?&FgoFdelF_qC$>F^qFj_n>FipF^`n>Fiu-F[v6%-F_s
6$-Fbs6$FbalFcdlF`vFgbnFav>Fcv-F[v6%-F_s6$-Fbs6$FdsFcdlF`vFgbnFav>F_oFicn>Fc\nF
[dn>F]]n-Fao6#-F^x6$Fa]nFcdn>F\zFedn>FcrFidn>Fjx7$F^z7$F`z"#6@$/FjnF]anC0>FdelF
cbn>F^qFjp>FipFjp?&FgoFdelF_qC$>F^qFj_n>FipF^`n>FiuFc[o>FcvFj[o>F_oFbin>Fc\nF[d
n>F]]n&-F^x6$Fi]nFgbnFdw>F\z-Fcu6$,$FgdnFguF\z>F`t,$-Fao6#-F^x6$-Fax6$FicnF\z7#
Ficn#FFFgr>FcrFidn>Fjx7%F^z7$F`z"#7F_hn@$-Fiw6#2F`tFgsC$>Fcr7&-Fcu6$Fgu&FcrFL&F
crF]r&FcrFdw-Fcu6$Fgu&FcrFe_l>Fjx7%F^zF`_o7#,$F`tFgu>F^w-Ffel6#-F[fl6#-Fbs6%F^x
FcrFZ-F`fl6#7$-FdflFe_mFjxF$F$F$%.classify_4dimGRF$F$F$F$C&@$FDFG@$-%.check_abe
lianGF$-F`fl6#7$-&Fjo6#%%diagG6&FFFFFFFF7$F^z7$F`zFgs@$-%0check_nilpotentGF$-F`
fl6#-%8classify_4dim_nilpotentGF$@$-%/check_solvableGF$-F`fl6#-%7classify_4dim_
solvableGF$F$F$F$%6central_decompositionGRF$65%)cent_algG%(der_algG%/central_fa
ctorG%"vG%%testGF=%"nG%#frG%*alg_basisG%.intersect_algGF,F>Fbdm%"BG%"CG%5comple
mentary_factorG%)alg_nameG%%ans1G%%ans2GF$F$C:@$FDFG>F_oFgn>Fiu%<_Vessiot_curre
nt_frame_nameG>F^q-F[fl6$F_oF_o>Fcv7#-FQ6$-FT6#(FYF\z/F\z;FFF_o@$%F_Vessiot_sho
w_intermediate_steps_flagG-%&printG6#%+Get~CenterG>FZF[bn@$Fceo-Feeo6#%4Get~Der
ived~AlgebraG>FN-F\o6$FcvFcv@$Fceo-Feeo6#%EIntersect~Center~and~Derived~Algebra
G>F^w-%Ccreate_intersection_of_subalgebrasG6#7$FZFN@$Fceo-Feeo6#%SFind~compleme
nt~C~of~intesection~algebra~in~centerG>Fjn-Fdo6$F^wFZ>Fjn-%9canonical_subspace_
basisG6$FjnFcv?(F\zFFFF-FdpFigmF_q?(FgoFFFFF_oF_q>&F^q6$FgoF\z-%+coeff_listG6$&
FjnF^jl7#Ffp@$Fceo-Feeo6#%@Find~complement~of~C~in~algebraG>Fip7$-FaoFigm-FaoFf
\n>F\q-Fdo6$FipFcv>F`t7$FfhoFhq>F`t-Fbgo6$F`tFcv?(F\zFFFF-FdpFhtF_q?(FgoFFFFF_o
F_q>&F^q6$Fgo,&F\zFFFegoFF-F[ho6$&F`tF^jlF^ho@%0FjnFjp>Ffw-%?subalgebra_to_Lie_
algebra_dataG6$Fjn-FW6$FiuFF>FfwFjp>Fiz-F]jo6$F`t-FW6$FiuFgr7%FfwFiz-Ff^m6#-Fdf
lF\`nF$F$F$%;classify_4dim_solvable_420GRF+6:F=F>F.F/F0F1F2F3F4F5F6F7F8F9F:F;F<
%#uvGFA%/list_of_coeffsG%1nilradical_checkG%$U_1G%$V_1GFc[nF$F$C4@$FDFG>FcrFO>F
jx-F\o6$FcrFcr>Ffw-Fdo6$FjxFcr-%'assumeG6$%#_aG%%realG-F]\p6$%#_bGF`\p>Fgo-Fa^n
6$-Fcu6$F_\pFhx-Fcu6$Fc\pFgx>Fc\n-F[fl6#-Fbs6%F^x-Fbs6%FaxFjxFgoFjx>Fc\n-%)simp
lifyG6#-%%map2G6%F^fm&7#-%&solveG6$-&Fjo6#%$detGF[en<$F_\pFc\pFLFb^p@%5/Fc\nFfs
0-Fbs6$FdsFc\nFfs>F\zQ(L(4,13)F$>F\zQ'L(4,7)F$-%)unassignG6#.F_\p-F^_p6#.Fc\p@$
0FianFgs>F\zQ'L(4,3)F$@$/F\zFg_pC.>Fiz-%,centralizerGFbil>F^w-FaoFjan>Ffw-Fdo6$
FijmFiz>Fiz7$F^w-Fao6#Ffw>Fcv-Fao6#-Fdo6$FizFcr>Fip-Fio6#-F_p6$FcvFiz?&F\qFipF_
q@$0&F\qFLFgsC%>Fgo-F[v6%-F_s6$-Fao6#&F\qF]rF`vFizFav>Fcv-Fcu6$*$FdapFguFcvF`u>
Fdel-Fao6#-Fdo6$7$F^wFgoFiz>F[[l-F^x6$-Fax6$FcvFdelFgbp>Fdel-Fcu6$*$F\^lFgu-F^^
n6$Fdel-Fcu6$FddlFgo>F_o7&F^wFdelFgoFdcl>Fiu7$F^z7$F`zF^r@$/F\zF\_pC,>Fiz-%8fin
d_nilradical_for_L47GF$>Fiz-Fe\nFi^m>FcvFh`p>F[[l-Fbs6%F^x-Fbs6%FaxFizFcv7$&Fiz
FL&FizFdw>Fcv-F^^n6$-Fa^n6$Fcv-Fcu6$&F[[lFijlFjdp-Fcu6$&F[[l6$FgrFF&FizF]r>Fgep
-F^^n6$-Fcu6$&F[[l6$F^rFgrFjdp-Fcu6$&F[[l6$FgrFgrFgep>Fidp-Fcu6$,$FafpFguFidp>F
cv-Fcu6$*$&F[[l6$FFFFFguFcv>F_o7$-FaoFi^mFcv>Fiu7$F^z7$F`z""(@$/F\zFj^pC0>Fe]n-
F^^n6$-Fcu6$-&Fjo6#%&traceG6#-F_p6$FgxFjxFhx-Fcu6$-F]hp6#-F_p6$FhxFjxFgx>Fip7#-
Fio6#-F_p6$Fe]nFjx>F^w-F[v6%-F_s6$-Fbs6$Fbal-Fao6#&Fip6$FFF^rF`vFjxFav>Fdel-F[v
6%-F_s6$-Fbs6$FdsFfipF`vFjxFav>8;-Fao6#-F^x6$-Fax6$F^wFe]nF`_n>Fe]n-Fcu6$,$*$Fb
jpFguFguFe]n>Fiz7$F^wFdel>F]]n-Fao6#-Fdo6$7#Fe]nFfw>F[[l-F^x6$-Fax6$FidpF]]nFiz
>F]]n-Fcu6$F`cp-Fa^n6$F]]n-Fcu6$FddlFe]n>F[[l-F^x6$-Fax6$F]]nFe]nFiz>Fe]n-F[v6%
7$FF-FaoF]\l7$Fe]nF_gpFav>F_o7%F_gpF]]nFe]n>Fiu7$F^z7$F`z"#8>F^q-Ffel6#-F[fl6#-
Fbs6%F^xF_oFcr-F`fl6#7$FjjoFiuF$F$F$%GLevi_decomposition_for_abelian_radicalGR6
#%$radG64%+parametersG%#npG%"pG%&basisGF2Fbdm%%tempG%#SSG%"SG%"qG%%varsGF=F>%"w
G%*general_SG%#eqG%&comboG%$ansGF$F$C5@$/FEFgr>FZ&FKFdw>FN-Fdp6#FZ>Fjn-Fdp6#Fah
m@$/FjnFgs-F`fl6#7$Fjp-%1frameBaseVectorsGF$>FcrF``q>Fcv-Fdo6$FahmFcr>F^wF]_m>F
del7#-FQ6$-FQ6$&%#_rG6$FgoF^q/Fgo;FFF^w/F^q;FFFjn>F\qFjp?(FgoFFFFF^wF_qC$>Fip-F
a^n6$&FcvFdim-F[v6%7#-FQ6$F]aqFbaqFahmFav>F\q7$FhqFip>Fjx7$-FaoFj_qFhq>F`tFjp?(
FgoFFFFF^wF_q?(F^q,&FgoFFFFFFFFF^wF_qC$>Ff[l-F^x6$-Fax6$&F\qFdim&F\qF\`nFjx>F`t
7$-FaoFht-Fao6#&Ff[l6#Fcaq>F\z-&Fjo6#%*genmatrixG6%F`tFdel%"bG>Ffw-&Fjo6#%)lins
olveG6&F\zF^dq.F\emFgs@$/Ffw%%NULLG-F[dm6#%[oVessiot~Error:~Linsolve~failed~to~
compute~semi-simple~subalgebraG>F_o-Fbs6$%)helmsimpG-F^fm6$<#-FQ6$/&FdelFdim&Ff
wFdim/Fgo;FF*&FjnFFF^wFFF\q>Fiu-Fg]p6%F^fm<#-FQ6$/&%#_tGFdim&FZFdim/Fgo;FF-FdpF
f\nF_o7$FahmFiuF$F$F$%8classify_indecomposableGRF$6#%(alg_dimGF$F$C)@$FDFG>FZFg
n@$/FZFF-F`fl6#7$-F[boFL7$F^z7$FFFF@$/FZFgr-F`fl6#-%.classify_2dimGF$@$/FZF^r-F
`fl6#-%.classify_3dimGF$@$/FZF`z-F`fl6#-F`aoF$-F[dm6#%UThe~dimension~of~the~Lie
~algebra~must~be~less~than~5GF$F$F$%3Levi_decompositionGRF$68%$RadG%"LG%#DSG%'s
ubalgGFbdmFg^q%#LDG%$PhiG%.masterAlgebraG%#LLGF`^q%$phiGFgco%+subalgDataG%"RGF.
%#DRGF[do%,subalg_pairG%*alg_name1G%*alg_name2G%*alg_name3GF$F$CE>FcvFedo>F^qFg
n@$/FceoF_qC$-Feeo6#%-Begin~Step~1G-Feeo6#-FW6$%6Computing~Radical~of~GFcv>FZ-%
(radicalGF$>FNF``q@$/Ff_qFgsC$@$FcjqC$-Feeo6#%7Algebra~is~semi-simpleG-Feeo6#%.
End~of~Step~1G-F`flFgfo@$/Ff_qFgnC$@$FcjqC$-Feeo6#%4Algebra~is~solvableGFi[r-F`
fl6#7$FZFjp@$FcjqC$-Feeo6#%/Begin~~Step~2~G-Feeo6#%9Computing~Derived~SeriesG>F
jn-%/derived_seriesGF$>Fcr-Fbgo6$&Fjn6#FguF``q>F[[l-FW6$Fcv%#_PG@$FcjqC%Fj\r-Fe
eo6#%gnThe~basis~for~the~terminal~algebra~in~the~derived~series~is:G-Feeo6#-%%S
howGF[[n>Fip-F]jo6%FcrF[[l%*inclusionG>FjxFjhm-%-Lie_alg_initG6#Fcdl@$FcjqC$-Fe
eo6#-FW6$%eoThe~terminal~algebra~in~the~~derived~series~is~initialized~as~Lie~a
lgebra~GF[[l-Feeo6#-%7Lie_bracket_mult_tableGF$@$Fcjq-Feeo6#-FW6$F\[rF[[l>F\qF^
[r>Fiu-%@matrix_to_Lie_algebra_transformG6%F[[lFcvFjx>Fdel-Fg]p6%%6Lie_algebra_
transformGFiuF\q@$Fcjq-Feeo6#%,End~~Step~2G?(F`tFgsFFF^qF_qC)@$Fcjq-Feeo6#%-Beg
in~Step~3G>Ff[l-F\o6$FdelFdel@$/-FdpFj`lFgsF`u@$Fcjq-Feeo6#-FW6%%BThe~~radical~
of~the~subalgebra~L_GF`t%0~is~not~AbelianG>F^w-%*Levi_loopG6%FcrFdelF`t>Fcr&F^w
FL>Fdel&F^wFdw@$FcjqC$-Feeo6#-FW6%%,Radical~of~GFc\n%1~is~~now~AbelianG-Feeo6#%
-Begin~Step~4G>F]]n-FW6$Fedo%$_S4G-%0change_frame_toGF^_m>Fiz-%Isubalgebra_pair
_to_Lie_algebra_data_pairG6&FcrF`_nF]]nFh^r>F\qFjdp>FfwFgep@$Fcjq-Feeo6#%QCompu
ting~Semi-Simple~Part~of~Levi~DecompositionG>Fgo-F``r6%F]]nFcvFfw>F^w-F]^qFiq>F
\z-Fg]p6%Fe`rFgoFdbr-%-remove_frameG6%F[[lFc\nF]]n@$Fcjq-Feeo6#%&Done!G7$FZF\zF
$F$F$%3matrix_centralizerGR6#'%1list_of_matricesG-F`v6#%&arrayG6-Fgco%#EQGFbdm%
"XGF__q%%newXGFj^q%"aGFeco%#IdGF\doF$F$C.-%%withG6$Fjo7$%(vectdimG%$rowG>FZ-Fcf
r6#-Fdfr6$F`hmFF>Fcr-F[fl6%FZFZR6$F=F>F$6$%)operatorG%&arrowGF$((%"xGFahmFafmF$
F$F$>F^w-Fger6%%)identityGFbhlFbhl>FN<"?&FjnFahmF_q>FN-%&unionG6$FN-F_s6$-Ff^m6
#,&-Fhbm6$FcrFjnFF-Fhbm6$FjnFcrFguFes>Fjx-F\^p6$FN-F_s6$FcrFes>F\z-Fg]p6%F^fmFj
xFcr>F_o-%*intersectG6$-F_s6$F\zFesF\ir>Fdel7#-Ff^mFe_m?&FcvF_oF_qC&>Fjn-Fbs6%%
%diffGF\zFcv>Fiu-F]hpFigm@$0FiuFgs>Fjn-Ff^m6#,&*&FZFFFjnFFFF*&FiuFFF^wFFFgu@$/-
F^x6$FjnFdelFjp>Fdel7$F`im-Ff^mFigmFdelF$F$F$%+decom_alg1GR6$%/set_cent_mod_JG%
"JG6'F__qFh^qF\doF.F=F$F$C(>FNFahm@$2FgrFgfq>FZF'@$/FgfqFgrC&>Fjn-F^x6$-Ff^m6#*
$&FNFdwFgr7$&FNFL-Fao6#Fafm@$/FegoFgs>FZF'@$32FgsFego2&FjnFLFgs>FZF_q@$3F_]s1Fg
sFa]s>FZF'@$/FgfqFF>FZF'@$/FZF'>FZ7#-FQ6$7$&FNFbil-&Fjo6#%(minpolyG6$Fa^sFiem/F
jx;FgrFgfqFZF$F$F$%6classify_decomposableGRF$6,F=F>%2original_alg_nameG%.decomp
ositionG%,factor_listG%$CBMG%4classification_infoG%+Lie_factorG%-factor_list2G%
0classify_matrixGF$F$C/@$FDFG>FjnFedo@$Feao-F`fl6#7$-F[bo6#-FQ6$FFFen7$F^z7$Fgn
Fgs>Fcr-%2totally_decomposeGF$>FjxF[`o>F\zF]`o>F_oFjp?&FiuFjxF_qC&@$-Fiw6#2F`z&
Fiu6%FFF^rFF-F[dm6#%enAt~least~one~factor~is~too~large~to~classify~at~this~time
.G-F[_r6%Fiu7#%.factor_vectorG7#%,factor_formG>F_o7$-FaoFifm-FifqF$-Fidr6#&F`]l
Fdw>Fcv7#-FQ6$&F_o6$FZFgr/FZ;FF-FdpFifm>F^w-F[bo6#-FQ6$&F_o6$FNFF/FNFfbs>F\z-Ff
^m6#-Fhbm6$F^wF\z-FHFigm-F`fl6#7$-Ff^mF^jlFcvF$F$F$Fc`sRF$65F=F>F.%#n1G%#n2G%#n
3GF\_s%(fr_nameG%'stack1G%'stack2GF__s%#D1G%/abelian_factorG%#B1G%#B2G%#B3G%#B4
G%/reorder_matrixG%-reorder_infoGF$F$C;@$FDFG>FcrFgn>F_oFedo>FcvFjp>F^wFjp>Fdel
-F[bo6#-FQ6$FF/FZ;FFFcr@$Feao-F`fl6#7$7$7%F`jmF_oFihl-%9structure_constants_lis
tGF$-Ff^mFgfl@$-%5check_indecomposableG6#F[fsFhes>FceoF'>Fgo-F_coF$>F^qFg_n@$/F
^qFjp>F^q-Fao6#Fjp>Fdel-Ff^m6#-Fhbm6$Fa`nFdel>Fcv7#&FgoFdw?(F$FFFFF$F`jlC&>FjxF
]_m>F\z&&FcvFbilF^as@%-Fbfs6#FjgsC)>F^w7$Fjgs-FaoFe_m>Fcv7#-FQ6$&FcvFf\n/FN;FF,
&FjxFFFguFF>Fip-F[fl6%,&FcrFFF\zFguF\zFgs>F\q-F[bo6#-FQ6$FF/Fjn;FFF]is>F`t-F[bo
6#-FQ6$FF/Fjn;FFF\z>Ff[l-F[fl6%F\zF]isFgs>Ffw-Ffel6#-&Fjo6#%,blockmatrixG6%FgrF
gr7&FipF\qF`tFf[lC*-F[_r6%FjgsFdasFfas>FiuFedo>Fgo-%2decompose_algebraGF$>Fcv7%
FdhsFg_nFdgs@%0FcrF\z>FipF_is>Fip-F[fl6$FgsFgs>F\qFa`n>Ffw-F[bo6$FipF\q-FidrF\]
l>Fdel-Ff^m6#-Fhbm6$FfwFdel>Fiz-%5reorder_Lie_algebrasGFe_m>F^w7$-Fao6#FidpF^q>
FipFjdp@%0F^qF[gs>F\q-F[bo6#-FQ6$FF/FN;FF&F^qF^as>F\qFe[t>FfwFi[t>FdelF]\t?(Fjn
FFFF-FdpFe_mF_q>&&F^wFigm6$FFFgr-FW6$FfhmFjn-FHFifm>FceoF_q-F`fl6#7$F^wF_fsF$6#
FceoF$Fe\nR6#'%-vector_spaceG-F`v6#%-Vessiot_vectG6&F=F>Fa[nFb[nF$F$C%@$0Fi_qF^
r-F[dm6#%Gthe~subalgebra~must~have~three~vectorsG?(FZFFFF,&Fi_qFFFguFFF_qC%>Fjn
&FahmFg_q?(FN,&FZFFFFFFFFFi_qF_qC$>Fcr&FahmFf\n@$0-%*coeff_setG6#-FaxFhhrFfsF`u
@$F]`tF`u-F`fl6#7%Fa`tFjnFcrF$F$F$%1Jacobson_radicalGR6#'%+matrix_algGFeer6.F\f
rF\do%"dG%"YGF=FierFbdmF__q%%newYGFj^qFe[sFecoF$F$C,>FNFahm>FjnFgfq>Fcr-Ff^m6#*
&&FZFLFFFg\sFF?(FjxFgrFFFjnF_q>Fcr-Ff^m6#,&FcrFF*&&FZFbilFFFa^sFFFF>F\zFjgr?&F_
oFNF_q>F\z-F^hr6$F\z<#-F]hp6#-Ff^m6#-Fhbm6$FcrF_o>Fiu-&Fjo6#%*nullspaceG6#-Fjcq
6$F\z7#-FQ6$F]bt/FjxFcaq>FdelFjp?&FgoFiuF_q>Fdel7$F`im-Ff^m6#-%$addG6$*&&FgoFbi
lFFFa^sFFFectFdelF$F$F$FbfsR6#'F,<$F`jm%%nameG6'Fh^qFe[s%&SmodJGF__qF\doF$F$C)>
FZ-F`er6#-%4list_of_ad_matricesGFj_q>FN-Ff`tFg_q>Fjn-Fdo6$FNFZ@$2FgrFego>FcrF'@
$/FegoFgrC&>Fjx-F^x6$-Ff^m6#*$&FjnFdwFgr7$Fa]sFfho@$/F[hmFgs>FcrF'@$32FgsF[hm2F
j[mFgs>FcrF_q@$3Ffft1FgsFj[m>FcrF'@$/FegoFF>FcrF_qFcrF$F$F$F_brR6%FbiqF`^qF.69%
+master_algGF`iqFjiq%#C1G%#C2G%"QG%'Q_dataG%%barSGFh^qFbdm%%barRG%%phi1GF=F>%*t
rans_eq1GF[jq%(new_radG%)temp_algG%+new_subalgG%&Q_algGF[do%%phi2GFiiqF$F$CB>FZ
Fedo>Ffw-FW6%FZ%%_S4_GFgem@$Fcjq-Feeo6#-FW6&%;Initializing~subalgebra~L_GFgem%1
~as~Lie~algebra~GFfw>F`t-Fhcr6&Fahm7#FafmFfwFh^r>Fe]n&F`tFdw>Fc\n&F`tF]r>F]]n-F
``r6%FfwFZFc\n>FNF``q>Fjn-F\o6$Fe]nFe]n@$FcjqC&-Feeo6#-FW6%%0The~radical~of~GFf
w%%~is:G-Feeo6#-Fd^r6#Fe]n-Feeo6#-FW6%%DThe~derived~algebra~of~~radical~of~GFfw
Fdjt-Feeo6#-Fd^rFigm>Fcr-Fbgo6$-Fdo6$FjnFe]nFN>Fjx-Fbgo6$-Fdo6$Fe]nFNFN>F\z-Fbg
o6$7$-FaoF[[nFailFN@$FcjqC$-Feeo6#%LThe~~complementary~basis~to~the~radical~is:
G-Feeo6#-Fd^rF^jl>F[[l-FW6$Ffw%#_QG>F_o-%8create_quotient_algebraG6%FjnF\zF[[l-
F[_rFifm@$FcjqC$-Feeo6#-FW6&%7The~Factor~Algebra~of~GFfw%0~by~Radical~is~GF[[lF
d_r>F^w-Ffel6#-F[fl6#-%.linear_combo2G6$F\zFN>F\q-FQ6$/(FYF^q-F]dt6$*&&F^w6$F^q
FipFF(FYFipFF/Fip;FF-&Fjo6#%'coldimGFe_m/F^q;FF-&Fjo6#%'rowdimGFe_m>Fgo-%*trans
formG6%F[[lFfw7#F\q>FdelF^[r@$Fcjq-Feeo6#%MCompute~Levi~decomposition~of~Factor
~AlgebraG>Fiu&-F]^qFgflFdw@$FcjqC&-Feeo6#-FW6%%1The~~radical~of~GF[[lFdjt-Feeo6
#-Fd^rFgfl-Feeo6#-FW6%%?The~semi-simple~subalgebra~of~GF[[lFdjt-Feeo6#-Fd^rF\]l
>Fcv-Fg]p6%Fe`rFgoFiu-FidrF]\l>Fiz-Fg]p6%Fe`rF]]n7$Feho-FaoF^_m>Ff[l-Fg]p6%Fe`r
F]]nFjn-FidrFf`p@$FcjqC$-Feeo6#-FW6%%2The~subalgebra~L_G,&FgemFFFFFF%$~isG-Feeo
6#-Fd^rFi^m7$FizFf[lF$F$F$%8classify_3dim_nilpotentGRF+6+F=F>F.F[glF/Fa[nFb[nF2
F8F$F$C)@$FDFggl>FcrFgn>Fjx7#-FQ6$-FT6#(FYFZFees?&F\zFjxF_qC$?&F_oFjxF_q@$0-F_`
t6#-Fax6$F\zF_oFfsF`u@$FbduF`u>Fiu7%FeduF\zF_o>Fcv-Ffel6#-F[fl6#-Fbs6%F^xFiuFjx
-F`fl6#7$-Ff^m6#-FdflF^_m7$F^z7$F^rFFF$F$F$F]hqRF+6)%$_dvG%#_HGF_glF:F5F6F8F$F$
C,@$FDFggl@$Feao-F`fl6#7$-F[fl6#7$7$FFFgs7$FgsFF7$F^z7$FgrFgs>FZ-Fax6$-FT6#%#z1
G-FT6#%#z2G>FN-Fao6#-Fdo6$7#FZ7$F]guF`gu>Fjn-Fio6#-F_p6$FN7$FNFZ?&FcrFjnF_q@$0F
[`oFgsF`u>Fjx-F[v6%-F_s6$&FcrFijlF`vF_huFav>F\z-Fcu6$,$*$F[`oFguFguFN>F_o-Ffel6
#-F[fl6#-Fbs6%F^x7$FjxF\zFigu-F`fl6#7$-Ff^m6#-FdflFifm7$F^z7$FgrFFF$F$F$FchqRF$
6%%$LCSG%,derived_algG%*name_listGF$F$C'@$FDFG@$Feao-F`fl6#7$-F[fl6#7%7%FFFgsFg
sF]cm7%FgsFgsFF7$F^z7$F^rFgs@$Fbbo-F`fl6#-FacuF$@$Fibo-F`fl6#-%7classify_3dim_s
olvableGF$@$-%2check_semi_simpleGF$-F`fl6#-FhflF$F$F$F$Fc\tR6#'%1full_factor_li
stG%)listlistG63F=F>F.%8number_of_total_factorsG%'colnumG%+final_listGF^_sFfhm%
5number_of_subfactorsG%(max_dimGFb_s%(factor2GF\gq%+factor_dimG%6cononical_form
_matrixG%$_ZMG%.matrix_blocksGF$F$C/>FcrFi_q>FjxFgs>F_oFjp?&FiuFahmF_qC$>F_o7$F
jas7$Fiu,&FjxFFFFFF>Fjx,&FjxFFF]asFF>F^qFjx>F\zFjp?(FZFFFFFcrF_qC)>FcvFgbs>F^w-
Fjt6#-FQ6$&&F_oFf\n6&FFFFF^rFF/FN;FFFcv?&FiuF_oF_q@$/&Fiu6%FFFFF^rFijmF`u>FdelF
jp?&FgoF_oF_q@$0FdgsFa]l>Fdel7$F`imFgo>F\z7$FfgmFiu>F_oFdel>FjxFgs>Ff[lFjp?(FZF
FFFFcrF_qC'>Fip&&F\zFg_qFg^v>&F`tFg_q-F[fl6%FipF^qFgs>Fjx&F^`vFdw-FQ6$-%'assign
G6#/.&F``v6$FN,(FjxFFFguFFFNFFFF/FN;FFFip>Ff[l7$-FaoFj`lF``v>F\q-&Fjo6#%,stackm
atrixG6#Fcav>F\z7#-FQ6$&F\z6$FZFFFees-F`fl6#7$F\z-Ff^mFiqF$F$F$FgboRF+6-F=F/F0F
5F6F_[nFa[pFb[p%-gamma_coeffsGF2F8F$F$C0@$FDFG>FNFO>FjnF[o>FcrF_`p>Fjx-Fao6#-Fd
o6$FihlFjn>F\z-Fdo6$7$FcrFjxFN>Fcv-Fbs6$Fao-Fbs6%F^x7$-Fax6$Fjx&F\zFdw-Fax6$Fjx
&F\zFLFihl>F_o-F^^n6$-Fcu6$Fi\mF`dv-Fcu6$&FcvFdwF]dv>Fiu-Fao6#-Fdo6$7#F_oF\z>Fi
u-Fcu6$*$-Fao6#-F^x6$-Fax6$FjxFiuFihlFguFiu>Fiu-Fcu6$*$&-F^x6$Ff^lFdcvFdwFguFiu
>F^w7&FcrFf^lF_oFiu>Fdel-Ffel6#-F[fl6#-Fbs6%F^xF^wFN-F`fl6#7$Fcfl7$F^z7$F`zFFF$
F$F$Fh[vRF+61F[glF/F0F1F2F3F5F6F4F8F9F:F;F<F=F$F$C7@$FDFggl>FZFgn>FN7#-FQ6$-FT6
#(FYF\q/F\qFbhl>FjnF[o>FcvF`o>Fdel-Fio6#-F_p6$Fcv7$FcvFeho>F^qFjp>FipFjp>FcrFjp
?&FgoFdelF_q@$Ff_nC$>F^qFj_n>FipF^`n@$/-FdpF\`nFgr>FcrQ)L(3,2,x)F$@$0-Fd]p6#Ff`
nFfs>FcrQ)L(3,4,x)F$@$3/F]ivFFF_an>FcrQ)L(3,2,1)F$@$3Fhiv/FatFF>FcrQ'L(3,3)F$@$
/FcrFjpC&-FW6$%1derived~algebra:G-%%showGFigm-FW6$%.Complement~H:G-FgjvF^_m-FW6
$%MEigenvectors~~of~~ad(H)~on~~derived~algebra:GF^q-F[dm6#%;Unable~to~classify~
algebraG@$/FcrFjivC&>F_o-F[v6%7#-FQ6$&FjhmFc\mFd\m7$Fa]sF_ftFav>Fiu-F[v6%Fh[w7%
FcvFa]sF_ftFav>Fjx7%F_oFiu-Fcu6$,$*$FeqFguFguFcv>F\z6$7$F^z7$F^rFgrFbcm@$/FcrF_
iv@%-Fiw6#2-Fet6#F\jm-Fet6#FeqC&>F_o-F[v6%FacnF`\wFav>Fiu-F[v6%FfcnF`\wFav>F\z6
$Fi\w7#*&F\jmFFFeqFgu>FjxFb\wC&>F_oFj]w>FiuFg]w>F\z6$Fi\w7#*&FeqFFF\jmFgu>Fjx7%
F_oFiu-Fcu6$,$*$F\jmFguFguFcv@$/FcrF_jvC'>F_o-F[v6%7#-FQ6$&FcdlFg[m/F>;FgrFZFjn
Fav>Fiu-Fao6#-Fdo6$F^evFjn>Fiu-Fcu6$*&FeqFF&-F^x6$-Fax6$FcvFiuF]dnFLFguFiu>FjxF
b\w>F\z7$F^z7$F^rF^r@$/FcrFeivC)@$-Fiw6#2-FdsFd]wFgsC$>FeqF\jm>FcdlFjhm>F_o-F[v
6%Fe[oFbhvFav>Fiu-F[v6%F\\oFbhvFav>Fcv-Fcu6$,$*$FcawFguFguFcv@$-Fiw6#1Fgs-FbalF
d]w>FjxF]dn@$-Fiw6#2FfbwFgs>Fjx7%FiuF_oFdcl>F\z6$7$F^z7$F^rF`z7#-Fet6#*&FfbwFFF
cawFgu>F^wFeel7$-Ff^m6#F_aoF\zF$F$F$F^coRF+6#FaiqF$F$C*@$FDFG>FceoF'>FZ-Fbs6$Fd
pFa]r>FceoF_q@$/FZ7%F`zF^rFgs-F`fl6#-F)F$@$/FZ7&F`zF^rFFFgs-F`fl6#-F\[nF$@$/FZ7
%F`zFgrFgs-F`fl6#-F[[pF$-F[dm6$%?The~algebra~is~not~in~the~listGFZF$Fc^tF$%>che
ck_absolute_indecomposableGRFbdt6)%(Adj_repGF\doFd^qF__qF=FbdmF[doF$F$C)>FZF\et
>FN-F`erFg_q>FN7#-FQ6$Fa^sFg^s>FjnFgfq>FcrF_q?&F\zFNF_q?&F_oFNF_q@$0-F]hp6#-Ff^
m6#-FhbmFfduFgs>FcrF'FcrF$F$F$F_dpRF$6*F=F/%#g1G%1complement_of_g1GF_\pFc\p%-th
ird_vectorG%%_eqsGF$F$C+>FNFO>FjnF[o>FcrFco>Fjx-%&evalnGFbil>F\z-F\hwF^jl>F_o-F
a^n6$-Fcu6$FjxF[`o-Fcu6$F\zF]`o>Fiu-F_s6$-Fbs6$Fao-Fbs6%F^x-Fbs6%Fax-Fbs6%FaxFj
nF_oF_oFjnFes>F_o-F_eq6#-F^fm6$<$/FjxFF/F\zFF-F^fm6$-F\^p6$Fiu<$FjxF\zF_o-F`fl6
#7$FehoF_oF$F$F$%5classify_Lie_algebraGRF$F$F$F$C%@$FDFG@$-Fbfs6#7$7%F`jm%)init
_algG7#FgnF]fs-F`fl6#F[bs-F`fl6#-Fi^sF$F$F$F$F^[tR6#%(lie_algG6,Fh^qFgcoFe[sFgd
tF[doF__qF,F^doFhco%-lie_alg_dataGF$F$C2@$Fceo-Feeo6#%QWarning:~Be~sure~to~remo
ve~central~factors~firstG@$/FEFgs>F^w-%ELie_algebra_name_to_Lie_algebra_dataGF$
@$3FD-%%typeG6$FJFedtC$FG>F^wF^\x@$3FD-Fc\x6$FJF`jm>F^wFahm@$/FEF`zC$-F^dm6&F^w
Fd_q&FKF]r&FKFe_l-F`fl6#%"%G@$Fceo-Feeo6#%inComputing~Set~Centralizer~and~Jacob
som~radical~for~ad~matricesG>FZ-F`er6#-F]etFe_m>FjnF_et>Fcr-Fdo6$FjnFZ>F\z-Fa[s
Ffhr@$/F\zF_q-F`flFifm@$Fceo-Feeo6#%ioFactoring~characteristic~polynomials~for~
~set~centralizer~mod~Jacobson~radicalG>Fjx-%+decom_alg2GF^jl@$/FjxFjpC&-Feeo6#%
8Koszul~needs~~you~help.G-Feeo6#%^oCompute~prime~factors~of~one~~of~the~~follow
ingminimal~polynomials.G-Feeo6#%`pRerun~~decompose_algebra~with~matrix~and~prim
e~factors~~ofcharacteristic~~polynomial~G-F`flF^jl@$Fceo-Feeo6#%OComputing~chan
ge~of~basis~to~decompose~algebraG-F^dm6&F^wFj[mF[\mFc_lF$F$F$F]_xR6#%&SpolyG6(F
=%"mGFgco%"fGFcdmFddmF$F$C(/FjnFF?(FZFFFFFi_qF_qC$>Fcr-%(factorsG6#&Ff_tFdw@$2F
F-Fdp6#F]`oC%>FNFZ/FjnFgrF`u@$Fh`x-F`flF\gs>Fjx)&&F]`oFLFL&F\bxFdw>F\z-Fcfm6#-F
fhm6#*&&F[`tFdwFFFjxFgu7%&F[`tFLFjxF\zF$F$F$F$,
I$mfl%Ee:\mathware\vessiot\programs\mfiles\GF$