notebook.community

Edit and run

Load pywicked, setup orbital spacesm, and create WickTheorem object





In [1]:



    


import pywicked as pw
from IPython.display import display, Math, Latex

pw.osi.reset()
pw.osi.add_space("v",pw.rdmtype.unoccupied,["a","b"])
pw.osi.add_space("o",pw.rdmtype.occupied,["i","j"])
wdt = pw.WDiagTheorem()

Create operators





In [2]:



    


T = pw.make_operator("T",["o->v","oo->vv"])
F = pw.make_operator("F", ["o->o","o->v","v->o","v->v"]);
V = pw.make_operator("V", ["oo->oo", "oo->ov", "oo->vv",
                           "ov->oo", "ov->ov", "ov->vv",
                           "vv->oo", "vv->ov", "vv->vv"]);
H = pw.make_operator("v", ["o->o"  , "o->v"  , "v->o"  , "v->v",
                           "oo->oo", "oo->ov", "oo->vv",
                           "ov->oo", "ov->ov", "ov->vv",
                           "vv->oo", "vv->ov", "vv->vv"]);

Setup the similarity-transformed Hamiltonian and compute expectation value





In [3]:



    


Hbar = pw.bch_series(H,T,4)
sum = wdt.contract_sum(pw.rational(1,1),Hbar, 0, 0);
display(Math("E = " + sum.latex()))



    


















    







$$E = \frac{1}{2} T^{i}_{a} T^{j}_{b} v^{a b}_{i j} \\ 
 T^{i}_{a} v^{a}_{i} \\ 
\frac{1}{4} T^{i j}_{a b} v^{a b}_{i j}$$

Create ambit equations





In [4]:



    


eqs = sum.to_manybody_equation("R")

print "Generated %d equations" % len(eqs)

for eq in eqs:
    print eq.ambit()



    


















    






Generated 3 equations
R += 0.500000 * T[o0,v0] * T[o1,v1] * v[v0,v1,o0,o1];
R += 1.000000 * T[o0,v0] * v[v0,o0];
R += 0.250000 * T[o0,o1,v0,v1] * v[v0,v1,o0,o1];

Compute the CCSD similarity-transformed Hamiltonian





In [6]:



    


Hbar = pw.bch_series(H,T,4)
sum = wdt.contract_sum(pw.rational(1,1),Hbar, 2, 4)

print sum.str()

eqs = sum.to_manybody_equation("R")

print "Generated %d equations" % len(eqs)

n1 = 0
n2 = 0
for eq in eqs:
    ambit = eq.ambit()
    if ("R[o0,v0]" in ambit):
        n1 = n1 + 1
    if ("R[o1,o0,v0,v1]" in ambit):
        n2 = n2 + 1
        print ambit     
        
print "n1 = %d" % n1
print "n2 = %d" % n2



    


















    






1/2 T^{o1}_{v3} T^{o0,o2}_{v0,v1} v^{v2,v3}_{o1,o2} { a+(v0) a+(v1) a-(v2) a-(o0) }
-1/4 T^{o0}_{v3} T^{o1,o2}_{v0,v1} v^{v2,v3}_{o1,o2} { a+(v0) a+(v1) a-(v2) a-(o0) }
-1/2 T^{o0}_{v3} v^{v2,v3}_{v0,v1} { a+(v0) a+(v1) a-(v2) a-(o0) }
-1/2 T^{o2}_{v2} T^{o0,o1}_{v0,v3} v^{v2,v3}_{v1,o2} { a+(v0) a+(v1) a-(o1) a-(o0) }
-1/2 T^{o2}_{v2} T^{o0,o3}_{v0,v1} v^{v2,o1}_{o2,o3} { a+(v0) a+(v1) a-(o1) a-(o0) }
-1/2 T^{o0}_{v2} T^{o2}_{v3} T^{o1,o3}_{v0,v1} v^{v2,v3}_{o2,o3} { a+(v0) a+(v1) a-(o1) a-(o0) }
1/8 T^{o0}_{v2} T^{o1}_{v3} T^{o2,o3}_{v0,v1} v^{v2,v3}_{o2,o3} { a+(v0) a+(v1) a-(o1) a-(o0) }
1/4 T^{o0}_{v2} T^{o1}_{v3} v^{v2,v3}_{v0,v1} { a+(v0) a+(v1) a-(o1) a-(o0) }
- T^{o0}_{v2} T^{o1,o2}_{v0,v3} v^{v2,v3}_{v1,o2} { a+(v0) a+(v1) a-(o1) a-(o0) }
1/4 T^{o0}_{v2} T^{o2,o3}_{v0,v1} v^{v2,o1}_{o2,o3} { a+(v0) a+(v1) a-(o1) a-(o0) }
1/2 T^{o0}_{v2} T^{o1,o2}_{v0,v1} v^{v2}_{o2} { a+(v0) a+(v1) a-(o1) a-(o0) }
- T^{o0}_{v2} v^{v1,v2}_{v0,o1} { a+(v0) a+(o1) a-(v1) a-(o0) }
 T^{o0}_{v2} v^{v1,v2}_{v0,o0} { a+(v0) a-(v1) }
1/2 T^{o0}_{v2} v^{v2,o1}_{v0,v1} { a+(v0) a+(v1) a-(o1) a-(o0) }
1/8 T^{o0,o1}_{v2,v3} v^{v2,v3}_{v0,v1} { a+(v0) a+(v1) a-(o1) a-(o0) }
-1/2 T^{o3}_{v1} T^{o0,o1}_{v0,v2} v^{v1,v2}_{o2,o3} { a+(v0) a+(o2) a-(o1) a-(o0) }
 T^{o1}_{v1} T^{o0,o2}_{v0,v2} v^{v1,v2}_{o1,o2} { a+(v0) a-(o0) }
 T^{o1}_{v1} v^{v0,v1}_{o0,o1} { a+(o0) a-(v0) }
- T^{o1}_{v1} v^{v1,o0}_{v0,o1} { a+(v0) a-(o0) }
1/2 T^{o0}_{v1} T^{o1}_{v2} v^{v1,v2}_{v0,o2} { a+(v0) a+(o2) a-(o1) a-(o0) }
 T^{o0}_{v1} T^{o1}_{v2} v^{v1,v2}_{v0,o1} { a+(v0) a-(o0) }
- T^{o0}_{v1} T^{o1,o3}_{v0,v2} v^{v1,v2}_{o2,o3} { a+(v0) a+(o2) a-(o1) a-(o0) }
-1/2 T^{o0}_{v1} T^{o1,o2}_{v0,v2} v^{v1,v2}_{o1,o2} { a+(v0) a-(o0) }
-1/2 T^{o0}_{v1} v^{v0,v1}_{o1,o2} { a+(o1) a+(o2) a-(v0) a-(o0) }
 T^{o0}_{v1} v^{v1}_{v0} { a+(v0) a-(o0) }
 T^{o0}_{v1} v^{v1,o2}_{v0,o1} { a+(v0) a+(o1) a-(o2) a-(o0) }
1/4 T^{o0,o1}_{v1,v2} v^{v1,v2}_{v0,o2} { a+(v0) a+(o2) a-(o1) a-(o0) }
1/2 T^{o0,o1}_{v1,v2} v^{v1,v2}_{v0,o1} { a+(v0) a-(o0) }
1/2 T^{o3}_{v0} T^{o0}_{v1} T^{o1}_{v2} v^{v1,v2}_{o2,o3} { a+(v0) a+(o2) a-(o1) a-(o0) }
 T^{o3}_{v0} T^{o0}_{v1} v^{v1,o2}_{o1,o3} { a+(v0) a+(o1) a-(o2) a-(o0) }
1/4 T^{o3}_{v0} T^{o0,o1}_{v1,v2} v^{v1,v2}_{o2,o3} { a+(v0) a+(o2) a-(o1) a-(o0) }
1/2 T^{o3}_{v0} v^{o1,o2}_{o0,o3} { a+(v0) a+(o0) a-(o2) a-(o1) }
-1/2 T^{o2}_{v0} T^{o3}_{v2} T^{o0,o1}_{v1,v3} v^{v2,v3}_{o2,o3} { a+(v0) a+(v1) a-(o1) a-(o0) }
1/2 T^{o2}_{v0} T^{o0}_{v2} T^{o1}_{v3} v^{v2,v3}_{v1,o2} { a+(v0) a+(v1) a-(o1) a-(o0) }
- T^{o2}_{v0} T^{o0}_{v2} T^{o1,o3}_{v1,v3} v^{v2,v3}_{o2,o3} { a+(v0) a+(v1) a-(o1) a-(o0) }
- T^{o2}_{v0} T^{o0}_{v2} v^{v1,v2}_{o1,o2} { a+(v0) a+(o1) a-(v1) a-(o0) }
 T^{o2}_{v0} T^{o0}_{v2} v^{v2,o1}_{v1,o2} { a+(v0) a+(v1) a-(o1) a-(o0) }
