template
|
template
|
struct X {
|
struct X {
|
double & f (const unsigned int i,
|
double & f (const unsigned int i,
|
const unsigned int j);
|
const unsigned int j);
|
void g (X &M);
|
void g (X &M);
|
};
|
};
|
|
|
template
|
template
|
void X::g (X &x) {
|
void X::g (X &x) {
|
double t14 = x.f(0,0)*x.f(1,1);
|
double t14 = x.f(0,0)*x.f(1,1);
|
double t15 = x.f(2,2)*x.f(3,3);
|
double t15 = x.f(2,2)*x.f(3,3);
|
double t17 = x.f(2,3)*x.f(3,2);
|
double t17 = x.f(2,3)*x.f(3,2);
|
double t19 = x.f(0,0)*x.f(2,1);
|
double t19 = x.f(0,0)*x.f(2,1);
|
double t20 = x.f(1,2)*x.f(3,3);
|
double t20 = x.f(1,2)*x.f(3,3);
|
double t22 = x.f(1,3)*x.f(3,2);
|
double t22 = x.f(1,3)*x.f(3,2);
|
double t24 = x.f(0,0)*x.f(3,1);
|
double t24 = x.f(0,0)*x.f(3,1);
|
double t25 = x.f(1,2)*x.f(2,3);
|
double t25 = x.f(1,2)*x.f(2,3);
|
double t27 = x.f(1,3)*x.f(2,2);
|
double t27 = x.f(1,3)*x.f(2,2);
|
double t29 = x.f(1,0)*x.f(0,1);
|
double t29 = x.f(1,0)*x.f(0,1);
|
double t32 = x.f(1,0)*x.f(2,1);
|
double t32 = x.f(1,0)*x.f(2,1);
|
double t33 = x.f(0,2)*x.f(3,3);
|
double t33 = x.f(0,2)*x.f(3,3);
|
double t35 = x.f(0,3)*x.f(3,2);
|
double t35 = x.f(0,3)*x.f(3,2);
|
double t37 = x.f(1,0)*x.f(3,1);
|
double t37 = x.f(1,0)*x.f(3,1);
|
double t38 = x.f(0,2)*x.f(2,3);
|
double t38 = x.f(0,2)*x.f(2,3);
|
double t40 = x.f(0,3)*x.f(2,2);
|
double t40 = x.f(0,3)*x.f(2,2);
|
double t42 = t14*t15-t14*t17-t19*t20+t19*t22+
|
double t42 = t14*t15-t14*t17-t19*t20+t19*t22+
|
t24*t25-t24*t27-t29*t15+t29*t17+
|
t24*t25-t24*t27-t29*t15+t29*t17+
|
t32*t33-t32*t35-t37*t38+t37*t40;
|
t32*t33-t32*t35-t37*t38+t37*t40;
|
double t43 = x.f(2,0)*x.f(0,1);
|
double t43 = x.f(2,0)*x.f(0,1);
|
double t46 = x.f(2,0)*x.f(1,1);
|
double t46 = x.f(2,0)*x.f(1,1);
|
double t49 = x.f(2,0)*x.f(3,1);
|
double t49 = x.f(2,0)*x.f(3,1);
|
double t50 = x.f(0,2)*x.f(1,3);
|
double t50 = x.f(0,2)*x.f(1,3);
|
double t52 = x.f(0,3)*x.f(1,2);
|
double t52 = x.f(0,3)*x.f(1,2);
|
double t54 = x.f(3,0)*x.f(0,1);
|
double t54 = x.f(3,0)*x.f(0,1);
|
double t57 = x.f(3,0)*x.f(1,1);
|
double t57 = x.f(3,0)*x.f(1,1);
|
double t60 = x.f(3,0)*x.f(2,1);
|
double t60 = x.f(3,0)*x.f(2,1);
|
double t63 = t43*t20-t43*t22-t46*t33+t46*t35+
|
double t63 = t43*t20-t43*t22-t46*t33+t46*t35+
|
t49*t50-t49*t52-t54*t25+t54*t27+
|
t49*t50-t49*t52-t54*t25+t54*t27+
|
t57*t38-t57*t40-t60*t50+t60*t52;
|
t57*t38-t57*t40-t60*t50+t60*t52;
|
double t65 = 1/(t42+t63);
|
double t65 = 1/(t42+t63);
|
double t71 = x.f(0,2)*x.f(2,1);
|
double t71 = x.f(0,2)*x.f(2,1);
|
double t73 = x.f(0,3)*x.f(2,1);
|
double t73 = x.f(0,3)*x.f(2,1);
|
double t75 = x.f(0,2)*x.f(3,1);
|
double t75 = x.f(0,2)*x.f(3,1);
|
double t77 = x.f(0,3)*x.f(3,1);
|
double t77 = x.f(0,3)*x.f(3,1);
|
double t81 = x.f(0,1)*x.f(1,2);
|
double t81 = x.f(0,1)*x.f(1,2);
|
double t83 = x.f(0,1)*x.f(1,3);
|
double t83 = x.f(0,1)*x.f(1,3);
|
double t85 = x.f(0,2)*x.f(1,1);
|
double t85 = x.f(0,2)*x.f(1,1);
|
double t87 = x.f(0,3)*x.f(1,1);
|
double t87 = x.f(0,3)*x.f(1,1);
|
double t101 = x.f(1,0)*x.f(2,2);
|
double t101 = x.f(1,0)*x.f(2,2);
|
double t103 = x.f(1,0)*x.f(2,3);
|
double t103 = x.f(1,0)*x.f(2,3);
|
double t105 = x.f(2,0)*x.f(1,2);
|
double t105 = x.f(2,0)*x.f(1,2);
|
double t107 = x.f(2,0)*x.f(1,3);
|
double t107 = x.f(2,0)*x.f(1,3);
|
double t109 = x.f(3,0)*x.f(1,2);
|
double t109 = x.f(3,0)*x.f(1,2);
|
double t111 = x.f(3,0)*x.f(1,3);
|
double t111 = x.f(3,0)*x.f(1,3);
|
double t115 = x.f(0,0)*x.f(2,2);
|
double t115 = x.f(0,0)*x.f(2,2);
|
double t117 = x.f(0,0)*x.f(2,3);
|
double t117 = x.f(0,0)*x.f(2,3);
|
double t119 = x.f(2,0)*x.f(0,2);
|
double t119 = x.f(2,0)*x.f(0,2);
|
double t121 = x.f(2,0)*x.f(0,3);
|
double t121 = x.f(2,0)*x.f(0,3);
|
}
|
}
|
|
|
template void X::g (X&);
|
template void X::g (X&);
|
|
|