template<typename T>
struct X {
    double & f (const unsigned int i,
                 const unsigned int j);
    void g (X &M);
};

template <typename T>
void X<T>::g (X &x) {
        double t14 = x.f(0,0)*x.f(1,1);
        double t15 = x.f(2,2)*x.f(3,3);
        double t17 = x.f(2,3)*x.f(3,2);
        double t19 = x.f(0,0)*x.f(2,1);
        double t20 = x.f(1,2)*x.f(3,3);
        double t22 = x.f(1,3)*x.f(3,2);
        double t24 = x.f(0,0)*x.f(3,1);
        double t25 = x.f(1,2)*x.f(2,3);
        double t27 = x.f(1,3)*x.f(2,2);
        double t29 = x.f(1,0)*x.f(0,1);
        double t32 = x.f(1,0)*x.f(2,1);
        double t33 = x.f(0,2)*x.f(3,3);
        double t35 = x.f(0,3)*x.f(3,2);
        double t37 = x.f(1,0)*x.f(3,1);
        double t38 = x.f(0,2)*x.f(2,3);
        double t40 = x.f(0,3)*x.f(2,2);
        double t42 = t14*t15-t14*t17-t19*t20+t19*t22+
                           t24*t25-t24*t27-t29*t15+t29*t17+
                           t32*t33-t32*t35-t37*t38+t37*t40;
        double t43 = x.f(2,0)*x.f(0,1);
        double t46 = x.f(2,0)*x.f(1,1);
        double t49 = x.f(2,0)*x.f(3,1);
        double t50 = x.f(0,2)*x.f(1,3);
        double t52 = x.f(0,3)*x.f(1,2);
        double t54 = x.f(3,0)*x.f(0,1);
        double t57 = x.f(3,0)*x.f(1,1);
        double t60 = x.f(3,0)*x.f(2,1);
        double t63 = t43*t20-t43*t22-t46*t33+t46*t35+
                           t49*t50-t49*t52-t54*t25+t54*t27+
                           t57*t38-t57*t40-t60*t50+t60*t52;
        double t65 = 1/(t42+t63);
        double t71 = x.f(0,2)*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 t77 = x.f(0,3)*x.f(3,1);
        double t81 = x.f(0,1)*x.f(1,2);
        double t83 = x.f(0,1)*x.f(1,3);
        double t85 = x.f(0,2)*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 t103 = x.f(1,0)*x.f(2,3);
        double t105 = x.f(2,0)*x.f(1,2);
        double t107 = x.f(2,0)*x.f(1,3);
        double t109 = x.f(3,0)*x.f(1,2);
        double t111 = x.f(3,0)*x.f(1,3);
        double t115 = x.f(0,0)*x.f(2,2);
        double t117 = x.f(0,0)*x.f(2,3);
        double t119 = x.f(2,0)*x.f(0,2);
        double t121 = x.f(2,0)*x.f(0,3);
}

template void X<double>::g (X<double>&);