1/4 T^{o2}_{v0} T^{o0,o1}_{v2,v3} v^{v2,v3}_{v1,o2} { a+(v0) a+(v1) a-(o1) a-(o0) }
1/4 T^{o2}_{v0} T^{o3}_{v1} T^{o0}_{v2} T^{o1}_{v3} v^{v2,v3}_{o2,o3} { a+(v0) a+(v1) a-(o1) a-(o0) }
1/2 T^{o2}_{v0} T^{o3}_{v1} T^{o0}_{v2} v^{v2,o1}_{o2,o3} { a+(v0) a+(v1) a-(o1) a-(o0) }
1/8 T^{o2}_{v0} T^{o3}_{v1} T^{o0,o1}_{v2,v3} v^{v2,v3}_{o2,o3} { a+(v0) a+(v1) a-(o1) a-(o0) }
1/4 T^{o2}_{v0} T^{o3}_{v1} v^{o0,o1}_{o2,o3} { a+(v0) a+(v1) a-(o1) a-(o0) }
- T^{o2}_{v0} T^{o0,o3}_{v1,v2} v^{v2,o1}_{o2,o3} { a+(v0) a+(v1) a-(o1) a-(o0) }
1/2 T^{o2}_{v0} T^{o0,o1}_{v1,v2} v^{v2}_{o2} { a+(v0) a+(v1) a-(o1) a-(o0) }
- T^{o2}_{v0} v^{v1,o1}_{o0,o2} { a+(v0) a+(o0) a-(v1) a-(o1) }
- T^{o2}_{v0} v^{v0,o1}_{o0,o2} { a+(o0) a-(o1) }
1/2 T^{o2}_{v0} v^{o0,o1}_{v1,o2} { a+(v0) a+(v1) a-(o1) a-(o0) }
- T^{o1}_{v0} T^{o0}_{v3} v^{v2,v3}_{v1,o1} { a+(v0) a+(v1) a-(v2) a-(o0) }
-1/2 T^{o1}_{v0} T^{o2}_{v1} T^{o0}_{v3} v^{v2,v3}_{o1,o2} { a+(v0) a+(v1) a-(v2) a-(o0) }
-1/2 T^{o1}_{v0} T^{o2}_{v1} v^{v2,o0}_{o1,o2} { a+(v0) a+(v1) a-(v2) a-(o0) }
 T^{o1}_{v0} T^{o2}_{v1} v^{v1,o0}_{o1,o2} { a+(v0) a-(o0) }
- T^{o1}_{v0} T^{o0}_{v1} T^{o2}_{v2} v^{v1,v2}_{o1,o2} { a+(v0) a-(o0) }
- T^{o1}_{v0} T^{o0}_{v1} v^{v1}_{o1} { a+(v0) a-(o0) }
 T^{o1}_{v0} T^{o0,o2}_{v1,v3} v^{v2,v3}_{o1,o2} { a+(v0) a+(v1) a-(v2) a-(o0) }
-1/2 T^{o1}_{v0} T^{o0,o2}_{v1,v2} v^{v1,v2}_{o1,o2} { a+(v0) a-(o0) }
- T^{o1}_{v0} v^{o0}_{o1} { a+(v0) a-(o0) }
1/2 T^{o1}_{v0} v^{v1,v2}_{o0,o1} { a+(v0) a+(o0) a-(v2) a-(v1) }
- T^{o1}_{v0} v^{v2,o0}_{v1,o1} { a+(v0) a+(v1) a-(v2) a-(o0) }
- T^{o0}_{v0} T^{o1}_{v2} v^{v1,v2}_{o0,o1} { a+(v0) a-(v1) }
 T^{o0}_{v0} T^{o2}_{v1} v^{v0,v1}_{o1,o2} { a+(o1) a-(o0) }
1/4 T^{o0}_{v0} T^{o1}_{v1} v^{v0,v1}_{o2,o3} { a+(o2) a+(o3) a-(o1) a-(o0) }
1/4 T^{o0}_{v0} T^{o1}_{v1} v^{v2,v3}_{o0,o1} { a+(v0) a+(v1) a-(v3) a-(v2) }
 T^{o0}_{v0} v^{v0}_{o1} { a+(o1) a-(o0) }
1/2 T^{o0}_{v0} v^{v0,o3}_{o1,o2} { a+(o1) a+(o2) a-(o3) a-(o0) }
- T^{o0}_{v0} v^{v1}_{o0} { a+(v0) a-(v1) }
1/2 T^{o0}_{v0} v^{v2,v3}_{v1,o0} { a+(v0) a+(v1) a-(v3) a-(v2) }
 T^{o0,o1}_{v0,v3} v^{v2,v3}_{v1,o1} { a+(v0) a+(v1) a-(v2) a-(o0) }
1/2 T^{o0,o2}_{v0,v2} T^{o1,o3}_{v1,v3} v^{v2,v3}_{o2,o3} { a+(v0) a+(v1) a-(o1) a-(o0) }
 T^{o0,o2}_{v0,v2} v^{v1,v2}_{o1,o2} { a+(v0) a+(o1) a-(v1) a-(o0) }
- T^{o0,o2}_{v0,v2} v^{v2,o1}_{v1,o2} { a+(v0) a+(v1) a-(o1) a-(o0) }
-1/4 T^{o0,o1}_{v0,v2} T^{o2,o3}_{v1,v3} v^{v2,v3}_{o2,o3} { a+(v0) a+(v1) a-(o1) a-(o0) }
-1/2 T^{o0,o1}_{v0,v2} v^{v1,v2}_{o0,o1} { a+(v0) a-(v1) }
1/2 T^{o0,o1}_{v0,v2} v^{v2}_{v1} { a+(v0) a+(v1) a-(o1) a-(o0) }
1/16 T^{o2,o3}_{v0,v1} T^{o0,o1}_{v2,v3} v^{v2,v3}_{o2,o3} { a+(v0) a+(v1) a-(o1) a-(o0) }
1/8 T^{o2,o3}_{v0,v1} v^{o0,o1}_{o2,o3} { a+(v0) a+(v1) a-(o1) a-(o0) }
-1/4 T^{o1,o2}_{v0,v1} v^{v2,o0}_{o1,o2} { a+(v0) a+(v1) a-(v2) a-(o0) }
1/2 T^{o1,o2}_{v0,v1} v^{v1,o0}_{o1,o2} { a+(v0) a-(o0) }
- T^{o0,o3}_{v0,v1} v^{v1,o2}_{o1,o3} { a+(v0) a+(o1) a-(o2) a-(o0) }
-1/4 T^{o0,o2}_{v0,v1} T^{o1,o3}_{v2,v3} v^{v2,v3}_{o2,o3} { a+(v0) a+(v1) a-(o1) a-(o0) }
-1/2 T^{o0,o2}_{v0,v1} v^{o1}_{o2} { a+(v0) a+(v1) a-(o1) a-(o0) }
1/2 T^{o0,o2}_{v0,v1} v^{v0,v1}_{o1,o2} { a+(o1) a-(o0) }
1/2 T^{o0,o1}_{v0,v1} v^{v1}_{o2} { a+(v0) a+(o2) a-(o1) a-(o0) }
1/8 T^{o0,o1}_{v0,v1} v^{v0,v1}_{o2,o3} { a+(o2) a+(o3) a-(o1) a-(o0) }
-1/2 T^{o0,o1}_{v0,v1} v^{v2}_{o1} { a+(v0) a+(v1) a-(v2) a-(o0) }
 T^{o0,o1}_{v0,v1} v^{v1}_{o1} { a+(v0) a-(o0) }
1/8 T^{o0,o1}_{v0,v1} v^{v2,v3}_{o0,o1} { a+(v0) a+(v1) a-(v3) a-(v2) }
 v^{o1}_{o0} { a+(o0) a-(o1) }
 v^{v0}_{o0} { a+(o0) a-(v0) }
1/4 v^{o2,o3}_{o0,o1} { a+(o0) a+(o1) a-(o3) a-(o2) }
-1/2 v^{v0,o2}_{o0,o1} { a+(o0) a+(o1) a-(v0) a-(o2) }
1/4 v^{v0,v1}_{o0,o1} { a+(o0) a+(o1) a-(v1) a-(v0) }
 v^{o0}_{v0} { a+(v0) a-(o0) }
 v^{v1}_{v0} { a+(v0) a-(v1) }
