#!/usr/bin/env cadabra2
{\mu,\nu,\rho}::Indices(curved, position=fixed).
{m,n,p,q,r,s,t,u,v}::Indices(flat, position=independent).
{m,n,p,q,r,s,t,u,v}::Integer(0..10).
T^{#{\mu}}::AntiSymmetric.
\psi_{\mu}::SelfAntiCommuting.
\psi_{\mu}::Spinor(dimension=11, type=Majorana).
\theta::Spinor(dimension=11, type=Majorana).
\epsilon::Spinor(dimension=11, type=Majorana).
{\theta,\epsilon,\psi_{\mu}}::AntiCommuting.
\bar{#}::DiracBar.
\delta^{m n}::KroneckerDelta.
\Gamma^{#{m}}::GammaMatrix(metric=\delta).
obj:= T^{\mu\nu\rho} e_{\nu}^{s}
     \bar{\theta} \Gamma^{r s} \psi_{\rho}
     \bar{\psi_{\mu}} \Gamma^{r} \epsilon;
fierz(_, $\theta, \epsilon, \psi_{\mu}, \psi_{\rho}$ );
converge(obj):
   join_gamma(_)
   distribute(_)
   eliminate_kronecker(_)
canonicalise(_)
rename_dummies(_);
tst:=1/4 T^{\mu \nu \rho} e_{\mu}^{m} \bar{\theta} \Gamma^{m n}
\epsilon \bar{\psi_{\nu}} \Gamma_{n} \psi_{\rho} 
+ 5/16 T^{\mu \nu \rho} e_{\mu}^{m} \bar{\theta} \epsilon
\bar{\psi_{\nu}} \Gamma_{m} \psi_{\rho} 
+ 3/32 T^{\mu \nu \rho} e_{\mu}^{m} \bar{\theta} \Gamma^{m n p} \epsilon
\bar{\psi_{\nu}} \Gamma_{n p} \psi_{\rho} 
+ 1/4 T^{\mu \nu \rho} e_{\mu}^{m} \bar{\theta} \Gamma^{n}
\epsilon \bar{\psi_{\nu}} \Gamma_{m n} \psi_{\rho} 
+ 1/384 T^{\mu \nu \rho} e_{\mu}^{m} \bar{\theta}
\Gamma^{n p q r} \epsilon \bar{\psi_{\nu}} \Gamma_{m n p q r} \psi_{\rho};