BF_C_H_Q4_3D.h

// ***********************************************************************
// Q4 element, conforming, 3D
// ***********************************************************************

static void C_H_Q4_3D_Funct(double xi, double eta, double zeta,
                          double *values)
{
  double t1, t2, t3, t4, t5, t6, t7, t8, t9, t10, t11, t12, t13, t14, t15;
  double t16, t17, t18, t19, t20, t21, t22, t23, t24, t25, t27, t28, t29;
  double t30, t31, t32, t36, t37, t39, t40, t42, t43, t45, t46, t47, t48;
  double t49, t50, t52, t54, t56, t58, t62, t63, t65, t67, t69, t71, t73;
  double t74, t76, t78, t80, t82, t84, t85, t87, t89, t91, t93, t95, t96;
  double t97, t98, t99, t127, t153, t179;

  t1 = xi/6.0;
  t2 = xi*xi;
  t3 = t2/6.0;
  t4 = t2*xi;
  t5 = 2.0/3.0*t4;
  t6 = t2*t2;
  t7 = 2.0/3.0*t6;
  t8 = t1-t3-t5+t7;
  t9 = eta/6.0;
  t10 = eta*eta;
  t11 = t10/6.0;
  t12 = t10*eta;
  t13 = 2.0/3.0*t12;
  t14 = t10*t10;
  t15 = 2.0/3.0*t14;
  t16 = t9-t11-t13+t15;
  t17 = t8*t16;
  t18 = zeta/6.0;
  t19 = zeta*zeta;
  t20 = t19/6.0;
  t21 = t19*zeta;
  t22 = 2.0/3.0*t21;
  t23 = t19*t19;
  t24 = 2.0/3.0*t23;
  t25 = t18-t20-t22+t24;
  t27 = 4.0/3.0*xi;
  t28 = 8.0/3.0*t2;
  t29 = 4.0/3.0*t4;
  t30 = 8.0/3.0*t6;
  t31 = -t27+t28+t29-t30;
  t32 = t31*t16;
  t36 = 1.0-5.0*t2+4.0*t6;
  t37 = t36*t16;
  t39 = t27+t28-t29-t30;
  t40 = t39*t16;
  t42 = -t1-t3+t5+t7;
  t43 = t42*t16;
  t45 = 4.0/3.0*eta;
  t46 = 8.0/3.0*t10;
  t47 = 4.0/3.0*t12;
  t48 = 8.0/3.0*t14;
  t49 = -t45+t46+t47-t48;
  t50 = t8*t49;
  t52 = t31*t49;
  t54 = t36*t49;
  t56 = t39*t49;
  t58 = t42*t49;
  t62 = 1.0-5.0*t10+4.0*t14;
  t63 = t8*t62;
  t65 = t31*t62;
  t67 = t36*t62;
  t69 = t39*t62;
  t71 = t42*t62;
  t73 = t45+t46-t47-t48;
  t74 = t8*t73;
  t76 = t31*t73;
  t78 = t36*t73;
  t80 = t39*t73;
  t82 = t42*t73;
  t84 = -t9-t11+t13+t15;
  t85 = t8*t84;
  t87 = t31*t84;
  t89 = t36*t84;
  t91 = t39*t84;
  t93 = t42*t84;
  t95 = 4.0/3.0*zeta;
  t96 = 8.0/3.0*t19;
  t97 = 4.0/3.0*t21;
  t98 = 8.0/3.0*t23;
  t99 = -t95+t96+t97-t98;
  t127 = 1.0-5.0*t19+4.0*t23;
  t153 = t95+t96-t97-t98;
  t179 = -t18-t20+t22+t24;

  values[0] = t17*t25;
  values[1] = t32*t25;
  values[2] = t37*t25;
  values[3] = t40*t25;
  values[4] = t43*t25;
  values[5] = t50*t25;
  values[6] = t52*t25;
  values[7] = t54*t25;
  values[8] = t56*t25;
  values[9] = t58*t25;
  values[10] = t63*t25;
  values[11] = t65*t25;
  values[12] = t67*t25;
  values[13] = t69*t25;
  values[14] = t71*t25;
  values[15] = t74*t25;
  values[16] = t76*t25;
  values[17] = t78*t25;
  values[18] = t80*t25;
  values[19] = t82*t25;
  values[20] = t85*t25;
  values[21] = t87*t25;
  values[22] = t89*t25;
  values[23] = t91*t25;
  values[24] = t93*t25;
  values[25] = t17*t99;
  values[26] = t32*t99;
  values[27] = t37*t99;
  values[28] = t40*t99;
  values[29] = t43*t99;
  values[30] = t50*t99;
  values[31] = t52*t99;
  values[32] = t54*t99;
  values[33] = t56*t99;
  values[34] = t58*t99;
  values[35] = t63*t99;
  values[36] = t65*t99;
  values[37] = t67*t99;
  values[38] = t69*t99;
  values[39] = t71*t99;
  values[40] = t74*t99;
  values[41] = t76*t99;
  values[42] = t78*t99;
  values[43] = t80*t99;
  values[44] = t82*t99;
  values[45] = t85*t99;
  values[46] = t87*t99;
  values[47] = t89*t99;
  values[48] = t91*t99;
  values[49] = t93*t99;
  values[50] = t17*t127;
  values[51] = t32*t127;
  values[52] = t37*t127;
  values[53] = t40*t127;
  values[54] = t43*t127;
  values[55] = t50*t127;
  values[56] = t52*t127;
  values[57] = t54*t127;
  values[58] = t56*t127;
  values[59] = t58*t127;
  values[60] = t63*t127;
  values[61] = t65*t127;
  values[62] = t67*t127;
  values[63] = t69*t127;
  values[64] = t71*t127;
  values[65] = t74*t127;
  values[66] = t76*t127;
  values[67] = t78*t127;
  values[68] = t80*t127;
  values[69] = t82*t127;
  values[70] = t85*t127;
  values[71] = t87*t127;
  values[72] = t89*t127;
  values[73] = t91*t127;
  values[74] = t93*t127;
  values[75] = t17*t153;
  values[76] = t32*t153;
  values[77] = t37*t153;
  values[78] = t40*t153;
  values[79] = t43*t153;
  values[80] = t50*t153;
  values[81] = t52*t153;
  values[82] = t54*t153;
  values[83] = t56*t153;
  values[84] = t58*t153;
  values[85] = t63*t153;
  values[86] = t65*t153;
  values[87] = t67*t153;
  values[88] = t69*t153;
  values[89] = t71*t153;
  values[90] = t74*t153;
  values[91] = t76*t153;
  values[92] = t78*t153;
  values[93] = t80*t153;
  values[94] = t82*t153;
  values[95] = t85*t153;
  values[96] = t87*t153;
  values[97] = t89*t153;
  values[98] = t91*t153;
  values[99] = t93*t153;
  values[100] = t17*t179;
  values[101] = t32*t179;
  values[102] = t37*t179;
  values[103] = t40*t179;
  values[104] = t43*t179;
  values[105] = t50*t179;
  values[106] = t52*t179;
  values[107] = t54*t179;
  values[108] = t56*t179;
  values[109] = t58*t179;
  values[110] = t63*t179;
  values[111] = t65*t179;
  values[112] = t67*t179;
  values[113] = t69*t179;
  values[114] = t71*t179;
  values[115] = t74*t179;
  values[116] = t76*t179;
  values[117] = t78*t179;
  values[118] = t80*t179;
  values[119] = t82*t179;
  values[120] = t85*t179;
  values[121] = t87*t179;
  values[122] = t89*t179;
  values[123] = t91*t179;
  values[124] = t93*t179;
}