1/2 v^{o1,o2}_{v0,o0} { a+(v0) a+(o0) a-(o2) a-(o1) }
- v^{v1,o1}_{v0,o0} { a+(v0) a+(o0) a-(v1) a-(o1) }
1/2 v^{v1,v2}_{v0,o0} { a+(v0) a+(o0) a-(v2) a-(v1) }
1/4 v^{o0,o1}_{v0,v1} { a+(v0) a+(v1) a-(o1) a-(o0) }
-1/2 v^{v2,o0}_{v0,v1} { a+(v0) a+(v1) a-(v2) a-(o0) }
1/4 v^{v2,v3}_{v0,v1} { a+(v0) a+(v1) a-(v3) a-(v2) }
Generated 99 equations
R[o1,o0,v0,v1] += -2.000000 * T[o2,v2] * T[o0,o1,v0,v3] * v[v2,v3,v1,o2];
R[o1,o0,v0,v1] += -2.000000 * T[o2,v2] * T[o0,o3,v0,v1] * v[v2,o1,o2,o3];
R[o1,o0,v0,v1] += -2.000000 * T[o0,v2] * T[o2,v3] * T[o1,o3,v0,v1] * v[v2,v3,o2,o3];
R[o1,o0,v0,v1] += 0.500000 * T[o0,v2] * T[o1,v3] * T[o2,o3,v0,v1] * v[v2,v3,o2,o3];
R[o1,o0,v0,v1] += 1.000000 * T[o0,v2] * T[o1,v3] * v[v2,v3,v0,v1];
R[o1,o0,v0,v1] += -4.000000 * T[o0,v2] * T[o1,o2,v0,v3] * v[v2,v3,v1,o2];
R[o1,o0,v0,v1] += 1.000000 * T[o0,v2] * T[o2,o3,v0,v1] * v[v2,o1,o2,o3];
R[o1,o0,v0,v1] += 2.000000 * T[o0,v2] * T[o1,o2,v0,v1] * v[v2,o2];
R[o1,o0,v0,v1] += 2.000000 * T[o0,v2] * v[v2,o1,v0,v1];
R[o1,o0,v0,v1] += 0.500000 * T[o0,o1,v2,v3] * v[v2,v3,v0,v1];
R[o1,o0,v0,v1] += -2.000000 * T[o2,v0] * T[o3,v2] * T[o0,o1,v1,v3] * v[v2,v3,o2,o3];
R[o1,o0,v0,v1] += 2.000000 * T[o2,v0] * T[o0,v2] * T[o1,v3] * v[v2,v3,v1,o2];
R[o1,o0,v0,v1] += -4.000000 * T[o2,v0] * T[o0,v2] * T[o1,o3,v1,v3] * v[v2,v3,o2,o3];
R[o1,o0,v0,v1] += 4.000000 * T[o2,v0] * T[o0,v2] * v[v2,o1,v1,o2];
R[o1,o0,v0,v1] += 1.000000 * T[o2,v0] * T[o0,o1,v2,v3] * v[v2,v3,v1,o2];
R[o1,o0,v0,v1] += 1.000000 * T[o2,v0] * T[o3,v1] * T[o0,v2] * T[o1,v3] * v[v2,v3,o2,o3];
R[o1,o0,v0,v1] += 2.000000 * T[o2,v0] * T[o3,v1] * T[o0,v2] * v[v2,o1,o2,o3];
R[o1,o0,v0,v1] += 0.500000 * T[o2,v0] * T[o3,v1] * T[o0,o1,v2,v3] * v[v2,v3,o2,o3];
R[o1,o0,v0,v1] += 1.000000 * T[o2,v0] * T[o3,v1] * v[o0,o1,o2,o3];
R[o1,o0,v0,v1] += -4.000000 * T[o2,v0] * T[o0,o3,v1,v2] * v[v2,o1,o2,o3];
R[o1,o0,v0,v1] += 2.000000 * T[o2,v0] * T[o0,o1,v1,v2] * v[v2,o2];
R[o1,o0,v0,v1] += 2.000000 * T[o2,v0] * v[o0,o1,v1,o2];
R[o1,o0,v0,v1] += 2.000000 * T[o0,o2,v0,v2] * T[o1,o3,v1,v3] * v[v2,v3,o2,o3];
R[o1,o0,v0,v1] += -4.000000 * T[o0,o2,v0,v2] * v[v2,o1,v1,o2];
R[o1,o0,v0,v1] += -1.000000 * T[o0,o1,v0,v2] * T[o2,o3,v1,v3] * v[v2,v3,o2,o3];
R[o1,o0,v0,v1] += 2.000000 * T[o0,o1,v0,v2] * v[v2,v1];
R[o1,o0,v0,v1] += 0.250000 * T[o2,o3,v0,v1] * T[o0,o1,v2,v3] * v[v2,v3,o2,o3];
R[o1,o0,v0,v1] += 0.500000 * T[o2,o3,v0,v1] * v[o0,o1,o2,o3];
R[o1,o0,v0,v1] += -1.000000 * T[o0,o2,v0,v1] * T[o1,o3,v2,v3] * v[v2,v3,o2,o3];
R[o1,o0,v0,v1] += -2.000000 * T[o0,o2,v0,v1] * v[o1,o2];
R[o1,o0,v0,v1] += 1.000000 * v[o0,o1,v0,v1];
n1 = 14
n2 = 31

Compute the CCSDTQ similarity-transformed Hamiltonian





In [7]:



    


T = pw.make_operator("T",["o->v","oo->vv","ooo->vvv","oooo->vvvv"])

Hbar = pw.bch_series(H,T,4)
sum = wdt.contract_sum(pw.rational(1,1),Hbar, 2, 8)

eqs = sum.to_manybody_equation("R")

print "Generated %d equations" % len(eqs)

numr1 = 0
numr2 = 0
for eq in eqs:
    print eq.ambit()



    


















    






