4 #define SET_GAMMAI_DER_ZERO(I) \ 5 d##I##gammai11_a(i, j, k) = 0.0; \ 6 d##I##gammai22_a(i, j, k) = 0.0; \ 7 d##I##gammai33_a(i, j, k) = 0.0; \ 8 d##I##gammai12_a(i, j, k) = 0.0; \ 9 d##I##gammai13_a(i, j, k) = 0.0; \ 10 d##I##gammai23_a(i, j, k) = 0.0 13 #define SET_GAMMAI_DER(I) \ 14 d##I##gammai11_a(i, j, k) = -4.0*derivative(i, j, k, I, DIFFphi_a)*gammai11 \ 15 + std::exp(-4.0*DIFFphi_a(i, j, k))*(derivative(i, j, k, I, DIFFgamma22_a) + derivative(i, j, k, I, DIFFgamma33_a) - 2.0*DIFFgamma23_a(i, j, k)*derivative(i, j, k, I, DIFFgamma23_a) + derivative(i, j, k, I, DIFFgamma22_a)*DIFFgamma33_a(i, j, k) + DIFFgamma22_a(i, j, k)*derivative(i, j, k, I, DIFFgamma33_a)); \ 17 d##I##gammai22_a(i, j, k) = -4.0*derivative(i, j, k, I, DIFFphi_a)*gammai22 \ 18 + std::exp(-4.0*DIFFphi_a(i, j, k))*(derivative(i, j, k, I, DIFFgamma11_a) + derivative(i, j, k, I, DIFFgamma33_a) - 2.0*DIFFgamma13_a(i, j, k)*derivative(i, j, k, I, DIFFgamma13_a) + derivative(i, j, k, I, DIFFgamma11_a)*DIFFgamma33_a(i, j, k) + DIFFgamma11_a(i, j, k)*derivative(i, j, k, I, DIFFgamma33_a)); \ 20 d##I##gammai33_a(i, j, k) = -4.0*derivative(i, j, k, I, DIFFphi_a)*gammai33 \ 21 + std::exp(-4.0*DIFFphi_a(i, j, k))*(derivative(i, j, k, I, DIFFgamma11_a) + derivative(i, j, k, I, DIFFgamma22_a) - 2.0*DIFFgamma12_a(i,j,k)*derivative(i, j, k, I, DIFFgamma12_a) + derivative(i, j, k, I, DIFFgamma11_a)*DIFFgamma22_a(i, j, k) + DIFFgamma11_a(i, j, k)*derivative(i, j, k, I, DIFFgamma22_a)); \ 23 d##I##gammai12_a(i, j, k) = -4.0*derivative(i, j, k, I, DIFFphi_a)*gammai12 \ 24 + std::exp(-4.0*DIFFphi_a(i, j, k))*(derivative(i, j, k, I, DIFFgamma13_a)*DIFFgamma23_a(i, j, k) + DIFFgamma13_a(i, j, k)*derivative(i, j, k, I, DIFFgamma23_a) - derivative(i, j, k, I, DIFFgamma12_a)*(1.0 + DIFFgamma33_a(i, j, k)) - DIFFgamma12_a(i, j, k)*derivative(i, j, k, I, DIFFgamma33_a)); \ 26 d##I##gammai13_a(i, j, k) = -4.0*derivative(i, j, k, I, DIFFphi_a)*gammai13 \ 27 + std::exp(-4.0*DIFFphi_a(i, j, k))*(derivative(i, j, k, I, DIFFgamma12_a)*DIFFgamma23_a(i, j, k) + DIFFgamma12_a(i, j, k)*derivative(i, j, k, I, DIFFgamma23_a) - derivative(i, j, k, I, DIFFgamma13_a)*(1.0 + DIFFgamma22_a(i, j, k)) - DIFFgamma13_a(i, j, k)*derivative(i, j, k, I, DIFFgamma22_a)); \ 29 d##I##gammai23_a(i, j, k) = -4.0*derivative(i, j, k, I, DIFFphi_a)*gammai23 \ 30 + std::exp(-4.0*DIFFphi_a(i, j, k))*(derivative(i, j, k, I, DIFFgamma12_a)*DIFFgamma13_a(i, j, k) + DIFFgamma12_a(i, j, k)*derivative(i, j, k, I, DIFFgamma13_a) - derivative(i, j, k, I, DIFFgamma23_a)*(1.0 + DIFFgamma11_a(i, j, k)) - DIFFgamma23_a(i, j, k)*derivative(i, j, k, I, DIFFgamma11_a)) 