static void C_H_Q4_3D_DeriveXi(double xi, double eta, double zeta,
                             double *values)
{
  double t1, t2, t3, t4, t5, t6, t7, t8, t9, t10, t11, t12, t13, t14, t15;
  double t16, t17, t18, t19, t20, t21, t22, t23, t25, t26, t27, t28, t29;
  double t33, t34, t36, t37, t39, t40, t42, t43, t44, t45, t46, t47, t49;
  double t51, t53, t55, t59, t60, t62, t64, t66, t68, t70, t71, t73, t75;
  double t77, t79, t81, t82, t84, t86, t88, t90, t92, t93, t94, t95, t96;
  double t124, t150, t176;

  t1 = xi/3.0;
  t2 = xi*xi;
  t3 = 2.0*t2;
  t4 = t2*xi;
  t5 = 8.0/3.0*t4;
  t6 = 1.0/6.0-t1-t3+t5;
  t7 = eta/6.0;
  t8 = eta*eta;
  t9 = t8/6.0;
  t10 = t8*eta;
  t11 = 2.0/3.0*t10;
  t12 = t8*t8;
  t13 = 2.0/3.0*t12;
  t14 = t7-t9-t11+t13;
  t15 = t6*t14;
  t16 = zeta/6.0;
  t17 = zeta*zeta;
  t18 = t17/6.0;
  t19 = t17*zeta;
  t20 = 2.0/3.0*t19;
  t21 = t17*t17;
  t22 = 2.0/3.0*t21;
  t23 = t16-t18-t20+t22;
  t25 = 16.0/3.0*xi;
  t26 = 4.0*t2;
  t27 = 32.0/3.0*t4;
  t28 = -4.0/3.0+t25+t26-t27;
  t29 = t28*t14;
  t33 = -10.0*xi+16.0*t4;
  t34 = t33*t14;
  t36 = 4.0/3.0+t25-t26-t27;
  t37 = t36*t14;
  t39 = -1.0/6.0-t1+t3+t5;
  t40 = t39*t14;
  t42 = 4.0/3.0*eta;
  t43 = 8.0/3.0*t8;
  t44 = 4.0/3.0*t10;
  t45 = 8.0/3.0*t12;
  t46 = -t42+t43+t44-t45;
  t47 = t6*t46;
  t49 = t28*t46;
  t51 = t33*t46;
  t53 = t36*t46;
  t55 = t39*t46;
  t59 = 1.0-5.0*t8+4.0*t12;
  t60 = t6*t59;
  t62 = t28*t59;
  t64 = t33*t59;
  t66 = t36*t59;
  t68 = t39*t59;
  t70 = t42+t43-t44-t45;
  t71 = t6*t70;
  t73 = t28*t70;
  t75 = t33*t70;
  t77 = t36*t70;
  t79 = t39*t70;
  t81 = -t7-t9+t11+t13;
  t82 = t6*t81;
  t84 = t28*t81;
  t86 = t33*t81;
  t88 = t36*t81;
  t90 = t39*t81;
  t92 = 4.0/3.0*zeta;
  t93 = 8.0/3.0*t17;
  t94 = 4.0/3.0*t19;
  t95 = 8.0/3.0*t21;
  t96 = -t92+t93+t94-t95;
  t124 = 1.0-5.0*t17+4.0*t21;
  t150 = t92+t93-t94-t95;
  t176 = -t16-t18+t20+t22;

  values[0] = t15*t23;
  values[1] = t29*t23;
  values[2] = t34*t23;
  values[3] = t37*t23;
  values[4] = t40*t23;
  values[5] = t47*t23;
  values[6] = t49*t23;
  values[7] = t51*t23;
  values[8] = t53*t23;
  values[9] = t55*t23;
  values[10] = t60*t23;
  values[11] = t62*t23;
  values[12] = t64*t23;
  values[13] = t66*t23;
  values[14] = t68*t23;
  values[15] = t71*t23;
  values[16] = t73*t23;
  values[17] = t75*t23;
  values[18] = t77*t23;
  values[19] = t79*t23;
  values[20] = t82*t23;
  values[21] = t84*t23;
  values[22] = t86*t23;
  values[23] = t88*t23;
  values[24] = t90*t23;
  values[25] = t15*t96;
  values[26] = t29*t96;
  values[27] = t34*t96;
  values[28] = t37*t96;
  values[29] = t40*t96;
  values[30] = t47*t96;
  values[31] = t49*t96;
  values[32] = t51*t96;
  values[33] = t53*t96;
  values[34] = t55*t96;
  values[35] = t60*t96;
  values[36] = t62*t96;
  values[37] = t64*t96;
  values[38] = t66*t96;
  values[39] = t68*t96;
  values[40] = t71*t96;
  values[41] = t73*t96;
  values[42] = t75*t96;
  values[43] = t77*t96;
  values[44] = t79*t96;
  values[45] = t82*t96;
  values[46] = t84*t96;
  values[47] = t86*t96;
  values[48] = t88*t96;
  values[49] = t90*t96;
  values[50] = t15*t124;
  values[51] = t29*t124;
  values[52] = t34*t124;
  values[53] = t37*t124;
  values[54] = t40*t124;
  values[55] = t47*t124;
  values[56] = t49*t124;
  values[57] = t51*t124;
  values[58] = t53*t124;
  values[59] = t55*t124;
  values[60] = t60*t124;
  values[61] = t62*t124;
  values[62] = t64*t124;
  values[63] = t66*t124;
  values[64] = t68*t124;
  values[65] = t71*t124;
  values[66] = t73*t124;
  values[67] = t75*t124;
  values[68] = t77*t124;
  values[69] = t79*t124;
  values[70] = t82*t124;
  values[71] = t84*t124;
  values[72] = t86*t124;
  values[73] = t88*t124;
  values[74] = t90*t124;
  values[75] = t15*t150;
  values[76] = t29*t150;
  values[77] = t34*t150;
  values[78] = t37*t150;
  values[79] = t40*t150;
  values[80] = t47*t150;
  values[81] = t49*t150;
  values[82] = t51*t150;
  values[83] = t53*t150;
  values[84] = t55*t150;
  values[85] = t60*t150;
  values[86] = t62*t150;
  values[87] = t64*t150;
  values[88] = t66*t150;
  values[89] = t68*t150;
  values[90] = t71*t150;
  values[91] = t73*t150;
  values[92] = t75*t150;
  values[93] = t77*t150;
  values[94] = t79*t150;
  values[95] = t82*t150;
  values[96] = t84*t150;
  values[97] = t86*t150;
  values[98] = t88*t150;
  values[99] = t90*t150;
  values[100] = t15*t176;
  values[101] = t29*t176;
  values[102] = t34*t176;
  values[103] = t37*t176;
  values[104] = t40*t176;
  values[105] = t47*t176;
  values[106] = t49*t176;
  values[107] = t51*t176;
  values[108] = t53*t176;
  values[109] = t55*t176;
  values[110] = t60*t176;
  values[111] = t62*t176;
  values[112] = t64*t176;
  values[113] = t66*t176;
  values[114] = t68*t176;
  values[115] = t71*t176;
  values[116] = t73*t176;
  values[117] = t75*t176;
  values[118] = t77*t176;
  values[119] = t79*t176;
  values[120] = t82*t176;
  values[121] = t84*t176;
  values[122] = t86*t176;
  values[123] = t88*t176;
  values[124] = t90*t176;
}