Generated 339 equations
R[v4,o2,o1,o0,v0,v1,v2,v3] += 1.000000 * T[o3,v5] * T[o0,o1,o2,o4,v0,v1,v2,v3] * v[v4,v5,o3,o4];
R[v4,o2,o1,o0,v0,v1,v2,v3] += 18.000000 * T[o0,v5] * T[o1,o3,v0,v1] * T[o2,o4,v2,v3] * v[v4,v5,o3,o4];
R[v4,o2,o1,o0,v0,v1,v2,v3] += -12.000000 * T[o0,v5] * T[o1,o2,o3,v0,v1,v2] * v[v4,v5,v3,o3];
R[v4,o2,o1,o0,v0,v1,v2,v3] += -1.500000 * T[o0,v5] * T[o1,o2,o3,o4,v0,v1,v2,v3] * v[v4,v5,o3,o4];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -24.000000 * T[o4,v4] * T[o0,o1,v0,v5] * T[o2,o3,o5,v1,v2,v3] * v[v4,v5,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -24.000000 * T[o4,v4] * T[o0,o5,v0,v1] * T[o1,o2,o3,v2,v3,v5] * v[v4,v5,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -4.000000 * T[o4,v4] * T[o0,o1,o2,o3,v0,v1,v2,v5] * v[v4,v5,v3,o4];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -4.000000 * T[o4,v4] * T[o0,o1,o2,o5,v0,v1,v2,v3] * v[v4,o3,o4,o5];
R[v3,o1,o0,v0,v1,v2] += 1.000000 * T[o2,v4] * T[o0,o1,o3,v0,v1,v2] * v[v3,v4,o2,o3];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -4.000000 * T[o0,v4] * T[o4,v5] * T[o1,o2,o3,o5,v0,v1,v2,v3] * v[v4,v5,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -36.000000 * T[o0,v4] * T[o1,v5] * T[o2,o4,v0,v1] * T[o3,o5,v2,v3] * v[v4,v5,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 24.000000 * T[o0,v4] * T[o1,v5] * T[o2,o3,o4,v0,v1,v2] * v[v4,v5,v3,o4];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 3.000000 * T[o0,v4] * T[o1,v5] * T[o2,o3,o4,o5,v0,v1,v2,v3] * v[v4,v5,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -48.000000 * T[o0,v4] * T[o1,o4,v0,v5] * T[o2,o3,o5,v1,v2,v3] * v[v4,v5,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -24.000000 * T[o0,v4] * T[o1,o2,v0,v5] * T[o3,o4,o5,v1,v2,v3] * v[v4,v5,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 12.000000 * T[o0,v4] * T[o4,o5,v0,v1] * T[o1,o2,o3,v2,v3,v5] * v[v4,v5,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -144.000000 * T[o0,v4] * T[o1,o4,v0,v1] * T[o2,o3,v2,v5] * v[v4,v5,v3,o4];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -72.000000 * T[o0,v4] * T[o1,o4,v0,v1] * T[o2,o5,v2,v3] * v[v4,o3,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 72.000000 * T[o0,v4] * T[o1,o4,v0,v1] * T[o2,o3,o5,v2,v3,v5] * v[v4,v5,o4,o5];
R[v3,o1,o0,v0,v1,v2] += 6.000000 * T[o0,v4] * T[o1,o2,v0,v1] * v[v3,v4,v2,o2];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 24.000000 * T[o0,v4] * T[o1,o2,o3,v0,v1,v5] * v[v4,v5,v2,v3];
R[v3,o2,o1,o0,v0,v1,v2,o3] += -3.000000 * T[o0,v4] * T[o1,o2,o4,v0,v1,v2] * v[v3,v4,o3,o4];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 48.000000 * T[o0,v4] * T[o1,o2,o4,v0,v1,v2] * v[v4,o3,v3,o4];
R[v3,o1,o0,v0,v1,v2] += 1.000000 * T[o0,v4] * T[o1,o2,o3,v0,v1,v2] * v[v3,v4,o2,o3];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -16.000000 * T[o0,v4] * T[o1,o2,o3,o4,v0,v1,v2,v5] * v[v4,v5,v3,o4];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 6.000000 * T[o0,v4] * T[o1,o2,o4,o5,v0,v1,v2,v3] * v[v4,o3,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 4.000000 * T[o0,v4] * T[o1,o2,o3,o4,v0,v1,v2,v3] * v[v4,o4];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -2.000000 * T[o0,o4,v4,v5] * T[o1,o2,o3,o5,v0,v1,v2,v3] * v[v4,v5,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 12.000000 * T[o0,o1,v4,v5] * T[o2,o3,o4,v0,v1,v2] * v[v4,v5,v3,o4];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 1.500000 * T[o0,o1,v4,v5] * T[o2,o3,o4,o5,v0,v1,v2,v3] * v[v4,v5,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,o4] += -1.000000 * T[o5,v3] * T[o0,o1,o2,o3,v0,v1,v2,v4] * v[v3,v4,o4,o5];
R[o2,o1,o0,v0,v1,v2] += 9.000000 * T[o3,v3] * T[o0,o4,v0,v1] * T[o1,o2,v2,v4] * v[v3,v4,o3,o4];
R[o2,o1,o0,v0,v1,v2] += -3.000000 * T[o3,v3] * T[o0,o1,o2,v0,v1,v4] * v[v3,v4,v2,o3];
R[o2,o1,o0,v0,v1,v2] += -3.000000 * T[o3,v3] * T[o0,o1,o4,v0,v1,v2] * v[v3,o2,o3,o4];
R[o2,o1,o0,v0,v1,v2] += 1.000000 * T[o3,v3] * T[o0,o1,o2,o4,v0,v1,v2,v4] * v[v3,v4,o3,o4];
R[v2,o0,v0,v1] += 1.000000 * T[o1,v3] * T[o0,o2,v0,v1] * v[v2,v3,o1,o2];
R[o2,o1,o0,v0,v1,v2] += 3.000000 * T[o0,v3] * T[o3,v4] * T[o1,o2,o4,v0,v1,v2] * v[v3,v4,o3,o4];
R[o2,o1,o0,v0,v1,v2] += 9.000000 * T[o0,v3] * T[o1,v4] * T[o2,o3,v0,v1] * v[v3,v4,v2,o3];
R[o3,o2,o1,o0,v0,v1,v2,o4] += 6.000000 * T[o0,v3] * T[o1,v4] * T[o2,o3,o5,v0,v1,v2] * v[v3,v4,o4,o5];
R[o2,o1,o0,v0,v1,v2] += 1.500000 * T[o0,v3] * T[o1,v4] * T[o2,o3,o4,v0,v1,v2] * v[v3,v4,o3,o4];
R[o2,o1,o0,v0,v1,v2] += -9.000000 * T[o0,v3] * T[o1,o2,v0,v4] * v[v3,v4,v1,v2];
R[o2,o1,o0,v0,v1,v2] += -4.500000 * T[o0,v3] * T[o3,o4,v0,v1] * T[o1,o2,v2,v4] * v[v3,v4,o3,o4];
R[o3,o2,o1,o0,v0,v1,v2,o4] += -36.000000 * T[o0,v3] * T[o1,o5,v0,v1] * T[o2,o3,v2,v4] * v[v3,v4,o4,o5];
R[o2,o1,o0,v0,v1,v2] += 18.000000 * T[o0,v3] * T[o1,o3,v0,v1] * T[o2,o4,v2,v4] * v[v3,v4,o3,o4];
R[v2,o1,o0,v0,v1,o2] += 2.000000 * T[o0,v3] * T[o1,o3,v0,v1] * v[v2,v3,o2,o3];
R[o2,o1,o0,v0,v1,v2] += -18.000000 * T[o0,v3] * T[o1,o3,v0,v1] * v[v3,o2,v2,o3];
R[v2,o0,v0,v1] += -0.500000 * T[o0,v3] * T[o1,o2,v0,v1] * v[v2,v3,o1,o2];
R[o3,o2,o1,o0,v0,v1,v2,o4] += 12.000000 * T[o0,v3] * T[o1,o2,o3,v0,v1,v4] * v[v3,v4,v2,o4];
R[o2,o1,o0,v0,v1,v2] += 9.000000 * T[o0,v3] * T[o1,o2,o3,v0,v1,v4] * v[v3,v4,v2,o3];
R[o2,o1,o0,v0,v1,v2] += -3.000000 * T[o0,v3] * T[o1,o3,o4,v0,v1,v2] * v[v3,o2,o3,o4];
R[o4,o2,o1,o0,v0,v1,v2,o3] += 12.000000 * T[o0,v3] * T[o1,o2,o5,v0,v1,v2] * v[v3,o4,o3,o5];
R[o2,o1,o0,v0,v1,v2] += -3.000000 * T[o0,v3] * T[o1,o2,o3,v0,v1,v2] * v[v3,o3];
R[o3,o2,o1,o0,v0,v1,v2,o4] += -4.000000 * T[o0,v3] * T[o1,o2,o3,o5,v0,v1,v2,v4] * v[v3,v4,o4,o5];
R[o2,o1,o0,v0,v1,v2] += -1.500000 * T[o0,v3] * T[o1,o2,o3,o4,v0,v1,v2,v4] * v[v3,v4,o3,o4];
R[v2,o0,v0,v1] += -1.000000 * T[o0,v3] * v[v2,v3,v0,v1];
R[o2,o1,o0,v0,v1,v2] += 1.500000 * T[o0,o3,v3,v4] * T[o1,o2,o4,v0,v1,v2] * v[v3,v4,o3,o4];
R[o3,o2,o1,o0,v0,v1,v2,o4] += 3.000000 * T[o0,o1,v3,v4] * T[o2,o3,o5,v0,v1,v2] * v[v3,v4,o4,o5];
R[o2,o1,o0,v0,v1,v2] += 0.750000 * T[o0,o1,v3,v4] * T[o2,o3,o4,v0,v1,v2] * v[v3,v4,o3,o4];
R[o2,o1,o0,v0,v1,o3] += -1.000000 * T[o4,v2] * T[o0,o1,o2,v0,v1,v3] * v[v2,v3,o3,o4];
R[o1,o0,v0,v1] += -2.000000 * T[o2,v2] * T[o0,o1,v0,v3] * v[v2,v3,v1,o2];
R[o1,o0,v0,v1] += -2.000000 * T[o2,v2] * T[o0,o3,v0,v1] * v[v2,o1,o2,o3];
R[o1,o0,v0,v1] += 1.000000 * T[o2,v2] * T[o0,o1,o3,v0,v1,v3] * v[v2,v3,o2,o3];
R[o1,o0,v0,v1] += -2.000000 * T[o0,v2] * T[o2,v3] * T[o1,o3,v0,v1] * v[v2,v3,o2,o3];
R[o2,o1,o0,v0,v1,o3] += 3.000000 * T[o0,v2] * T[o1,v3] * T[o2,o4,v0,v1] * v[v2,v3,o3,o4];
R[o1,o0,v0,v1] += 0.500000 * T[o0,v2] * T[o1,v3] * T[o2,o3,v0,v1] * v[v2,v3,o2,o3];
R[o1,o0,v0,v1] += 1.000000 * T[o0,v2] * T[o1,v3] * v[v2,v3,v0,v1];
R[o2,o1,o0,v0,v1,o3] += -6.000000 * T[o0,v2] * T[o1,o2,v0,v3] * v[v2,v3,v1,o3];
R[o1,o0,v0,v1] += -4.000000 * T[o0,v2] * T[o1,o2,v0,v3] * v[v2,v3,v1,o2];
R[o1,o0,v0,v1] += 1.000000 * T[o0,v2] * T[o2,o3,v0,v1] * v[v2,o1,o2,o3];
R[o3,o1,o0,v0,v1,o2] += -6.000000 * T[o0,v2] * T[o1,o4,v0,v1] * v[v2,o3,o2,o4];
R[o1,o0,v0,v1] += 2.000000 * T[o0,v2] * T[o1,o2,v0,v1] * v[v2,o2];
R[o2,o1,o0,v0,v1,o3] += 3.000000 * T[o0,v2] * T[o1,o2,o4,v0,v1,v3] * v[v2,v3,o3,o4];
R[o3,o2,o1,o0,v0,v1,o4,o5] += 4.000000 * T[o0,v2] * T[o1,o2,o3,v0,v1,v3] * v[v2,v3,o4,o5];
R[o1,o0,v0,v1] += 1.000000 * T[o0,v2] * T[o1,o2,o3,v0,v1,v3] * v[v2,v3,o2,o3];
R[v1,o0,v0,o1] += -1.000000 * T[o0,v2] * v[v1,v2,v0,o1];
R[v1,v0] += 1.000000 * T[o0,v2] * v[v1,v2,v0,o0];
R[o1,o0,v0,v1] += 2.000000 * T[o0,v2] * v[v2,o1,v0,v1];
R[o1,o0,v0,v1] += 0.500000 * T[o0,o1,v2,v3] * v[v2,v3,v0,v1];
R[o1,o0,v0,o2] += -1.000000 * T[o3,v1] * T[o0,o1,v0,v2] * v[v1,v2,o2,o3];
R[o0,v0] += 1.000000 * T[o1,v1] * T[o0,o2,v0,v2] * v[v1,v2,o1,o2];
R[v0,o0] += 1.000000 * T[o1,v1] * v[v0,v1,o0,o1];
R[o0,v0] += -1.000000 * T[o1,v1] * v[v1,o0,v0,o1];
R[o1,o0,v0,o2] += 1.000000 * T[o0,v1] * T[o1,v2] * v[v1,v2,v0,o2];
R[o0,v0] += 1.000000 * T[o0,v1] * T[o1,v2] * v[v1,v2,v0,o1];
R[o1,o0,v0,o2] += -2.000000 * T[o0,v1] * T[o1,o3,v0,v2] * v[v1,v2,o2,o3];
R[o2,o1,o0,v0,o3,o4] += -3.000000 * T[o0,v1] * T[o1,o2,v0,v2] * v[v1,v2,o3,o4];
R[o0,v0] += -0.500000 * T[o0,v1] * T[o1,o2,v0,v2] * v[v1,v2,o1,o2];
R[v0,o0,o1,o2] += -1.000000 * T[o0,v1] * v[v0,v1,o1,o2];
R[o0,v0] += 1.000000 * T[o0,v1] * v[v1,v0];
R[o2,o0,v0,o1] += 2.000000 * T[o0,v1] * v[v1,o2,v0,o1];
R[o1,o0,v0,o2] += 0.500000 * T[o0,o1,v1,v2] * v[v1,v2,v0,o2];
R[o0,v0] += 0.500000 * T[o0,o1,v1,v2] * v[v1,v2,v0,o1];
R[o3,o2,o1,o0,v0,v1,v2,o4] += 12.000000 * T[o5,v0] * T[o0,v3] * T[o1,o2,o3,v1,v2,v4] * v[v3,v4,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,o4] += -18.000000 * T[o5,v0] * T[o0,o1,v1,v3] * T[o2,o3,v2,v4] * v[v3,v4,o4,o5];
R[o4,o2,o1,o0,v0,v1,v2,o3] += 12.000000 * T[o5,v0] * T[o0,o1,o2,v1,v2,v3] * v[v3,o4,o3,o5];
R[o3,o2,o1,o0,v0,v1,v2,o4] += 1.500000 * T[o5,v0] * T[o0,o1,o2,o3,v1,v2,v3,v4] * v[v3,v4,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -4.000000 * T[o4,v0] * T[o5,v4] * T[o0,o1,o2,o3,v1,v2,v3,v5] * v[v4,v5,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 24.000000 * T[o4,v0] * T[o0,v4] * T[o1,v5] * T[o2,o3,o5,v1,v2,v3] * v[v4,v5,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -144.000000 * T[o4,v0] * T[o0,v4] * T[o1,o5,v1,v2] * T[o2,o3,v3,v5] * v[v4,v5,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 48.000000 * T[o4,v0] * T[o0,v4] * T[o1,o2,o3,v1,v2,v5] * v[v4,v5,v3,o4];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 48.000000 * T[o4,v0] * T[o0,v4] * T[o1,o2,o5,v1,v2,v3] * v[v4,o3,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -16.000000 * T[o4,v0] * T[o0,v4] * T[o1,o2,o3,o5,v1,v2,v3,v5] * v[v4,v5,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 12.000000 * T[o4,v0] * T[o0,o1,v4,v5] * T[o2,o3,o5,v1,v2,v3] * v[v4,v5,o4,o5];
R[o2,o1,o0,v0,v1,o3] += 6.000000 * T[o4,v0] * T[o0,v2] * T[o1,o2,v1,v3] * v[v2,v3,o3,o4];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 24.000000 * T[o4,v0] * T[o5,v1] * T[o0,v4] * T[o1,o2,o3,v2,v3,v5] * v[v4,v5,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -36.000000 * T[o4,v0] * T[o5,v1] * T[o0,o1,v2,v4] * T[o2,o3,v3,v5] * v[v4,v5,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 24.000000 * T[o4,v0] * T[o5,v1] * T[o0,o1,o2,v2,v3,v4] * v[v4,o3,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 3.000000 * T[o4,v0] * T[o5,v1] * T[o0,o1,o2,o3,v2,v3,v4,v5] * v[v4,v5,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -48.000000 * T[o4,v0] * T[o0,o5,v1,v4] * T[o1,o2,o3,v2,v3,v5] * v[v4,v5,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -72.000000 * T[o4,v0] * T[o0,o1,v1,v4] * T[o2,o3,v2,v5] * v[v4,v5,v3,o4];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 72.000000 * T[o4,v0] * T[o0,o1,v1,v4] * T[o2,o3,o5,v2,v3,v5] * v[v4,v5,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -144.000000 * T[o4,v0] * T[o0,o5,v1,v2] * T[o1,o2,v3,v4] * v[v4,o3,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -24.000000 * T[o4,v0] * T[o0,o5,v1,v2] * T[o1,o2,o3,v3,v4,v5] * v[v4,v5,o4,o5];
R[o3,o1,o0,v0,v1,o2] += -6.000000 * T[o4,v0] * T[o0,o1,v1,v2] * v[v2,o3,o2,o4];
R[v3,o2,o1,o0,v0,v1,v2,o3] += -3.000000 * T[o4,v0] * T[o0,o1,o2,v1,v2,v4] * v[v3,v4,o3,o4];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 48.000000 * T[o4,v0] * T[o0,o1,o2,v1,v2,v4] * v[v4,o3,v3,o4];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 6.000000 * T[o4,v0] * T[o0,o1,o2,o3,v1,v2,v4,v5] * v[v4,v5,v3,o4];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 24.000000 * T[o4,v0] * T[o0,o1,o5,v1,v2,v3] * v[o2,o3,o4,o5];
R[o2,o1,o0,v0,v1,o3] += -1.000000 * T[o4,v0] * T[o0,o1,o2,v1,v2,v3] * v[v2,v3,o3,o4];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -16.000000 * T[o4,v0] * T[o0,o1,o2,o5,v1,v2,v3,v4] * v[v4,o3,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 4.000000 * T[o4,v0] * T[o0,o1,o2,o3,v1,v2,v3,v4] * v[v4,o4];
R[v4,o2,o1,o0,v0,v1,v2,v3] += -12.000000 * T[o3,v0] * T[o0,v5] * T[o1,o2,o4,v1,v2,v3] * v[v4,v5,o3,o4];
R[o2,o1,o0,v0,v1,v2] += 3.000000 * T[o3,v0] * T[o4,v3] * T[o0,o1,o2,v1,v2,v4] * v[v3,v4,o3,o4];
R[o2,o1,o0,v0,v1,v2] += -9.000000 * T[o3,v0] * T[o0,v3] * T[o1,v4] * T[o2,o4,v1,v2] * v[v3,v4,o3,o4];
R[o2,o1,o0,v0,v1,v2] += 18.000000 * T[o3,v0] * T[o0,v3] * T[o1,o2,v1,v4] * v[v3,v4,v2,o3];
R[o2,o1,o0,v0,v1,v2] += 18.000000 * T[o3,v0] * T[o0,v3] * T[o1,o4,v1,v2] * v[v3,o2,o3,o4];
R[o2,o1,o0,v0,v1,v2] += -9.000000 * T[o3,v0] * T[o0,v3] * T[o1,o2,o4,v1,v2,v4] * v[v3,v4,o3,o4];
R[o2,o1,o0,v0,v1,v2] += -9.000000 * T[o3,v0] * T[o4,v1] * T[o0,v3] * T[o1,o2,v2,v4] * v[v3,v4,o3,o4];
R[o2,o1,o0,v0,v1,v2] += 9.000000 * T[o3,v0] * T[o4,v1] * T[o0,o1,v2,v3] * v[v3,o2,o3,o4];
R[v4,o2,o1,o0,v0,v1,v2,v3] += -6.000000 * T[o3,v0] * T[o4,v1] * T[o0,o1,o2,v2,v3,v5] * v[v4,v5,o3,o4];
R[o2,o1,o0,v0,v1,v2] += 1.500000 * T[o3,v0] * T[o4,v1] * T[o0,o1,o2,v2,v3,v4] * v[v3,v4,o3,o4];
R[o1,o0,v0,o2] += 1.000000 * T[o3,v0] * T[o0,v1] * T[o1,v2] * v[v1,v2,o2,o3];
R[o2,o0,v0,o1] += 2.000000 * T[o3,v0] * T[o0,v1] * v[v1,o2,o1,o3];
R[o2,o1,o0,v0,v1,v2] += 18.000000 * T[o3,v0] * T[o0,o1,v1,v3] * T[o2,o4,v2,v4] * v[v3,v4,o3,o4];
R[v2,o1,o0,v0,v1,o2] += 2.000000 * T[o3,v0] * T[o0,o1,v1,v3] * v[v2,v3,o2,o3];
R[o2,o1,o0,v0,v1,v2] += -18.000000 * T[o3,v0] * T[o0,o1,v1,v3] * v[v3,o2,v2,o3];
R[o2,o1,o0,v0,v1,v2] += -3.000000 * T[o3,v0] * T[o0,o1,o2,v1,v3,v4] * v[v3,v4,v2,o3];
R[v4,o2,o1,o0,v0,v1,v2,v3] += 36.000000 * T[o3,v0] * T[o0,o4,v1,v2] * T[o1,o2,v3,v5] * v[v4,v5,o3,o4];
R[o2,o1,o0,v0,v1,v2] += -4.500000 * T[o3,v0] * T[o0,o4,v1,v2] * T[o1,o2,v3,v4] * v[v3,v4,o3,o4];
R[o2,o1,o0,v0,v1,v2] += -9.000000 * T[o3,v0] * T[o0,o4,v1,v2] * v[o1,o2,o3,o4];
R[o1,o0,v0,o2] += 0.500000 * T[o3,v0] * T[o0,o1,v1,v2] * v[v1,v2,o2,o3];
R[v4,o2,o1,o0,v0,v1,v2,v3] += -12.000000 * T[o3,v0] * T[o0,o1,o2,v1,v2,v5] * v[v4,v5,v3,o3];
R[v4,o2,o1,o0,v0,v1,v2,v3] += -12.000000 * T[o3,v0] * T[o0,o1,o4,v1,v2,v3] * v[v4,o2,o3,o4];
R[o2,o1,o0,v0,v1,v2] += 9.000000 * T[o3,v0] * T[o0,o1,o4,v1,v2,v3] * v[v3,o2,o3,o4];
R[o2,o1,o0,v0,v1,v2] += -3.000000 * T[o3,v0] * T[o0,o1,o2,v1,v2,v3] * v[v3,o3];
R[v4,o2,o1,o0,v0,v1,v2,v3] += 4.000000 * T[o3,v0] * T[o0,o1,o2,o4,v1,v2,v3,v5] * v[v4,v5,o3,o4];
R[o2,o1,o0,v0,v1,v2] += -1.500000 * T[o3,v0] * T[o0,o1,o2,o4,v1,v2,v3,v4] * v[v3,v4,o3,o4];
R[o2,o1,v0,o0] += 1.000000 * T[o3,v0] * v[o1,o2,o0,o3];
R[v3,o1,o0,v0,v1,v2] += -6.000000 * T[o2,v0] * T[o0,v4] * T[o1,o3,v1,v2] * v[v3,v4,o2,o3];
R[o1,o0,v0,v1] += -2.000000 * T[o2,v0] * T[o3,v2] * T[o0,o1,v1,v3] * v[v2,v3,o2,o3];
R[o1,o0,v0,v1] += 2.000000 * T[o2,v0] * T[o0,v2] * T[o1,v3] * v[v2,v3,v1,o2];
R[o1,o0,v0,v1] += -4.000000 * T[o2,v0] * T[o0,v2] * T[o1,o3,v1,v3] * v[v2,v3,o2,o3];
R[v1,o0,v0,o1] += -1.000000 * T[o2,v0] * T[o0,v2] * v[v1,v2,o1,o2];
R[o1,o0,v0,v1] += 4.000000 * T[o2,v0] * T[o0,v2] * v[v2,o1,v1,o2];
R[o1,o0,v0,v1] += 1.000000 * T[o2,v0] * T[o0,o1,v2,v3] * v[v2,v3,v1,o2];
R[o1,o0,v0,v1] += 1.000000 * T[o2,v0] * T[o3,v1] * T[o0,v2] * T[o1,v3] * v[v2,v3,o2,o3];
R[o1,o0,v0,v1] += 2.000000 * T[o2,v0] * T[o3,v1] * T[o0,v2] * v[v2,o1,o2,o3];
R[v3,o1,o0,v0,v1,v2] += -3.000000 * T[o2,v0] * T[o3,v1] * T[o0,o1,v2,v4] * v[v3,v4,o2,o3];
R[o1,o0,v0,v1] += 0.500000 * T[o2,v0] * T[o3,v1] * T[o0,o1,v2,v3] * v[v2,v3,o2,o3];
R[o1,o0,v0,v1] += 1.000000 * T[o2,v0] * T[o3,v1] * v[o0,o1,o2,o3];
R[v3,o1,o0,v0,v1,v2] += 6.000000 * T[o2,v0] * T[o0,o1,v1,v4] * v[v3,v4,v2,o2];
R[v3,o1,o0,v0,v1,v2] += 6.000000 * T[o2,v0] * T[o0,o3,v1,v2] * v[v3,o1,o2,o3];
R[o1,o0,v0,v1] += -4.000000 * T[o2,v0] * T[o0,o3,v1,v2] * v[v2,o1,o2,o3];
R[o1,o0,v0,v1] += 2.000000 * T[o2,v0] * T[o0,o1,v1,v2] * v[v2,o2];
R[v3,o1,o0,v0,v1,v2] += -3.000000 * T[o2,v0] * T[o0,o1,o3,v1,v2,v4] * v[v3,v4,o2,o3];
R[v5,v4,o1,o0,v0,v1,v2,v3] += 4.000000 * T[o2,v0] * T[o0,o1,o3,v1,v2,v3] * v[v4,v5,o2,o3];
R[o1,o0,v0,v1] += 1.000000 * T[o2,v0] * T[o0,o1,o3,v1,v2,v3] * v[v2,v3,o2,o3];
R[v1,o1,v0,o0] += -1.000000 * T[o2,v0] * v[v1,o1,o0,o2];
R[o1,o0] += -1.000000 * T[o2,v0] * v[v0,o1,o0,o2];
R[o1,o0,v0,v1] += 2.000000 * T[o2,v0] * v[o0,o1,v1,o2];
R[v2,o0,v0,v1] += -2.000000 * T[o1,v0] * T[o0,v3] * v[v2,v3,v1,o1];
R[v2,o0,v0,v1] += -1.000000 * T[o1,v0] * T[o2,v1] * T[o0,v3] * v[v2,v3,o1,o2];
R[v2,o0,v0,v1] += -1.000000 * T[o1,v0] * T[o2,v1] * v[v2,o0,o1,o2];
R[o0,v0] += 1.000000 * T[o1,v0] * T[o2,v1] * v[v1,o0,o1,o2];
R[o0,v0] += -1.000000 * T[o1,v0] * T[o0,v1] * T[o2,v2] * v[v1,v2,o1,o2];
R[o0,v0] += -1.000000 * T[o1,v0] * T[o0,v1] * v[v1,o1];
R[v2,o0,v0,v1] += 2.000000 * T[o1,v0] * T[o0,o2,v1,v3] * v[v2,v3,o1,o2];
R[v4,v3,o0,v0,v1,v2] += -3.000000 * T[o1,v0] * T[o0,o2,v1,v2] * v[v3,v4,o1,o2];
R[o0,v0] += -0.500000 * T[o1,v0] * T[o0,o2,v1,v2] * v[v1,v2,o1,o2];
R[o0,v0] += -1.000000 * T[o1,v0] * v[o0,o1];
R[v2,v1,v0,o0] += 1.000000 * T[o1,v0] * v[v1,v2,o0,o1];
R[v2,o0,v0,v1] += -2.000000 * T[o1,v0] * v[v2,o0,v1,o1];
R[v1,v0] += -1.000000 * T[o0,v0] * T[o1,v2] * v[v1,v2,o0,o1];
R[o0,o1] += 1.000000 * T[o0,v0] * T[o2,v1] * v[v0,v1,o1,o2];
R[o1,o0,o2,o3] += 1.000000 * T[o0,v0] * T[o1,v1] * v[v0,v1,o2,o3];
R[v3,v2,v0,v1] += 1.000000 * T[o0,v0] * T[o1,v1] * v[v2,v3,o0,o1];
R[o0,o1] += 1.000000 * T[o0,v0] * v[v0,o1];
R[o3,o0,o1,o2] += 2.000000 * T[o0,v0] * v[v0,o3,o1,o2];
R[v1,v0] += -1.000000 * T[o0,v0] * v[v1,o0];
R[v3,v2,v0,v1] += 2.000000 * T[o0,v0] * v[v2,v3,v1,o0];
R[v4,o2,o1,o0,v0,v1,v2,v3] += 12.000000 * T[o0,o3,v0,v5] * T[o1,o2,o4,v1,v2,v3] * v[v4,v5,o3,o4];
R[v4,o2,o1,o0,v0,v1,v2,v3] += 6.000000 * T[o0,o1,v0,v5] * T[o2,o3,o4,v1,v2,v3] * v[v4,v5,o3,o4];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -2.000000 * T[o4,o5,v0,v4] * T[o0,o1,o2,o3,v1,v2,v3,v5] * v[v4,v5,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -48.000000 * T[o0,o4,v0,v4] * T[o1,o2,o3,v1,v2,v5] * v[v4,v5,v3,o4];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -48.000000 * T[o0,o4,v0,v4] * T[o1,o2,o5,v1,v2,v3] * v[v4,o3,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 16.000000 * T[o0,o4,v0,v4] * T[o1,o2,o3,o5,v1,v2,v3,v5] * v[v4,v5,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -36.000000 * T[o0,o1,v0,v4] * T[o2,o3,v1,v5] * v[v4,v5,v2,v3];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 72.000000 * T[o0,o1,v0,v4] * T[o2,o3,o4,v1,v2,v5] * v[v4,v5,v3,o4];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -24.000000 * T[o0,o1,v0,v4] * T[o2,o4,o5,v1,v2,v3] * v[v4,o3,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -24.000000 * T[o0,o1,v0,v4] * T[o2,o3,o4,v1,v2,v3] * v[v4,o4];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -12.000000 * T[o0,o1,v0,v4] * T[o2,o3,o4,o5,v1,v2,v3,v5] * v[v4,v5,o4,o5];
R[v3,o1,o0,v0,v1,v2] += -3.000000 * T[o0,o1,v0,v4] * v[v3,v4,v1,v2];
R[o2,o1,o0,v0,v1,v2] += 1.500000 * T[o3,o4,v0,v3] * T[o0,o1,o2,v1,v2,v4] * v[v3,v4,o3,o4];
R[o3,o2,o1,o0,v0,v1,v2,o4] += -12.000000 * T[o0,o5,v0,v3] * T[o1,o2,o3,v1,v2,v4] * v[v3,v4,o4,o5];
R[o2,o1,o0,v0,v1,v2] += 9.000000 * T[o0,o3,v0,v3] * T[o1,o2,o4,v1,v2,v4] * v[v3,v4,o3,o4];
R[o3,o2,o1,o0,v0,v1,v2,o4] += -18.000000 * T[o0,o1,v0,v3] * T[o2,o3,v1,v4] * v[v3,v4,v2,o4];
R[o2,o1,o0,v0,v1,v2] += -18.000000 * T[o0,o1,v0,v3] * T[o2,o3,v1,v4] * v[v3,v4,v2,o3];
R[o3,o2,o1,o0,v0,v1,v2,o4] += 18.000000 * T[o0,o1,v0,v3] * T[o2,o3,o5,v1,v2,v4] * v[v3,v4,o4,o5];
R[o2,o1,o0,v0,v1,v2] += 4.500000 * T[o0,o1,v0,v3] * T[o2,o3,o4,v1,v2,v4] * v[v3,v4,o3,o4];
R[v2,o1,o0,v0,v1,o2] += -2.000000 * T[o0,o1,v0,v3] * v[v2,v3,v1,o2];
R[v2,o0,v0,v1] += 2.000000 * T[o0,o1,v0,v3] * v[v2,v3,v1,o1];
R[o2,o1,o0,v0,v1,v2] += 9.000000 * T[o0,o1,v0,v3] * v[v3,o2,v1,v2];
R[o2,o1,o0,v0,v1,v2] += 1.500000 * T[o0,o1,o2,v0,v3,v4] * v[v3,v4,v1,v2];
R[o1,o0,v0,v1] += 2.000000 * T[o0,o2,v0,v2] * T[o1,o3,v1,v3] * v[v2,v3,o2,o3];
R[v1,o0,v0,o1] += 1.000000 * T[o0,o2,v0,v2] * v[v1,v2,o1,o2];
R[o1,o0,v0,v1] += -4.000000 * T[o0,o2,v0,v2] * v[v2,o1,v1,o2];
R[o2,o1,o0,v0,v1,o3] += -6.000000 * T[o0,o1,v0,v2] * T[o2,o4,v1,v3] * v[v2,v3,o3,o4];
R[o3,o2,o1,o0,v0,v1,o4,o5] += -6.000000 * T[o0,o1,v0,v2] * T[o2,o3,v1,v3] * v[v2,v3,o4,o5];
R[o1,o0,v0,v1] += -1.000000 * T[o0,o1,v0,v2] * T[o2,o3,v1,v3] * v[v2,v3,o2,o3];
R[v1,o1,o0,v0,o2,o3] += -1.000000 * T[o0,o1,v0,v2] * v[v1,v2,o2,o3];
R[v1,v0] += -0.500000 * T[o0,o1,v0,v2] * v[v1,v2,o0,o1];
R[o1,o0,v0,v1] += 2.000000 * T[o0,o1,v0,v2] * v[v2,v1];
R[o3,o1,o0,v0,v1,o2] += 6.000000 * T[o0,o1,v0,v2] * v[v2,o3,v1,o2];
R[o2,o1,o0,v0,v1,o3] += 1.000000 * T[o0,o1,o2,v0,v2,v3] * v[v2,v3,v1,o3];
R[o1,o0,v0,v1] += 1.000000 * T[o0,o1,o2,v0,v2,v3] * v[v2,v3,v1,o2];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -18.000000 * T[o4,o5,v0,v1] * T[o0,o1,v2,v4] * T[o2,o3,v3,v5] * v[v4,v5,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 12.000000 * T[o4,o5,v0,v1] * T[o0,o1,o2,v2,v3,v4] * v[v4,o3,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 1.500000 * T[o4,o5,v0,v1] * T[o0,o1,o2,o3,v2,v3,v4,v5] * v[v4,v5,o4,o5];
R[o2,o1,o0,v0,v1,v2] += 4.500000 * T[o3,o4,v0,v1] * T[o0,o1,v2,v3] * v[v3,o2,o3,o4];
R[v4,o2,o1,o0,v0,v1,v2,v3] += -3.000000 * T[o3,o4,v0,v1] * T[o0,o1,o2,v2,v3,v5] * v[v4,v5,o3,o4];
R[o2,o1,o0,v0,v1,v2] += 0.750000 * T[o3,o4,v0,v1] * T[o0,o1,o2,v2,v3,v4] * v[v3,v4,o3,o4];
R[v3,o1,o0,v0,v1,v2] += -1.500000 * T[o2,o3,v0,v1] * T[o0,o1,v2,v4] * v[v3,v4,o2,o3];
R[o1,o0,v0,v1] += 0.250000 * T[o2,o3,v0,v1] * T[o0,o1,v2,v3] * v[v2,v3,o2,o3];
R[o1,o0,v0,v1] += 0.500000 * T[o2,o3,v0,v1] * v[o0,o1,o2,o3];
R[v2,o0,v0,v1] += -0.500000 * T[o1,o2,v0,v1] * v[v2,o0,o1,o2];
R[o0,v0] += 0.500000 * T[o1,o2,v0,v1] * v[v1,o0,o1,o2];
R[o4,o2,o1,o0,v0,v1,v2,o3] += -36.000000 * T[o0,o5,v0,v1] * T[o1,o2,v2,v3] * v[v3,o4,o3,o5];
R[o3,o2,o1,o0,v0,v1,v2,o4] += -6.000000 * T[o0,o5,v0,v1] * T[o1,o2,o3,v2,v3,v4] * v[v3,v4,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 144.000000 * T[o0,o4,v0,v1] * T[o1,o2,v2,v4] * T[o3,o5,v3,v5] * v[v4,v5,o4,o5];
R[v3,o2,o1,o0,v0,v1,v2,o3] += 9.000000 * T[o0,o4,v0,v1] * T[o1,o2,v2,v4] * v[v3,v4,o3,o4];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -144.000000 * T[o0,o4,v0,v1] * T[o1,o2,v2,v4] * v[v4,o3,v3,o4];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -24.000000 * T[o0,o4,v0,v1] * T[o1,o2,o3,v2,v4,v5] * v[v4,v5,v3,o4];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -18.000000 * T[o0,o4,v0,v1] * T[o1,o5,v2,v3] * T[o2,o3,v4,v5] * v[v4,v5,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -36.000000 * T[o0,o4,v0,v1] * T[o1,o5,v2,v3] * v[o2,o3,o4,o5];
R[o2,o1,o0,v0,v1,o3] += 1.500000 * T[o0,o4,v0,v1] * T[o1,o2,v2,v3] * v[v2,v3,o3,o4];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 72.000000 * T[o0,o4,v0,v1] * T[o1,o2,o5,v2,v3,v4] * v[v4,o3,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -24.000000 * T[o0,o4,v0,v1] * T[o1,o2,o3,v2,v3,v4] * v[v4,o4];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -12.000000 * T[o0,o4,v0,v1] * T[o1,o2,o3,o5,v2,v3,v4,v5] * v[v4,v5,o4,o5];
R[o3,o2,o0,v0,v1,o1] += 3.000000 * T[o0,o4,v0,v1] * v[o2,o3,o1,o4];
R[o2,o1,o0,v0,v1,v2] += 4.500000 * T[o0,o3,v0,v1] * T[o1,o2,v3,v4] * v[v3,v4,v2,o3];
R[v4,o2,o1,o0,v0,v1,v2,v3] += 36.000000 * T[o0,o3,v0,v1] * T[o1,o2,v2,v5] * v[v4,v5,v3,o3];
R[v4,o2,o1,o0,v0,v1,v2,v3] += 18.000000 * T[o0,o3,v0,v1] * T[o1,o4,v2,v3] * v[v4,o2,o3,o4];
R[o2,o1,o0,v0,v1,v2] += -18.000000 * T[o0,o3,v0,v1] * T[o1,o4,v2,v3] * v[v3,o2,o3,o4];
R[o2,o1,o0,v0,v1,v2] += 9.000000 * T[o0,o3,v0,v1] * T[o1,o2,v2,v3] * v[v3,o3];
R[v4,o2,o1,o0,v0,v1,v2,v3] += -18.000000 * T[o0,o3,v0,v1] * T[o1,o2,o4,v2,v3,v5] * v[v4,v5,o3,o4];
R[o2,o1,o0,v0,v1,v2] += 4.500000 * T[o0,o3,v0,v1] * T[o1,o2,o4,v2,v3,v4] * v[v3,v4,o3,o4];
R[v2,o2,o0,v0,v1,o1] += -2.000000 * T[o0,o3,v0,v1] * v[v2,o2,o1,o3];
R[o2,o0,v0,o1] += -2.000000 * T[o0,o3,v0,v1] * v[v1,o2,o1,o3];
R[o2,o1,o0,v0,v1,v2] += 9.000000 * T[o0,o3,v0,v1] * v[o1,o2,v2,o3];
R[v3,o1,o0,v0,v1,v2] += 6.000000 * T[o0,o2,v0,v1] * T[o1,o3,v2,v4] * v[v3,v4,o2,o3];
R[v5,v4,o1,o0,v0,v1,v2,v3] += -6.000000 * T[o0,o2,v0,v1] * T[o1,o3,v2,v3] * v[v4,v5,o2,o3];
R[o1,o0,v0,v1] += -1.000000 * T[o0,o2,v0,v1] * T[o1,o3,v2,v3] * v[v2,v3,o2,o3];
R[o1,o0,v0,v1] += -2.000000 * T[o0,o2,v0,v1] * v[o1,o2];
R[v3,v2,o0,v0,v1,o1] += 1.000000 * T[o0,o2,v0,v1] * v[v2,v3,o1,o2];
R[o0,o1] += 0.500000 * T[o0,o2,v0,v1] * v[v0,v1,o1,o2];
R[v3,o1,o0,v0,v1,v2] += -6.000000 * T[o0,o2,v0,v1] * v[v3,o1,v2,o2];
R[o1,o0,v0,o2] += 1.000000 * T[o0,o1,v0,v1] * v[v1,o2];
R[o4,o1,o0,v0,o2,o3] += 3.000000 * T[o0,o1,v0,v1] * v[v1,o4,o2,o3];
R[o1,o0,o2,o3] += 0.500000 * T[o0,o1,v0,v1] * v[v0,v1,o2,o3];
R[v2,o0,v0,v1] += -1.000000 * T[o0,o1,v0,v1] * v[v2,o1];
R[o0,v0] += 1.000000 * T[o0,o1,v0,v1] * v[v1,o1];
R[v3,v2,v0,v1] += 0.500000 * T[o0,o1,v0,v1] * v[v2,v3,o0,o1];
R[v4,v3,o0,v0,v1,v2] += 3.000000 * T[o0,o1,v0,v1] * v[v3,v4,v2,o1];
R[v4,o2,o1,o0,v0,v1,v2,v3] += -6.000000 * T[o0,o1,o2,v0,v1,v5] * v[v4,v5,v2,v3];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 18.000000 * T[o0,o1,o4,v0,v1,v4] * T[o2,o3,o5,v2,v3,v5] * v[v4,v5,o4,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 12.000000 * T[o0,o1,o2,v0,v1,v4] * T[o3,o4,o5,v2,v3,v5] * v[v4,v5,o4,o5];
R[v3,o2,o1,o0,v0,v1,v2,o3] += -3.000000 * T[o0,o1,o2,v0,v1,v4] * v[v3,v4,v2,o3];
R[v3,o1,o0,v0,v1,v2] += 3.000000 * T[o0,o1,o2,v0,v1,v4] * v[v3,v4,v2,o2];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 24.000000 * T[o0,o1,o2,v0,v1,v4] * v[v4,o3,v2,v3];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 3.000000 * T[o0,o1,o2,o3,v0,v1,v4,v5] * v[v4,v5,v2,v3];
R[v2,o1,o0,v0,v1,o2] += 1.000000 * T[o0,o1,o3,v0,v1,v3] * v[v2,v3,o2,o3];
R[o2,o1,o0,v0,v1,v2] += -9.000000 * T[o0,o1,o3,v0,v1,v3] * v[v3,o2,v2,o3];
R[v2,o2,o1,o0,v0,v1,o3,o4] += -1.000000 * T[o0,o1,o2,v0,v1,v3] * v[v2,v3,o3,o4];
R[v2,o0,v0,v1] += -0.500000 * T[o0,o1,o2,v0,v1,v3] * v[v2,v3,o1,o2];
R[o2,o1,o0,v0,v1,v2] += 3.000000 * T[o0,o1,o2,v0,v1,v3] * v[v3,v2];
R[o4,o2,o1,o0,v0,v1,v2,o3] += 12.000000 * T[o0,o1,o2,v0,v1,v3] * v[v3,o4,v2,o3];
R[o3,o2,o1,o0,v0,v1,v2,o4] += 1.500000 * T[o0,o1,o2,o3,v0,v1,v3,v4] * v[v3,v4,v2,o4];
R[o2,o1,o0,v0,v1,v2] += 1.500000 * T[o0,o1,o2,o3,v0,v1,v3,v4] * v[v3,v4,v2,o3];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 4.000000 * T[o0,o4,o5,v0,v1,v2] * T[o1,o2,o3,v3,v4,v5] * v[v4,v5,o4,o5];
R[o2,o1,o0,v0,v1,v2] += 1.500000 * T[o0,o3,o4,v0,v1,v2] * v[o1,o2,o3,o4];
R[v3,o1,o0,v0,v1,v2] += -1.000000 * T[o0,o2,o3,v0,v1,v2] * v[v3,o1,o2,o3];
R[o1,o0,v0,v1] += 1.000000 * T[o0,o2,o3,v0,v1,v2] * v[v2,o1,o2,o3];
R[o4,o3,o1,o0,v0,v1,v2,o2] += 6.000000 * T[o0,o1,o5,v0,v1,v2] * v[o3,o4,o2,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 12.000000 * T[o0,o1,o4,v0,v1,v2] * T[o2,o3,o5,v3,v4,v5] * v[v4,v5,o4,o5];
R[v3,o3,o1,o0,v0,v1,v2,o2] += -3.000000 * T[o0,o1,o4,v0,v1,v2] * v[v3,o3,o2,o4];
R[o3,o1,o0,v0,v1,o2] += -3.000000 * T[o0,o1,o4,v0,v1,v2] * v[v2,o3,o2,o4];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 24.000000 * T[o0,o1,o4,v0,v1,v2] * v[o2,o3,v3,o4];
R[o2,o1,o0,v0,v1,v2] += -3.000000 * T[o0,o1,o3,v0,v1,v2] * v[o2,o3];
R[v4,v3,o1,o0,v0,v1,v2,o2] += 1.000000 * T[o0,o1,o3,v0,v1,v2] * v[v3,v4,o2,o3];
R[o1,o0,v0,o2] += 0.500000 * T[o0,o1,o3,v0,v1,v2] * v[v1,v2,o2,o3];
R[v4,o2,o1,o0,v0,v1,v2,v3] += -12.000000 * T[o0,o1,o3,v0,v1,v2] * v[v4,o2,v3,o3];
R[o2,o1,o0,v0,v1,o3] += 1.000000 * T[o0,o1,o2,v0,v1,v2] * v[v2,o3];
R[o5,o2,o1,o0,v0,v1,o3,o4] += 4.000000 * T[o0,o1,o2,v0,v1,v2] * v[v2,o5,o3,o4];
R[o2,o1,o0,v0,o3,o4] += 0.500000 * T[o0,o1,o2,v0,v1,v2] * v[v1,v2,o3,o4];
R[v3,o1,o0,v0,v1,v2] += -1.000000 * T[o0,o1,o2,v0,v1,v2] * v[v3,o2];
R[o1,o0,v0,v1] += 1.000000 * T[o0,o1,o2,v0,v1,v2] * v[v2,o2];
R[v4,v3,o0,v0,v1,v2] += 0.500000 * T[o0,o1,o2,v0,v1,v2] * v[v3,v4,o1,o2];
R[o0,v0] += 0.250000 * T[o0,o1,o2,v0,v1,v2] * v[v1,v2,o1,o2];
R[v5,v4,o1,o0,v0,v1,v2,v3] += 4.000000 * T[o0,o1,o2,v0,v1,v2] * v[v4,v5,v3,o2];
R[v4,o2,o1,o0,v0,v1,v2,v3] += 4.000000 * T[o0,o1,o2,o3,v0,v1,v2,v5] * v[v4,v5,v3,o3];
R[v3,o2,o1,o0,v0,v1,v2,o3] += 1.000000 * T[o0,o1,o2,o4,v0,v1,v2,v4] * v[v3,v4,o3,o4];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -16.000000 * T[o0,o1,o2,o4,v0,v1,v2,v4] * v[v4,o3,v3,o4];
R[v3,o1,o0,v0,v1,v2] += -0.500000 * T[o0,o1,o2,o3,v0,v1,v2,v4] * v[v3,v4,o2,o3];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 4.000000 * T[o0,o1,o2,o3,v0,v1,v2,v4] * v[v4,v3];
R[o3,o2,o1,o0,v0,v1,v2,v3] += 3.000000 * T[o0,o1,o4,o5,v0,v1,v2,v3] * v[o2,o3,o4,o5];
R[v4,o2,o1,o0,v0,v1,v2,v3] += -1.500000 * T[o0,o1,o3,o4,v0,v1,v2,v3] * v[v4,o2,o3,o4];
R[o2,o1,o0,v0,v1,v2] += 1.500000 * T[o0,o1,o3,o4,v0,v1,v2,v3] * v[v3,o2,o3,o4];
R[o4,o2,o1,o0,v0,v1,v2,o3] += -4.000000 * T[o0,o1,o2,o5,v0,v1,v2,v3] * v[v3,o4,o3,o5];
R[o3,o2,o1,o0,v0,v1,v2,v3] += -4.000000 * T[o0,o1,o2,o4,v0,v1,v2,v3] * v[o3,o4];
R[o2,o1,o0,v0,v1,o3] += 0.500000 * T[o0,o1,o2,o4,v0,v1,v2,v3] * v[v2,v3,o3,o4];
R[o3,o2,o1,o0,v0,v1,v2,o4] += 1.000000 * T[o0,o1,o2,o3,v0,v1,v2,v3] * v[v3,o4];
R[o3,o2,o1,o0,v0,v1,o4,o5] += 0.500000 * T[o0,o1,o2,o3,v0,v1,v2,v3] * v[v2,v3,o4,o5];
R[v4,o2,o1,o0,v0,v1,v2,v3] += -1.000000 * T[o0,o1,o2,o3,v0,v1,v2,v3] * v[v4,o3];
R[o2,o1,o0,v0,v1,v2] += 1.000000 * T[o0,o1,o2,o3,v0,v1,v2,v3] * v[v3,o3];
R[v5,v4,o1,o0,v0,v1,v2,v3] += 0.500000 * T[o0,o1,o2,o3,v0,v1,v2,v3] * v[v4,v5,o2,o3];
R[o1,o0,v0,v1] += 0.250000 * T[o0,o1,o2,o3,v0,v1,v2,v3] * v[v2,v3,o2,o3];
R[o1,o0] += 1.000000 * v[o1,o0];
R[v0,o0] += 1.000000 * v[v0,o0];
R[o3,o2,o0,o1] += 1.000000 * v[o2,o3,o0,o1];
R[v0,o2,o0,o1] += -1.000000 * v[v0,o2,o0,o1];
R[v1,v0,o0,o1] += 1.000000 * v[v0,v1,o0,o1];
R[o0,v0] += 1.000000 * v[o0,v0];
R[v1,v0] += 1.000000 * v[v1,v0];
R[o2,o1,v0,o0] += 1.000000 * v[o1,o2,v0,o0];
R[v1,o1,v0,o0] += -1.000000 * v[v1,o1,v0,o0];
R[v2,v1,v0,o0] += 1.000000 * v[v1,v2,v0,o0];
R[o1,o0,v0,v1] += 1.000000 * v[o0,o1,v0,v1];
R[v2,o0,v0,v1] += -1.000000 * v[v2,o0,v0,v1];
R[v3,v2,v0,v1] += 1.000000 * v[v2,v3,v0,v1];

















In [8]:



    


wdt.set_print(pw.WDiagPrint.summary)
Hthis = pw.make_operator("H", ["ov->oo"]);
Hbar2 = pw.bch_series(Hthis,T,4)
sum = wdt.contract_sum(pw.rational(1,1),Hbar2, 2, 2)

print sum.str()

eqs = sum.to_manybody_equation("R")



    


















    






-1/2 H^{v1,o0}_{o1,o2} T^{o2}_{v0} T^{o1}_{v1} { a+(v0) a-(o0) }
1/2 H^{v1,o0}_{o1,o2} T^{o1}_{v0} T^{o2}_{v1} { a+(v0) a-(o0) }
1/2 H^{v1,o0}_{o1,o2} T^{o1,o2}_{v0,v1} { a+(v0) a-(o0) }
-1 H^{v0,o1}_{o0,o2} T^{o2}_{v0} { a+(o0) a-(o1) }

















In [ ]:

Content source: fevangelista/pyWicked

Similar notebooks:

notebook.community | gallery | about