Hammer  1.0.0
Helicity Amplitude Module for Matrix Element Reweighting
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
RateLbLcLepNu.cc
Go to the documentation of this file.
1 ///
2 /// @file RateLbLcLepNu.cc
3 /// @brief \f$ \Lambda_b \rightarrow \Lambda_c \tau\nu \f$ total rate
4 ///
5 
6 //**** This file is a part of the HAMMER library
7 //**** Copyright (C) 2016 - 2020 The HAMMER Collaboration
8 //**** HAMMER is licensed under version 3 of the GPL; see COPYING for details
9 //**** Please note the MCnet academic guidelines; see GUIDELINES for details
10 
11 // -*- C++ -*-
12 #include "Hammer/Rates/RateLbLcLepNu.hh"
13 #include "Hammer/IndexLabels.hh"
14 #include "Hammer/Math/Constants.hh"
15 #include "Hammer/Particle.hh"
16 #include "Hammer/Tools/Pdg.hh"
17 #include "Hammer/Math/MultiDim/SparseContainer.hh"
18 #include <cmath>
19 
20 using namespace std;
21 
22 namespace Hammer {
23 
24  namespace MD = MultiDimensional;
25 
26  RateLbLcLepNu::RateLbLcLepNu() {
27  // Create tensor rank and dimensions
28  vector<IndexType> dims{{11, 12, 11, 12, _nPoints}};
29  string name{"RateLbLcLepNuQ2"};
30  auto& pdg = PID::instance();
31 
32  addProcessSignature(-PID::LAMBDAB, {PID::LAMBDACMINUS, PID::NU_TAU, PID::ANTITAU});
33  addIntegrationBoundaries({PS::makeQ2Function(pdg.getMass(PID::ANTITAU), pdg.getMass(PID::LAMBDAB)-pdg.getMass(PID::LAMBDACMINUS))});
34  addTensor(Tensor{name, MD::makeEmptySparse(dims, {WILSON_BCTAUNU, FF_LBLC, WILSON_BCTAUNU_HC, FF_LBLC_HC, INTEGRATION_INDEX})});
35 
36  addProcessSignature(-PID::LAMBDAB, {PID::LAMBDACMINUS, PID::NU_MU, PID::ANTIMUON});
37  addIntegrationBoundaries({PS::makeQ2Function(pdg.getMass(PID::ANTIMUON), pdg.getMass(PID::LAMBDAB)-pdg.getMass(PID::LAMBDACMINUS))});
38  addTensor(Tensor{name, MD::makeEmptySparse(dims, {WILSON_BCTAUNU, FF_LBLC, WILSON_BCTAUNU_HC, FF_LBLC_HC, INTEGRATION_INDEX})});
39 
40  addProcessSignature(-PID::LAMBDAB, {PID::LAMBDACMINUS, PID::NU_E, PID::POSITRON});
41  addIntegrationBoundaries({PS::makeQ2Function(pdg.getMass(PID::POSITRON), pdg.getMass(PID::LAMBDAB)-pdg.getMass(PID::LAMBDACMINUS))});
42  addTensor(Tensor{name, MD::makeEmptySparse(dims, {WILSON_BCTAUNU, FF_LBLC, WILSON_BCTAUNU_HC, FF_LBLC_HC, INTEGRATION_INDEX})});
43 
44  setSignatureIndex();
45  }
46 
47  Tensor RateLbLcLepNu::evalAtPSPoint(const vector<double>& point) {
48  auto labs = getTensor().labels();
49  labs.pop_back();
50  auto dimensions = getTensor().dims();
51  dimensions.pop_back();
52  Tensor result{"RateLbLcLepNu", MD::makeEmptySparse(dimensions, labs)};
53 
54  const double Mb = masses()[0];
55  const double Mc = masses()[1];
56  const double Mt = masses()[3];
57 
58  // kinematic objects
59  const double Sqq = point[0];
60  const double mSqq = Sqq/(Mb*Mb);
61  const double w = (Mb*Mb + Mc*Mc - Sqq)/(2.*Mb*Mc);
62  const double rC = Mc/Mb;
63  const double rt = Mt/Mb;
64 
65  const double mSqq2 = pow(mSqq,2.);
66  const double wSq = w*w;
67  const double rCSq = rC*rC;
68  const double rtSq = pow(rt,2.);
69 
70  const double RateNorm = (pow(GFermi,2.)*pow(Mb,5.)*pow(rC,3.)*pow(mSqq - rtSq,2.)*sqrt(-1 + wSq))/(16.*mSqq*pi3) ;
71 
72  // set non-zero tensor elements
73  result.element({0,2,0,2}) = (2.*(2.*mSqq2*(-1 + w) + mSqq*(1 + 2.*rC + rCSq + rtSq)*(-1 + w) + rtSq*(1 + 2.*w - 2.*rC*(2 + w) + rCSq*(1 + 2.*w))))/(3.*mSqq2) ;
74  result.element({0,2,0,3}) = ((1 + w)*(2.*mSqq*rC*(1 + rC)*(-1 + w) + rtSq*(3 - 2.*rC*(2 + w) + rCSq*(-1 + 4.*w))))/(3.*mSqq2) ;
75  result.element({0,2,0,4}) = ((1 + w)*(2.*mSqq*(1 + rC)*(-1 + w) + rtSq*(-1 + 3.*rCSq + 4.*w - 2.*rC*(2 + w))))/(3.*mSqq2) ;
76  result.element({0,2,1,0}) = ((-1 + rC)*rt*(1 + w))/mSqq ;
77  result.element({0,2,3,0}) = ((-1 + rC)*rt*(1 + w))/mSqq ;
78  result.element({0,2,5,2}) = (2.*(2.*mSqq2*(-1 + w) + mSqq*(1 + 2.*rC + rCSq + rtSq)*(-1 + w) + rtSq*(1 + 2.*w - 2.*rC*(2 + w) + rCSq*(1 + 2.*w))))/(3.*mSqq2) ;
79  result.element({0,2,5,3}) = ((1 + w)*(2.*mSqq*rC*(1 + rC)*(-1 + w) + rtSq*(3 - 2.*rC*(2 + w) + rCSq*(-1 + 4.*w))))/(3.*mSqq2) ;
80  result.element({0,2,5,4}) = ((1 + w)*(2.*mSqq*(1 + rC)*(-1 + w) + rtSq*(-1 + 3.*rCSq + 4.*w - 2.*rC*(2 + w))))/(3.*mSqq2) ;
81  result.element({0,2,7,2}) = (2.*(2.*mSqq2*(-1 + w) + mSqq*(1 + 2.*rC + rCSq + rtSq)*(-1 + w) + rtSq*(1 + 2.*w - 2.*rC*(2 + w) + rCSq*(1 + 2.*w))))/(3.*mSqq2) ;
82  result.element({0,2,7,3}) = ((1 + w)*(2.*mSqq*rC*(1 + rC)*(-1 + w) + rtSq*(3 - 2.*rC*(2 + w) + rCSq*(-1 + 4.*w))))/(3.*mSqq2) ;
83  result.element({0,2,7,4}) = ((1 + w)*(2.*mSqq*(1 + rC)*(-1 + w) + rtSq*(-1 + 3.*rCSq + 4.*w - 2.*rC*(2 + w))))/(3.*mSqq2) ;
84  result.element({0,2,9,8}) = (-12.*(1 + rC)*rt*(-1 + w))/mSqq ;
85  result.element({0,2,9,9}) = (-4.*rt*(-1 + w)*(-3 + rC*(-1 + 2.*w)))/mSqq ;
86  result.element({0,2,9,10}) = (-4.*rt*(1 + 3.*rC - 2.*w)*(-1 + w))/mSqq ;
87  result.element({0,2,9,11}) = (4.*(1 + rC)*rt*(-1 + wSq))/mSqq ;
88  result.element({0,3,0,2}) = ((1 + w)*(2.*mSqq*rC*(1 + rC)*(-1 + w) + rtSq*(3 - 2.*rC*(2 + w) + rCSq*(-1 + 4.*w))))/(3.*mSqq2) ;
89  result.element({0,3,0,3}) = ((1 + w)*(2.*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6.*rC*w + rCSq*(-1 + 4.*wSq))))/(3.*mSqq2) ;
90  result.element({0,3,0,4}) = ((1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
91  result.element({0,3,1,0}) = (rt*(1 + w)*(-1 + rC*w))/mSqq ;
92  result.element({0,3,3,0}) = (rt*(1 + w)*(-1 + rC*w))/mSqq ;
93  result.element({0,3,5,2}) = ((1 + w)*(2.*mSqq*rC*(1 + rC)*(-1 + w) + rtSq*(3 - 2.*rC*(2 + w) + rCSq*(-1 + 4.*w))))/(3.*mSqq2) ;
94  result.element({0,3,5,3}) = ((1 + w)*(2.*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6.*rC*w + rCSq*(-1 + 4.*wSq))))/(3.*mSqq2) ;
95  result.element({0,3,5,4}) = ((1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
96  result.element({0,3,7,2}) = ((1 + w)*(2.*mSqq*rC*(1 + rC)*(-1 + w) + rtSq*(3 - 2.*rC*(2 + w) + rCSq*(-1 + 4.*w))))/(3.*mSqq2) ;
97  result.element({0,3,7,3}) = ((1 + w)*(2.*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6.*rC*w + rCSq*(-1 + 4.*wSq))))/(3.*mSqq2) ;
98  result.element({0,3,7,4}) = ((1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
99  result.element({0,3,9,8}) = (-4.*rC*rt*(-1 + wSq))/mSqq ;
100  result.element({0,3,9,9}) = (4.*rC*rt*(-1 + wSq))/mSqq ;
101  result.element({0,3,9,10}) = (-4.*rC*rt*(-1 + wSq))/mSqq ;
102  result.element({0,3,9,11}) = (4.*rC*rt*(-1 + w)*pow(1 + w,2.))/mSqq ;
103  result.element({0,4,0,2}) = ((1 + w)*(2.*mSqq*(1 + rC)*(-1 + w) + rtSq*(-1 + 3.*rCSq + 4.*w - 2.*rC*(2 + w))))/(3.*mSqq2) ;
104  result.element({0,4,0,3}) = ((1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
105  result.element({0,4,0,4}) = ((1 + w)*(2.*mSqq*(-1 + wSq) + rtSq*(-1 + 3.*rCSq - 6.*rC*w + 4.*wSq)))/(3.*mSqq2) ;
106  result.element({0,4,1,0}) = (rt*(rC - w)*(1 + w))/mSqq ;
107  result.element({0,4,3,0}) = (rt*(rC - w)*(1 + w))/mSqq ;
108  result.element({0,4,5,2}) = ((1 + w)*(2.*mSqq*(1 + rC)*(-1 + w) + rtSq*(-1 + 3.*rCSq + 4.*w - 2.*rC*(2 + w))))/(3.*mSqq2) ;
109  result.element({0,4,5,3}) = ((1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
110  result.element({0,4,5,4}) = ((1 + w)*(2.*mSqq*(-1 + wSq) + rtSq*(-1 + 3.*rCSq - 6.*rC*w + 4.*wSq)))/(3.*mSqq2) ;
111  result.element({0,4,7,2}) = ((1 + w)*(2.*mSqq*(1 + rC)*(-1 + w) + rtSq*(-1 + 3.*rCSq + 4.*w - 2.*rC*(2 + w))))/(3.*mSqq2) ;
112  result.element({0,4,7,3}) = ((1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
113  result.element({0,4,7,4}) = ((1 + w)*(2.*mSqq*(-1 + wSq) + rtSq*(-1 + 3.*rCSq - 6.*rC*w + 4.*wSq)))/(3.*mSqq2) ;
114  result.element({0,4,9,8}) = (-4.*rt*(-1 + wSq))/mSqq ;
115  result.element({0,4,9,9}) = (4.*rt*(-1 + wSq))/mSqq ;
116  result.element({0,4,9,10}) = (-4.*rt*(-1 + wSq))/mSqq ;
117  result.element({0,4,9,11}) = (4.*rt*(-1 + w)*pow(1 + w,2.))/mSqq ;
118  result.element({0,5,0,5}) = (2.*(2.*mSqq2*(1 + w) + mSqq*(1 - 2.*rC + rCSq + rtSq)*(1 + w) + rtSq*(-1 + 2.*rC*(-2 + w) + 2.*w + rCSq*(-1 + 2.*w))))/(3.*mSqq2) ;
119  result.element({0,5,0,6}) = ((-1 + w)*(2.*mSqq*(-1 + rC)*rC*(1 + w) + rtSq*(-3 + 2.*rC*(-2 + w) + rCSq*(1 + 4.*w))))/(3.*mSqq2) ;
120  result.element({0,5,0,7}) = ((-1 + w)*(rtSq*(-1 + 3.*rCSq - 2.*rC*(-2 + w) - 4.*w) + 2.*mSqq*(-1 + rC)*(1 + w)))/(3.*mSqq2) ;
121  result.element({0,5,1,1}) = ((1 + rC)*rt*(-1 + w))/mSqq ;
122  result.element({0,5,3,1}) = -(((1 + rC)*rt*(-1 + w))/mSqq) ;
123  result.element({0,5,5,5}) = (-2.*(2.*mSqq2*(1 + w) + mSqq*(1 - 2.*rC + rCSq + rtSq)*(1 + w) + rtSq*(-1 + 2.*rC*(-2 + w) + 2.*w + rCSq*(-1 + 2.*w))))/(3.*mSqq2) ;
124  result.element({0,5,5,6}) = -((-1 + w)*(2.*mSqq*(-1 + rC)*rC*(1 + w) + rtSq*(-3 + 2.*rC*(-2 + w) + rCSq*(1 + 4.*w))))/(3.*mSqq2) ;
125  result.element({0,5,5,7}) = -((-1 + w)*(rtSq*(-1 + 3.*rCSq - 2.*rC*(-2 + w) - 4.*w) + 2.*mSqq*(-1 + rC)*(1 + w)))/(3.*mSqq2) ;
126  result.element({0,5,7,5}) = (2.*(2.*mSqq2*(1 + w) + mSqq*(1 - 2.*rC + rCSq + rtSq)*(1 + w) + rtSq*(-1 + 2.*rC*(-2 + w) + 2.*w + rCSq*(-1 + 2.*w))))/(3.*mSqq2) ;
127  result.element({0,5,7,6}) = ((-1 + w)*(2.*mSqq*(-1 + rC)*rC*(1 + w) + rtSq*(-3 + 2.*rC*(-2 + w) + rCSq*(1 + 4.*w))))/(3.*mSqq2) ;
128  result.element({0,5,7,7}) = ((-1 + w)*(rtSq*(-1 + 3.*rCSq - 2.*rC*(-2 + w) - 4.*w) + 2.*mSqq*(-1 + rC)*(1 + w)))/(3.*mSqq2) ;
129  result.element({0,5,9,8}) = (-12.*(-1 + rC)*rt*(1 + w))/mSqq ;
130  result.element({0,5,9,9}) = (-8.*rC*rt*(-1 + wSq))/mSqq ;
131  result.element({0,5,9,10}) = (-8.*rt*(-1 + wSq))/mSqq ;
132  result.element({0,6,0,5}) = ((-1 + w)*(2.*mSqq*(-1 + rC)*rC*(1 + w) + rtSq*(-3 + 2.*rC*(-2 + w) + rCSq*(1 + 4.*w))))/(3.*mSqq2) ;
133  result.element({0,6,0,6}) = ((-1 + w)*(2.*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6.*rC*w + rCSq*(-1 + 4.*wSq))))/(3.*mSqq2) ;
134  result.element({0,6,0,7}) = ((-1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
135  result.element({0,6,1,1}) = (rt*(-1 + w)*(-1 + rC*w))/mSqq ;
136  result.element({0,6,3,1}) = -((rt*(-1 + w)*(-1 + rC*w))/mSqq) ;
137  result.element({0,6,5,5}) = -((-1 + w)*(2.*mSqq*(-1 + rC)*rC*(1 + w) + rtSq*(-3 + 2.*rC*(-2 + w) + rCSq*(1 + 4.*w))))/(3.*mSqq2) ;
138  result.element({0,6,5,6}) = -((-1 + w)*(2.*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6.*rC*w + rCSq*(-1 + 4.*wSq))))/(3.*mSqq2) ;
139  result.element({0,6,5,7}) = -((-1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
140  result.element({0,6,7,5}) = ((-1 + w)*(2.*mSqq*(-1 + rC)*rC*(1 + w) + rtSq*(-3 + 2.*rC*(-2 + w) + rCSq*(1 + 4.*w))))/(3.*mSqq2) ;
141  result.element({0,6,7,6}) = ((-1 + w)*(2.*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6.*rC*w + rCSq*(-1 + 4.*wSq))))/(3.*mSqq2) ;
142  result.element({0,6,7,7}) = ((-1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
143  result.element({0,6,9,8}) = (-4.*rC*rt*(-1 + wSq))/mSqq ;
144  result.element({0,7,0,5}) = ((-1 + w)*(rtSq*(-1 + 3.*rCSq - 2.*rC*(-2 + w) - 4.*w) + 2.*mSqq*(-1 + rC)*(1 + w)))/(3.*mSqq2) ;
145  result.element({0,7,0,6}) = ((-1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
146  result.element({0,7,0,7}) = ((-1 + w)*(2.*mSqq*(-1 + wSq) + rtSq*(-1 + 3.*rCSq - 6.*rC*w + 4.*wSq)))/(3.*mSqq2) ;
147  result.element({0,7,1,1}) = (rt*(rC - w)*(-1 + w))/mSqq ;
148  result.element({0,7,3,1}) = -((rt*(rC - w)*(-1 + w))/mSqq) ;
149  result.element({0,7,5,5}) = -((-1 + w)*(rtSq*(-1 + 3.*rCSq - 2.*rC*(-2 + w) - 4.*w) + 2.*mSqq*(-1 + rC)*(1 + w)))/(3.*mSqq2) ;
150  result.element({0,7,5,6}) = -((-1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
151  result.element({0,7,5,7}) = -((-1 + w)*(2.*mSqq*(-1 + wSq) + rtSq*(-1 + 3.*rCSq - 6.*rC*w + 4.*wSq)))/(3.*mSqq2) ;
152  result.element({0,7,7,5}) = ((-1 + w)*(rtSq*(-1 + 3.*rCSq - 2.*rC*(-2 + w) - 4.*w) + 2.*mSqq*(-1 + rC)*(1 + w)))/(3.*mSqq2) ;
153  result.element({0,7,7,6}) = ((-1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
154  result.element({0,7,7,7}) = ((-1 + w)*(2.*mSqq*(-1 + wSq) + rtSq*(-1 + 3.*rCSq - 6.*rC*w + 4.*wSq)))/(3.*mSqq2) ;
155  result.element({0,7,9,8}) = (-4.*rt*(-1 + wSq))/mSqq ;
156  result.element({1,0,0,2}) = ((-1 + rC)*rt*(1 + w))/mSqq ;
157  result.element({1,0,0,3}) = (rt*(1 + w)*(-1 + rC*w))/mSqq ;
158  result.element({1,0,0,4}) = (rt*(rC - w)*(1 + w))/mSqq ;
159  result.element({1,0,1,0}) = 1 + w ;
160  result.element({1,0,3,0}) = 1 + w ;
161  result.element({1,0,5,2}) = ((-1 + rC)*rt*(1 + w))/mSqq ;
162  result.element({1,0,5,3}) = (rt*(1 + w)*(-1 + rC*w))/mSqq ;
163  result.element({1,0,5,4}) = (rt*(rC - w)*(1 + w))/mSqq ;
164  result.element({1,0,7,2}) = ((-1 + rC)*rt*(1 + w))/mSqq ;
165  result.element({1,0,7,3}) = (rt*(1 + w)*(-1 + rC*w))/mSqq ;
166  result.element({1,0,7,4}) = (rt*(rC - w)*(1 + w))/mSqq ;
167  result.element({1,1,0,5}) = ((1 + rC)*rt*(-1 + w))/mSqq ;
168  result.element({1,1,0,6}) = (rt*(-1 + w)*(-1 + rC*w))/mSqq ;
169  result.element({1,1,0,7}) = (rt*(rC - w)*(-1 + w))/mSqq ;
170  result.element({1,1,1,1}) = -1 + w ;
171  result.element({1,1,3,1}) = 1 - w ;
172  result.element({1,1,5,5}) = -(((1 + rC)*rt*(-1 + w))/mSqq) ;
173  result.element({1,1,5,6}) = -((rt*(-1 + w)*(-1 + rC*w))/mSqq) ;
174  result.element({1,1,5,7}) = -((rt*(rC - w)*(-1 + w))/mSqq) ;
175  result.element({1,1,7,5}) = ((1 + rC)*rt*(-1 + w))/mSqq ;
176  result.element({1,1,7,6}) = (rt*(-1 + w)*(-1 + rC*w))/mSqq ;
177  result.element({1,1,7,7}) = (rt*(rC - w)*(-1 + w))/mSqq ;
178  result.element({2,0,2,0}) = 1 + w ;
179  result.element({2,0,4,0}) = 1 + w ;
180  result.element({2,0,6,2}) = ((-1 + rC)*rt*(1 + w))/mSqq ;
181  result.element({2,0,6,3}) = (rt*(1 + w)*(-1 + rC*w))/mSqq ;
182  result.element({2,0,6,4}) = (rt*(rC - w)*(1 + w))/mSqq ;
183  result.element({2,0,8,2}) = ((-1 + rC)*rt*(1 + w))/mSqq ;
184  result.element({2,0,8,3}) = (rt*(1 + w)*(-1 + rC*w))/mSqq ;
185  result.element({2,0,8,4}) = (rt*(rC - w)*(1 + w))/mSqq ;
186  result.element({2,1,2,1}) = -1 + w ;
187  result.element({2,1,4,1}) = 1 - w ;
188  result.element({2,1,6,5}) = -(((1 + rC)*rt*(-1 + w))/mSqq) ;
189  result.element({2,1,6,6}) = -((rt*(-1 + w)*(-1 + rC*w))/mSqq) ;
190  result.element({2,1,6,7}) = -((rt*(rC - w)*(-1 + w))/mSqq) ;
191  result.element({2,1,8,5}) = ((1 + rC)*rt*(-1 + w))/mSqq ;
192  result.element({2,1,8,6}) = (rt*(-1 + w)*(-1 + rC*w))/mSqq ;
193  result.element({2,1,8,7}) = (rt*(rC - w)*(-1 + w))/mSqq ;
194  result.element({3,0,0,2}) = ((-1 + rC)*rt*(1 + w))/mSqq ;
195  result.element({3,0,0,3}) = (rt*(1 + w)*(-1 + rC*w))/mSqq ;
196  result.element({3,0,0,4}) = (rt*(rC - w)*(1 + w))/mSqq ;
197  result.element({3,0,1,0}) = 1 + w ;
198  result.element({3,0,3,0}) = 1 + w ;
199  result.element({3,0,5,2}) = ((-1 + rC)*rt*(1 + w))/mSqq ;
200  result.element({3,0,5,3}) = (rt*(1 + w)*(-1 + rC*w))/mSqq ;
201  result.element({3,0,5,4}) = (rt*(rC - w)*(1 + w))/mSqq ;
202  result.element({3,0,7,2}) = ((-1 + rC)*rt*(1 + w))/mSqq ;
203  result.element({3,0,7,3}) = (rt*(1 + w)*(-1 + rC*w))/mSqq ;
204  result.element({3,0,7,4}) = (rt*(rC - w)*(1 + w))/mSqq ;
205  result.element({3,1,0,5}) = -(((1 + rC)*rt*(-1 + w))/mSqq) ;
206  result.element({3,1,0,6}) = -((rt*(-1 + w)*(-1 + rC*w))/mSqq) ;
207  result.element({3,1,0,7}) = -((rt*(rC - w)*(-1 + w))/mSqq) ;
208  result.element({3,1,1,1}) = 1 - w ;
209  result.element({3,1,3,1}) = -1 + w ;
210  result.element({3,1,5,5}) = ((1 + rC)*rt*(-1 + w))/mSqq ;
211  result.element({3,1,5,6}) = (rt*(-1 + w)*(-1 + rC*w))/mSqq ;
212  result.element({3,1,5,7}) = (rt*(rC - w)*(-1 + w))/mSqq ;
213  result.element({3,1,7,5}) = -(((1 + rC)*rt*(-1 + w))/mSqq) ;
214  result.element({3,1,7,6}) = -((rt*(-1 + w)*(-1 + rC*w))/mSqq) ;
215  result.element({3,1,7,7}) = -((rt*(rC - w)*(-1 + w))/mSqq) ;
216  result.element({4,0,2,0}) = 1 + w ;
217  result.element({4,0,4,0}) = 1 + w ;
218  result.element({4,0,6,2}) = ((-1 + rC)*rt*(1 + w))/mSqq ;
219  result.element({4,0,6,3}) = (rt*(1 + w)*(-1 + rC*w))/mSqq ;
220  result.element({4,0,6,4}) = (rt*(rC - w)*(1 + w))/mSqq ;
221  result.element({4,0,8,2}) = ((-1 + rC)*rt*(1 + w))/mSqq ;
222  result.element({4,0,8,3}) = (rt*(1 + w)*(-1 + rC*w))/mSqq ;
223  result.element({4,0,8,4}) = (rt*(rC - w)*(1 + w))/mSqq ;
224  result.element({4,1,2,1}) = 1 - w ;
225  result.element({4,1,4,1}) = -1 + w ;
226  result.element({4,1,6,5}) = ((1 + rC)*rt*(-1 + w))/mSqq ;
227  result.element({4,1,6,6}) = (rt*(-1 + w)*(-1 + rC*w))/mSqq ;
228  result.element({4,1,6,7}) = (rt*(rC - w)*(-1 + w))/mSqq ;
229  result.element({4,1,8,5}) = -(((1 + rC)*rt*(-1 + w))/mSqq) ;
230  result.element({4,1,8,6}) = -((rt*(-1 + w)*(-1 + rC*w))/mSqq) ;
231  result.element({4,1,8,7}) = -((rt*(rC - w)*(-1 + w))/mSqq) ;
232  result.element({5,2,0,2}) = (2.*(2.*mSqq2*(-1 + w) + mSqq*(1 + 2.*rC + rCSq + rtSq)*(-1 + w) + rtSq*(1 + 2.*w - 2.*rC*(2 + w) + rCSq*(1 + 2.*w))))/(3.*mSqq2) ;
233  result.element({5,2,0,3}) = ((1 + w)*(2.*mSqq*rC*(1 + rC)*(-1 + w) + rtSq*(3 - 2.*rC*(2 + w) + rCSq*(-1 + 4.*w))))/(3.*mSqq2) ;
234  result.element({5,2,0,4}) = ((1 + w)*(2.*mSqq*(1 + rC)*(-1 + w) + rtSq*(-1 + 3.*rCSq + 4.*w - 2.*rC*(2 + w))))/(3.*mSqq2) ;
235  result.element({5,2,1,0}) = ((-1 + rC)*rt*(1 + w))/mSqq ;
236  result.element({5,2,3,0}) = ((-1 + rC)*rt*(1 + w))/mSqq ;
237  result.element({5,2,5,2}) = (2.*(2.*mSqq2*(-1 + w) + mSqq*(1 + 2.*rC + rCSq + rtSq)*(-1 + w) + rtSq*(1 + 2.*w - 2.*rC*(2 + w) + rCSq*(1 + 2.*w))))/(3.*mSqq2) ;
238  result.element({5,2,5,3}) = ((1 + w)*(2.*mSqq*rC*(1 + rC)*(-1 + w) + rtSq*(3 - 2.*rC*(2 + w) + rCSq*(-1 + 4.*w))))/(3.*mSqq2) ;
239  result.element({5,2,5,4}) = ((1 + w)*(2.*mSqq*(1 + rC)*(-1 + w) + rtSq*(-1 + 3.*rCSq + 4.*w - 2.*rC*(2 + w))))/(3.*mSqq2) ;
240  result.element({5,2,7,2}) = (2.*(2.*mSqq2*(-1 + w) + mSqq*(1 + 2.*rC + rCSq + rtSq)*(-1 + w) + rtSq*(1 + 2.*w - 2.*rC*(2 + w) + rCSq*(1 + 2.*w))))/(3.*mSqq2) ;
241  result.element({5,2,7,3}) = ((1 + w)*(2.*mSqq*rC*(1 + rC)*(-1 + w) + rtSq*(3 - 2.*rC*(2 + w) + rCSq*(-1 + 4.*w))))/(3.*mSqq2) ;
242  result.element({5,2,7,4}) = ((1 + w)*(2.*mSqq*(1 + rC)*(-1 + w) + rtSq*(-1 + 3.*rCSq + 4.*w - 2.*rC*(2 + w))))/(3.*mSqq2) ;
243  result.element({5,2,9,8}) = (-12.*(1 + rC)*rt*(-1 + w))/mSqq ;
244  result.element({5,2,9,9}) = (-4.*rt*(-1 + w)*(-3 + rC*(-1 + 2.*w)))/mSqq ;
245  result.element({5,2,9,10}) = (-4.*rt*(1 + 3.*rC - 2.*w)*(-1 + w))/mSqq ;
246  result.element({5,2,9,11}) = (4.*(1 + rC)*rt*(-1 + wSq))/mSqq ;
247  result.element({5,3,0,2}) = ((1 + w)*(2.*mSqq*rC*(1 + rC)*(-1 + w) + rtSq*(3 - 2.*rC*(2 + w) + rCSq*(-1 + 4.*w))))/(3.*mSqq2) ;
248  result.element({5,3,0,3}) = ((1 + w)*(2.*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6.*rC*w + rCSq*(-1 + 4.*wSq))))/(3.*mSqq2) ;
249  result.element({5,3,0,4}) = ((1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
250  result.element({5,3,1,0}) = (rt*(1 + w)*(-1 + rC*w))/mSqq ;
251  result.element({5,3,3,0}) = (rt*(1 + w)*(-1 + rC*w))/mSqq ;
252  result.element({5,3,5,2}) = ((1 + w)*(2.*mSqq*rC*(1 + rC)*(-1 + w) + rtSq*(3 - 2.*rC*(2 + w) + rCSq*(-1 + 4.*w))))/(3.*mSqq2) ;
253  result.element({5,3,5,3}) = ((1 + w)*(2.*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6.*rC*w + rCSq*(-1 + 4.*wSq))))/(3.*mSqq2) ;
254  result.element({5,3,5,4}) = ((1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
255  result.element({5,3,7,2}) = ((1 + w)*(2.*mSqq*rC*(1 + rC)*(-1 + w) + rtSq*(3 - 2.*rC*(2 + w) + rCSq*(-1 + 4.*w))))/(3.*mSqq2) ;
256  result.element({5,3,7,3}) = ((1 + w)*(2.*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6.*rC*w + rCSq*(-1 + 4.*wSq))))/(3.*mSqq2) ;
257  result.element({5,3,7,4}) = ((1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
258  result.element({5,3,9,8}) = (-4.*rC*rt*(-1 + wSq))/mSqq ;
259  result.element({5,3,9,9}) = (4.*rC*rt*(-1 + wSq))/mSqq ;
260  result.element({5,3,9,10}) = (-4.*rC*rt*(-1 + wSq))/mSqq ;
261  result.element({5,3,9,11}) = (4.*rC*rt*(-1 + w)*pow(1 + w,2.))/mSqq ;
262  result.element({5,4,0,2}) = ((1 + w)*(2.*mSqq*(1 + rC)*(-1 + w) + rtSq*(-1 + 3.*rCSq + 4.*w - 2.*rC*(2 + w))))/(3.*mSqq2) ;
263  result.element({5,4,0,3}) = ((1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
264  result.element({5,4,0,4}) = ((1 + w)*(2.*mSqq*(-1 + wSq) + rtSq*(-1 + 3.*rCSq - 6.*rC*w + 4.*wSq)))/(3.*mSqq2) ;
265  result.element({5,4,1,0}) = (rt*(rC - w)*(1 + w))/mSqq ;
266  result.element({5,4,3,0}) = (rt*(rC - w)*(1 + w))/mSqq ;
267  result.element({5,4,5,2}) = ((1 + w)*(2.*mSqq*(1 + rC)*(-1 + w) + rtSq*(-1 + 3.*rCSq + 4.*w - 2.*rC*(2 + w))))/(3.*mSqq2) ;
268  result.element({5,4,5,3}) = ((1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
269  result.element({5,4,5,4}) = ((1 + w)*(2.*mSqq*(-1 + wSq) + rtSq*(-1 + 3.*rCSq - 6.*rC*w + 4.*wSq)))/(3.*mSqq2) ;
270  result.element({5,4,7,2}) = ((1 + w)*(2.*mSqq*(1 + rC)*(-1 + w) + rtSq*(-1 + 3.*rCSq + 4.*w - 2.*rC*(2 + w))))/(3.*mSqq2) ;
271  result.element({5,4,7,3}) = ((1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
272  result.element({5,4,7,4}) = ((1 + w)*(2.*mSqq*(-1 + wSq) + rtSq*(-1 + 3.*rCSq - 6.*rC*w + 4.*wSq)))/(3.*mSqq2) ;
273  result.element({5,4,9,8}) = (-4.*rt*(-1 + wSq))/mSqq ;
274  result.element({5,4,9,9}) = (4.*rt*(-1 + wSq))/mSqq ;
275  result.element({5,4,9,10}) = (-4.*rt*(-1 + wSq))/mSqq ;
276  result.element({5,4,9,11}) = (4.*rt*(-1 + w)*pow(1 + w,2.))/mSqq ;
277  result.element({5,5,0,5}) = (-2.*(2.*mSqq2*(1 + w) + mSqq*(1 - 2.*rC + rCSq + rtSq)*(1 + w) + rtSq*(-1 + 2.*rC*(-2 + w) + 2.*w + rCSq*(-1 + 2.*w))))/(3.*mSqq2) ;
278  result.element({5,5,0,6}) = -((-1 + w)*(2.*mSqq*(-1 + rC)*rC*(1 + w) + rtSq*(-3 + 2.*rC*(-2 + w) + rCSq*(1 + 4.*w))))/(3.*mSqq2) ;
279  result.element({5,5,0,7}) = -((-1 + w)*(rtSq*(-1 + 3.*rCSq - 2.*rC*(-2 + w) - 4.*w) + 2.*mSqq*(-1 + rC)*(1 + w)))/(3.*mSqq2) ;
280  result.element({5,5,1,1}) = -(((1 + rC)*rt*(-1 + w))/mSqq) ;
281  result.element({5,5,3,1}) = ((1 + rC)*rt*(-1 + w))/mSqq ;
282  result.element({5,5,5,5}) = (2.*(2.*mSqq2*(1 + w) + mSqq*(1 - 2.*rC + rCSq + rtSq)*(1 + w) + rtSq*(-1 + 2.*rC*(-2 + w) + 2.*w + rCSq*(-1 + 2.*w))))/(3.*mSqq2) ;
283  result.element({5,5,5,6}) = ((-1 + w)*(2.*mSqq*(-1 + rC)*rC*(1 + w) + rtSq*(-3 + 2.*rC*(-2 + w) + rCSq*(1 + 4.*w))))/(3.*mSqq2) ;
284  result.element({5,5,5,7}) = ((-1 + w)*(rtSq*(-1 + 3.*rCSq - 2.*rC*(-2 + w) - 4.*w) + 2.*mSqq*(-1 + rC)*(1 + w)))/(3.*mSqq2) ;
285  result.element({5,5,7,5}) = (-2.*(2.*mSqq2*(1 + w) + mSqq*(1 - 2.*rC + rCSq + rtSq)*(1 + w) + rtSq*(-1 + 2.*rC*(-2 + w) + 2.*w + rCSq*(-1 + 2.*w))))/(3.*mSqq2) ;
286  result.element({5,5,7,6}) = -((-1 + w)*(2.*mSqq*(-1 + rC)*rC*(1 + w) + rtSq*(-3 + 2.*rC*(-2 + w) + rCSq*(1 + 4.*w))))/(3.*mSqq2) ;
287  result.element({5,5,7,7}) = -((-1 + w)*(rtSq*(-1 + 3.*rCSq - 2.*rC*(-2 + w) - 4.*w) + 2.*mSqq*(-1 + rC)*(1 + w)))/(3.*mSqq2) ;
288  result.element({5,5,9,8}) = (12.*(-1 + rC)*rt*(1 + w))/mSqq ;
289  result.element({5,5,9,9}) = (8.*rC*rt*(-1 + wSq))/mSqq ;
290  result.element({5,5,9,10}) = (8.*rt*(-1 + wSq))/mSqq ;
291  result.element({5,6,0,5}) = -((-1 + w)*(2.*mSqq*(-1 + rC)*rC*(1 + w) + rtSq*(-3 + 2.*rC*(-2 + w) + rCSq*(1 + 4.*w))))/(3.*mSqq2) ;
292  result.element({5,6,0,6}) = -((-1 + w)*(2.*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6.*rC*w + rCSq*(-1 + 4.*wSq))))/(3.*mSqq2) ;
293  result.element({5,6,0,7}) = -((-1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
294  result.element({5,6,1,1}) = -((rt*(-1 + w)*(-1 + rC*w))/mSqq) ;
295  result.element({5,6,3,1}) = (rt*(-1 + w)*(-1 + rC*w))/mSqq ;
296  result.element({5,6,5,5}) = ((-1 + w)*(2.*mSqq*(-1 + rC)*rC*(1 + w) + rtSq*(-3 + 2.*rC*(-2 + w) + rCSq*(1 + 4.*w))))/(3.*mSqq2) ;
297  result.element({5,6,5,6}) = ((-1 + w)*(2.*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6.*rC*w + rCSq*(-1 + 4.*wSq))))/(3.*mSqq2) ;
298  result.element({5,6,5,7}) = ((-1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
299  result.element({5,6,7,5}) = -((-1 + w)*(2.*mSqq*(-1 + rC)*rC*(1 + w) + rtSq*(-3 + 2.*rC*(-2 + w) + rCSq*(1 + 4.*w))))/(3.*mSqq2) ;
300  result.element({5,6,7,6}) = -((-1 + w)*(2.*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6.*rC*w + rCSq*(-1 + 4.*wSq))))/(3.*mSqq2) ;
301  result.element({5,6,7,7}) = -((-1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
302  result.element({5,6,9,8}) = (4.*rC*rt*(-1 + wSq))/mSqq ;
303  result.element({5,7,0,5}) = -((-1 + w)*(rtSq*(-1 + 3.*rCSq - 2.*rC*(-2 + w) - 4.*w) + 2.*mSqq*(-1 + rC)*(1 + w)))/(3.*mSqq2) ;
304  result.element({5,7,0,6}) = -((-1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
305  result.element({5,7,0,7}) = -((-1 + w)*(2.*mSqq*(-1 + wSq) + rtSq*(-1 + 3.*rCSq - 6.*rC*w + 4.*wSq)))/(3.*mSqq2) ;
306  result.element({5,7,1,1}) = -((rt*(rC - w)*(-1 + w))/mSqq) ;
307  result.element({5,7,3,1}) = (rt*(rC - w)*(-1 + w))/mSqq ;
308  result.element({5,7,5,5}) = ((-1 + w)*(rtSq*(-1 + 3.*rCSq - 2.*rC*(-2 + w) - 4.*w) + 2.*mSqq*(-1 + rC)*(1 + w)))/(3.*mSqq2) ;
309  result.element({5,7,5,6}) = ((-1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
310  result.element({5,7,5,7}) = ((-1 + w)*(2.*mSqq*(-1 + wSq) + rtSq*(-1 + 3.*rCSq - 6.*rC*w + 4.*wSq)))/(3.*mSqq2) ;
311  result.element({5,7,7,5}) = -((-1 + w)*(rtSq*(-1 + 3.*rCSq - 2.*rC*(-2 + w) - 4.*w) + 2.*mSqq*(-1 + rC)*(1 + w)))/(3.*mSqq2) ;
312  result.element({5,7,7,6}) = -((-1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
313  result.element({5,7,7,7}) = -((-1 + w)*(2.*mSqq*(-1 + wSq) + rtSq*(-1 + 3.*rCSq - 6.*rC*w + 4.*wSq)))/(3.*mSqq2) ;
314  result.element({5,7,9,8}) = (4.*rt*(-1 + wSq))/mSqq ;
315  result.element({6,2,2,0}) = ((-1 + rC)*rt*(1 + w))/mSqq ;
316  result.element({6,2,4,0}) = ((-1 + rC)*rt*(1 + w))/mSqq ;
317  result.element({6,2,6,2}) = (2.*(2.*mSqq2*(-1 + w) + mSqq*(1 + 2.*rC + rCSq + rtSq)*(-1 + w) + rtSq*(1 + 2.*w - 2.*rC*(2 + w) + rCSq*(1 + 2.*w))))/(3.*mSqq2) ;
318  result.element({6,2,6,3}) = ((1 + w)*(2.*mSqq*rC*(1 + rC)*(-1 + w) + rtSq*(3 - 2.*rC*(2 + w) + rCSq*(-1 + 4.*w))))/(3.*mSqq2) ;
319  result.element({6,2,6,4}) = ((1 + w)*(2.*mSqq*(1 + rC)*(-1 + w) + rtSq*(-1 + 3.*rCSq + 4.*w - 2.*rC*(2 + w))))/(3.*mSqq2) ;
320  result.element({6,2,8,2}) = (2.*(2.*mSqq2*(-1 + w) + mSqq*(1 + 2.*rC + rCSq + rtSq)*(-1 + w) + rtSq*(1 + 2.*w - 2.*rC*(2 + w) + rCSq*(1 + 2.*w))))/(3.*mSqq2) ;
321  result.element({6,2,8,3}) = ((1 + w)*(2.*mSqq*rC*(1 + rC)*(-1 + w) + rtSq*(3 - 2.*rC*(2 + w) + rCSq*(-1 + 4.*w))))/(3.*mSqq2) ;
322  result.element({6,2,8,4}) = ((1 + w)*(2.*mSqq*(1 + rC)*(-1 + w) + rtSq*(-1 + 3.*rCSq + 4.*w - 2.*rC*(2 + w))))/(3.*mSqq2) ;
323  result.element({6,2,10,8}) = (-12.*(1 + rC)*rt*(-1 + w))/mSqq ;
324  result.element({6,2,10,9}) = (-4.*rt*(-1 + w)*(-3 + rC*(-1 + 2.*w)))/mSqq ;
325  result.element({6,2,10,10}) = (-4.*rt*(1 + 3.*rC - 2.*w)*(-1 + w))/mSqq ;
326  result.element({6,2,10,11}) = (4.*(1 + rC)*rt*(-1 + wSq))/mSqq ;
327  result.element({6,3,2,0}) = (rt*(1 + w)*(-1 + rC*w))/mSqq ;
328  result.element({6,3,4,0}) = (rt*(1 + w)*(-1 + rC*w))/mSqq ;
329  result.element({6,3,6,2}) = ((1 + w)*(2.*mSqq*rC*(1 + rC)*(-1 + w) + rtSq*(3 - 2.*rC*(2 + w) + rCSq*(-1 + 4.*w))))/(3.*mSqq2) ;
330  result.element({6,3,6,3}) = ((1 + w)*(2.*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6.*rC*w + rCSq*(-1 + 4.*wSq))))/(3.*mSqq2) ;
331  result.element({6,3,6,4}) = ((1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
332  result.element({6,3,8,2}) = ((1 + w)*(2.*mSqq*rC*(1 + rC)*(-1 + w) + rtSq*(3 - 2.*rC*(2 + w) + rCSq*(-1 + 4.*w))))/(3.*mSqq2) ;
333  result.element({6,3,8,3}) = ((1 + w)*(2.*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6.*rC*w + rCSq*(-1 + 4.*wSq))))/(3.*mSqq2) ;
334  result.element({6,3,8,4}) = ((1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
335  result.element({6,3,10,8}) = (-4.*rC*rt*(-1 + wSq))/mSqq ;
336  result.element({6,3,10,9}) = (4.*rC*rt*(-1 + wSq))/mSqq ;
337  result.element({6,3,10,10}) = (-4.*rC*rt*(-1 + wSq))/mSqq ;
338  result.element({6,3,10,11}) = (4.*rC*rt*(-1 + w)*pow(1 + w,2.))/mSqq ;
339  result.element({6,4,2,0}) = (rt*(rC - w)*(1 + w))/mSqq ;
340  result.element({6,4,4,0}) = (rt*(rC - w)*(1 + w))/mSqq ;
341  result.element({6,4,6,2}) = ((1 + w)*(2.*mSqq*(1 + rC)*(-1 + w) + rtSq*(-1 + 3.*rCSq + 4.*w - 2.*rC*(2 + w))))/(3.*mSqq2) ;
342  result.element({6,4,6,3}) = ((1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
343  result.element({6,4,6,4}) = ((1 + w)*(2.*mSqq*(-1 + wSq) + rtSq*(-1 + 3.*rCSq - 6.*rC*w + 4.*wSq)))/(3.*mSqq2) ;
344  result.element({6,4,8,2}) = ((1 + w)*(2.*mSqq*(1 + rC)*(-1 + w) + rtSq*(-1 + 3.*rCSq + 4.*w - 2.*rC*(2 + w))))/(3.*mSqq2) ;
345  result.element({6,4,8,3}) = ((1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
346  result.element({6,4,8,4}) = ((1 + w)*(2.*mSqq*(-1 + wSq) + rtSq*(-1 + 3.*rCSq - 6.*rC*w + 4.*wSq)))/(3.*mSqq2) ;
347  result.element({6,4,10,8}) = (-4.*rt*(-1 + wSq))/mSqq ;
348  result.element({6,4,10,9}) = (4.*rt*(-1 + wSq))/mSqq ;
349  result.element({6,4,10,10}) = (-4.*rt*(-1 + wSq))/mSqq ;
350  result.element({6,4,10,11}) = (4.*rt*(-1 + w)*pow(1 + w,2.))/mSqq ;
351  result.element({6,5,2,1}) = -(((1 + rC)*rt*(-1 + w))/mSqq) ;
352  result.element({6,5,4,1}) = ((1 + rC)*rt*(-1 + w))/mSqq ;
353  result.element({6,5,6,5}) = (2.*(2.*mSqq2*(1 + w) + mSqq*(1 - 2.*rC + rCSq + rtSq)*(1 + w) + rtSq*(-1 + 2.*rC*(-2 + w) + 2.*w + rCSq*(-1 + 2.*w))))/(3.*mSqq2) ;
354  result.element({6,5,6,6}) = ((-1 + w)*(2.*mSqq*(-1 + rC)*rC*(1 + w) + rtSq*(-3 + 2.*rC*(-2 + w) + rCSq*(1 + 4.*w))))/(3.*mSqq2) ;
355  result.element({6,5,6,7}) = ((-1 + w)*(rtSq*(-1 + 3.*rCSq - 2.*rC*(-2 + w) - 4.*w) + 2.*mSqq*(-1 + rC)*(1 + w)))/(3.*mSqq2) ;
356  result.element({6,5,8,5}) = (-2.*(2.*mSqq2*(1 + w) + mSqq*(1 - 2.*rC + rCSq + rtSq)*(1 + w) + rtSq*(-1 + 2.*rC*(-2 + w) + 2.*w + rCSq*(-1 + 2.*w))))/(3.*mSqq2) ;
357  result.element({6,5,8,6}) = -((-1 + w)*(2.*mSqq*(-1 + rC)*rC*(1 + w) + rtSq*(-3 + 2.*rC*(-2 + w) + rCSq*(1 + 4.*w))))/(3.*mSqq2) ;
358  result.element({6,5,8,7}) = -((-1 + w)*(rtSq*(-1 + 3.*rCSq - 2.*rC*(-2 + w) - 4.*w) + 2.*mSqq*(-1 + rC)*(1 + w)))/(3.*mSqq2) ;
359  result.element({6,5,10,8}) = (-12.*(-1 + rC)*rt*(1 + w))/mSqq ;
360  result.element({6,5,10,9}) = (-8.*rC*rt*(-1 + wSq))/mSqq ;
361  result.element({6,5,10,10}) = (-8.*rt*(-1 + wSq))/mSqq ;
362  result.element({6,6,2,1}) = -((rt*(-1 + w)*(-1 + rC*w))/mSqq) ;
363  result.element({6,6,4,1}) = (rt*(-1 + w)*(-1 + rC*w))/mSqq ;
364  result.element({6,6,6,5}) = ((-1 + w)*(2.*mSqq*(-1 + rC)*rC*(1 + w) + rtSq*(-3 + 2.*rC*(-2 + w) + rCSq*(1 + 4.*w))))/(3.*mSqq2) ;
365  result.element({6,6,6,6}) = ((-1 + w)*(2.*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6.*rC*w + rCSq*(-1 + 4.*wSq))))/(3.*mSqq2) ;
366  result.element({6,6,6,7}) = ((-1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
367  result.element({6,6,8,5}) = -((-1 + w)*(2.*mSqq*(-1 + rC)*rC*(1 + w) + rtSq*(-3 + 2.*rC*(-2 + w) + rCSq*(1 + 4.*w))))/(3.*mSqq2) ;
368  result.element({6,6,8,6}) = -((-1 + w)*(2.*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6.*rC*w + rCSq*(-1 + 4.*wSq))))/(3.*mSqq2) ;
369  result.element({6,6,8,7}) = -((-1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
370  result.element({6,6,10,8}) = (-4.*rC*rt*(-1 + wSq))/mSqq ;
371  result.element({6,7,2,1}) = -((rt*(rC - w)*(-1 + w))/mSqq) ;
372  result.element({6,7,4,1}) = (rt*(rC - w)*(-1 + w))/mSqq ;
373  result.element({6,7,6,5}) = ((-1 + w)*(rtSq*(-1 + 3.*rCSq - 2.*rC*(-2 + w) - 4.*w) + 2.*mSqq*(-1 + rC)*(1 + w)))/(3.*mSqq2) ;
374  result.element({6,7,6,6}) = ((-1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
375  result.element({6,7,6,7}) = ((-1 + w)*(2.*mSqq*(-1 + wSq) + rtSq*(-1 + 3.*rCSq - 6.*rC*w + 4.*wSq)))/(3.*mSqq2) ;
376  result.element({6,7,8,5}) = -((-1 + w)*(rtSq*(-1 + 3.*rCSq - 2.*rC*(-2 + w) - 4.*w) + 2.*mSqq*(-1 + rC)*(1 + w)))/(3.*mSqq2) ;
377  result.element({6,7,8,6}) = -((-1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
378  result.element({6,7,8,7}) = -((-1 + w)*(2.*mSqq*(-1 + wSq) + rtSq*(-1 + 3.*rCSq - 6.*rC*w + 4.*wSq)))/(3.*mSqq2) ;
379  result.element({6,7,10,8}) = (-4.*rt*(-1 + wSq))/mSqq ;
380  result.element({7,2,0,2}) = (2.*(2.*mSqq2*(-1 + w) + mSqq*(1 + 2.*rC + rCSq + rtSq)*(-1 + w) + rtSq*(1 + 2.*w - 2.*rC*(2 + w) + rCSq*(1 + 2.*w))))/(3.*mSqq2) ;
381  result.element({7,2,0,3}) = ((1 + w)*(2.*mSqq*rC*(1 + rC)*(-1 + w) + rtSq*(3 - 2.*rC*(2 + w) + rCSq*(-1 + 4.*w))))/(3.*mSqq2) ;
382  result.element({7,2,0,4}) = ((1 + w)*(2.*mSqq*(1 + rC)*(-1 + w) + rtSq*(-1 + 3.*rCSq + 4.*w - 2.*rC*(2 + w))))/(3.*mSqq2) ;
383  result.element({7,2,1,0}) = ((-1 + rC)*rt*(1 + w))/mSqq ;
384  result.element({7,2,3,0}) = ((-1 + rC)*rt*(1 + w))/mSqq ;
385  result.element({7,2,5,2}) = (2.*(2.*mSqq2*(-1 + w) + mSqq*(1 + 2.*rC + rCSq + rtSq)*(-1 + w) + rtSq*(1 + 2.*w - 2.*rC*(2 + w) + rCSq*(1 + 2.*w))))/(3.*mSqq2) ;
386  result.element({7,2,5,3}) = ((1 + w)*(2.*mSqq*rC*(1 + rC)*(-1 + w) + rtSq*(3 - 2.*rC*(2 + w) + rCSq*(-1 + 4.*w))))/(3.*mSqq2) ;
387  result.element({7,2,5,4}) = ((1 + w)*(2.*mSqq*(1 + rC)*(-1 + w) + rtSq*(-1 + 3.*rCSq + 4.*w - 2.*rC*(2 + w))))/(3.*mSqq2) ;
388  result.element({7,2,7,2}) = (2.*(2.*mSqq2*(-1 + w) + mSqq*(1 + 2.*rC + rCSq + rtSq)*(-1 + w) + rtSq*(1 + 2.*w - 2.*rC*(2 + w) + rCSq*(1 + 2.*w))))/(3.*mSqq2) ;
389  result.element({7,2,7,3}) = ((1 + w)*(2.*mSqq*rC*(1 + rC)*(-1 + w) + rtSq*(3 - 2.*rC*(2 + w) + rCSq*(-1 + 4.*w))))/(3.*mSqq2) ;
390  result.element({7,2,7,4}) = ((1 + w)*(2.*mSqq*(1 + rC)*(-1 + w) + rtSq*(-1 + 3.*rCSq + 4.*w - 2.*rC*(2 + w))))/(3.*mSqq2) ;
391  result.element({7,2,9,8}) = (-12.*(1 + rC)*rt*(-1 + w))/mSqq ;
392  result.element({7,2,9,9}) = (-4.*rt*(-1 + w)*(-3 + rC*(-1 + 2.*w)))/mSqq ;
393  result.element({7,2,9,10}) = (-4.*rt*(1 + 3.*rC - 2.*w)*(-1 + w))/mSqq ;
394  result.element({7,2,9,11}) = (4.*(1 + rC)*rt*(-1 + wSq))/mSqq ;
395  result.element({7,3,0,2}) = ((1 + w)*(2.*mSqq*rC*(1 + rC)*(-1 + w) + rtSq*(3 - 2.*rC*(2 + w) + rCSq*(-1 + 4.*w))))/(3.*mSqq2) ;
396  result.element({7,3,0,3}) = ((1 + w)*(2.*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6.*rC*w + rCSq*(-1 + 4.*wSq))))/(3.*mSqq2) ;
397  result.element({7,3,0,4}) = ((1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
398  result.element({7,3,1,0}) = (rt*(1 + w)*(-1 + rC*w))/mSqq ;
399  result.element({7,3,3,0}) = (rt*(1 + w)*(-1 + rC*w))/mSqq ;
400  result.element({7,3,5,2}) = ((1 + w)*(2.*mSqq*rC*(1 + rC)*(-1 + w) + rtSq*(3 - 2.*rC*(2 + w) + rCSq*(-1 + 4.*w))))/(3.*mSqq2) ;
401  result.element({7,3,5,3}) = ((1 + w)*(2.*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6.*rC*w + rCSq*(-1 + 4.*wSq))))/(3.*mSqq2) ;
402  result.element({7,3,5,4}) = ((1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
403  result.element({7,3,7,2}) = ((1 + w)*(2.*mSqq*rC*(1 + rC)*(-1 + w) + rtSq*(3 - 2.*rC*(2 + w) + rCSq*(-1 + 4.*w))))/(3.*mSqq2) ;
404  result.element({7,3,7,3}) = ((1 + w)*(2.*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6.*rC*w + rCSq*(-1 + 4.*wSq))))/(3.*mSqq2) ;
405  result.element({7,3,7,4}) = ((1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
406  result.element({7,3,9,8}) = (-4.*rC*rt*(-1 + wSq))/mSqq ;
407  result.element({7,3,9,9}) = (4.*rC*rt*(-1 + wSq))/mSqq ;
408  result.element({7,3,9,10}) = (-4.*rC*rt*(-1 + wSq))/mSqq ;
409  result.element({7,3,9,11}) = (4.*rC*rt*(-1 + w)*pow(1 + w,2.))/mSqq ;
410  result.element({7,4,0,2}) = ((1 + w)*(2.*mSqq*(1 + rC)*(-1 + w) + rtSq*(-1 + 3.*rCSq + 4.*w - 2.*rC*(2 + w))))/(3.*mSqq2) ;
411  result.element({7,4,0,3}) = ((1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
412  result.element({7,4,0,4}) = ((1 + w)*(2.*mSqq*(-1 + wSq) + rtSq*(-1 + 3.*rCSq - 6.*rC*w + 4.*wSq)))/(3.*mSqq2) ;
413  result.element({7,4,1,0}) = (rt*(rC - w)*(1 + w))/mSqq ;
414  result.element({7,4,3,0}) = (rt*(rC - w)*(1 + w))/mSqq ;
415  result.element({7,4,5,2}) = ((1 + w)*(2.*mSqq*(1 + rC)*(-1 + w) + rtSq*(-1 + 3.*rCSq + 4.*w - 2.*rC*(2 + w))))/(3.*mSqq2) ;
416  result.element({7,4,5,3}) = ((1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
417  result.element({7,4,5,4}) = ((1 + w)*(2.*mSqq*(-1 + wSq) + rtSq*(-1 + 3.*rCSq - 6.*rC*w + 4.*wSq)))/(3.*mSqq2) ;
418  result.element({7,4,7,2}) = ((1 + w)*(2.*mSqq*(1 + rC)*(-1 + w) + rtSq*(-1 + 3.*rCSq + 4.*w - 2.*rC*(2 + w))))/(3.*mSqq2) ;
419  result.element({7,4,7,3}) = ((1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
420  result.element({7,4,7,4}) = ((1 + w)*(2.*mSqq*(-1 + wSq) + rtSq*(-1 + 3.*rCSq - 6.*rC*w + 4.*wSq)))/(3.*mSqq2) ;
421  result.element({7,4,9,8}) = (-4.*rt*(-1 + wSq))/mSqq ;
422  result.element({7,4,9,9}) = (4.*rt*(-1 + wSq))/mSqq ;
423  result.element({7,4,9,10}) = (-4.*rt*(-1 + wSq))/mSqq ;
424  result.element({7,4,9,11}) = (4.*rt*(-1 + w)*pow(1 + w,2.))/mSqq ;
425  result.element({7,5,0,5}) = (2.*(2.*mSqq2*(1 + w) + mSqq*(1 - 2.*rC + rCSq + rtSq)*(1 + w) + rtSq*(-1 + 2.*rC*(-2 + w) + 2.*w + rCSq*(-1 + 2.*w))))/(3.*mSqq2) ;
426  result.element({7,5,0,6}) = ((-1 + w)*(2.*mSqq*(-1 + rC)*rC*(1 + w) + rtSq*(-3 + 2.*rC*(-2 + w) + rCSq*(1 + 4.*w))))/(3.*mSqq2) ;
427  result.element({7,5,0,7}) = ((-1 + w)*(rtSq*(-1 + 3.*rCSq - 2.*rC*(-2 + w) - 4.*w) + 2.*mSqq*(-1 + rC)*(1 + w)))/(3.*mSqq2) ;
428  result.element({7,5,1,1}) = ((1 + rC)*rt*(-1 + w))/mSqq ;
429  result.element({7,5,3,1}) = -(((1 + rC)*rt*(-1 + w))/mSqq) ;
430  result.element({7,5,5,5}) = (-2.*(2.*mSqq2*(1 + w) + mSqq*(1 - 2.*rC + rCSq + rtSq)*(1 + w) + rtSq*(-1 + 2.*rC*(-2 + w) + 2.*w + rCSq*(-1 + 2.*w))))/(3.*mSqq2) ;
431  result.element({7,5,5,6}) = -((-1 + w)*(2.*mSqq*(-1 + rC)*rC*(1 + w) + rtSq*(-3 + 2.*rC*(-2 + w) + rCSq*(1 + 4.*w))))/(3.*mSqq2) ;
432  result.element({7,5,5,7}) = -((-1 + w)*(rtSq*(-1 + 3.*rCSq - 2.*rC*(-2 + w) - 4.*w) + 2.*mSqq*(-1 + rC)*(1 + w)))/(3.*mSqq2) ;
433  result.element({7,5,7,5}) = (2.*(2.*mSqq2*(1 + w) + mSqq*(1 - 2.*rC + rCSq + rtSq)*(1 + w) + rtSq*(-1 + 2.*rC*(-2 + w) + 2.*w + rCSq*(-1 + 2.*w))))/(3.*mSqq2) ;
434  result.element({7,5,7,6}) = ((-1 + w)*(2.*mSqq*(-1 + rC)*rC*(1 + w) + rtSq*(-3 + 2.*rC*(-2 + w) + rCSq*(1 + 4.*w))))/(3.*mSqq2) ;
435  result.element({7,5,7,7}) = ((-1 + w)*(rtSq*(-1 + 3.*rCSq - 2.*rC*(-2 + w) - 4.*w) + 2.*mSqq*(-1 + rC)*(1 + w)))/(3.*mSqq2) ;
436  result.element({7,5,9,8}) = (-12.*(-1 + rC)*rt*(1 + w))/mSqq ;
437  result.element({7,5,9,9}) = (-8.*rC*rt*(-1 + wSq))/mSqq ;
438  result.element({7,5,9,10}) = (-8.*rt*(-1 + wSq))/mSqq ;
439  result.element({7,6,0,5}) = ((-1 + w)*(2.*mSqq*(-1 + rC)*rC*(1 + w) + rtSq*(-3 + 2.*rC*(-2 + w) + rCSq*(1 + 4.*w))))/(3.*mSqq2) ;
440  result.element({7,6,0,6}) = ((-1 + w)*(2.*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6.*rC*w + rCSq*(-1 + 4.*wSq))))/(3.*mSqq2) ;
441  result.element({7,6,0,7}) = ((-1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
442  result.element({7,6,1,1}) = (rt*(-1 + w)*(-1 + rC*w))/mSqq ;
443  result.element({7,6,3,1}) = -((rt*(-1 + w)*(-1 + rC*w))/mSqq) ;
444  result.element({7,6,5,5}) = -((-1 + w)*(2.*mSqq*(-1 + rC)*rC*(1 + w) + rtSq*(-3 + 2.*rC*(-2 + w) + rCSq*(1 + 4.*w))))/(3.*mSqq2) ;
445  result.element({7,6,5,6}) = -((-1 + w)*(2.*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6.*rC*w + rCSq*(-1 + 4.*wSq))))/(3.*mSqq2) ;
446  result.element({7,6,5,7}) = -((-1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
447  result.element({7,6,7,5}) = ((-1 + w)*(2.*mSqq*(-1 + rC)*rC*(1 + w) + rtSq*(-3 + 2.*rC*(-2 + w) + rCSq*(1 + 4.*w))))/(3.*mSqq2) ;
448  result.element({7,6,7,6}) = ((-1 + w)*(2.*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6.*rC*w + rCSq*(-1 + 4.*wSq))))/(3.*mSqq2) ;
449  result.element({7,6,7,7}) = ((-1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
450  result.element({7,6,9,8}) = (-4.*rC*rt*(-1 + wSq))/mSqq ;
451  result.element({7,7,0,5}) = ((-1 + w)*(rtSq*(-1 + 3.*rCSq - 2.*rC*(-2 + w) - 4.*w) + 2.*mSqq*(-1 + rC)*(1 + w)))/(3.*mSqq2) ;
452  result.element({7,7,0,6}) = ((-1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
453  result.element({7,7,0,7}) = ((-1 + w)*(2.*mSqq*(-1 + wSq) + rtSq*(-1 + 3.*rCSq - 6.*rC*w + 4.*wSq)))/(3.*mSqq2) ;
454  result.element({7,7,1,1}) = (rt*(rC - w)*(-1 + w))/mSqq ;
455  result.element({7,7,3,1}) = -((rt*(rC - w)*(-1 + w))/mSqq) ;
456  result.element({7,7,5,5}) = -((-1 + w)*(rtSq*(-1 + 3.*rCSq - 2.*rC*(-2 + w) - 4.*w) + 2.*mSqq*(-1 + rC)*(1 + w)))/(3.*mSqq2) ;
457  result.element({7,7,5,6}) = -((-1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
458  result.element({7,7,5,7}) = -((-1 + w)*(2.*mSqq*(-1 + wSq) + rtSq*(-1 + 3.*rCSq - 6.*rC*w + 4.*wSq)))/(3.*mSqq2) ;
459  result.element({7,7,7,5}) = ((-1 + w)*(rtSq*(-1 + 3.*rCSq - 2.*rC*(-2 + w) - 4.*w) + 2.*mSqq*(-1 + rC)*(1 + w)))/(3.*mSqq2) ;
460  result.element({7,7,7,6}) = ((-1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
461  result.element({7,7,7,7}) = ((-1 + w)*(2.*mSqq*(-1 + wSq) + rtSq*(-1 + 3.*rCSq - 6.*rC*w + 4.*wSq)))/(3.*mSqq2) ;
462  result.element({7,7,9,8}) = (-4.*rt*(-1 + wSq))/mSqq ;
463  result.element({8,2,2,0}) = ((-1 + rC)*rt*(1 + w))/mSqq ;
464  result.element({8,2,4,0}) = ((-1 + rC)*rt*(1 + w))/mSqq ;
465  result.element({8,2,6,2}) = (2.*(2.*mSqq2*(-1 + w) + mSqq*(1 + 2.*rC + rCSq + rtSq)*(-1 + w) + rtSq*(1 + 2.*w - 2.*rC*(2 + w) + rCSq*(1 + 2.*w))))/(3.*mSqq2) ;
466  result.element({8,2,6,3}) = ((1 + w)*(2.*mSqq*rC*(1 + rC)*(-1 + w) + rtSq*(3 - 2.*rC*(2 + w) + rCSq*(-1 + 4.*w))))/(3.*mSqq2) ;
467  result.element({8,2,6,4}) = ((1 + w)*(2.*mSqq*(1 + rC)*(-1 + w) + rtSq*(-1 + 3.*rCSq + 4.*w - 2.*rC*(2 + w))))/(3.*mSqq2) ;
468  result.element({8,2,8,2}) = (2.*(2.*mSqq2*(-1 + w) + mSqq*(1 + 2.*rC + rCSq + rtSq)*(-1 + w) + rtSq*(1 + 2.*w - 2.*rC*(2 + w) + rCSq*(1 + 2.*w))))/(3.*mSqq2) ;
469  result.element({8,2,8,3}) = ((1 + w)*(2.*mSqq*rC*(1 + rC)*(-1 + w) + rtSq*(3 - 2.*rC*(2 + w) + rCSq*(-1 + 4.*w))))/(3.*mSqq2) ;
470  result.element({8,2,8,4}) = ((1 + w)*(2.*mSqq*(1 + rC)*(-1 + w) + rtSq*(-1 + 3.*rCSq + 4.*w - 2.*rC*(2 + w))))/(3.*mSqq2) ;
471  result.element({8,2,10,8}) = (-12.*(1 + rC)*rt*(-1 + w))/mSqq ;
472  result.element({8,2,10,9}) = (-4.*rt*(-1 + w)*(-3 + rC*(-1 + 2.*w)))/mSqq ;
473  result.element({8,2,10,10}) = (-4.*rt*(1 + 3.*rC - 2.*w)*(-1 + w))/mSqq ;
474  result.element({8,2,10,11}) = (4.*(1 + rC)*rt*(-1 + wSq))/mSqq ;
475  result.element({8,3,2,0}) = (rt*(1 + w)*(-1 + rC*w))/mSqq ;
476  result.element({8,3,4,0}) = (rt*(1 + w)*(-1 + rC*w))/mSqq ;
477  result.element({8,3,6,2}) = ((1 + w)*(2.*mSqq*rC*(1 + rC)*(-1 + w) + rtSq*(3 - 2.*rC*(2 + w) + rCSq*(-1 + 4.*w))))/(3.*mSqq2) ;
478  result.element({8,3,6,3}) = ((1 + w)*(2.*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6.*rC*w + rCSq*(-1 + 4.*wSq))))/(3.*mSqq2) ;
479  result.element({8,3,6,4}) = ((1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
480  result.element({8,3,8,2}) = ((1 + w)*(2.*mSqq*rC*(1 + rC)*(-1 + w) + rtSq*(3 - 2.*rC*(2 + w) + rCSq*(-1 + 4.*w))))/(3.*mSqq2) ;
481  result.element({8,3,8,3}) = ((1 + w)*(2.*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6.*rC*w + rCSq*(-1 + 4.*wSq))))/(3.*mSqq2) ;
482  result.element({8,3,8,4}) = ((1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
483  result.element({8,3,10,8}) = (-4.*rC*rt*(-1 + wSq))/mSqq ;
484  result.element({8,3,10,9}) = (4.*rC*rt*(-1 + wSq))/mSqq ;
485  result.element({8,3,10,10}) = (-4.*rC*rt*(-1 + wSq))/mSqq ;
486  result.element({8,3,10,11}) = (4.*rC*rt*(-1 + w)*pow(1 + w,2.))/mSqq ;
487  result.element({8,4,2,0}) = (rt*(rC - w)*(1 + w))/mSqq ;
488  result.element({8,4,4,0}) = (rt*(rC - w)*(1 + w))/mSqq ;
489  result.element({8,4,6,2}) = ((1 + w)*(2.*mSqq*(1 + rC)*(-1 + w) + rtSq*(-1 + 3.*rCSq + 4.*w - 2.*rC*(2 + w))))/(3.*mSqq2) ;
490  result.element({8,4,6,3}) = ((1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
491  result.element({8,4,6,4}) = ((1 + w)*(2.*mSqq*(-1 + wSq) + rtSq*(-1 + 3.*rCSq - 6.*rC*w + 4.*wSq)))/(3.*mSqq2) ;
492  result.element({8,4,8,2}) = ((1 + w)*(2.*mSqq*(1 + rC)*(-1 + w) + rtSq*(-1 + 3.*rCSq + 4.*w - 2.*rC*(2 + w))))/(3.*mSqq2) ;
493  result.element({8,4,8,3}) = ((1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
494  result.element({8,4,8,4}) = ((1 + w)*(2.*mSqq*(-1 + wSq) + rtSq*(-1 + 3.*rCSq - 6.*rC*w + 4.*wSq)))/(3.*mSqq2) ;
495  result.element({8,4,10,8}) = (-4.*rt*(-1 + wSq))/mSqq ;
496  result.element({8,4,10,9}) = (4.*rt*(-1 + wSq))/mSqq ;
497  result.element({8,4,10,10}) = (-4.*rt*(-1 + wSq))/mSqq ;
498  result.element({8,4,10,11}) = (4.*rt*(-1 + w)*pow(1 + w,2.))/mSqq ;
499  result.element({8,5,2,1}) = ((1 + rC)*rt*(-1 + w))/mSqq ;
500  result.element({8,5,4,1}) = -(((1 + rC)*rt*(-1 + w))/mSqq) ;
501  result.element({8,5,6,5}) = (-2.*(2.*mSqq2*(1 + w) + mSqq*(1 - 2.*rC + rCSq + rtSq)*(1 + w) + rtSq*(-1 + 2.*rC*(-2 + w) + 2.*w + rCSq*(-1 + 2.*w))))/(3.*mSqq2) ;
502  result.element({8,5,6,6}) = -((-1 + w)*(2.*mSqq*(-1 + rC)*rC*(1 + w) + rtSq*(-3 + 2.*rC*(-2 + w) + rCSq*(1 + 4.*w))))/(3.*mSqq2) ;
503  result.element({8,5,6,7}) = -((-1 + w)*(rtSq*(-1 + 3.*rCSq - 2.*rC*(-2 + w) - 4.*w) + 2.*mSqq*(-1 + rC)*(1 + w)))/(3.*mSqq2) ;
504  result.element({8,5,8,5}) = (2.*(2.*mSqq2*(1 + w) + mSqq*(1 - 2.*rC + rCSq + rtSq)*(1 + w) + rtSq*(-1 + 2.*rC*(-2 + w) + 2.*w + rCSq*(-1 + 2.*w))))/(3.*mSqq2) ;
505  result.element({8,5,8,6}) = ((-1 + w)*(2.*mSqq*(-1 + rC)*rC*(1 + w) + rtSq*(-3 + 2.*rC*(-2 + w) + rCSq*(1 + 4.*w))))/(3.*mSqq2) ;
506  result.element({8,5,8,7}) = ((-1 + w)*(rtSq*(-1 + 3.*rCSq - 2.*rC*(-2 + w) - 4.*w) + 2.*mSqq*(-1 + rC)*(1 + w)))/(3.*mSqq2) ;
507  result.element({8,5,10,8}) = (12.*(-1 + rC)*rt*(1 + w))/mSqq ;
508  result.element({8,5,10,9}) = (8.*rC*rt*(-1 + wSq))/mSqq ;
509  result.element({8,5,10,10}) = (8.*rt*(-1 + wSq))/mSqq ;
510  result.element({8,6,2,1}) = (rt*(-1 + w)*(-1 + rC*w))/mSqq ;
511  result.element({8,6,4,1}) = -((rt*(-1 + w)*(-1 + rC*w))/mSqq) ;
512  result.element({8,6,6,5}) = -((-1 + w)*(2.*mSqq*(-1 + rC)*rC*(1 + w) + rtSq*(-3 + 2.*rC*(-2 + w) + rCSq*(1 + 4.*w))))/(3.*mSqq2) ;
513  result.element({8,6,6,6}) = -((-1 + w)*(2.*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6.*rC*w + rCSq*(-1 + 4.*wSq))))/(3.*mSqq2) ;
514  result.element({8,6,6,7}) = -((-1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
515  result.element({8,6,8,5}) = ((-1 + w)*(2.*mSqq*(-1 + rC)*rC*(1 + w) + rtSq*(-3 + 2.*rC*(-2 + w) + rCSq*(1 + 4.*w))))/(3.*mSqq2) ;
516  result.element({8,6,8,6}) = ((-1 + w)*(2.*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6.*rC*w + rCSq*(-1 + 4.*wSq))))/(3.*mSqq2) ;
517  result.element({8,6,8,7}) = ((-1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
518  result.element({8,6,10,8}) = (4.*rC*rt*(-1 + wSq))/mSqq ;
519  result.element({8,7,2,1}) = (rt*(rC - w)*(-1 + w))/mSqq ;
520  result.element({8,7,4,1}) = -((rt*(rC - w)*(-1 + w))/mSqq) ;
521  result.element({8,7,6,5}) = -((-1 + w)*(rtSq*(-1 + 3.*rCSq - 2.*rC*(-2 + w) - 4.*w) + 2.*mSqq*(-1 + rC)*(1 + w)))/(3.*mSqq2) ;
522  result.element({8,7,6,6}) = -((-1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
523  result.element({8,7,6,7}) = -((-1 + w)*(2.*mSqq*(-1 + wSq) + rtSq*(-1 + 3.*rCSq - 6.*rC*w + 4.*wSq)))/(3.*mSqq2) ;
524  result.element({8,7,8,5}) = ((-1 + w)*(rtSq*(-1 + 3.*rCSq - 2.*rC*(-2 + w) - 4.*w) + 2.*mSqq*(-1 + rC)*(1 + w)))/(3.*mSqq2) ;
525  result.element({8,7,8,6}) = ((-1 + w)*(2.*mSqq*rC*(-1 + wSq) + rtSq*(3.*w + 3.*rCSq*w - 2.*rC*(2 + wSq))))/(3.*mSqq2) ;
526  result.element({8,7,8,7}) = ((-1 + w)*(2.*mSqq*(-1 + wSq) + rtSq*(-1 + 3.*rCSq - 6.*rC*w + 4.*wSq)))/(3.*mSqq2) ;
527  result.element({8,7,10,8}) = (4.*rt*(-1 + wSq))/mSqq ;
528  result.element({9,8,0,2}) = (-12.*(1 + rC)*rt*(-1 + w))/mSqq ;
529  result.element({9,8,0,3}) = (-4.*rC*rt*(-1 + wSq))/mSqq ;
530  result.element({9,8,0,4}) = (-4.*rt*(-1 + wSq))/mSqq ;
531  result.element({9,8,0,5}) = (-12.*(-1 + rC)*rt*(1 + w))/mSqq ;
532  result.element({9,8,0,6}) = (-4.*rC*rt*(-1 + wSq))/mSqq ;
533  result.element({9,8,0,7}) = (-4.*rt*(-1 + wSq))/mSqq ;
534  result.element({9,8,5,2}) = (-12.*(1 + rC)*rt*(-1 + w))/mSqq ;
535  result.element({9,8,5,3}) = (-4.*rC*rt*(-1 + wSq))/mSqq ;
536  result.element({9,8,5,4}) = (-4.*rt*(-1 + wSq))/mSqq ;
537  result.element({9,8,5,5}) = (12.*(-1 + rC)*rt*(1 + w))/mSqq ;
538  result.element({9,8,5,6}) = (4.*rC*rt*(-1 + wSq))/mSqq ;
539  result.element({9,8,5,7}) = (4.*rt*(-1 + wSq))/mSqq ;
540  result.element({9,8,7,2}) = (-12.*(1 + rC)*rt*(-1 + w))/mSqq ;
541  result.element({9,8,7,3}) = (-4.*rC*rt*(-1 + wSq))/mSqq ;
542  result.element({9,8,7,4}) = (-4.*rt*(-1 + wSq))/mSqq ;
543  result.element({9,8,7,5}) = (-12.*(-1 + rC)*rt*(1 + w))/mSqq ;
544  result.element({9,8,7,6}) = (-4.*rC*rt*(-1 + wSq))/mSqq ;
545  result.element({9,8,7,7}) = (-4.*rt*(-1 + wSq))/mSqq ;
546  result.element({9,8,9,8}) = (32.*(mSqq + 2.*rtSq)*(-4.*rC + (2 + mSqq)*w + 2.*rCSq*w))/(3.*mSqq2) ;
547  result.element({9,8,9,9}) = (-16.*(mSqq + 2.*rtSq)*(-1 + w)*(2 + mSqq + 4.*rC - 2.*rCSq*(1 + 2.*w)))/(3.*mSqq2) ;
548  result.element({9,8,9,10}) = (16.*(mSqq + 2.*rtSq)*(-2 + mSqq + 4.*rC + 2.*rCSq - 4.*w)*(-1 + w))/(3.*mSqq2) ;
549  result.element({9,8,9,11}) = (-16.*(mSqq + 2.*rtSq)*(-1 + wSq))/(3.*mSqq) ;
550  result.element({9,9,0,2}) = (-4.*rt*(-1 + w)*(-3 + rC*(-1 + 2.*w)))/mSqq ;
551  result.element({9,9,0,3}) = (4.*rC*rt*(-1 + wSq))/mSqq ;
552  result.element({9,9,0,4}) = (4.*rt*(-1 + wSq))/mSqq ;
553  result.element({9,9,0,5}) = (-8.*rC*rt*(-1 + wSq))/mSqq ;
554  result.element({9,9,5,2}) = (-4.*rt*(-1 + w)*(-3 + rC*(-1 + 2.*w)))/mSqq ;
555  result.element({9,9,5,3}) = (4.*rC*rt*(-1 + wSq))/mSqq ;
556  result.element({9,9,5,4}) = (4.*rt*(-1 + wSq))/mSqq ;
557  result.element({9,9,5,5}) = (8.*rC*rt*(-1 + wSq))/mSqq ;
558  result.element({9,9,7,2}) = (-4.*rt*(-1 + w)*(-3 + rC*(-1 + 2.*w)))/mSqq ;
559  result.element({9,9,7,3}) = (4.*rC*rt*(-1 + wSq))/mSqq ;
560  result.element({9,9,7,4}) = (4.*rt*(-1 + wSq))/mSqq ;
561  result.element({9,9,7,5}) = (-8.*rC*rt*(-1 + wSq))/mSqq ;
562  result.element({9,9,9,8}) = (-16.*(mSqq + 2.*rtSq)*(-1 + w)*(2 + mSqq + 4.*rC - 2.*rCSq*(1 + 2.*w)))/(3.*mSqq2) ;
563  result.element({9,9,9,9}) = (16.*(mSqq + 2.*rtSq)*(-1 + w)*(2 + mSqq - 4.*rC*w + rCSq*(-2 + 4.*wSq)))/(3.*mSqq2) ;
564  result.element({9,9,9,10}) = (-16.*(mSqq + 2.*rtSq)*(-1 + w)*(mSqq - 2.*(-2.*rC + w + rCSq*w)))/(3.*mSqq2) ;
565  result.element({9,9,9,11}) = (16.*(mSqq + 2.*rtSq)*(-1 + wSq))/(3.*mSqq) ;
566  result.element({9,10,0,2}) = (-4.*rt*(1 + 3.*rC - 2.*w)*(-1 + w))/mSqq ;
567  result.element({9,10,0,3}) = (-4.*rC*rt*(-1 + wSq))/mSqq ;
568  result.element({9,10,0,4}) = (-4.*rt*(-1 + wSq))/mSqq ;
569  result.element({9,10,0,5}) = (-8.*rt*(-1 + wSq))/mSqq ;
570  result.element({9,10,5,2}) = (-4.*rt*(1 + 3.*rC - 2.*w)*(-1 + w))/mSqq ;
571  result.element({9,10,5,3}) = (-4.*rC*rt*(-1 + wSq))/mSqq ;
572  result.element({9,10,5,4}) = (-4.*rt*(-1 + wSq))/mSqq ;
573  result.element({9,10,5,5}) = (8.*rt*(-1 + wSq))/mSqq ;
574  result.element({9,10,7,2}) = (-4.*rt*(1 + 3.*rC - 2.*w)*(-1 + w))/mSqq ;
575  result.element({9,10,7,3}) = (-4.*rC*rt*(-1 + wSq))/mSqq ;
576  result.element({9,10,7,4}) = (-4.*rt*(-1 + wSq))/mSqq ;
577  result.element({9,10,7,5}) = (-8.*rt*(-1 + wSq))/mSqq ;
578  result.element({9,10,9,8}) = (16.*(mSqq + 2.*rtSq)*(-2 + mSqq + 4.*rC + 2.*rCSq - 4.*w)*(-1 + w))/(3.*mSqq2) ;
579  result.element({9,10,9,9}) = (-16.*(mSqq + 2.*rtSq)*(-1 + w)*(mSqq - 2.*(-2.*rC + w + rCSq*w)))/(3.*mSqq2) ;
580  result.element({9,10,9,10}) = (16.*(mSqq + 2.*rtSq)*(-1 + w)*(-2 + mSqq + 2.*rCSq - 4.*rC*w + 4.*wSq))/(3.*mSqq2) ;
581  result.element({9,10,9,11}) = (-16.*(mSqq + 2.*rtSq)*(-1 + wSq))/(3.*mSqq) ;
582  result.element({9,11,0,2}) = (4.*(1 + rC)*rt*(-1 + wSq))/mSqq ;
583  result.element({9,11,0,3}) = (4.*rC*rt*(-1 + w)*pow(1 + w,2.))/mSqq ;
584  result.element({9,11,0,4}) = (4.*rt*(-1 + w)*pow(1 + w,2.))/mSqq ;
585  result.element({9,11,5,2}) = (4.*(1 + rC)*rt*(-1 + wSq))/mSqq ;
586  result.element({9,11,5,3}) = (4.*rC*rt*(-1 + w)*pow(1 + w,2.))/mSqq ;
587  result.element({9,11,5,4}) = (4.*rt*(-1 + w)*pow(1 + w,2.))/mSqq ;
588  result.element({9,11,7,2}) = (4.*(1 + rC)*rt*(-1 + wSq))/mSqq ;
589  result.element({9,11,7,3}) = (4.*rC*rt*(-1 + w)*pow(1 + w,2.))/mSqq ;
590  result.element({9,11,7,4}) = (4.*rt*(-1 + w)*pow(1 + w,2.))/mSqq ;
591  result.element({9,11,9,8}) = (-16.*(mSqq + 2.*rtSq)*(-1 + wSq))/(3.*mSqq) ;
592  result.element({9,11,9,9}) = (16.*(mSqq + 2.*rtSq)*(-1 + wSq))/(3.*mSqq) ;
593  result.element({9,11,9,10}) = (-16.*(mSqq + 2.*rtSq)*(-1 + wSq))/(3.*mSqq) ;
594  result.element({9,11,9,11}) = (16.*(mSqq + 2.*rtSq)*(-1 + w)*pow(1 + w,2.))/(3.*mSqq) ;
595  result.element({10,8,6,2}) = (-12.*(1 + rC)*rt*(-1 + w))/mSqq ;
596  result.element({10,8,6,3}) = (-4.*rC*rt*(-1 + wSq))/mSqq ;
597  result.element({10,8,6,4}) = (-4.*rt*(-1 + wSq))/mSqq ;
598  result.element({10,8,6,5}) = (-12.*(-1 + rC)*rt*(1 + w))/mSqq ;
599  result.element({10,8,6,6}) = (-4.*rC*rt*(-1 + wSq))/mSqq ;
600  result.element({10,8,6,7}) = (-4.*rt*(-1 + wSq))/mSqq ;
601  result.element({10,8,8,2}) = (-12.*(1 + rC)*rt*(-1 + w))/mSqq ;
602  result.element({10,8,8,3}) = (-4.*rC*rt*(-1 + wSq))/mSqq ;
603  result.element({10,8,8,4}) = (-4.*rt*(-1 + wSq))/mSqq ;
604  result.element({10,8,8,5}) = (12.*(-1 + rC)*rt*(1 + w))/mSqq ;
605  result.element({10,8,8,6}) = (4.*rC*rt*(-1 + wSq))/mSqq ;
606  result.element({10,8,8,7}) = (4.*rt*(-1 + wSq))/mSqq ;
607  result.element({10,8,10,8}) = (32.*(mSqq + 2.*rtSq)*(-4.*rC + (2 + mSqq)*w + 2.*rCSq*w))/(3.*mSqq2) ;
608  result.element({10,8,10,9}) = (-16.*(mSqq + 2.*rtSq)*(-1 + w)*(2 + mSqq + 4.*rC - 2.*rCSq*(1 + 2.*w)))/(3.*mSqq2) ;
609  result.element({10,8,10,10}) = (16.*(mSqq + 2.*rtSq)*(-2 + mSqq + 4.*rC + 2.*rCSq - 4.*w)*(-1 + w))/(3.*mSqq2) ;
610  result.element({10,8,10,11}) = (-16.*(mSqq + 2.*rtSq)*(-1 + wSq))/(3.*mSqq) ;
611  result.element({10,9,6,2}) = (-4.*rt*(-1 + w)*(-3 + rC*(-1 + 2.*w)))/mSqq ;
612  result.element({10,9,6,3}) = (4.*rC*rt*(-1 + wSq))/mSqq ;
613  result.element({10,9,6,4}) = (4.*rt*(-1 + wSq))/mSqq ;
614  result.element({10,9,6,5}) = (-8.*rC*rt*(-1 + wSq))/mSqq ;
615  result.element({10,9,8,2}) = (-4.*rt*(-1 + w)*(-3 + rC*(-1 + 2.*w)))/mSqq ;
616  result.element({10,9,8,3}) = (4.*rC*rt*(-1 + wSq))/mSqq ;
617  result.element({10,9,8,4}) = (4.*rt*(-1 + wSq))/mSqq ;
618  result.element({10,9,8,5}) = (8.*rC*rt*(-1 + wSq))/mSqq ;
619  result.element({10,9,10,8}) = (-16.*(mSqq + 2.*rtSq)*(-1 + w)*(2 + mSqq + 4.*rC - 2.*rCSq*(1 + 2.*w)))/(3.*mSqq2) ;
620  result.element({10,9,10,9}) = (16.*(mSqq + 2.*rtSq)*(-1 + w)*(2 + mSqq - 4.*rC*w + rCSq*(-2 + 4.*wSq)))/(3.*mSqq2) ;
621  result.element({10,9,10,10}) = (-16.*(mSqq + 2.*rtSq)*(-1 + w)*(mSqq - 2.*(-2.*rC + w + rCSq*w)))/(3.*mSqq2) ;
622  result.element({10,9,10,11}) = (16.*(mSqq + 2.*rtSq)*(-1 + wSq))/(3.*mSqq) ;
623  result.element({10,10,6,2}) = (-4.*rt*(1 + 3.*rC - 2.*w)*(-1 + w))/mSqq ;
624  result.element({10,10,6,3}) = (-4.*rC*rt*(-1 + wSq))/mSqq ;
625  result.element({10,10,6,4}) = (-4.*rt*(-1 + wSq))/mSqq ;
626  result.element({10,10,6,5}) = (-8.*rt*(-1 + wSq))/mSqq ;
627  result.element({10,10,8,2}) = (-4.*rt*(1 + 3.*rC - 2.*w)*(-1 + w))/mSqq ;
628  result.element({10,10,8,3}) = (-4.*rC*rt*(-1 + wSq))/mSqq ;
629  result.element({10,10,8,4}) = (-4.*rt*(-1 + wSq))/mSqq ;
630  result.element({10,10,8,5}) = (8.*rt*(-1 + wSq))/mSqq ;
631  result.element({10,10,10,8}) = (16.*(mSqq + 2.*rtSq)*(-2 + mSqq + 4.*rC + 2.*rCSq - 4.*w)*(-1 + w))/(3.*mSqq2) ;
632  result.element({10,10,10,9}) = (-16.*(mSqq + 2.*rtSq)*(-1 + w)*(mSqq - 2.*(-2.*rC + w + rCSq*w)))/(3.*mSqq2) ;
633  result.element({10,10,10,10}) = (16.*(mSqq + 2.*rtSq)*(-1 + w)*(-2 + mSqq + 2.*rCSq - 4.*rC*w + 4.*wSq))/(3.*mSqq2) ;
634  result.element({10,10,10,11}) = (-16.*(mSqq + 2.*rtSq)*(-1 + wSq))/(3.*mSqq) ;
635  result.element({10,11,6,2}) = (4.*(1 + rC)*rt*(-1 + wSq))/mSqq ;
636  result.element({10,11,6,3}) = (4.*rC*rt*(-1 + w)*pow(1 + w,2.))/mSqq ;
637  result.element({10,11,6,4}) = (4.*rt*(-1 + w)*pow(1 + w,2.))/mSqq ;
638  result.element({10,11,8,2}) = (4.*(1 + rC)*rt*(-1 + wSq))/mSqq ;
639  result.element({10,11,8,3}) = (4.*rC*rt*(-1 + w)*pow(1 + w,2.))/mSqq ;
640  result.element({10,11,8,4}) = (4.*rt*(-1 + w)*pow(1 + w,2.))/mSqq ;
641  result.element({10,11,10,8}) = (-16.*(mSqq + 2.*rtSq)*(-1 + wSq))/(3.*mSqq) ;
642  result.element({10,11,10,9}) = (16.*(mSqq + 2.*rtSq)*(-1 + wSq))/(3.*mSqq) ;
643  result.element({10,11,10,10}) = (-16.*(mSqq + 2.*rtSq)*(-1 + wSq))/(3.*mSqq) ;
644  result.element({10,11,10,11}) = (16.*(mSqq + 2.*rtSq)*(-1 + w)*pow(1 + w,2.))/(3.*mSqq) ;
645 
646  result *= (RateNorm);
647 
648  return result;
649  }
650 
651 } // namespace Hammer
Hammer::MultiDimensional::makeEmptySparse
TensorData makeEmptySparse(const IndexList &dimensions, const LabelsList &labels)
Definition: SparseContainer.cc:287
Hammer::Tensor::labels
LabelsList labels() const
get the labels of all the indices at once
Definition: Tensor.cc:83
Hammer::GFermi
static constexpr double GFermi
Definition: Constants.hh:29
IndexLabels.hh
Tensor indices label definitions.
Hammer::pi3
static constexpr double pi3
Definition: Constants.hh:25
SparseContainer.hh
Sparse tensor data container.
Hammer::Tensor
Multidimensional tensor class with complex numbers as elements.
Definition: Tensor.hh:33
Constants.hh
Various numerical constants.
Pdg.hh
Hammer particle data class.
RateLbLcLepNu.hh
total rate
Hammer::PS::makeQ2Function
BoundaryFunction makeQ2Function(double qmin, double qmax)
Definition: RateBase.hh:122
Particle.hh
Hammer particle class.