static void C_H_Q4_3D_DeriveEta(double xi, double eta, double zeta,
                             double *values)
{
  double t1, t2, t3, t4, t5, t6, t7, t8, t9, t10, t11, t12, t13, t14, t15;
  double t16, t17, t18, t19, t20, t21, t22, t23, t25, t26, t27, t28, t29;
  double t30, t34, t35, t37, t38, t40, t41, t43, t44, t45, t46, t47, t49;
  double t51, t53, t55, t59, t60, t62, t64, t66, t68, t70, t71, t73, t75;
  double t77, t79, t81, t82, t84, t86, t88, t90, t92, t93, t94, t95, t96;
  double t124, t150, t176;

  t1 = xi/6.0;
  t2 = xi*xi;
  t3 = t2/6.0;
  t4 = t2*xi;
  t5 = 2.0/3.0*t4;
  t6 = t2*t2;
  t7 = 2.0/3.0*t6;
  t8 = t1-t3-t5+t7;
  t9 = eta/3.0;
  t10 = eta*eta;
  t11 = 2.0*t10;
  t12 = t10*eta;
  t13 = 8.0/3.0*t12;
  t14 = 1.0/6.0-t9-t11+t13;
  t15 = t8*t14;
  t16 = zeta/6.0;
  t17 = zeta*zeta;
  t18 = t17/6.0;
  t19 = t17*zeta;
  t20 = 2.0/3.0*t19;
  t21 = t17*t17;
  t22 = 2.0/3.0*t21;
  t23 = t16-t18-t20+t22;
  t25 = 4.0/3.0*xi;
  t26 = 8.0/3.0*t2;
  t27 = 4.0/3.0*t4;
  t28 = 8.0/3.0*t6;
  t29 = -t25+t26+t27-t28;
  t30 = t29*t14;
  t34 = 1.0-5.0*t2+4.0*t6;
  t35 = t34*t14;
  t37 = t25+t26-t27-t28;
  t38 = t37*t14;
  t40 = -t1-t3+t5+t7;
  t41 = t40*t14;
  t43 = 16.0/3.0*eta;
  t44 = 4.0*t10;
  t45 = 32.0/3.0*t12;
  t46 = -4.0/3.0+t43+t44-t45;
  t47 = t8*t46;
  t49 = t29*t46;
  t51 = t34*t46;
  t53 = t37*t46;
  t55 = t40*t46;
  t59 = -10.0*eta+16.0*t12;
  t60 = t8*t59;
  t62 = t29*t59;
  t64 = t34*t59;
  t66 = t37*t59;
  t68 = t40*t59;
  t70 = 4.0/3.0+t43-t44-t45;
  t71 = t8*t70;
  t73 = t29*t70;
  t75 = t34*t70;
  t77 = t37*t70;
  t79 = t40*t70;
  t81 = -1.0/6.0-t9+t11+t13;
  t82 = t8*t81;
  t84 = t29*t81;
  t86 = t34*t81;
  t88 = t37*t81;
  t90 = t40*t81;
  t92 = 4.0/3.0*zeta;
  t93 = 8.0/3.0*t17;
  t94 = 4.0/3.0*t19;
  t95 = 8.0/3.0*t21;
  t96 = -t92+t93+t94-t95;
  t124 = 1.0-5.0*t17+4.0*t21;
  t150 = t92+t93-t94-t95;
  t176 = -t16-t18+t20+t22;

  values[0] = t15*t23;
  values[1] = t30*t23;
  values[2] = t35*t23;
  values[3] = t38*t23;
  values[4] = t41*t23;
  values[5] = t47*t23;
  values[6] = t49*t23;
  values[7] = t51*t23;
  values[8] = t53*t23;
  values[9] = t55*t23;
  values[10] = t60*t23;
  values[11] = t62*t23;
  values[12] = t64*t23;
  values[13] = t66*t23;
  values[14] = t68*t23;
  values[15] = t71*t23;
  values[16] = t73*t23;
  values[17] = t75*t23;
  values[18] = t77*t23;
  values[19] = t79*t23;
  values[20] = t82*t23;
  values[21] = t84*t23;
  values[22] = t86*t23;
  values[23] = t88*t23;
  values[24] = t90*t23;
  values[25] = t15*t96;
  values[26] = t30*t96;
  values[27] = t35*t96;
  values[28] = t38*t96;
  values[29] = t41*t96;
  values[30] = t47*t96;
  values[31] = t49*t96;
  values[32] = t51*t96;
  values[33] = t53*t96;
  values[34] = t55*t96;
  values[35] = t60*t96;
  values[36] = t62*t96;
  values[37] = t64*t96;
  values[38] = t66*t96;
  values[39] = t68*t96;
  values[40] = t71*t96;
  values[41] = t73*t96;
  values[42] = t75*t96;
  values[43] = t77*t96;
  values[44] = t79*t96;
  values[45] = t82*t96;
  values[46] = t84*t96;
  values[47] = t86*t96;
  values[48] = t88*t96;
  values[49] = t90*t96;
  values[50] = t15*t124;
  values[51] = t30*t124;
  values[52] = t35*t124;
  values[53] = t38*t124;
  values[54] = t41*t124;
  values[55] = t47*t124;
  values[56] = t49*t124;
  values[57] = t51*t124;
  values[58] = t53*t124;
  values[59] = t55*t124;
  values[60] = t60*t124;
  values[61] = t62*t124;
  values[62] = t64*t124;
  values[63] = t66*t124;
  values[64] = t68*t124;
  values[65] = t71*t124;
  values[66] = t73*t124;
  values[67] = t75*t124;
  values[68] = t77*t124;
  values[69] = t79*t124;
  values[70] = t82*t124;
  values[71] = t84*t124;
  values[72] = t86*t124;
  values[73] = t88*t124;
  values[74] = t90*t124;
  values[75] = t15*t150;
  values[76] = t30*t150;
  values[77] = t35*t150;
  values[78] = t38*t150;
  values[79] = t41*t150;
  values[80] = t47*t150;
  values[81] = t49*t150;
  values[82] = t51*t150;
  values[83] = t53*t150;
  values[84] = t55*t150;
  values[85] = t60*t150;
  values[86] = t62*t150;
  values[87] = t64*t150;
  values[88] = t66*t150;
  values[89] = t68*t150;
  values[90] = t71*t150;
  values[91] = t73*t150;
  values[92] = t75*t150;
  values[93] = t77*t150;
  values[94] = t79*t150;
  values[95] = t82*t150;
  values[96] = t84*t150;
  values[97] = t86*t150;
  values[98] = t88*t150;
  values[99] = t90*t150;
  values[100] = t15*t176;
  values[101] = t30*t176;
  values[102] = t35*t176;
  values[103] = t38*t176;
  values[104] = t41*t176;
  values[105] = t47*t176;
  values[106] = t49*t176;
  values[107] = t51*t176;
  values[108] = t53*t176;
  values[109] = t55*t176;
  values[110] = t60*t176;
  values[111] = t62*t176;
  values[112] = t64*t176;
  values[113] = t66*t176;
  values[114] = t68*t176;
  values[115] = t71*t176;
  values[116] = t73*t176;
  values[117] = t75*t176;
  values[118] = t77*t176;
  values[119] = t79*t176;
  values[120] = t82*t176;
  values[121] = t84*t176;
  values[122] = t86*t176;
  values[123] = t88*t176;
  values[124] = t90*t176;
}

static void C_H_Q4_3D_DeriveZeta(double xi, double eta, double zeta,
                             double *values)
{
  double t1, t2, t3, t4, t5, t6, t7, t8, t9, t10, t11, t12, t13, t14, t15;
  double t16, t17, t18, t19, t20, t21, t22, t23, t25, t26, t27, t28, t29;
  double t30, t34, t35, t37, t38, t40, t41, t43, t44, t45, t46, t47, t48;
  double t50, t52, t54, t56, t60, t61, t63, t65, t67, t69, t71, t72, t74;
  double t76, t78, t80, t82, t83, t85, t87, t89, t91, t93, t94, t95, t96;
  double t124, t150, t176;

  t1 = xi/6.0;
  t2 = xi*xi;
  t3 = t2/6.0;
  t4 = t2*xi;
  t5 = 2.0/3.0*t4;
  t6 = t2*t2;
  t7 = 2.0/3.0*t6;
  t8 = t1-t3-t5+t7;
  t9 = eta/6.0;
  t10 = eta*eta;
  t11 = t10/6.0;
  t12 = t10*eta;
  t13 = 2.0/3.0*t12;
  t14 = t10*t10;
  t15 = 2.0/3.0*t14;
  t16 = t9-t11-t13+t15;
  t17 = t8*t16;
  t18 = zeta/3.0;
  t19 = zeta*zeta;
  t20 = 2.0*t19;
  t21 = t19*zeta;
  t22 = 8.0/3.0*t21;
  t23 = 1.0/6.0-t18-t20+t22;
  t25 = 4.0/3.0*xi;
  t26 = 8.0/3.0*t2;
  t27 = 4.0/3.0*t4;
  t28 = 8.0/3.0*t6;
  t29 = -t25+t26+t27-t28;
  t30 = t29*t16;
  t34 = 1.0-5.0*t2+4.0*t6;
  t35 = t34*t16;
  t37 = t25+t26-t27-t28;
  t38 = t37*t16;
  t40 = -t1-t3+t5+t7;
  t41 = t40*t16;
  t43 = 4.0/3.0*eta;
  t44 = 8.0/3.0*t10;
  t45 = 4.0/3.0*t12;
  t46 = 8.0/3.0*t14;
  t47 = -t43+t44+t45-t46;
  t48 = t8*t47;
  t50 = t29*t47;
  t52 = t34*t47;
  t54 = t37*t47;
  t56 = t40*t47;
  t60 = 1.0-5.0*t10+4.0*t14;
  t61 = t8*t60;
  t63 = t29*t60;
  t65 = t34*t60;
  t67 = t37*t60;
  t69 = t40*t60;
  t71 = t43+t44-t45-t46;
  t72 = t8*t71;
  t74 = t29*t71;
  t76 = t34*t71;
  t78 = t37*t71;
  t80 = t40*t71;
  t82 = -t9-t11+t13+t15;
  t83 = t8*t82;
  t85 = t29*t82;
  t87 = t34*t82;
  t89 = t37*t82;
  t91 = t40*t82;
  t93 = 16.0/3.0*zeta;
  t94 = 4.0*t19;
  t95 = 32.0/3.0*t21;
  t96 = -4.0/3.0+t93+t94-t95;
  t124 = -10.0*zeta+16.0*t21;
  t150 = 4.0/3.0+t93-t94-t95;
  t176 = -1.0/6.0-t18+t20+t22;

  values[0] = t17*t23;
  values[1] = t30*t23;
  values[2] = t35*t23;
  values[3] = t38*t23;
  values[4] = t41*t23;
  values[5] = t48*t23;
  values[6] = t50*t23;
  values[7] = t52*t23;
  values[8] = t54*t23;
  values[9] = t56*t23;
  values[10] = t61*t23;
  values[11] = t63*t23;
  values[12] = t65*t23;
  values[13] = t67*t23;
  values[14] = t69*t23;
  values[15] = t72*t23;
  values[16] = t74*t23;
  values[17] = t76*t23;
  values[18] = t78*t23;
  values[19] = t80*t23;
  values[20] = t83*t23;
  values[21] = t85*t23;
  values[22] = t87*t23;
  values[23] = t89*t23;
  values[24] = t91*t23;
  values[25] = t17*t96;
  values[26] = t30*t96;
  values[27] = t35*t96;
  values[28] = t38*t96;
  values[29] = t41*t96;
  values[30] = t48*t96;
  values[31] = t50*t96;
  values[32] = t52*t96;
  values[33] = t54*t96;
  values[34] = t56*t96;
  values[35] = t61*t96;
  values[36] = t63*t96;
  values[37] = t65*t96;
  values[38] = t67*t96;
  values[39] = t69*t96;
  values[40] = t72*t96;
  values[41] = t74*t96;
  values[42] = t76*t96;
  values[43] = t78*t96;
  values[44] = t80*t96;
  values[45] = t83*t96;
  values[46] = t85*t96;
  values[47] = t87*t96;
  values[48] = t89*t96;
  values[49] = t91*t96;
  values[50] = t17*t124;
  values[51] = t30*t124;
  values[52] = t35*t124;
  values[53] = t38*t124;
  values[54] = t41*t124;
  values[55] = t48*t124;
  values[56] = t50*t124;
  values[57] = t52*t124;
  values[58] = t54*t124;
  values[59] = t56*t124;
  values[60] = t61*t124;
  values[61] = t63*t124;
  values[62] = t65*t124;
  values[63] = t67*t124;
  values[64] = t69*t124;
  values[65] = t72*t124;
  values[66] = t74*t124;
  values[67] = t76*t124;
  values[68] = t78*t124;
  values[69] = t80*t124;
  values[70] = t83*t124;
  values[71] = t85*t124;
  values[72] = t87*t124;
  values[73] = t89*t124;
  values[74] = t91*t124;
  values[75] = t17*t150;
  values[76] = t30*t150;
  values[77] = t35*t150;
  values[78] = t38*t150;
  values[79] = t41*t150;
  values[80] = t48*t150;
  values[81] = t50*t150;
  values[82] = t52*t150;
  values[83] = t54*t150;
  values[84] = t56*t150;
  values[85] = t61*t150;
  values[86] = t63*t150;
  values[87] = t65*t150;
  values[88] = t67*t150;
  values[89] = t69*t150;
  values[90] = t72*t150;
  values[91] = t74*t150;
  values[92] = t76*t150;
  values[93] = t78*t150;
  values[94] = t80*t150;
  values[95] = t83*t150;
  values[96] = t85*t150;
  values[97] = t87*t150;
  values[98] = t89*t150;
  values[99] = t91*t150;
  values[100] = t17*t176;
  values[101] = t30*t176;
  values[102] = t35*t176;
  values[103] = t38*t176;
  values[104] = t41*t176;
  values[105] = t48*t176;
  values[106] = t50*t176;
  values[107] = t52*t176;
  values[108] = t54*t176;
  values[109] = t56*t176;
  values[110] = t61*t176;
  values[111] = t63*t176;
  values[112] = t65*t176;
  values[113] = t67*t176;
  values[114] = t69*t176;
  values[115] = t72*t176;
  values[116] = t74*t176;
  values[117] = t76*t176;
  values[118] = t78*t176;
  values[119] = t80*t176;
  values[120] = t83*t176;
  values[121] = t85*t176;
  values[122] = t87*t176;
  values[123] = t89*t176;
  values[124] = t91*t176;
}

static void C_H_Q4_3D_DeriveXiXi(double xi, double eta, double zeta,
                             double *values)
{
  double t1, t2, t3, t4, t5, t6, t7, t8, t9, t10, t11, t12, t13, t14, t15;
  double t16, t17, t18, t19, t20, t21, t23, t24, t25, t26, t29, t30, t32;
  double t33, t35, t36, t38, t39, t40, t41, t42, t43, t45, t47, t49, t51;
  double t55, t56, t58, t60, t62, t64, t66, t67, t69, t71, t73, t75, t77;
  double t78, t80, t82, t84, t86, t88, t89, t90, t91, t92, t120, t146, t172;

  t1 = 4.0*xi;
  t2 = xi*xi;
  t3 = 8.0*t2;
  t4 = -1.0/3.0-t1+t3;
  t5 = eta/6.0;
  t6 = eta*eta;
  t7 = t6/6.0;
  t8 = t6*eta;
  t9 = 2.0/3.0*t8;
  t10 = t6*t6;
  t11 = 2.0/3.0*t10;
  t12 = t5-t7-t9+t11;
  t13 = t4*t12;
  t14 = zeta/6.0;
  t15 = zeta*zeta;
  t16 = t15/6.0;
  t17 = t15*zeta;
  t18 = 2.0/3.0*t17;
  t19 = t15*t15;
  t20 = 2.0/3.0*t19;
  t21 = t14-t16-t18+t20;
  t23 = 8.0*xi;
  t24 = 32.0*t2;
  t25 = 16.0/3.0+t23-t24;
  t26 = t25*t12;
  t29 = -10.0+48.0*t2;
  t30 = t29*t12;
  t32 = 16.0/3.0-t23-t24;
  t33 = t32*t12;
  t35 = -1.0/3.0+t1+t3;
  t36 = t35*t12;
  t38 = 4.0/3.0*eta;
  t39 = 8.0/3.0*t6;
  t40 = 4.0/3.0*t8;
  t41 = 8.0/3.0*t10;
  t42 = -t38+t39+t40-t41;
  t43 = t4*t42;
  t45 = t25*t42;
  t47 = t29*t42;
  t49 = t32*t42;
  t51 = t35*t42;
  t55 = 1.0-5.0*t6+4.0*t10;
  t56 = t4*t55;
  t58 = t25*t55;
  t60 = t29*t55;
  t62 = t32*t55;
  t64 = t35*t55;
  t66 = t38+t39-t40-t41;
  t67 = t4*t66;
  t69 = t25*t66;
  t71 = t29*t66;
  t73 = t32*t66;
  t75 = t35*t66;
  t77 = -t5-t7+t9+t11;
  t78 = t4*t77;
  t80 = t25*t77;
  t82 = t29*t77;
  t84 = t32*t77;
  t86 = t35*t77;
  t88 = 4.0/3.0*zeta;
  t89 = 8.0/3.0*t15;
  t90 = 4.0/3.0*t17;
  t91 = 8.0/3.0*t19;
  t92 = -t88+t89+t90-t91;
  t120 = 1.0-5.0*t15+4.0*t19;
  t146 = t88+t89-t90-t91;
  t172 = -t14-t16+t18+t20;

  values[0] = t13*t21;
  values[1] = t26*t21;
  values[2] = t30*t21;
  values[3] = t33*t21;
  values[4] = t36*t21;
  values[5] = t43*t21;
  values[6] = t45*t21;
  values[7] = t47*t21;
  values[8] = t49*t21;
  values[9] = t51*t21;
  values[10] = t56*t21;
  values[11] = t58*t21;
  values[12] = t60*t21;
  values[13] = t62*t21;
  values[14] = t64*t21;
  values[15] = t67*t21;
  values[16] = t69*t21;
  values[17] = t71*t21;
  values[18] = t73*t21;
  values[19] = t75*t21;
  values[20] = t78*t21;
  values[21] = t80*t21;
  values[22] = t82*t21;
  values[23] = t84*t21;
  values[24] = t86*t21;
  values[25] = t13*t92;
  values[26] = t26*t92;
  values[27] = t30*t92;
  values[28] = t33*t92;
  values[29] = t36*t92;
  values[30] = t43*t92;
  values[31] = t45*t92;
  values[32] = t47*t92;
  values[33] = t49*t92;
  values[34] = t51*t92;
  values[35] = t56*t92;
  values[36] = t58*t92;
  values[37] = t60*t92;
  values[38] = t62*t92;
  values[39] = t64*t92;
  values[40] = t67*t92;
  values[41] = t69*t92;
  values[42] = t71*t92;
  values[43] = t73*t92;
  values[44] = t75*t92;
  values[45] = t78*t92;
  values[46] = t80*t92;
  values[47] = t82*t92;
  values[48] = t84*t92;
  values[49] = t86*t92;
  values[50] = t13*t120;
  values[51] = t26*t120;
  values[52] = t30*t120;
  values[53] = t33*t120;
  values[54] = t36*t120;
  values[55] = t43*t120;
  values[56] = t45*t120;
  values[57] = t47*t120;
  values[58] = t49*t120;
  values[59] = t51*t120;
  values[60] = t56*t120;
  values[61] = t58*t120;
  values[62] = t60*t120;
  values[63] = t62*t120;
  values[64] = t64*t120;
  values[65] = t67*t120;
  values[66] = t69*t120;
  values[67] = t71*t120;
  values[68] = t73*t120;
  values[69] = t75*t120;
  values[70] = t78*t120;
  values[71] = t80*t120;
  values[72] = t82*t120;
  values[73] = t84*t120;
  values[74] = t86*t120;
  values[75] = t13*t146;
  values[76] = t26*t146;
  values[77] = t30*t146;
  values[78] = t33*t146;
  values[79] = t36*t146;
  values[80] = t43*t146;
  values[81] = t45*t146;
  values[82] = t47*t146;
  values[83] = t49*t146;
  values[84] = t51*t146;
  values[85] = t56*t146;
  values[86] = t58*t146;
  values[87] = t60*t146;
  values[88] = t62*t146;
  values[89] = t64*t146;
  values[90] = t67*t146;
  values[91] = t69*t146;
  values[92] = t71*t146;
  values[93] = t73*t146;
  values[94] = t75*t146;
  values[95] = t78*t146;
  values[96] = t80*t146;
  values[97] = t82*t146;
  values[98] = t84*t146;
  values[99] = t86*t146;
  values[100] = t13*t172;
  values[101] = t26*t172;
  values[102] = t30*t172;
  values[103] = t33*t172;
  values[104] = t36*t172;
  values[105] = t43*t172;
  values[106] = t45*t172;
  values[107] = t47*t172;
  values[108] = t49*t172;
  values[109] = t51*t172;
  values[110] = t56*t172;
  values[111] = t58*t172;
  values[112] = t60*t172;
  values[113] = t62*t172;
  values[114] = t64*t172;
  values[115] = t67*t172;
  values[116] = t69*t172;
  values[117] = t71*t172;
  values[118] = t73*t172;
  values[119] = t75*t172;
  values[120] = t78*t172;
  values[121] = t80*t172;
  values[122] = t82*t172;
  values[123] = t84*t172;
  values[124] = t86*t172;
}

static void C_H_Q4_3D_DeriveXiEta(double xi, double eta, double zeta,
                             double *values)
{
  double t1, t2, t3, t4, t5, t6, t7, t8, t9, t10, t11, t12, t13, t14, t15;
  double t16, t17, t18, t19, t20, t21, t23, t24, t25, t26, t27, t31, t32;
  double t34, t35, t37, t38, t40, t41, t42, t43, t44, t46, t48, t50, t52;
  double t56, t57, t59, t61, t63, t65, t67, t68, t70, t72, t74, t76, t78;
  double t79, t81, t83, t85, t87, t89, t90, t91, t92, t93, t121, t147, t173;

  t1 = xi/3.0;
  t2 = xi*xi;
  t3 = 2.0*t2;
  t4 = t2*xi;
  t5 = 8.0/3.0*t4;
  t6 = 1.0/6.0-t1-t3+t5;
  t7 = eta/3.0;
  t8 = eta*eta;
  t9 = 2.0*t8;
  t10 = t8*eta;
  t11 = 8.0/3.0*t10;
  t12 = 1.0/6.0-t7-t9+t11;
  t13 = t6*t12;
  t14 = zeta/6.0;
  t15 = zeta*zeta;
  t16 = t15/6.0;
  t17 = t15*zeta;
  t18 = 2.0/3.0*t17;
  t19 = t15*t15;
  t20 = 2.0/3.0*t19;
  t21 = t14-t16-t18+t20;
  t23 = 16.0/3.0*xi;
  t24 = 4.0*t2;
  t25 = 32.0/3.0*t4;
  t26 = -4.0/3.0+t23+t24-t25;
  t27 = t26*t12;
  t31 = -10.0*xi+16.0*t4;
  t32 = t31*t12;
  t34 = 4.0/3.0+t23-t24-t25;
  t35 = t34*t12;
  t37 = -1.0/6.0-t1+t3+t5;
  t38 = t37*t12;
  t40 = 16.0/3.0*eta;
  t41 = 4.0*t8;
  t42 = 32.0/3.0*t10;
  t43 = -4.0/3.0+t40+t41-t42;
  t44 = t6*t43;
  t46 = t26*t43;
  t48 = t31*t43;
  t50 = t34*t43;
  t52 = t37*t43;
  t56 = -10.0*eta+16.0*t10;
  t57 = t6*t56;
  t59 = t26*t56;
  t61 = t31*t56;
  t63 = t34*t56;
  t65 = t37*t56;
  t67 = 4.0/3.0+t40-t41-t42;
  t68 = t6*t67;
  t70 = t26*t67;
  t72 = t31*t67;
  t74 = t34*t67;
  t76 = t37*t67;
  t78 = -1.0/6.0-t7+t9+t11;
  t79 = t6*t78;
  t81 = t26*t78;
  t83 = t31*t78;
  t85 = t34*t78;
  t87 = t37*t78;
  t89 = 4.0/3.0*zeta;
  t90 = 8.0/3.0*t15;
  t91 = 4.0/3.0*t17;
  t92 = 8.0/3.0*t19;
  t93 = -t89+t90+t91-t92;
  t121 = 1.0-5.0*t15+4.0*t19;
  t147 = t89+t90-t91-t92;
  t173 = -t14-t16+t18+t20;

  values[0] = t13*t21;
  values[1] = t27*t21;
  values[2] = t32*t21;
  values[3] = t35*t21;
  values[4] = t38*t21;
  values[5] = t44*t21;
  values[6] = t46*t21;
  values[7] = t48*t21;
  values[8] = t50*t21;
  values[9] = t52*t21;
  values[10] = t57*t21;
  values[11] = t59*t21;
  values[12] = t61*t21;
  values[13] = t63*t21;
  values[14] = t65*t21;
  values[15] = t68*t21;
  values[16] = t70*t21;
  values[17] = t72*t21;
  values[18] = t74*t21;
  values[19] = t76*t21;
  values[20] = t79*t21;
  values[21] = t81*t21;
  values[22] = t83*t21;
  values[23] = t85*t21;
  values[24] = t87*t21;
  values[25] = t13*t93;
  values[26] = t27*t93;
  values[27] = t32*t93;
  values[28] = t35*t93;
  values[29] = t38*t93;
  values[30] = t44*t93;
  values[31] = t46*t93;
  values[32] = t48*t93;
  values[33] = t50*t93;
  values[34] = t52*t93;
  values[35] = t57*t93;
  values[36] = t59*t93;
  values[37] = t61*t93;
  values[38] = t63*t93;
  values[39] = t65*t93;
  values[40] = t68*t93;
  values[41] = t70*t93;
  values[42] = t72*t93;
  values[43] = t74*t93;
  values[44] = t76*t93;
  values[45] = t79*t93;
  values[46] = t81*t93;
  values[47] = t83*t93;
  values[48] = t85*t93;
  values[49] = t87*t93;
  values[50] = t13*t121;
  values[51] = t27*t121;
  values[52] = t32*t121;
  values[53] = t35*t121;
  values[54] = t38*t121;
  values[55] = t44*t121;
  values[56] = t46*t121;
  values[57] = t48*t121;
  values[58] = t50*t121;
  values[59] = t52*t121;
  values[60] = t57*t121;
  values[61] = t59*t121;
  values[62] = t61*t121;
  values[63] = t63*t121;
  values[64] = t65*t121;
  values[65] = t68*t121;
  values[66] = t70*t121;
  values[67] = t72*t121;
  values[68] = t74*t121;
  values[69] = t76*t121;
  values[70] = t79*t121;
  values[71] = t81*t121;
  values[72] = t83*t121;
  values[73] = t85*t121;
  values[74] = t87*t121;
  values[75] = t13*t147;
  values[76] = t27*t147;
  values[77] = t32*t147;
  values[78] = t35*t147;
  values[79] = t38*t147;
  values[80] = t44*t147;
  values[81] = t46*t147;
  values[82] = t48*t147;
  values[83] = t50*t147;
  values[84] = t52*t147;
  values[85] = t57*t147;
  values[86] = t59*t147;
  values[87] = t61*t147;
  values[88] = t63*t147;
  values[89] = t65*t147;
  values[90] = t68*t147;
  values[91] = t70*t147;
  values[92] = t72*t147;
  values[93] = t74*t147;
  values[94] = t76*t147;
  values[95] = t79*t147;
  values[96] = t81*t147;
  values[97] = t83*t147;
  values[98] = t85*t147;
  values[99] = t87*t147;
  values[100] = t13*t173;
  values[101] = t27*t173;
  values[102] = t32*t173;
  values[103] = t35*t173;
  values[104] = t38*t173;
  values[105] = t44*t173;
  values[106] = t46*t173;
  values[107] = t48*t173;
  values[108] = t50*t173;
  values[109] = t52*t173;
  values[110] = t57*t173;
  values[111] = t59*t173;
  values[112] = t61*t173;
  values[113] = t63*t173;
  values[114] = t65*t173;
  values[115] = t68*t173;
  values[116] = t70*t173;
  values[117] = t72*t173;
  values[118] = t74*t173;
  values[119] = t76*t173;
  values[120] = t79*t173;
  values[121] = t81*t173;
  values[122] = t83*t173;
  values[123] = t85*t173;
  values[124] = t87*t173;
}

static void C_H_Q4_3D_DeriveXiZeta(double xi, double eta, double zeta,
                             double *values)
{
  double t1, t2, t3, t4, t5, t6, t7, t8, t9, t10, t11, t12, t13, t14, t15;
  double t16, t17, t18, t19, t20, t21, t23, t24, t25, t26, t27, t31, t32;
  double t34, t35, t37, t38, t40, t41, t42, t43, t44, t45, t47, t49, t51;
  double t53, t57, t58, t60, t62, t64, t66, t68, t69, t71, t73, t75, t77;
  double t79, t80, t82, t84, t86, t88, t90, t91, t92, t93, t121, t147, t173;

  t1 = xi/3.0;
  t2 = xi*xi;
  t3 = 2.0*t2;
  t4 = t2*xi;
  t5 = 8.0/3.0*t4;
  t6 = 1.0/6.0-t1-t3+t5;
  t7 = eta/6.0;
  t8 = eta*eta;
  t9 = t8/6.0;
  t10 = t8*eta;
  t11 = 2.0/3.0*t10;
  t12 = t8*t8;
  t13 = 2.0/3.0*t12;
  t14 = t7-t9-t11+t13;
  t15 = t6*t14;
  t16 = zeta/3.0;
  t17 = zeta*zeta;
  t18 = 2.0*t17;
  t19 = t17*zeta;
  t20 = 8.0/3.0*t19;
  t21 = 1.0/6.0-t16-t18+t20;
  t23 = 16.0/3.0*xi;
  t24 = 4.0*t2;
  t25 = 32.0/3.0*t4;
  t26 = -4.0/3.0+t23+t24-t25;
  t27 = t26*t14;
  t31 = -10.0*xi+16.0*t4;
  t32 = t31*t14;
  t34 = 4.0/3.0+t23-t24-t25;
  t35 = t34*t14;
  t37 = -1.0/6.0-t1+t3+t5;
  t38 = t37*t14;
  t40 = 4.0/3.0*eta;
  t41 = 8.0/3.0*t8;
  t42 = 4.0/3.0*t10;
  t43 = 8.0/3.0*t12;
  t44 = -t40+t41+t42-t43;
  t45 = t6*t44;
  t47 = t26*t44;
  t49 = t31*t44;
  t51 = t34*t44;
  t53 = t37*t44;
  t57 = 1.0-5.0*t8+4.0*t12;
  t58 = t6*t57;
  t60 = t26*t57;
  t62 = t31*t57;
  t64 = t34*t57;
  t66 = t37*t57;
  t68 = t40+t41-t42-t43;
  t69 = t6*t68;
  t71 = t26*t68;
  t73 = t31*t68;
  t75 = t34*t68;
  t77 = t37*t68;
  t79 = -t7-t9+t11+t13;
  t80 = t6*t79;
  t82 = t26*t79;
  t84 = t31*t79;
  t86 = t34*t79;
  t88 = t37*t79;
  t90 = 16.0/3.0*zeta;
  t91 = 4.0*t17;
  t92 = 32.0/3.0*t19;
  t93 = -4.0/3.0+t90+t91-t92;
  t121 = -10.0*zeta+16.0*t19;
  t147 = 4.0/3.0+t90-t91-t92;
  t173 = -1.0/6.0-t16+t18+t20;

  values[0] = t15*t21;
  values[1] = t27*t21;
  values[2] = t32*t21;
  values[3] = t35*t21;
  values[4] = t38*t21;
  values[5] = t45*t21;
  values[6] = t47*t21;
  values[7] = t49*t21;
  values[8] = t51*t21;
  values[9] = t53*t21;
  values[10] = t58*t21;
  values[11] = t60*t21;
  values[12] = t62*t21;
  values[13] = t64*t21;
  values[14] = t66*t21;
  values[15] = t69*t21;
  values[16] = t71*t21;
  values[17] = t73*t21;
  values[18] = t75*t21;
  values[19] = t77*t21;
  values[20] = t80*t21;
  values[21] = t82*t21;
  values[22] = t84*t21;
  values[23] = t86*t21;
  values[24] = t88*t21;
  values[25] = t15*t93;
  values[26] = t27*t93;
  values[27] = t32*t93;
  values[28] = t35*t93;
  values[29] = t38*t93;
  values[30] = t45*t93;
  values[31] = t47*t93;
  values[32] = t49*t93;
  values[33] = t51*t93;
  values[34] = t53*t93;
  values[35] = t58*t93;
  values[36] = t60*t93;
  values[37] = t62*t93;
  values[38] = t64*t93;
  values[39] = t66*t93;
  values[40] = t69*t93;
  values[41] = t71*t93;
  values[42] = t73*t93;
  values[43] = t75*t93;
  values[44] = t77*t93;
  values[45] = t80*t93;
  values[46] = t82*t93;
  values[47] = t84*t93;
  values[48] = t86*t93;
  values[49] = t88*t93;
  values[50] = t15*t121;
  values[51] = t27*t121;
  values[52] = t32*t121;
  values[53] = t35*t121;
  values[54] = t38*t121;
  values[55] = t45*t121;
  values[56] = t47*t121;
  values[57] = t49*t121;
  values[58] = t51*t121;
  values[59] = t53*t121;
  values[60] = t58*t121;
  values[61] = t60*t121;
  values[62] = t62*t121;
  values[63] = t64*t121;
  values[64] = t66*t121;
  values[65] = t69*t121;
  values[66] = t71*t121;
  values[67] = t73*t121;
  values[68] = t75*t121;
  values[69] = t77*t121;
  values[70] = t80*t121;
  values[71] = t82*t121;
  values[72] = t84*t121;
  values[73] = t86*t121;
  values[74] = t88*t121;
  values[75] = t15*t147;
  values[76] = t27*t147;
  values[77] = t32*t147;
  values[78] = t35*t147;
  values[79] = t38*t147;
  values[80] = t45*t147;
  values[81] = t47*t147;
  values[82] = t49*t147;
  values[83] = t51*t147;
  values[84] = t53*t147;
  values[85] = t58*t147;
  values[86] = t60*t147;
  values[87] = t62*t147;
  values[88] = t64*t147;
  values[89] = t66*t147;
  values[90] = t69*t147;
  values[91] = t71*t147;
  values[92] = t73*t147;
  values[93] = t75*t147;
  values[94] = t77*t147;
  values[95] = t80*t147;
  values[96] = t82*t147;
  values[97] = t84*t147;
  values[98] = t86*t147;
  values[99] = t88*t147;
  values[100] = t15*t173;
  values[101] = t27*t173;
  values[102] = t32*t173;
  values[103] = t35*t173;
  values[104] = t38*t173;
  values[105] = t45*t173;
  values[106] = t47*t173;
  values[107] = t49*t173;
  values[108] = t51*t173;
  values[109] = t53*t173;
  values[110] = t58*t173;
  values[111] = t60*t173;
  values[112] = t62*t173;
  values[113] = t64*t173;
  values[114] = t66*t173;
  values[115] = t69*t173;
  values[116] = t71*t173;
  values[117] = t73*t173;
  values[118] = t75*t173;
  values[119] = t77*t173;
  values[120] = t80*t173;
  values[121] = t82*t173;
  values[122] = t84*t173;
  values[123] = t86*t173;
  values[124] = t88*t173;
}

static void C_H_Q4_3D_DeriveEtaEta(double xi, double eta, double zeta,
                             double *values)
{
  double t1, t2, t3, t4, t5, t6, t7, t8, t9, t10, t11, t12, t13, t14, t15;
  double t16, t17, t18, t19, t20, t21, t23, t24, t25, t26, t27, t28, t32;
  double t33, t35, t36, t38, t39, t41, t42, t43, t44, t46, t48, t50, t52;
  double t55, t56, t58, t60, t62, t64, t66, t67, t69, t71, t73, t75, t77;
  double t78, t80, t82, t84, t86, t88, t89, t90, t91, t92, t120, t146, t172;

  t1 = xi/6.0;
  t2 = xi*xi;
  t3 = t2/6.0;
  t4 = t2*xi;
  t5 = 2.0/3.0*t4;
  t6 = t2*t2;
  t7 = 2.0/3.0*t6;
  t8 = t1-t3-t5+t7;
  t9 = 4.0*eta;
  t10 = eta*eta;
  t11 = 8.0*t10;
  t12 = -1.0/3.0-t9+t11;
  t13 = t8*t12;
  t14 = zeta/6.0;
  t15 = zeta*zeta;
  t16 = t15/6.0;
  t17 = t15*zeta;
  t18 = 2.0/3.0*t17;
  t19 = t15*t15;
  t20 = 2.0/3.0*t19;
  t21 = t14-t16-t18+t20;
  t23 = 4.0/3.0*xi;
  t24 = 8.0/3.0*t2;
  t25 = 4.0/3.0*t4;
  t26 = 8.0/3.0*t6;
  t27 = -t23+t24+t25-t26;
  t28 = t27*t12;
  t32 = 1.0-5.0*t2+4.0*t6;
  t33 = t32*t12;
  t35 = t23+t24-t25-t26;
  t36 = t35*t12;
  t38 = -t1-t3+t5+t7;
  t39 = t38*t12;
  t41 = 8.0*eta;
  t42 = 32.0*t10;
  t43 = 16.0/3.0+t41-t42;
  t44 = t8*t43;
  t46 = t27*t43;
  t48 = t32*t43;
  t50 = t35*t43;
  t52 = t38*t43;
  t55 = -10.0+48.0*t10;
  t56 = t8*t55;
  t58 = t27*t55;
  t60 = t32*t55;
  t62 = t35*t55;
  t64 = t38*t55;
  t66 = 16.0/3.0-t41-t42;
  t67 = t8*t66;
  t69 = t27*t66;
  t71 = t32*t66;
  t73 = t35*t66;
  t75 = t38*t66;
  t77 = -1.0/3.0+t9+t11;
  t78 = t8*t77;
  t80 = t27*t77;
  t82 = t32*t77;
  t84 = t35*t77;
  t86 = t38*t77;
  t88 = 4.0/3.0*zeta;
  t89 = 8.0/3.0*t15;
  t90 = 4.0/3.0*t17;
  t91 = 8.0/3.0*t19;
  t92 = -t88+t89+t90-t91;
  t120 = 1.0-5.0*t15+4.0*t19;
  t146 = t88+t89-t90-t91;
  t172 = -t14-t16+t18+t20;

  values[0] = t13*t21;
  values[1] = t28*t21;
  values[2] = t33*t21;
  values[3] = t36*t21;
  values[4] = t39*t21;
  values[5] = t44*t21;
  values[6] = t46*t21;
  values[7] = t48*t21;
  values[8] = t50*t21;
  values[9] = t52*t21;
  values[10] = t56*t21;
  values[11] = t58*t21;
  values[12] = t60*t21;
  values[13] = t62*t21;
  values[14] = t64*t21;
  values[15] = t67*t21;
  values[16] = t69*t21;
  values[17] = t71*t21;
  values[18] = t73*t21;
  values[19] = t75*t21;
  values[20] = t78*t21;
  values[21] = t80*t21;
  values[22] = t82*t21;
  values[23] = t84*t21;
  values[24] = t86*t21;
  values[25] = t13*t92;
  values[26] = t28*t92;
  values[27] = t33*t92;
  values[28] = t36*t92;
  values[29] = t39*t92;
  values[30] = t44*t92;
  values[31] = t46*t92;
  values[32] = t48*t92;
  values[33] = t50*t92;
  values[34] = t52*t92;
  values[35] = t56*t92;
  values[36] = t58*t92;
  values[37] = t60*t92;
  values[38] = t62*t92;
  values[39] = t64*t92;
  values[40] = t67*t92;
  values[41] = t69*t92;
  values[42] = t71*t92;
  values[43] = t73*t92;
  values[44] = t75*t92;
  values[45] = t78*t92;
  values[46] = t80*t92;
  values[47] = t82*t92;
  values[48] = t84*t92;
  values[49] = t86*t92;
  values[50] = t13*t120;
  values[51] = t28*t120;
  values[52] = t33*t120;
  values[53] = t36*t120;
  values[54] = t39*t120;
  values[55] = t44*t120;
  values[56] = t46*t120;
  values[57] = t48*t120;
  values[58] = t50*t120;
  values[59] = t52*t120;
  values[60] = t56*t120;
  values[61] = t58*t120;
  values[62] = t60*t120;
  values[63] = t62*t120;
  values[64] = t64*t120;
  values[65] = t67*t120;
  values[66] = t69*t120;
  values[67] = t71*t120;
  values[68] = t73*t120;
  values[69] = t75*t120;
  values[70] = t78*t120;
  values[71] = t80*t120;
  values[72] = t82*t120;
  values[73] = t84*t120;
  values[74] = t86*t120;
  values[75] = t13*t146;
  values[76] = t28*t146;
  values[77] = t33*t146;
  values[78] = t36*t146;
  values[79] = t39*t146;
  values[80] = t44*t146;
  values[81] = t46*t146;
  values[82] = t48*t146;
  values[83] = t50*t146;
  values[84] = t52*t146;
  values[85] = t56*t146;
  values[86] = t58*t146;
  values[87] = t60*t146;
  values[88] = t62*t146;
  values[89] = t64*t146;
  values[90] = t67*t146;
  values[91] = t69*t146;
  values[92] = t71*t146;
  values[93] = t73*t146;
  values[94] = t75*t146;
  values[95] = t78*t146;
  values[96] = t80*t146;
  values[97] = t82*t146;
  values[98] = t84*t146;
  values[99] = t86*t146;
  values[100] = t13*t172;
  values[101] = t28*t172;
  values[102] = t33*t172;
  values[103] = t36*t172;
  values[104] = t39*t172;
  values[105] = t44*t172;
  values[106] = t46*t172;
  values[107] = t48*t172;
  values[108] = t50*t172;
  values[109] = t52*t172;
  values[110] = t56*t172;
  values[111] = t58*t172;
  values[112] = t60*t172;
  values[113] = t62*t172;
  values[114] = t64*t172;
  values[115] = t67*t172;
  values[116] = t69*t172;
  values[117] = t71*t172;
  values[118] = t73*t172;
  values[119] = t75*t172;
  values[120] = t78*t172;
  values[121] = t80*t172;
  values[122] = t82*t172;
  values[123] = t84*t172;
  values[124] = t86*t172;
}

static void C_H_Q4_3D_DeriveEtaZeta(double xi, double eta, double zeta,
                             double *values)
{
  double t1, t2, t3, t4, t5, t6, t7, t8, t9, t10, t11, t12, t13, t14, t15;
  double t16, t17, t18, t19, t20, t21, t23, t24, t25, t26, t27, t28, t32;
  double t33, t35, t36, t38, t39, t41, t42, t43, t44, t45, t47, t49, t51;
  double t53, t57, t58, t60, t62, t64, t66, t68, t69, t71, t73, t75, t77;
  double t79, t80, t82, t84, t86, t88, t90, t91, t92, t93, t121, t147, t173;

  t1 = xi/6.0;
  t2 = xi*xi;
  t3 = t2/6.0;
  t4 = t2*xi;
  t5 = 2.0/3.0*t4;
  t6 = t2*t2;
  t7 = 2.0/3.0*t6;
  t8 = t1-t3-t5+t7;
  t9 = eta/3.0;
  t10 = eta*eta;
  t11 = 2.0*t10;
  t12 = t10*eta;
  t13 = 8.0/3.0*t12;
  t14 = 1.0/6.0-t9-t11+t13;
  t15 = t8*t14;
  t16 = zeta/3.0;
  t17 = zeta*zeta;
  t18 = 2.0*t17;
  t19 = t17*zeta;
  t20 = 8.0/3.0*t19;
  t21 = 1.0/6.0-t16-t18+t20;
  t23 = 4.0/3.0*xi;
  t24 = 8.0/3.0*t2;
  t25 = 4.0/3.0*t4;
  t26 = 8.0/3.0*t6;
  t27 = -t23+t24+t25-t26;
  t28 = t27*t14;
  t32 = 1.0-5.0*t2+4.0*t6;
  t33 = t32*t14;
  t35 = t23+t24-t25-t26;
  t36 = t35*t14;
  t38 = -t1-t3+t5+t7;
  t39 = t38*t14;
  t41 = 16.0/3.0*eta;
  t42 = 4.0*t10;
  t43 = 32.0/3.0*t12;
  t44 = -4.0/3.0+t41+t42-t43;
  t45 = t8*t44;
  t47 = t27*t44;
  t49 = t32*t44;
  t51 = t35*t44;
  t53 = t38*t44;
  t57 = -10.0*eta+16.0*t12;
  t58 = t8*t57;
  t60 = t27*t57;
  t62 = t32*t57;
  t64 = t35*t57;
  t66 = t38*t57;
  t68 = 4.0/3.0+t41-t42-t43;
  t69 = t8*t68;
  t71 = t27*t68;
  t73 = t32*t68;
  t75 = t35*t68;
  t77 = t38*t68;
  t79 = -1.0/6.0-t9+t11+t13;
  t80 = t8*t79;
  t82 = t27*t79;
  t84 = t32*t79;
  t86 = t35*t79;
  t88 = t38*t79;
  t90 = 16.0/3.0*zeta;
  t91 = 4.0*t17;
  t92 = 32.0/3.0*t19;
  t93 = -4.0/3.0+t90+t91-t92;
  t121 = -10.0*zeta+16.0*t19;
  t147 = 4.0/3.0+t90-t91-t92;
  t173 = -1.0/6.0-t16+t18+t20;

  values[0] = t15*t21;
  values[1] = t28*t21;
  values[2] = t33*t21;
  values[3] = t36*t21;
  values[4] = t39*t21;
  values[5] = t45*t21;
  values[6] = t47*t21;
  values[7] = t49*t21;
  values[8] = t51*t21;
  values[9] = t53*t21;
  values[10] = t58*t21;
  values[11] = t60*t21;
  values[12] = t62*t21;
  values[13] = t64*t21;
  values[14] = t66*t21;
  values[15] = t69*t21;
  values[16] = t71*t21;
  values[17] = t73*t21;
  values[18] = t75*t21;
  values[19] = t77*t21;
  values[20] = t80*t21;
  values[21] = t82*t21;
  values[22] = t84*t21;
  values[23] = t86*t21;
  values[24] = t88*t21;
  values[25] = t15*t93;
  values[26] = t28*t93;
  values[27] = t33*t93;
  values[28] = t36*t93;
  values[29] = t39*t93;
  values[30] = t45*t93;
  values[31] = t47*t93;
  values[32] = t49*t93;
  values[33] = t51*t93;
  values[34] = t53*t93;
  values[35] = t58*t93;
  values[36] = t60*t93;
  values[37] = t62*t93;
  values[38] = t64*t93;
  values[39] = t66*t93;
  values[40] = t69*t93;
  values[41] = t71*t93;
  values[42] = t73*t93;
  values[43] = t75*t93;
  values[44] = t77*t93;
  values[45] = t80*t93;
  values[46] = t82*t93;
  values[47] = t84*t93;
  values[48] = t86*t93;
  values[49] = t88*t93;
  values[50] = t15*t121;
  values[51] = t28*t121;
  values[52] = t33*t121;
  values[53] = t36*t121;
  values[54] = t39*t121;
  values[55] = t45*t121;
  values[56] = t47*t121;
  values[57] = t49*t121;
  values[58] = t51*t121;
  values[59] = t53*t121;
  values[60] = t58*t121;
  values[61] = t60*t121;
  values[62] = t62*t121;
  values[63] = t64*t121;
  values[64] = t66*t121;
  values[65] = t69*t121;
  values[66] = t71*t121;
  values[67] = t73*t121;
  values[68] = t75*t121;
  values[69] = t77*t121;
  values[70] = t80*t121;
  values[71] = t82*t121;
  values[72] = t84*t121;
  values[73] = t86*t121;
  values[74] = t88*t121;
  values[75] = t15*t147;
  values[76] = t28*t147;
  values[77] = t33*t147;
  values[78] = t36*t147;
  values[79] = t39*t147;
  values[80] = t45*t147;
  values[81] = t47*t147;
  values[82] = t49*t147;
  values[83] = t51*t147;
  values[84] = t53*t147;
  values[85] = t58*t147;
  values[86] = t60*t147;
  values[87] = t62*t147;
  values[88] = t64*t147;
  values[89] = t66*t147;
  values[90] = t69*t147;
  values[91] = t71*t147;
  values[92] = t73*t147;
  values[93] = t75*t147;
  values[94] = t77*t147;
  values[95] = t80*t147;
  values[96] = t82*t147;
  values[97] = t84*t147;
  values[98] = t86*t147;
  values[99] = t88*t147;
  values[100] = t15*t173;
  values[101] = t28*t173;
  values[102] = t33*t173;
  values[103] = t36*t173;
  values[104] = t39*t173;
  values[105] = t45*t173;
  values[106] = t47*t173;
  values[107] = t49*t173;
  values[108] = t51*t173;
  values[109] = t53*t173;
  values[110] = t58*t173;
  values[111] = t60*t173;
  values[112] = t62*t173;
  values[113] = t64*t173;
  values[114] = t66*t173;
  values[115] = t69*t173;
  values[116] = t71*t173;
  values[117] = t73*t173;
  values[118] = t75*t173;
  values[119] = t77*t173;
  values[120] = t80*t173;
  values[121] = t82*t173;
  values[122] = t84*t173;
  values[123] = t86*t173;
  values[124] = t88*t173;
}

static void C_H_Q4_3D_DeriveZetaZeta(double xi, double eta, double zeta,
                             double *values)
{
  double t1, t2, t3, t4, t5, t6, t7, t8, t9, t10, t11, t12, t13, t14, t15;
  double t16, t17, t18, t19, t20, t21, t23, t24, t25, t26, t27, t28, t32;
  double t33, t35, t36, t38, t39, t41, t42, t43, t44, t45, t46, t48, t50;
  double t52, t54, t58, t59, t61, t63, t65, t67, t69, t70, t72, t74, t76;
  double t78, t80, t81, t83, t85, t87, t89, t91, t92, t93, t120, t146, t172;

  t1 = xi/6.0;
  t2 = xi*xi;
  t3 = t2/6.0;
  t4 = t2*xi;
  t5 = 2.0/3.0*t4;
  t6 = t2*t2;
  t7 = 2.0/3.0*t6;
  t8 = t1-t3-t5+t7;
  t9 = eta/6.0;
  t10 = eta*eta;
  t11 = t10/6.0;
  t12 = t10*eta;
  t13 = 2.0/3.0*t12;
  t14 = t10*t10;
  t15 = 2.0/3.0*t14;
  t16 = t9-t11-t13+t15;
  t17 = t8*t16;
  t18 = 4.0*zeta;
  t19 = zeta*zeta;
  t20 = 8.0*t19;
  t21 = -1.0/3.0-t18+t20;
  t23 = 4.0/3.0*xi;
  t24 = 8.0/3.0*t2;
  t25 = 4.0/3.0*t4;
  t26 = 8.0/3.0*t6;
  t27 = -t23+t24+t25-t26;
  t28 = t27*t16;
  t32 = 1.0-5.0*t2+4.0*t6;
  t33 = t32*t16;
  t35 = t23+t24-t25-t26;
  t36 = t35*t16;
  t38 = -t1-t3+t5+t7;
  t39 = t38*t16;
  t41 = 4.0/3.0*eta;
  t42 = 8.0/3.0*t10;
  t43 = 4.0/3.0*t12;
  t44 = 8.0/3.0*t14;
  t45 = -t41+t42+t43-t44;
  t46 = t8*t45;
  t48 = t27*t45;
  t50 = t32*t45;
  t52 = t35*t45;
  t54 = t38*t45;
  t58 = 1.0-5.0*t10+4.0*t14;
  t59 = t8*t58;
  t61 = t27*t58;
  t63 = t32*t58;
  t65 = t35*t58;
  t67 = t38*t58;
  t69 = t41+t42-t43-t44;
  t70 = t8*t69;
  t72 = t27*t69;
  t74 = t32*t69;
  t76 = t35*t69;
  t78 = t38*t69;
  t80 = -t9-t11+t13+t15;
  t81 = t8*t80;
  t83 = t27*t80;
  t85 = t32*t80;
  t87 = t35*t80;
  t89 = t38*t80;
  t91 = 8.0*zeta;
  t92 = 32.0*t19;
  t93 = 16.0/3.0+t91-t92;
  t120 = -10.0+48.0*t19;
  t146 = 16.0/3.0-t91-t92;
  t172 = -1.0/3.0+t18+t20;

  values[0] = t17*t21;
  values[1] = t28*t21;
  values[2] = t33*t21;
  values[3] = t36*t21;
  values[4] = t39*t21;
  values[5] = t46*t21;
  values[6] = t48*t21;
  values[7] = t50*t21;
  values[8] = t52*t21;
  values[9] = t54*t21;
  values[10] = t59*t21;
  values[11] = t61*t21;
  values[12] = t63*t21;
  values[13] = t65*t21;
  values[14] = t67*t21;
  values[15] = t70*t21;
  values[16] = t72*t21;
  values[17] = t74*t21;
  values[18] = t76*t21;
  values[19] = t78*t21;
  values[20] = t81*t21;
  values[21] = t83*t21;
  values[22] = t85*t21;
  values[23] = t87*t21;
  values[24] = t89*t21;
  values[25] = t17*t93;
  values[26] = t28*t93;
  values[27] = t33*t93;
  values[28] = t36*t93;
  values[29] = t39*t93;
  values[30] = t46*t93;
  values[31] = t48*t93;
  values[32] = t50*t93;
  values[33] = t52*t93;
  values[34] = t54*t93;
  values[35] = t59*t93;
  values[36] = t61*t93;
  values[37] = t63*t93;
  values[38] = t65*t93;
  values[39] = t67*t93;
  values[40] = t70*t93;
  values[41] = t72*t93;
  values[42] = t74*t93;
  values[43] = t76*t93;
  values[44] = t78*t93;
  values[45] = t81*t93;
  values[46] = t83*t93;
  values[47] = t85*t93;
  values[48] = t87*t93;
  values[49] = t89*t93;
  values[50] = t17*t120;
  values[51] = t28*t120;
  values[52] = t33*t120;
  values[53] = t36*t120;
  values[54] = t39*t120;
  values[55] = t46*t120;
  values[56] = t48*t120;
  values[57] = t50*t120;
  values[58] = t52*t120;
  values[59] = t54*t120;
  values[60] = t59*t120;
  values[61] = t61*t120;
  values[62] = t63*t120;
  values[63] = t65*t120;
  values[64] = t67*t120;
  values[65] = t70*t120;
  values[66] = t72*t120;
  values[67] = t74*t120;
  values[68] = t76*t120;
  values[69] = t78*t120;
  values[70] = t81*t120;
  values[71] = t83*t120;
  values[72] = t85*t120;
  values[73] = t87*t120;
  values[74] = t89*t120;
  values[75] = t17*t146;
  values[76] = t28*t146;
  values[77] = t33*t146;
  values[78] = t36*t146;
  values[79] = t39*t146;
  values[80] = t46*t146;
  values[81] = t48*t146;
  values[82] = t50*t146;
  values[83] = t52*t146;
  values[84] = t54*t146;
  values[85] = t59*t146;
  values[86] = t61*t146;
  values[87] = t63*t146;
  values[88] = t65*t146;
  values[89] = t67*t146;
  values[90] = t70*t146;
  values[91] = t72*t146;
  values[92] = t74*t146;
  values[93] = t76*t146;
  values[94] = t78*t146;
  values[95] = t81*t146;
  values[96] = t83*t146;
  values[97] = t85*t146;
  values[98] = t87*t146;
  values[99] = t89*t146;
  values[100] = t17*t172;
  values[101] = t28*t172;
  values[102] = t33*t172;
  values[103] = t36*t172;
  values[104] = t39*t172;
  values[105] = t46*t172;
  values[106] = t48*t172;
  values[107] = t50*t172;
  values[108] = t52*t172;
  values[109] = t54*t172;
  values[110] = t59*t172;
  values[111] = t61*t172;
  values[112] = t63*t172;
  values[113] = t65*t172;
  values[114] = t67*t172;
  values[115] = t70*t172;
  values[116] = t72*t172;
  values[117] = t74*t172;
  values[118] = t76*t172;
  values[119] = t78*t172;
  values[120] = t81*t172;
  values[121] = t83*t172;
  values[122] = t85*t172;
  values[123] = t87*t172;
  values[124] = t89*t172;
}

TBaseFunct3D *BF_C_H_Q4_3D_Obj = 
new TBaseFunct3D(125, BF_C_H_Q4_3D, BFUnitHexahedron, 
                 C_H_Q4_3D_Funct, C_H_Q4_3D_DeriveXi,
                 C_H_Q4_3D_DeriveEta, C_H_Q4_3D_DeriveZeta,
                 C_H_Q4_3D_DeriveXiXi, C_H_Q4_3D_DeriveXiEta,
                 C_H_Q4_3D_DeriveXiZeta, C_H_Q4_3D_DeriveEtaEta,
                 C_H_Q4_3D_DeriveEtaZeta, C_H_Q4_3D_DeriveZetaZeta,
                 4, 4,
                 0, NULL);