BF_C_H_Q3_3D.h

// ***********************************************************************
// Q3 element, conforming, 3D
// ***********************************************************************

static void C_H_Q3_3D_Funct(double xi, double eta, double zeta,
                          double *values)
{
  double t1, t3, t4, t5, t6, t7, t9, t10, t11, t14, t15, t16, t19, t20;
  double t22, t24, t25, t26, t27, t30, t31, t34, t35, t40, t41, t42, t44;
  double t45, t46, t48, t50, t51, t53, t56, t57, t64, t65, t70, t75, t80;
  double t85, t88, t89, t90, t95, t100, t105, t112, t113, t118, t123, t128;

  t1 = xi-1.0;
  t3 = t1*(xi+1.0/3.0);
  t4 = xi-1.0/3.0;
  t5 = eta-1.0;
  t6 = t4*t5;
  t7 = t3*t6;
  t9 = eta-1.0/3.0;
  t10 = (eta+1.0/3.0)*t9;
  t11 = zeta-1.0;
  t14 = zeta-1.0/3.0;
  t15 = t11*(zeta+1.0/3.0)*t14;
  t16 = t10*t15;
  t19 = t1*(xi+1.0);
  t20 = t19*t6;
  t22 = -xi-1.0;
  t24 = t22*(-xi+1.0);
  t25 = -xi-1.0/3.0;
  t26 = t25*t5;
  t27 = t24*t26;
  t30 = t22*(-xi+1.0/3.0);
  t31 = t30*t26;
  t34 = (eta+1.0)*t9;
  t35 = t34*t15;
  t40 = -eta-1.0;
  t41 = t4*t40;
  t42 = t3*t41;
  t44 = -eta-1.0/3.0;
  t45 = (-eta+1.0)*t44;
  t46 = t45*t15;
  t48 = t19*t41;
  t50 = t25*t40;
  t51 = t24*t50;
  t53 = t30*t50;
  t56 = (-eta+1.0/3.0)*t44;
  t57 = t56*t15;
  t64 = t11*(zeta+1.0)*t14;
  t65 = t10*t64;
  t70 = t34*t64;
  t75 = t45*t64;
  t80 = t56*t64;
  t85 = -zeta-1.0;
  t88 = -zeta-1.0/3.0;
  t89 = t85*(-zeta+1.0)*t88;
  t90 = t10*t89;
  t95 = t34*t89;
  t100 = t45*t89;
  t105 = t56*t89;
  t112 = t85*(-zeta+1.0/3.0)*t88;
  t113 = t10*t112;
  t118 = t34*t112;
  t123 = t45*t112;
  t128 = t56*t112;

  values[0] = -729.0/4096.0*t7*t16;
  values[1] = 2187.0/4096.0*t20*t16;
  values[2] = 2187.0/4096.0*t27*t16;
  values[3] = -729.0/4096.0*t31*t16;
  values[4] = 2187.0/4096.0*t7*t35;
  values[5] = -6561.0/4096.0*t20*t35;
  values[6] = -6561.0/4096.0*t27*t35;
  values[7] = 2187.0/4096.0*t31*t35;
  values[8] = 2187.0/4096.0*t42*t46;
  values[9] = -6561.0/4096.0*t48*t46;
  values[10] = -6561.0/4096.0*t51*t46;
  values[11] = 2187.0/4096.0*t53*t46;
  values[12] = -729.0/4096.0*t42*t57;
  values[13] = 2187.0/4096.0*t48*t57;
  values[14] = 2187.0/4096.0*t51*t57;
  values[15] = -729.0/4096.0*t53*t57;
  values[16] = 2187.0/4096.0*t7*t65;
  values[17] = -6561.0/4096.0*t20*t65;
  values[18] = -6561.0/4096.0*t27*t65;
  values[19] = 2187.0/4096.0*t31*t65;
  values[20] = -6561.0/4096.0*t7*t70;
  values[21] = 19683.0/4096.0*t20*t70;
  values[22] = 19683.0/4096.0*t27*t70;
  values[23] = -6561.0/4096.0*t31*t70;
  values[24] = -6561.0/4096.0*t42*t75;
  values[25] = 19683.0/4096.0*t48*t75;
  values[26] = 19683.0/4096.0*t51*t75;
  values[27] = -6561.0/4096.0*t53*t75;
  values[28] = 2187.0/4096.0*t42*t80;
  values[29] = -6561.0/4096.0*t48*t80;
  values[30] = -6561.0/4096.0*t51*t80;
  values[31] = 2187.0/4096.0*t53*t80;
  values[32] = 2187.0/4096.0*t7*t90;
  values[33] = -6561.0/4096.0*t20*t90;
  values[34] = -6561.0/4096.0*t27*t90;
  values[35] = 2187.0/4096.0*t31*t90;
  values[36] = -6561.0/4096.0*t7*t95;
  values[37] = 19683.0/4096.0*t20*t95;
  values[38] = 19683.0/4096.0*t27*t95;
  values[39] = -6561.0/4096.0*t31*t95;
  values[40] = -6561.0/4096.0*t42*t100;
  values[41] = 19683.0/4096.0*t48*t100;
  values[42] = 19683.0/4096.0*t51*t100;
  values[43] = -6561.0/4096.0*t53*t100;
  values[44] = 2187.0/4096.0*t42*t105;
  values[45] = -6561.0/4096.0*t48*t105;
  values[46] = -6561.0/4096.0*t51*t105;
  values[47] = 2187.0/4096.0*t53*t105;
  values[48] = -729.0/4096.0*t7*t113;
  values[49] = 2187.0/4096.0*t20*t113;
  values[50] = 2187.0/4096.0*t27*t113;
  values[51] = -729.0/4096.0*t31*t113;
  values[52] = 2187.0/4096.0*t7*t118;
  values[53] = -6561.0/4096.0*t20*t118;
  values[54] = -6561.0/4096.0*t27*t118;
  values[55] = 2187.0/4096.0*t31*t118;
  values[56] = 2187.0/4096.0*t42*t123;
  values[57] = -6561.0/4096.0*t48*t123;
  values[58] = -6561.0/4096.0*t51*t123;
  values[59] = 2187.0/4096.0*t53*t123;
  values[60] = -729.0/4096.0*t42*t128;
  values[61] = 2187.0/4096.0*t48*t128;
  values[62] = 2187.0/4096.0*t51*t128;
  values[63] = -729.0/4096.0*t53*t128;
}

static void C_H_Q3_3D_DeriveXi(double xi, double eta, double zeta,
                             double *values)
{
  double t1, t2, t3, t4, t6, t7, t8, t9, t10, t12, t13, t14, t16, t17, t18;
  double t19, t20, t21, t24, t25, t26, t28, t29, t32, t33, t34, t35, t37;
  double t38, t39, t40, t41, t42, t45, t46, t47, t49, t50, t54, t55, t57;
  double t58, t59, t62, t64, t67, t69, t70, t71, t74, t76, t79, t81, t82;
  double t83, t84, t85, t87, t88, t89, t92, t94, t97, t99, t100, t101, t104;
  double t106, t110, t111, t113, t114, t115, t118, t120, t123, t125, t126;
  double t127, t130, t132, t136, t137, t139, t146, t153, t160, t166, t168;
  double t175, t182, t189, t195, t196, t198, t199, t200, t202, t209, t216;
  double t223, t229, t230, t232, t239, t246, t253, t260, t261, t263, t270;
  double t277, t284, t290, t292, t299, t306, t313;

  t1 = xi+1.0/3.0;
  t2 = xi-1.0/3.0;
  t3 = t1*t2;
  t4 = eta-1.0;
  t6 = t4*(eta+1.0/3.0);
  t7 = t3*t6;
  t8 = eta-1.0/3.0;
  t9 = zeta-1.0;
  t10 = t8*t9;
  t12 = zeta-1.0/3.0;
  t13 = (zeta+1.0/3.0)*t12;
  t14 = t10*t13;
  t16 = xi-1.0;
  t17 = t16*t2;
  t18 = t17*t6;
  t19 = t18*t14;
  t20 = t16*t1;
  t21 = t20*t6;
  t24 = xi+1.0;
  t25 = t24*t2;
  t26 = t25*t6;
  t28 = t16*t24;
  t29 = t28*t6;
  t32 = -xi+1.0;
  t33 = -xi-1.0/3.0;
  t34 = t32*t33;
  t35 = t34*t6;
  t37 = -xi-1.0;
  t38 = t37*t33;
  t39 = t38*t6;
  t40 = t39*t14;
  t41 = t37*t32;
  t42 = t41*t6;
  t45 = -xi+1.0/3.0;
  t46 = t45*t33;
  t47 = t46*t6;
  t49 = t37*t45;
  t50 = t49*t6;
  t54 = t4*(eta+1.0);
  t55 = t3*t54;
  t57 = t17*t54;
  t58 = t57*t14;
  t59 = t20*t54;
  t62 = t25*t54;
  t64 = t28*t54;
  t67 = t34*t54;
  t69 = t38*t54;
  t70 = t69*t14;
  t71 = t41*t54;
  t74 = t46*t54;
  t76 = t49*t54;
  t79 = -eta-1.0;
  t81 = t79*(-eta+1.0);
  t82 = t3*t81;
  t83 = -eta-1.0/3.0;
  t84 = t83*t9;
  t85 = t84*t13;
  t87 = t17*t81;
  t88 = t87*t85;
  t89 = t20*t81;
  t92 = t25*t81;
  t94 = t28*t81;
  t97 = t34*t81;
  t99 = t38*t81;
  t100 = t99*t85;
  t101 = t41*t81;
  t104 = t46*t81;
  t106 = t49*t81;
  t110 = t79*(-eta+1.0/3.0);
  t111 = t3*t110;
  t113 = t17*t110;
  t114 = t113*t85;
  t115 = t20*t110;
  t118 = t25*t110;
  t120 = t28*t110;
  t123 = t34*t110;
  t125 = t38*t110;
  t126 = t125*t85;
  t127 = t41*t110;
  t130 = t46*t110;
  t132 = t49*t110;
  t136 = (zeta+1.0)*t12;
  t137 = t10*t136;
  t139 = t18*t137;
  t146 = t39*t137;
  t153 = t57*t137;
  t160 = t69*t137;
  t166 = t84*t136;
  t168 = t87*t166;
  t175 = t99*t166;
  t182 = t113*t166;
  t189 = t125*t166;
  t195 = -zeta-1.0;
  t196 = t8*t195;
  t198 = -zeta-1.0/3.0;
  t199 = (-zeta+1.0)*t198;
  t200 = t196*t199;
  t202 = t18*t200;
  t209 = t39*t200;
  t216 = t57*t200;
  t223 = t69*t200;
  t229 = t83*t195;
  t230 = t229*t199;
  t232 = t87*t230;
  t239 = t99*t230;
  t246 = t113*t230;
  t253 = t125*t230;
  t260 = (-zeta+1.0/3.0)*t198;
  t261 = t196*t260;
  t263 = t18*t261;
  t270 = t39*t261;
  t277 = t57*t261;
  t284 = t69*t261;
  t290 = t229*t260;
  t292 = t87*t290;
  t299 = t99*t290;
  t306 = t113*t290;
  t313 = t125*t290;

  values[0] = -729.0/4096.0*t7*t14-729.0/4096.0*t19-729.0/4096.0*t21*t14;
  values[1] = 2187.0/4096.0*t26*t14+2187.0/4096.0*t19+2187.0/4096.0*t29*t14;
  values[2] = -2187.0/4096.0*t35*t14-2187.0/4096.0*t40-2187.0/4096.0*t42*t14;
  values[3] = 729.0/4096.0*t47*t14+729.0/4096.0*t40+729.0/4096.0*t50*t14;
  values[4] = 2187.0/4096.0*t55*t14+2187.0/4096.0*t58+2187.0/4096.0*t59*t14;
  values[5] = -6561.0/4096.0*t62*t14-6561.0/4096.0*t58-6561.0/4096.0*t64*t14;
  values[6] = 6561.0/4096.0*t67*t14+6561.0/4096.0*t70+6561.0/4096.0*t71*t14;
  values[7] = -2187.0/4096.0*t74*t14-2187.0/4096.0*t70-2187.0/4096.0*t76*t14;
  values[8] = 2187.0/4096.0*t82*t85+2187.0/4096.0*t88+2187.0/4096.0*t89*t85;
  values[9] = -6561.0/4096.0*t92*t85-6561.0/4096.0*t88-6561.0/4096.0*t94*t85;
  values[10] = 6561.0/4096.0*t97*t85+6561.0/4096.0*t100+6561.0/4096.0*t101*t85;
  values[11] = -2187.0/4096.0*t104*t85-2187.0/4096.0*t100-2187.0/4096.0*t106*t85;
  values[12] = -729.0/4096.0*t111*t85-729.0/4096.0*t114-729.0/4096.0*t115*t85;
  values[13] = 2187.0/4096.0*t118*t85+2187.0/4096.0*t114+2187.0/4096.0*t120*t85;
  values[14] = -2187.0/4096.0*t123*t85-2187.0/4096.0*t126-2187.0/4096.0*t127*t85;
  values[15] = 729.0/4096.0*t130*t85+729.0/4096.0*t126+729.0/4096.0*t132*t85;
  values[16] = 2187.0/4096.0*t7*t137+2187.0/4096.0*t139+2187.0/4096.0*t21*t137;
  values[17] = -6561.0/4096.0*t26*t137-6561.0/4096.0*t139-6561.0/4096.0*t29*t137;
  values[18] = 6561.0/4096.0*t35*t137+6561.0/4096.0*t146+6561.0/4096.0*t42*t137;
  values[19] = -2187.0/4096.0*t47*t137-2187.0/4096.0*t146-2187.0/4096.0*t50*t137;
  values[20] = -6561.0/4096.0*t55*t137-6561.0/4096.0*t153-6561.0/4096.0*t59*t137;
  values[21] = 19683.0/4096.0*t62*t137+19683.0/4096.0*t153+19683.0/4096.0*t64*t137;
  values[22] = -19683.0/4096.0*t67*t137-19683.0/4096.0*t160-19683.0/4096.0*t71*t137;
  values[23] = 6561.0/4096.0*t74*t137+6561.0/4096.0*t160+6561.0/4096.0*t76*t137;
  values[24] = -6561.0/4096.0*t82*t166-6561.0/4096.0*t168-6561.0/4096.0*t89*t166;
  values[25] = 19683.0/4096.0*t92*t166+19683.0/4096.0*t168+19683.0/4096.0*t94*t166;
  values[26] = -19683.0/4096.0*t97*t166-19683.0/4096.0*t175-19683.0/4096.0*t101*t166;
  values[27] = 6561.0/4096.0*t104*t166+6561.0/4096.0*t175+6561.0/4096.0*t106*t166;
  values[28] = 2187.0/4096.0*t111*t166+2187.0/4096.0*t182+2187.0/4096.0*t115*t166;
  values[29] = -6561.0/4096.0*t118*t166-6561.0/4096.0*t182-6561.0/4096.0*t120*t166;
  values[30] = 6561.0/4096.0*t123*t166+6561.0/4096.0*t189+6561.0/4096.0*t127*t166;
  values[31] = -2187.0/4096.0*t130*t166-2187.0/4096.0*t189-2187.0/4096.0*t132*t166;
  values[32] = 2187.0/4096.0*t7*t200+2187.0/4096.0*t202+2187.0/4096.0*t21*t200;
  values[33] = -6561.0/4096.0*t26*t200-6561.0/4096.0*t202-6561.0/4096.0*t29*t200;
  values[34] = 6561.0/4096.0*t35*t200+6561.0/4096.0*t209+6561.0/4096.0*t42*t200;
  values[35] = -2187.0/4096.0*t47*t200-2187.0/4096.0*t209-2187.0/4096.0*t50*t200;
  values[36] = -6561.0/4096.0*t55*t200-6561.0/4096.0*t216-6561.0/4096.0*t59*t200;
  values[37] = 19683.0/4096.0*t62*t200+19683.0/4096.0*t216+19683.0/4096.0*t64*t200;
  values[38] = -19683.0/4096.0*t67*t200-19683.0/4096.0*t223-19683.0/4096.0*t71*t200;
  values[39] = 6561.0/4096.0*t74*t200+6561.0/4096.0*t223+6561.0/4096.0*t76*t200;
  values[40] = -6561.0/4096.0*t82*t230-6561.0/4096.0*t232-6561.0/4096.0*t89*t230;
  values[41] = 19683.0/4096.0*t92*t230+19683.0/4096.0*t232+19683.0/4096.0*t94*t230;
  values[42] = -19683.0/4096.0*t97*t230-19683.0/4096.0*t239-19683.0/4096.0*t101*t230;
  values[43] = 6561.0/4096.0*t104*t230+6561.0/4096.0*t239+6561.0/4096.0*t106*t230;
  values[44] = 2187.0/4096.0*t111*t230+2187.0/4096.0*t246+2187.0/4096.0*t115*t230;
  values[45] = -6561.0/4096.0*t118*t230-6561.0/4096.0*t246-6561.0/4096.0*t120*t230;
  values[46] = 6561.0/4096.0*t123*t230+6561.0/4096.0*t253+6561.0/4096.0*t127*t230;
  values[47] = -2187.0/4096.0*t130*t230-2187.0/4096.0*t253-2187.0/4096.0*t132*t230;
  values[48] = -729.0/4096.0*t7*t261-729.0/4096.0*t263-729.0/4096.0*t21*t261;
  values[49] = 2187.0/4096.0*t26*t261+2187.0/4096.0*t263+2187.0/4096.0*t29*t261;
  values[50] = -2187.0/4096.0*t35*t261-2187.0/4096.0*t270-2187.0/4096.0*t42*t261;
  values[51] = 729.0/4096.0*t47*t261+729.0/4096.0*t270+729.0/4096.0*t50*t261;
  values[52] = 2187.0/4096.0*t55*t261+2187.0/4096.0*t277+2187.0/4096.0*t59*t261;
  values[53] = -6561.0/4096.0*t62*t261-6561.0/4096.0*t277-6561.0/4096.0*t64*t261;
  values[54] = 6561.0/4096.0*t67*t261+6561.0/4096.0*t284+6561.0/4096.0*t71*t261;
  values[55] = -2187.0/4096.0*t74*t261-2187.0/4096.0*t284-2187.0/4096.0*t76*t261;
  values[56] = 2187.0/4096.0*t82*t290+2187.0/4096.0*t292+2187.0/4096.0*t89*t290;
  values[57] = -6561.0/4096.0*t92*t290-6561.0/4096.0*t292-6561.0/4096.0*t94*t290;
  values[58] = 6561.0/4096.0*t97*t290+6561.0/4096.0*t299+6561.0/4096.0*t101*t290;
  values[59] = -2187.0/4096.0*t104*t290-2187.0/4096.0*t299-2187.0/4096.0*t106*t290;
  values[60] = -729.0/4096.0*t111*t290-729.0/4096.0*t306-729.0/4096.0*t115*t290;
  values[61] = 2187.0/4096.0*t118*t290+2187.0/4096.0*t306+2187.0/4096.0*t120*t290;
  values[62] = -2187.0/4096.0*t123*t290-2187.0/4096.0*t313-2187.0/4096.0*t127*t290;
  values[63] = 729.0/4096.0*t130*t290+729.0/4096.0*t313+729.0/4096.0*t132*t290;
}

static void C_H_Q3_3D_DeriveEta(double xi, double eta, double zeta,
                             double *values)
{
  double t1, t3, t4, t5, t6, t7, t8, t9, t10, t12, t13, t14, t16, t17;
  double t18, t19, t20, t21, t25, t26, t28, t29, t32, t34, t35, t36, t37;
  double t39, t40, t41, t45, t46, t48, t49, t52, t53, t54, t56, t57, t60;
  double t64, t65, t69, t73, t74, t75, t76, t77, t78, t80, t81, t82, t83;
  double t84, t85, t88, t90, t91, t94, t95, t97, t98, t99, t102, t104;
  double t105, t108, t109, t110, t112, t113, t116, t120, t121, t125, t130;
  double t131, t133, t134, t138, t142, t146, t150, t162, t164, t165, t169;
  double t173, t177, t181, t193, t194, t196, t197, t198, t200, t201, t202;
  double t206, t210, t214, t218, t219, t231, t232, t234, t235, t236, t240;
  double t244, t248, t252, t253, t266, t267, t269, t270, t274, t278, t282;
  double t286, t298, t300, t301, t305, t309, t313, t317;

  t1 = xi-1.0;
  t3 = t1*(xi+1.0/3.0);
  t4 = xi-1.0/3.0;
  t5 = eta+1.0/3.0;
  t6 = t4*t5;
  t7 = t3*t6;
  t8 = eta-1.0/3.0;
  t9 = zeta-1.0;
  t10 = t8*t9;
  t12 = zeta-1.0/3.0;
  t13 = (zeta+1.0/3.0)*t12;
  t14 = t10*t13;
  t16 = eta-1.0;
  t17 = t4*t16;
  t18 = t3*t17;
  t19 = t18*t14;
  t20 = t5*t9;
  t21 = t20*t13;
  t25 = t1*(xi+1.0);
  t26 = t25*t6;
  t28 = t25*t17;
  t29 = t28*t14;
  t32 = -xi-1.0;
  t34 = t32*(-xi+1.0);
  t35 = -xi-1.0/3.0;
  t36 = t35*t5;
  t37 = t34*t36;
  t39 = t35*t16;
  t40 = t34*t39;
  t41 = t40*t14;
  t45 = t32*(-xi+1.0/3.0);
  t46 = t45*t36;
  t48 = t45*t39;
  t49 = t48*t14;
  t52 = eta+1.0;
  t53 = t4*t52;
  t54 = t3*t53;
  t56 = t52*t9;
  t57 = t56*t13;
  t60 = t25*t53;
  t64 = t35*t52;
  t65 = t34*t64;
  t69 = t45*t64;
  t73 = -eta+1.0;
  t74 = t4*t73;
  t75 = t3*t74;
  t76 = -eta-1.0/3.0;
  t77 = t76*t9;
  t78 = t77*t13;
  t80 = -eta-1.0;
  t81 = t4*t80;
  t82 = t3*t81;
  t83 = t82*t78;
  t84 = t73*t9;
  t85 = t84*t13;
  t88 = t25*t74;
  t90 = t25*t81;
  t91 = t90*t78;
  t94 = t35*t73;
  t95 = t34*t94;
  t97 = t35*t80;
  t98 = t34*t97;
  t99 = t98*t78;
  t102 = t45*t94;
  t104 = t45*t97;
  t105 = t104*t78;
  t108 = -eta+1.0/3.0;
  t109 = t4*t108;
  t110 = t3*t109;
  t112 = t108*t9;
  t113 = t112*t13;
  t116 = t25*t109;
  t120 = t35*t108;
  t121 = t34*t120;
  t125 = t45*t120;
  t130 = (zeta+1.0)*t12;
  t131 = t10*t130;
  t133 = t18*t131;
  t134 = t20*t130;
  t138 = t28*t131;
  t142 = t40*t131;
  t146 = t48*t131;
  t150 = t56*t130;
  t162 = t77*t130;
  t164 = t82*t162;
  t165 = t84*t130;
  t169 = t90*t162;
  t173 = t98*t162;
  t177 = t104*t162;
  t181 = t112*t130;
  t193 = -zeta-1.0;
  t194 = t8*t193;
  t196 = -zeta-1.0/3.0;
  t197 = (-zeta+1.0)*t196;
  t198 = t194*t197;
  t200 = t18*t198;
  t201 = t5*t193;
  t202 = t201*t197;
  t206 = t28*t198;
  t210 = t40*t198;
  t214 = t48*t198;
  t218 = t52*t193;
  t219 = t218*t197;
  t231 = t76*t193;
  t232 = t231*t197;
  t234 = t82*t232;
  t235 = t73*t193;
  t236 = t235*t197;
  t240 = t90*t232;
  t244 = t98*t232;
  t248 = t104*t232;
  t252 = t108*t193;
  t253 = t252*t197;
  t266 = (-zeta+1.0/3.0)*t196;
  t267 = t194*t266;
  t269 = t18*t267;
  t270 = t201*t266;
  t274 = t28*t267;
  t278 = t40*t267;
  t282 = t48*t267;
  t286 = t218*t266;
  t298 = t231*t266;
  t300 = t82*t298;
  t301 = t235*t266;
  t305 = t90*t298;
  t309 = t98*t298;
  t313 = t104*t298;
  t317 = t252*t266;

  values[0] = -729.0/4096.0*t7*t14-729.0/4096.0*t19-729.0/4096.0*t18*t21;
  values[1] = 2187.0/4096.0*t26*t14+2187.0/4096.0*t29+2187.0/4096.0*t28*t21;
  values[2] = 2187.0/4096.0*t37*t14+2187.0/4096.0*t41+2187.0/4096.0*t40*t21;
  values[3] = -729.0/4096.0*t46*t14-729.0/4096.0*t49-729.0/4096.0*t48*t21;
  values[4] = 2187.0/4096.0*t54*t14+2187.0/4096.0*t19+2187.0/4096.0*t18*t57;
  values[5] = -6561.0/4096.0*t60*t14-6561.0/4096.0*t29-6561.0/4096.0*t28*t57;
  values[6] = -6561.0/4096.0*t65*t14-6561.0/4096.0*t41-6561.0/4096.0*t40*t57;
  values[7] = 2187.0/4096.0*t69*t14+2187.0/4096.0*t49+2187.0/4096.0*t48*t57;
  values[8] = -2187.0/4096.0*t75*t78-2187.0/4096.0*t83-2187.0/4096.0*t82*t85;
  values[9] = 6561.0/4096.0*t88*t78+6561.0/4096.0*t91+6561.0/4096.0*t90*t85;
  values[10] = 6561.0/4096.0*t95*t78+6561.0/4096.0*t99+6561.0/4096.0*t98*t85;
  values[11] = -2187.0/4096.0*t102*t78-2187.0/4096.0*t105-2187.0/4096.0*t104*t85;
  values[12] = 729.0/4096.0*t110*t78+729.0/4096.0*t83+729.0/4096.0*t82*t113;
  values[13] = -2187.0/4096.0*t116*t78-2187.0/4096.0*t91-2187.0/4096.0*t90*t113;
  values[14] = -2187.0/4096.0*t121*t78-2187.0/4096.0*t99-2187.0/4096.0*t98*t113;
  values[15] = 729.0/4096.0*t125*t78+729.0/4096.0*t105+729.0/4096.0*t104*t113;
  values[16] = 2187.0/4096.0*t7*t131+2187.0/4096.0*t133+2187.0/4096.0*t18*t134;
  values[17] = -6561.0/4096.0*t26*t131-6561.0/4096.0*t138-6561.0/4096.0*t28*t134;
  values[18] = -6561.0/4096.0*t37*t131-6561.0/4096.0*t142-6561.0/4096.0*t40*t134;
  values[19] = 2187.0/4096.0*t46*t131+2187.0/4096.0*t146+2187.0/4096.0*t48*t134;
  values[20] = -6561.0/4096.0*t54*t131-6561.0/4096.0*t133-6561.0/4096.0*t18*t150;
  values[21] = 19683.0/4096.0*t60*t131+19683.0/4096.0*t138+19683.0/4096.0*t28*t150;
  values[22] = 19683.0/4096.0*t65*t131+19683.0/4096.0*t142+19683.0/4096.0*t40*t150;
  values[23] = -6561.0/4096.0*t69*t131-6561.0/4096.0*t146-6561.0/4096.0*t48*t150;
  values[24] = 6561.0/4096.0*t75*t162+6561.0/4096.0*t164+6561.0/4096.0*t82*t165;
  values[25] = -19683.0/4096.0*t88*t162-19683.0/4096.0*t169-19683.0/4096.0*t90*t165;
  values[26] = -19683.0/4096.0*t95*t162-19683.0/4096.0*t173-19683.0/4096.0*t98*t165;
  values[27] = 6561.0/4096.0*t102*t162+6561.0/4096.0*t177+6561.0/4096.0*t104*t165;
  values[28] = -2187.0/4096.0*t110*t162-2187.0/4096.0*t164-2187.0/4096.0*t82*t181;
  values[29] = 6561.0/4096.0*t116*t162+6561.0/4096.0*t169+6561.0/4096.0*t90*t181;
  values[30] = 6561.0/4096.0*t121*t162+6561.0/4096.0*t173+6561.0/4096.0*t98*t181;
  values[31] = -2187.0/4096.0*t125*t162-2187.0/4096.0*t177-2187.0/4096.0*t104*t181;
  values[32] = 2187.0/4096.0*t7*t198+2187.0/4096.0*t200+2187.0/4096.0*t18*t202;
  values[33] = -6561.0/4096.0*t26*t198-6561.0/4096.0*t206-6561.0/4096.0*t28*t202;
  values[34] = -6561.0/4096.0*t37*t198-6561.0/4096.0*t210-6561.0/4096.0*t40*t202;
  values[35] = 2187.0/4096.0*t46*t198+2187.0/4096.0*t214+2187.0/4096.0*t48*t202;
  values[36] = -6561.0/4096.0*t54*t198-6561.0/4096.0*t200-6561.0/4096.0*t18*t219;
  values[37] = 19683.0/4096.0*t60*t198+19683.0/4096.0*t206+19683.0/4096.0*t28*t219;
  values[38] = 19683.0/4096.0*t65*t198+19683.0/4096.0*t210+19683.0/4096.0*t40*t219;
  values[39] = -6561.0/4096.0*t69*t198-6561.0/4096.0*t214-6561.0/4096.0*t48*t219;
  values[40] = 6561.0/4096.0*t75*t232+6561.0/4096.0*t234+6561.0/4096.0*t82*t236;
  values[41] = -19683.0/4096.0*t88*t232-19683.0/4096.0*t240-19683.0/4096.0*t90*t236;
  values[42] = -19683.0/4096.0*t95*t232-19683.0/4096.0*t244-19683.0/4096.0*t98*t236;
  values[43] = 6561.0/4096.0*t102*t232+6561.0/4096.0*t248+6561.0/4096.0*t104*t236;
  values[44] = -2187.0/4096.0*t110*t232-2187.0/4096.0*t234-2187.0/4096.0*t82*t253;
  values[45] = 6561.0/4096.0*t116*t232+6561.0/4096.0*t240+6561.0/4096.0*t90*t253;
  values[46] = 6561.0/4096.0*t121*t232+6561.0/4096.0*t244+6561.0/4096.0*t98*t253;
  values[47] = -2187.0/4096.0*t125*t232-2187.0/4096.0*t248-2187.0/4096.0*t104*t253;
  values[48] = -729.0/4096.0*t7*t267-729.0/4096.0*t269-729.0/4096.0*t18*t270;
  values[49] = 2187.0/4096.0*t26*t267+2187.0/4096.0*t274+2187.0/4096.0*t28*t270;
  values[50] = 2187.0/4096.0*t37*t267+2187.0/4096.0*t278+2187.0/4096.0*t40*t270;
  values[51] = -729.0/4096.0*t46*t267-729.0/4096.0*t282-729.0/4096.0*t48*t270;
  values[52] = 2187.0/4096.0*t54*t267+2187.0/4096.0*t269+2187.0/4096.0*t18*t286;
  values[53] = -6561.0/4096.0*t60*t267-6561.0/4096.0*t274-6561.0/4096.0*t28*t286;
  values[54] = -6561.0/4096.0*t65*t267-6561.0/4096.0*t278-6561.0/4096.0*t40*t286;
  values[55] = 2187.0/4096.0*t69*t267+2187.0/4096.0*t282+2187.0/4096.0*t48*t286;
  values[56] = -2187.0/4096.0*t75*t298-2187.0/4096.0*t300-2187.0/4096.0*t82*t301;
  values[57] = 6561.0/4096.0*t88*t298+6561.0/4096.0*t305+6561.0/4096.0*t90*t301;
  values[58] = 6561.0/4096.0*t95*t298+6561.0/4096.0*t309+6561.0/4096.0*t98*t301;
  values[59] = -2187.0/4096.0*t102*t298-2187.0/4096.0*t313-2187.0/4096.0*t104*t301;
  values[60] = 729.0/4096.0*t110*t298+729.0/4096.0*t300+729.0/4096.0*t82*t317;
  values[61] = -2187.0/4096.0*t116*t298-2187.0/4096.0*t305-2187.0/4096.0*t90*t317;
  values[62] = -2187.0/4096.0*t121*t298-2187.0/4096.0*t309-2187.0/4096.0*t98*t317;
  values[63] = 729.0/4096.0*t125*t298+729.0/4096.0*t313+729.0/4096.0*t104*t317;
}

static void C_H_Q3_3D_DeriveZeta(double xi, double eta, double zeta,
                             double *values)
{
  double t1, t3, t4, t5, t6, t7, t9, t10, t11, t12, t13, t14, t16, t17;
  double t18, t19, t20, t21, t25, t26, t28, t31, t33, t34, t35, t36, t38;
  double t42, t43, t45, t49, t50, t52, t53, t54, t58, t62, t66, t69, t70;
  double t71, t73, t74, t75, t77, t78, t79, t82, t84, t87, t88, t90, t93;
  double t95, t99, t100, t102, t103, t104, t108, t112, t116, t119, t120;
  double t121, t123, t124, t136, t138, t150, t152, t164, t166, t178, t179;
  double t180, t181, t183, t184, t185, t186, t187, t188, t192, t196, t200;
  double t203, t205, t206, t207, t211, t215, t219, t222, t224, t225, t226;
  double t230, t234, t238, t241, t243, t244, t245, t249, t253, t257, t260;
  double t261, t262, t264, t265, t277, t279, t291, t293, t305, t307;

  t1 = xi-1.0;
  t3 = t1*(xi+1.0/3.0);
  t4 = xi-1.0/3.0;
  t5 = eta-1.0;
  t6 = t4*t5;
  t7 = t3*t6;
  t9 = eta-1.0/3.0;
  t10 = (eta+1.0/3.0)*t9;
  t11 = zeta+1.0/3.0;
  t12 = zeta-1.0/3.0;
  t13 = t11*t12;
  t14 = t10*t13;
  t16 = zeta-1.0;
  t17 = t16*t12;
  t18 = t10*t17;
  t19 = t7*t18;
  t20 = t16*t11;
  t21 = t10*t20;
  t25 = t1*(xi+1.0);
  t26 = t25*t6;
  t28 = t26*t18;
  t31 = -xi-1.0;
  t33 = t31*(-xi+1.0);
  t34 = -xi-1.0/3.0;
  t35 = t34*t5;
  t36 = t33*t35;
  t38 = t36*t18;
  t42 = t31*(-xi+1.0/3.0);
  t43 = t42*t35;
  t45 = t43*t18;
  t49 = (eta+1.0)*t9;
  t50 = t49*t13;
  t52 = t49*t17;
  t53 = t7*t52;
  t54 = t49*t20;
  t58 = t26*t52;
  t62 = t36*t52;
  t66 = t43*t52;
  t69 = -eta-1.0;
  t70 = t4*t69;
  t71 = t3*t70;
  t73 = -eta-1.0/3.0;
  t74 = (-eta+1.0)*t73;
  t75 = t74*t13;
  t77 = t74*t17;
  t78 = t71*t77;
  t79 = t74*t20;
  t82 = t25*t70;
  t84 = t82*t77;
  t87 = t34*t69;
  t88 = t33*t87;
  t90 = t88*t77;
  t93 = t42*t87;
  t95 = t93*t77;
  t99 = (-eta+1.0/3.0)*t73;
  t100 = t99*t13;
  t102 = t99*t17;
  t103 = t71*t102;
  t104 = t99*t20;
  t108 = t82*t102;
  t112 = t88*t102;
  t116 = t93*t102;
  t119 = zeta+1.0;
  t120 = t119*t12;
  t121 = t10*t120;
  t123 = t16*t119;
  t124 = t10*t123;
  t136 = t49*t120;
  t138 = t49*t123;
  t150 = t74*t120;
  t152 = t74*t123;
  t164 = t99*t120;
  t166 = t99*t123;
  t178 = -zeta+1.0;
  t179 = -zeta-1.0/3.0;
  t180 = t178*t179;
  t181 = t10*t180;
  t183 = -zeta-1.0;
  t184 = t183*t179;
  t185 = t10*t184;
  t186 = t7*t185;
  t187 = t183*t178;
  t188 = t10*t187;
  t192 = t26*t185;
  t196 = t36*t185;
  t200 = t43*t185;
  t203 = t49*t180;
  t205 = t49*t184;
  t206 = t7*t205;
  t207 = t49*t187;
  t211 = t26*t205;
  t215 = t36*t205;
  t219 = t43*t205;
  t222 = t74*t180;
  t224 = t74*t184;
  t225 = t71*t224;
  t226 = t74*t187;
  t230 = t82*t224;
  t234 = t88*t224;
  t238 = t93*t224;
  t241 = t99*t180;
  t243 = t99*t184;
  t244 = t71*t243;
  t245 = t99*t187;
  t249 = t82*t243;
  t253 = t88*t243;
  t257 = t93*t243;
  t260 = -zeta+1.0/3.0;
  t261 = t260*t179;
  t262 = t10*t261;
  t264 = t183*t260;
  t265 = t10*t264;
  t277 = t49*t261;
  t279 = t49*t264;
  t291 = t74*t261;
  t293 = t74*t264;
  t305 = t99*t261;
  t307 = t99*t264;

  values[0] = -729.0/4096.0*t7*t14-729.0/4096.0*t19-729.0/4096.0*t7*t21;
  values[1] = 2187.0/4096.0*t26*t14+2187.0/4096.0*t28+2187.0/4096.0*t26*t21;
  values[2] = 2187.0/4096.0*t36*t14+2187.0/4096.0*t38+2187.0/4096.0*t36*t21;
  values[3] = -729.0/4096.0*t43*t14-729.0/4096.0*t45-729.0/4096.0*t43*t21;
  values[4] = 2187.0/4096.0*t7*t50+2187.0/4096.0*t53+2187.0/4096.0*t7*t54;
  values[5] = -6561.0/4096.0*t26*t50-6561.0/4096.0*t58-6561.0/4096.0*t26*t54;
  values[6] = -6561.0/4096.0*t36*t50-6561.0/4096.0*t62-6561.0/4096.0*t36*t54;
  values[7] = 2187.0/4096.0*t43*t50+2187.0/4096.0*t66+2187.0/4096.0*t43*t54;
  values[8] = 2187.0/4096.0*t71*t75+2187.0/4096.0*t78+2187.0/4096.0*t71*t79;
  values[9] = -6561.0/4096.0*t82*t75-6561.0/4096.0*t84-6561.0/4096.0*t82*t79;
  values[10] = -6561.0/4096.0*t88*t75-6561.0/4096.0*t90-6561.0/4096.0*t88*t79;
  values[11] = 2187.0/4096.0*t93*t75+2187.0/4096.0*t95+2187.0/4096.0*t93*t79;
  values[12] = -729.0/4096.0*t71*t100-729.0/4096.0*t103-729.0/4096.0*t71*t104;
  values[13] = 2187.0/4096.0*t82*t100+2187.0/4096.0*t108+2187.0/4096.0*t82*t104;
  values[14] = 2187.0/4096.0*t88*t100+2187.0/4096.0*t112+2187.0/4096.0*t88*t104;
  values[15] = -729.0/4096.0*t93*t100-729.0/4096.0*t116-729.0/4096.0*t93*t104;
  values[16] = 2187.0/4096.0*t7*t121+2187.0/4096.0*t19+2187.0/4096.0*t7*t124;
  values[17] = -6561.0/4096.0*t26*t121-6561.0/4096.0*t28-6561.0/4096.0*t26*t124;
  values[18] = -6561.0/4096.0*t36*t121-6561.0/4096.0*t38-6561.0/4096.0*t36*t124;
  values[19] = 2187.0/4096.0*t43*t121+2187.0/4096.0*t45+2187.0/4096.0*t43*t124;
  values[20] = -6561.0/4096.0*t7*t136-6561.0/4096.0*t53-6561.0/4096.0*t7*t138;
  values[21] = 19683.0/4096.0*t26*t136+19683.0/4096.0*t58+19683.0/4096.0*t26*t138;
  values[22] = 19683.0/4096.0*t36*t136+19683.0/4096.0*t62+19683.0/4096.0*t36*t138;
  values[23] = -6561.0/4096.0*t43*t136-6561.0/4096.0*t66-6561.0/4096.0*t43*t138;
  values[24] = -6561.0/4096.0*t71*t150-6561.0/4096.0*t78-6561.0/4096.0*t71*t152;
  values[25] = 19683.0/4096.0*t82*t150+19683.0/4096.0*t84+19683.0/4096.0*t82*t152;
  values[26] = 19683.0/4096.0*t88*t150+19683.0/4096.0*t90+19683.0/4096.0*t88*t152;
  values[27] = -6561.0/4096.0*t93*t150-6561.0/4096.0*t95-6561.0/4096.0*t93*t152;
  values[28] = 2187.0/4096.0*t71*t164+2187.0/4096.0*t103+2187.0/4096.0*t71*t166;
  values[29] = -6561.0/4096.0*t82*t164-6561.0/4096.0*t108-6561.0/4096.0*t82*t166;
  values[30] = -6561.0/4096.0*t88*t164-6561.0/4096.0*t112-6561.0/4096.0*t88*t166;
  values[31] = 2187.0/4096.0*t93*t164+2187.0/4096.0*t116+2187.0/4096.0*t93*t166;
  values[32] = -2187.0/4096.0*t7*t181-2187.0/4096.0*t186-2187.0/4096.0*t7*t188;
  values[33] = 6561.0/4096.0*t26*t181+6561.0/4096.0*t192+6561.0/4096.0*t26*t188;
  values[34] = 6561.0/4096.0*t36*t181+6561.0/4096.0*t196+6561.0/4096.0*t36*t188;
  values[35] = -2187.0/4096.0*t43*t181-2187.0/4096.0*t200-2187.0/4096.0*t43*t188;
  values[36] = 6561.0/4096.0*t7*t203+6561.0/4096.0*t206+6561.0/4096.0*t7*t207;
  values[37] = -19683.0/4096.0*t26*t203-19683.0/4096.0*t211-19683.0/4096.0*t26*t207;
  values[38] = -19683.0/4096.0*t36*t203-19683.0/4096.0*t215-19683.0/4096.0*t36*t207;
  values[39] = 6561.0/4096.0*t43*t203+6561.0/4096.0*t219+6561.0/4096.0*t43*t207;
  values[40] = 6561.0/4096.0*t71*t222+6561.0/4096.0*t225+6561.0/4096.0*t71*t226;
  values[41] = -19683.0/4096.0*t82*t222-19683.0/4096.0*t230-19683.0/4096.0*t82*t226;
  values[42] = -19683.0/4096.0*t88*t222-19683.0/4096.0*t234-19683.0/4096.0*t88*t226;
  values[43] = 6561.0/4096.0*t93*t222+6561.0/4096.0*t238+6561.0/4096.0*t93*t226;
  values[44] = -2187.0/4096.0*t71*t241-2187.0/4096.0*t244-2187.0/4096.0*t71*t245;
  values[45] = 6561.0/4096.0*t82*t241+6561.0/4096.0*t249+6561.0/4096.0*t82*t245;
  values[46] = 6561.0/4096.0*t88*t241+6561.0/4096.0*t253+6561.0/4096.0*t88*t245;
  values[47] = -2187.0/4096.0*t93*t241-2187.0/4096.0*t257-2187.0/4096.0*t93*t245;
  values[48] = 729.0/4096.0*t7*t262+729.0/4096.0*t186+729.0/4096.0*t7*t265;
  values[49] = -2187.0/4096.0*t26*t262-2187.0/4096.0*t192-2187.0/4096.0*t26*t265;
  values[50] = -2187.0/4096.0*t36*t262-2187.0/4096.0*t196-2187.0/4096.0*t36*t265;
  values[51] = 729.0/4096.0*t43*t262+729.0/4096.0*t200+729.0/4096.0*t43*t265;
  values[52] = -2187.0/4096.0*t7*t277-2187.0/4096.0*t206-2187.0/4096.0*t7*t279;
  values[53] = 6561.0/4096.0*t26*t277+6561.0/4096.0*t211+6561.0/4096.0*t26*t279;
  values[54] = 6561.0/4096.0*t36*t277+6561.0/4096.0*t215+6561.0/4096.0*t36*t279;
  values[55] = -2187.0/4096.0*t43*t277-2187.0/4096.0*t219-2187.0/4096.0*t43*t279;
  values[56] = -2187.0/4096.0*t71*t291-2187.0/4096.0*t225-2187.0/4096.0*t71*t293;
  values[57] = 6561.0/4096.0*t82*t291+6561.0/4096.0*t230+6561.0/4096.0*t82*t293;
  values[58] = 6561.0/4096.0*t88*t291+6561.0/4096.0*t234+6561.0/4096.0*t88*t293;
  values[59] = -2187.0/4096.0*t93*t291-2187.0/4096.0*t238-2187.0/4096.0*t93*t293;
  values[60] = 729.0/4096.0*t71*t305+729.0/4096.0*t244+729.0/4096.0*t71*t307;
  values[61] = -2187.0/4096.0*t82*t305-2187.0/4096.0*t249-2187.0/4096.0*t82*t307;
  values[62] = -2187.0/4096.0*t88*t305-2187.0/4096.0*t253-2187.0/4096.0*t88*t307;
  values[63] = 729.0/4096.0*t93*t305+729.0/4096.0*t257+729.0/4096.0*t93*t307;
}

static void C_H_Q3_3D_DeriveXiXi(double xi, double eta, double zeta,
                             double *values)
{
  double t1, t2, t3, t4, t5, t6, t7, t8, t10, t11, t12, t13, t14, t15;
  double t16, t18, t19, t20, t21, t23, t24, t25, t28, t29, t30, t31, t32;
  double t33, t34, t36, t37, t38, t39, t41, t42, t43, t46, t47, t48, t49;
  double t51, t52, t54, t57, t58, t59, t61, t62, t64, t67, t68, t69, t70;
  double t71, t72, t73, t74, t75, t76, t78, t79, t80, t82, t83, t86, t87;
  double t88, t89, t90, t92, t93, t94, t96, t97, t100, t101, t102, t103;
  double t105, t106, t108, t111, t112, t113, t115, t116, t118, t122, t123;
  double t124, t126, t130, t132, t136, t138, t142, t144, t148, t149, t151;
  double t155, t157, t161, t163, t167, t169, t173, t174, t176, t177, t178;
  double t179, t181, t185, t187, t191, t193, t197, t199, t203, t204, t205;
  double t207, t211, t213, t217, t219, t223, t225, t230, t231, t232, t234;
  double t238, t240, t244, t246, t250, t252, t256, t257, t259, t263, t265;
  double t269, t271, t275, t277;

  t1 = xi-1.0/3.0;
  t2 = eta-1.0;
  t3 = t1*t2;
  t4 = eta+1.0/3.0;
  t5 = t3*t4;
  t6 = eta-1.0/3.0;
  t7 = zeta-1.0;
  t8 = t6*t7;
  t10 = zeta-1.0/3.0;
  t11 = (zeta+1.0/3.0)*t10;
  t12 = t8*t11;
  t13 = t5*t12;
  t14 = xi+1.0/3.0;
  t15 = t14*t2;
  t16 = t15*t4;
  t18 = xi-1.0;
  t19 = t18*t2;
  t20 = t19*t4;
  t21 = t20*t12;
  t23 = xi+1.0;
  t24 = t23*t2;
  t25 = t24*t4;
  t28 = -xi-1.0/3.0;
  t29 = t28*t2;
  t30 = t29*t4;
  t31 = t30*t12;
  t32 = -xi+1.0;
  t33 = t32*t2;
  t34 = t33*t4;
  t36 = -xi-1.0;
  t37 = t36*t2;
  t38 = t37*t4;
  t39 = t38*t12;
  t41 = -xi+1.0/3.0;
  t42 = t41*t2;
  t43 = t42*t4;
  t46 = eta+1.0;
  t47 = t3*t46;
  t48 = t47*t12;
  t49 = t15*t46;
  t51 = t19*t46;
  t52 = t51*t12;
  t54 = t24*t46;
  t57 = t29*t46;
  t58 = t57*t12;
  t59 = t33*t46;
  t61 = t37*t46;
  t62 = t61*t12;
  t64 = t42*t46;
  t67 = -eta-1.0;
  t68 = t1*t67;
  t69 = -eta+1.0;
  t70 = t68*t69;
  t71 = -eta-1.0/3.0;
  t72 = t71*t7;
  t73 = t72*t11;
  t74 = t70*t73;
  t75 = t14*t67;
  t76 = t75*t69;
  t78 = t18*t67;
  t79 = t78*t69;
  t80 = t79*t73;
  t82 = t23*t67;
  t83 = t82*t69;
  t86 = t28*t67;
  t87 = t86*t69;
  t88 = t87*t73;
  t89 = t32*t67;
  t90 = t89*t69;
  t92 = t36*t67;
  t93 = t92*t69;
  t94 = t93*t73;
  t96 = t41*t67;
  t97 = t96*t69;
  t100 = -eta+1.0/3.0;
  t101 = t68*t100;
  t102 = t101*t73;
  t103 = t75*t100;
  t105 = t78*t100;
  t106 = t105*t73;
  t108 = t82*t100;
  t111 = t86*t100;
  t112 = t111*t73;
  t113 = t89*t100;
  t115 = t92*t100;
  t116 = t115*t73;
  t118 = t96*t100;
  t122 = (zeta+1.0)*t10;
  t123 = t8*t122;
  t124 = t5*t123;
  t126 = t20*t123;
  t130 = t30*t123;
  t132 = t38*t123;
  t136 = t47*t123;
  t138 = t51*t123;
  t142 = t57*t123;
  t144 = t61*t123;
  t148 = t72*t122;
  t149 = t70*t148;
  t151 = t79*t148;
  t155 = t87*t148;
  t157 = t93*t148;
  t161 = t101*t148;
  t163 = t105*t148;
  t167 = t111*t148;
  t169 = t115*t148;
  t173 = -zeta-1.0;
  t174 = t6*t173;
  t176 = -zeta-1.0/3.0;
  t177 = (-zeta+1.0)*t176;
  t178 = t174*t177;
  t179 = t5*t178;
  t181 = t20*t178;
  t185 = t30*t178;
  t187 = t38*t178;
  t191 = t47*t178;
  t193 = t51*t178;
  t197 = t57*t178;
  t199 = t61*t178;
  t203 = t71*t173;
  t204 = t203*t177;
  t205 = t70*t204;
  t207 = t79*t204;
  t211 = t87*t204;
  t213 = t93*t204;
  t217 = t101*t204;
  t219 = t105*t204;
  t223 = t111*t204;
  t225 = t115*t204;
  t230 = (-zeta+1.0/3.0)*t176;
  t231 = t174*t230;
  t232 = t5*t231;
  t234 = t20*t231;
  t238 = t30*t231;
  t240 = t38*t231;
  t244 = t47*t231;
  t246 = t51*t231;
  t250 = t57*t231;
  t252 = t61*t231;
  t256 = t203*t230;
  t257 = t70*t256;
  t259 = t79*t256;
  t263 = t87*t256;
  t265 = t93*t256;
  t269 = t101*t256;
  t271 = t105*t256;
  t275 = t111*t256;
  t277 = t115*t256;

  values[0] = -729.0/2048.0*t13-729.0/2048.0*t16*t12-729.0/2048.0*t21;
  values[1] = 2187.0/2048.0*t13+2187.0/2048.0*t25*t12+2187.0/2048.0*t21;
  values[2] = 2187.0/2048.0*t31+2187.0/2048.0*t34*t12+2187.0/2048.0*t39;
  values[3] = -729.0/2048.0*t31-729.0/2048.0*t43*t12-729.0/2048.0*t39;
  values[4] = 2187.0/2048.0*t48+2187.0/2048.0*t49*t12+2187.0/2048.0*t52;
  values[5] = -6561.0/2048.0*t48-6561.0/2048.0*t54*t12-6561.0/2048.0*t52;
  values[6] = -6561.0/2048.0*t58-6561.0/2048.0*t59*t12-6561.0/2048.0*t62;
  values[7] = 2187.0/2048.0*t58+2187.0/2048.0*t64*t12+2187.0/2048.0*t62;
  values[8] = 2187.0/2048.0*t74+2187.0/2048.0*t76*t73+2187.0/2048.0*t80;
  values[9] = -6561.0/2048.0*t74-6561.0/2048.0*t83*t73-6561.0/2048.0*t80;
  values[10] = -6561.0/2048.0*t88-6561.0/2048.0*t90*t73-6561.0/2048.0*t94;
  values[11] = 2187.0/2048.0*t88+2187.0/2048.0*t97*t73+2187.0/2048.0*t94;
  values[12] = -729.0/2048.0*t102-729.0/2048.0*t103*t73-729.0/2048.0*t106;
  values[13] = 2187.0/2048.0*t102+2187.0/2048.0*t108*t73+2187.0/2048.0*t106;
  values[14] = 2187.0/2048.0*t112+2187.0/2048.0*t113*t73+2187.0/2048.0*t116;
  values[15] = -729.0/2048.0*t112-729.0/2048.0*t118*t73-729.0/2048.0*t116;
  values[16] = 2187.0/2048.0*t124+2187.0/2048.0*t16*t123+2187.0/2048.0*t126;
  values[17] = -6561.0/2048.0*t124-6561.0/2048.0*t25*t123-6561.0/2048.0*t126;
  values[18] = -6561.0/2048.0*t130-6561.0/2048.0*t34*t123-6561.0/2048.0*t132;
  values[19] = 2187.0/2048.0*t130+2187.0/2048.0*t43*t123+2187.0/2048.0*t132;
  values[20] = -6561.0/2048.0*t136-6561.0/2048.0*t49*t123-6561.0/2048.0*t138;
  values[21] = 19683.0/2048.0*t136+19683.0/2048.0*t54*t123+19683.0/2048.0*t138;
  values[22] = 19683.0/2048.0*t142+19683.0/2048.0*t59*t123+19683.0/2048.0*t144;
  values[23] = -6561.0/2048.0*t142-6561.0/2048.0*t64*t123-6561.0/2048.0*t144;
  values[24] = -6561.0/2048.0*t149-6561.0/2048.0*t76*t148-6561.0/2048.0*t151;
  values[25] = 19683.0/2048.0*t149+19683.0/2048.0*t83*t148+19683.0/2048.0*t151;
  values[26] = 19683.0/2048.0*t155+19683.0/2048.0*t90*t148+19683.0/2048.0*t157;
  values[27] = -6561.0/2048.0*t155-6561.0/2048.0*t97*t148-6561.0/2048.0*t157;
  values[28] = 2187.0/2048.0*t161+2187.0/2048.0*t103*t148+2187.0/2048.0*t163;
  values[29] = -6561.0/2048.0*t161-6561.0/2048.0*t108*t148-6561.0/2048.0*t163;
  values[30] = -6561.0/2048.0*t167-6561.0/2048.0*t113*t148-6561.0/2048.0*t169;
  values[31] = 2187.0/2048.0*t167+2187.0/2048.0*t118*t148+2187.0/2048.0*t169;
  values[32] = 2187.0/2048.0*t179+2187.0/2048.0*t16*t178+2187.0/2048.0*t181;
  values[33] = -6561.0/2048.0*t179-6561.0/2048.0*t25*t178-6561.0/2048.0*t181;
  values[34] = -6561.0/2048.0*t185-6561.0/2048.0*t34*t178-6561.0/2048.0*t187;
  values[35] = 2187.0/2048.0*t185+2187.0/2048.0*t43*t178+2187.0/2048.0*t187;
  values[36] = -6561.0/2048.0*t191-6561.0/2048.0*t49*t178-6561.0/2048.0*t193;
  values[37] = 19683.0/2048.0*t191+19683.0/2048.0*t54*t178+19683.0/2048.0*t193;
  values[38] = 19683.0/2048.0*t197+19683.0/2048.0*t59*t178+19683.0/2048.0*t199;
  values[39] = -6561.0/2048.0*t197-6561.0/2048.0*t64*t178-6561.0/2048.0*t199;
  values[40] = -6561.0/2048.0*t205-6561.0/2048.0*t76*t204-6561.0/2048.0*t207;
  values[41] = 19683.0/2048.0*t205+19683.0/2048.0*t83*t204+19683.0/2048.0*t207;
  values[42] = 19683.0/2048.0*t211+19683.0/2048.0*t90*t204+19683.0/2048.0*t213;
  values[43] = -6561.0/2048.0*t211-6561.0/2048.0*t97*t204-6561.0/2048.0*t213;
  values[44] = 2187.0/2048.0*t217+2187.0/2048.0*t103*t204+2187.0/2048.0*t219;
  values[45] = -6561.0/2048.0*t217-6561.0/2048.0*t108*t204-6561.0/2048.0*t219;
  values[46] = -6561.0/2048.0*t223-6561.0/2048.0*t113*t204-6561.0/2048.0*t225;
  values[47] = 2187.0/2048.0*t223+2187.0/2048.0*t118*t204+2187.0/2048.0*t225;
  values[48] = -729.0/2048.0*t232-729.0/2048.0*t16*t231-729.0/2048.0*t234;
  values[49] = 2187.0/2048.0*t232+2187.0/2048.0*t25*t231+2187.0/2048.0*t234;
  values[50] = 2187.0/2048.0*t238+2187.0/2048.0*t34*t231+2187.0/2048.0*t240;
  values[51] = -729.0/2048.0*t238-729.0/2048.0*t43*t231-729.0/2048.0*t240;
  values[52] = 2187.0/2048.0*t244+2187.0/2048.0*t49*t231+2187.0/2048.0*t246;
  values[53] = -6561.0/2048.0*t244-6561.0/2048.0*t54*t231-6561.0/2048.0*t246;
  values[54] = -6561.0/2048.0*t250-6561.0/2048.0*t59*t231-6561.0/2048.0*t252;
  values[55] = 2187.0/2048.0*t250+2187.0/2048.0*t64*t231+2187.0/2048.0*t252;
  values[56] = 2187.0/2048.0*t257+2187.0/2048.0*t76*t256+2187.0/2048.0*t259;
  values[57] = -6561.0/2048.0*t257-6561.0/2048.0*t83*t256-6561.0/2048.0*t259;
  values[58] = -6561.0/2048.0*t263-6561.0/2048.0*t90*t256-6561.0/2048.0*t265;
  values[59] = 2187.0/2048.0*t263+2187.0/2048.0*t97*t256+2187.0/2048.0*t265;
  values[60] = -729.0/2048.0*t269-729.0/2048.0*t103*t256-729.0/2048.0*t271;
  values[61] = 2187.0/2048.0*t269+2187.0/2048.0*t108*t256+2187.0/2048.0*t271;
  values[62] = 2187.0/2048.0*t275+2187.0/2048.0*t113*t256+2187.0/2048.0*t277;
  values[63] = -729.0/2048.0*t275-729.0/2048.0*t118*t256-729.0/2048.0*t277;
}

static void C_H_Q3_3D_DeriveXiEta(double xi, double eta, double zeta,
                             double *values)
{
  double t1, t2, t3, t4, t5, t6, t7, t8, t10, t11, t12, t14, t15, t16;
  double t17, t18, t20, t21, t22, t23, t24, t25, t26, t27, t28, t30;
  double t31, t34, t35, t36, t38, t39, t41, t42, t44, t45, t48, t49;
  double t50, t51, t53, t54, t56, t57, t58, t59, t60, t61, t62, t63;
  double t64, t66, t67, t70, t71, t72, t74, t75, t77, t78, t80, t81;
  double t84, t85, t87, t88, t90, t91, t92, t93, t97, t100, t104, t107;
  double t108, t109, t110, t114, t117, t121, t122, t123, t124, t125;
  double t127, t128, t129, t130, t131, t133, t134, t135, t136, t137;
  double t138, t140, t141, t144, t146, t147, t149, t151, t152, t155;
  double t157, t158, t160, t161, t162, t163, t164, t165, t167, t168;
  double t171, t173, t174, t176, t178, t179, t182, t183, t185, t186;
  double t188, t189, t190, t191, t195, t198, t202, t205, t206, t207;
  double t208, t212, t215, t220, t221, t223, t224, t226, t227, t228;
  double t230, t234, t237, t241, t243, t244, t245, t247, t251, t254;
  double t258, t260, t261, t272, t273, t282, t284, t285, t287, t288;
  double t289, t291, t295, t298, t302, t304, t305, t306, t308, t312;
  double t315, t319, t321, t322, t333, t334, t343, t344, t346, t347;
  double t348, t350, t351, t352, t354, t355, t356, t358, t362, t365;
  double t369, t371, t372, t373, t375, t379, t382, t386, t387, t389;
  double t390, t401, t402, t411, t412, t414, t415, t416, t418, t419;
  double t420, t422, t426, t429, t433, t435, t436, t437, t439, t443;
  double t446, t450, t451, t453, t454, t465, t466, t476, t477, t479;
  double t480, t482, t483, t484, t486, t490, t493, t497, t499, t500;
  double t501, t503, t507, t510, t514, t516, t517, t528, t529, t538;
  double t540, t541, t543, t544, t545, t547, t551, t554, t558, t560;
  double t561, t562, t564, t568, t571, t575, t577, t578, t589, t590;

  t1 = xi+1.0/3.0;
  t2 = xi-1.0/3.0;
  t3 = t1*t2;
  t4 = eta+1.0/3.0;
  t5 = t3*t4;
  t6 = eta-1.0/3.0;
  t7 = zeta-1.0;
  t8 = t6*t7;
  t10 = zeta-1.0/3.0;
  t11 = (zeta+1.0/3.0)*t10;
  t12 = t8*t11;
  t14 = eta-1.0;
  t15 = t3*t14;
  t16 = t15*t12;
  t17 = t4*t7;
  t18 = t17*t11;
  t20 = xi-1.0;
  t21 = t20*t2;
  t22 = t21*t4;
  t23 = t22*t12;
  t24 = t21*t14;
  t25 = t24*t12;
  t26 = t24*t18;
  t27 = t20*t1;
  t28 = t27*t4;
  t30 = t27*t14;
  t31 = t30*t12;
  t34 = xi+1.0;
  t35 = t34*t2;
  t36 = t35*t4;
  t38 = t35*t14;
  t39 = t38*t12;
  t41 = t20*t34;
  t42 = t41*t4;
  t44 = t41*t14;
  t45 = t44*t12;
  t48 = -xi+1.0;
  t49 = -xi-1.0/3.0;
  t50 = t48*t49;
  t51 = t50*t4;
  t53 = t50*t14;
  t54 = t53*t12;
  t56 = -xi-1.0;
  t57 = t56*t49;
  t58 = t57*t4;
  t59 = t58*t12;
  t60 = t57*t14;
  t61 = t60*t12;
  t62 = t60*t18;
  t63 = t56*t48;
  t64 = t63*t4;
  t66 = t63*t14;
  t67 = t66*t12;
  t70 = -xi+1.0/3.0;
  t71 = t70*t49;
  t72 = t71*t4;
  t74 = t71*t14;
  t75 = t74*t12;
  t77 = t56*t70;
  t78 = t77*t4;
  t80 = t77*t14;
  t81 = t80*t12;
  t84 = eta+1.0;
  t85 = t3*t84;
  t87 = t84*t7;
  t88 = t87*t11;
  t90 = t21*t84;
  t91 = t90*t12;
  t92 = t24*t88;
  t93 = t27*t84;
  t97 = t35*t84;
  t100 = t41*t84;
  t104 = t50*t84;
  t107 = t57*t84;
  t108 = t107*t12;
  t109 = t60*t88;
  t110 = t63*t84;
  t114 = t71*t84;
  t117 = t77*t84;
  t121 = -eta+1.0;
  t122 = t3*t121;
  t123 = -eta-1.0/3.0;
  t124 = t123*t7;
  t125 = t124*t11;
  t127 = -eta-1.0;
  t128 = t3*t127;
  t129 = t128*t125;
  t130 = t121*t7;
  t131 = t130*t11;
  t133 = t21*t121;
  t134 = t133*t125;
  t135 = t21*t127;
  t136 = t135*t125;
  t137 = t135*t131;
  t138 = t27*t121;
  t140 = t27*t127;
  t141 = t140*t125;
  t144 = t35*t121;
  t146 = t35*t127;
  t147 = t146*t125;
  t149 = t41*t121;
  t151 = t41*t127;
  t152 = t151*t125;
  t155 = t50*t121;
  t157 = t50*t127;
  t158 = t157*t125;
  t160 = t57*t121;
  t161 = t160*t125;
  t162 = t57*t127;
  t163 = t162*t125;
  t164 = t162*t131;
  t165 = t63*t121;
  t167 = t63*t127;
  t168 = t167*t125;
  t171 = t71*t121;
  t173 = t71*t127;
  t174 = t173*t125;
  t176 = t77*t121;
  t178 = t77*t127;
  t179 = t178*t125;
  t182 = -eta+1.0/3.0;
  t183 = t3*t182;
  t185 = t182*t7;
  t186 = t185*t11;
  t188 = t21*t182;
  t189 = t188*t125;
  t190 = t135*t186;
  t191 = t27*t182;
  t195 = t35*t182;
  t198 = t41*t182;
  t202 = t50*t182;
  t205 = t57*t182;
  t206 = t205*t125;
  t207 = t162*t186;
  t208 = t63*t182;
  t212 = t71*t182;
  t215 = t77*t182;
  t220 = (zeta+1.0)*t10;
  t221 = t8*t220;
  t223 = t15*t221;
  t224 = t17*t220;
  t226 = t22*t221;
  t227 = t24*t221;
  t228 = t24*t224;
  t230 = t30*t221;
  t234 = t38*t221;
  t237 = t44*t221;
  t241 = t53*t221;
  t243 = t58*t221;
  t244 = t60*t221;
  t245 = t60*t224;
  t247 = t66*t221;
  t251 = t74*t221;
  t254 = t80*t221;
  t258 = t87*t220;
  t260 = t90*t221;
  t261 = t24*t258;
  t272 = t107*t221;
  t273 = t60*t258;
  t282 = t124*t220;
  t284 = t128*t282;
  t285 = t130*t220;
  t287 = t133*t282;
  t288 = t135*t282;
  t289 = t135*t285;
  t291 = t140*t282;
  t295 = t146*t282;
  t298 = t151*t282;
  t302 = t157*t282;
  t304 = t160*t282;
  t305 = t162*t282;
  t306 = t162*t285;
  t308 = t167*t282;
  t312 = t173*t282;
  t315 = t178*t282;
  t319 = t185*t220;
  t321 = t188*t282;
  t322 = t135*t319;
  t333 = t205*t282;
  t334 = t162*t319;
  t343 = -zeta-1.0;
  t344 = t6*t343;
  t346 = -zeta-1.0/3.0;
  t347 = (-zeta+1.0)*t346;
  t348 = t344*t347;
  t350 = t15*t348;
  t351 = t4*t343;
  t352 = t351*t347;
  t354 = t22*t348;
  t355 = t24*t348;
  t356 = t24*t352;
  t358 = t30*t348;
  t362 = t38*t348;
  t365 = t44*t348;
  t369 = t53*t348;
  t371 = t58*t348;
  t372 = t60*t348;
  t373 = t60*t352;
  t375 = t66*t348;
  t379 = t74*t348;
  t382 = t80*t348;
  t386 = t84*t343;
  t387 = t386*t347;
  t389 = t90*t348;
  t390 = t24*t387;
  t401 = t107*t348;
  t402 = t60*t387;
  t411 = t123*t343;
  t412 = t411*t347;
  t414 = t128*t412;
  t415 = t121*t343;
  t416 = t415*t347;
  t418 = t133*t412;
  t419 = t135*t412;
  t420 = t135*t416;
  t422 = t140*t412;
  t426 = t146*t412;
  t429 = t151*t412;
  t433 = t157*t412;
  t435 = t160*t412;
  t436 = t162*t412;
  t437 = t162*t416;
  t439 = t167*t412;
  t443 = t173*t412;
  t446 = t178*t412;
  t450 = t182*t343;
  t451 = t450*t347;
  t453 = t188*t412;
  t454 = t135*t451;
  t465 = t205*t412;
  t466 = t162*t451;
  t476 = (-zeta+1.0/3.0)*t346;
  t477 = t344*t476;
  t479 = t15*t477;
  t480 = t351*t476;
  t482 = t22*t477;
  t483 = t24*t477;
  t484 = t24*t480;
  t486 = t30*t477;
  t490 = t38*t477;
  t493 = t44*t477;
  t497 = t53*t477;
  t499 = t58*t477;
  t500 = t60*t477;
  t501 = t60*t480;
  t503 = t66*t477;
  t507 = t74*t477;
  t510 = t80*t477;
  t514 = t386*t476;
  t516 = t90*t477;
  t517 = t24*t514;
  t528 = t107*t477;
  t529 = t60*t514;
  t538 = t411*t476;
  t540 = t128*t538;
  t541 = t415*t476;
  t543 = t133*t538;
  t544 = t135*t538;
  t545 = t135*t541;
  t547 = t140*t538;
  t551 = t146*t538;
  t554 = t151*t538;
  t558 = t157*t538;
  t560 = t160*t538;
  t561 = t162*t538;
  t562 = t162*t541;
  t564 = t167*t538;
  t568 = t173*t538;
  t571 = t178*t538;
  t575 = t450*t476;
  t577 = t188*t538;
  t578 = t135*t575;
  t589 = t205*t538;
  t590 = t162*t575;

  values[0] = -729.0/4096.0*t5*t12-729.0/4096.0*t16-729.0/4096.0*t15*t18-729.0/4096.0*t23-729.0/4096.0*t25-729.0/4096.0*t26-729.0/4096.0*t28*t12-729.0/4096.0*t31-729.0/4096.0*t30*t18;
  values[1] = 2187.0/4096.0*t36*t12+2187.0/4096.0*t39+2187.0/4096.0*t38*t18+2187.0/4096.0*t23+2187.0/4096.0*t25+2187.0/4096.0*t26+2187.0/4096.0*t42*t12+2187.0/4096.0*t45+2187.0/4096.0*t44*t18;
  values[2] = -2187.0/4096.0*t51*t12-2187.0/4096.0*t54-2187.0/4096.0*t53*t18-2187.0/4096.0*t59-2187.0/4096.0*t61-2187.0/4096.0*t62-2187.0/4096.0*t64*t12-2187.0/4096.0*t67-2187.0/4096.0*t66*t18;
  values[3] = 729.0/4096.0*t72*t12+729.0/4096.0*t75+729.0/4096.0*t74*t18+729.0/4096.0*t59+729.0/4096.0*t61+729.0/4096.0*t62+729.0/4096.0*t78*t12+729.0/4096.0*t81+729.0/4096.0*t80*t18;
  values[4] = 2187.0/4096.0*t85*t12+2187.0/4096.0*t16+2187.0/4096.0*t15*t88+2187.0/4096.0*t91+2187.0/4096.0*t25+2187.0/4096.0*t92+2187.0/4096.0*t93*t12+2187.0/4096.0*t31+2187.0/4096.0*t30*t88;
  values[5] = -6561.0/4096.0*t97*t12-6561.0/4096.0*t39-6561.0/4096.0*t38*t88-6561.0/4096.0*t91-6561.0/4096.0*t25-6561.0/4096.0*t92-6561.0/4096.0*t100*t12-6561.0/4096.0*t45-6561.0/4096.0*t44*t88;
  values[6] = 6561.0/4096.0*t104*t12+6561.0/4096.0*t54+6561.0/4096.0*t53*t88+6561.0/4096.0*t108+6561.0/4096.0*t61+6561.0/4096.0*t109+6561.0/4096.0*t110*t12+6561.0/4096.0*t67+6561.0/4096.0*t66*t88;
  values[7] = -2187.0/4096.0*t114*t12-2187.0/4096.0*t75-2187.0/4096.0*t74*t88-2187.0/4096.0*t108-2187.0/4096.0*t61-2187.0/4096.0*t109-2187.0/4096.0*t117*t12-2187.0/4096.0*t81-2187.0/4096.0*t80*t88;
  values[8] = -2187.0/4096.0*t122*t125-2187.0/4096.0*t129-2187.0/4096.0*t128*t131-2187.0/4096.0*t134-2187.0/4096.0*t136-2187.0/4096.0*t137-2187.0/4096.0*t138*t125-2187.0/4096.0*t141-2187.0/4096.0*t140*t131;
  values[9] = 6561.0/4096.0*t144*t125+6561.0/4096.0*t147+6561.0/4096.0*t146*t131+6561.0/4096.0*t134+6561.0/4096.0*t136+6561.0/4096.0*t137+6561.0/4096.0*t149*t125+6561.0/4096.0*t152+6561.0/4096.0*t151*t131;
  values[10] = -6561.0/4096.0*t155*t125-6561.0/4096.0*t158-6561.0/4096.0*t157*t131-6561.0/4096.0*t161-6561.0/4096.0*t163-6561.0/4096.0*t164-6561.0/4096.0*t165*t125-6561.0/4096.0*t168-6561.0/4096.0*t167*t131;
  values[11] = 2187.0/4096.0*t171*t125+2187.0/4096.0*t174+2187.0/4096.0*t173*t131+2187.0/4096.0*t161+2187.0/4096.0*t163+2187.0/4096.0*t164+2187.0/4096.0*t176*t125+2187.0/4096.0*t179+2187.0/4096.0*t178*t131;
  values[12] = 729.0/4096.0*t183*t125+729.0/4096.0*t129+729.0/4096.0*t128*t186+729.0/4096.0*t189+729.0/4096.0*t136+729.0/4096.0*t190+729.0/4096.0*t191*t125+729.0/4096.0*t141+729.0/4096.0*t140*t186;
  values[13] = -2187.0/4096.0*t195*t125-2187.0/4096.0*t147-2187.0/4096.0*t146*t186-2187.0/4096.0*t189-2187.0/4096.0*t136-2187.0/4096.0*t190-2187.0/4096.0*t198*t125-2187.0/4096.0*t152-2187.0/4096.0*t151*t186;
  values[14] = 2187.0/4096.0*t202*t125+2187.0/4096.0*t158+2187.0/4096.0*t157*t186+2187.0/4096.0*t206+2187.0/4096.0*t163+2187.0/4096.0*t207+2187.0/4096.0*t208*t125+2187.0/4096.0*t168+2187.0/4096.0*t167*t186;
  values[15] = -729.0/4096.0*t212*t125-729.0/4096.0*t174-729.0/4096.0*t173*t186-729.0/4096.0*t206-729.0/4096.0*t163-729.0/4096.0*t207-729.0/4096.0*t215*t125-729.0/4096.0*t179-729.0/4096.0*t178*t186;
  values[16] = 2187.0/4096.0*t5*t221+2187.0/4096.0*t223+2187.0/4096.0*t15*t224+2187.0/4096.0*t226+2187.0/4096.0*t227+2187.0/4096.0*t228+2187.0/4096.0*t28*t221+2187.0/4096.0*t230+2187.0/4096.0*t30*t224;
  values[17] = -6561.0/4096.0*t36*t221-6561.0/4096.0*t234-6561.0/4096.0*t38*t224-6561.0/4096.0*t226-6561.0/4096.0*t227-6561.0/4096.0*t228-6561.0/4096.0*t42*t221-6561.0/4096.0*t237-6561.0/4096.0*t44*t224;
  values[18] = 6561.0/4096.0*t51*t221+6561.0/4096.0*t241+6561.0/4096.0*t53*t224+6561.0/4096.0*t243+6561.0/4096.0*t244+6561.0/4096.0*t245+6561.0/4096.0*t64*t221+6561.0/4096.0*t247+6561.0/4096.0*t66*t224;
  values[19] = -2187.0/4096.0*t72*t221-2187.0/4096.0*t251-2187.0/4096.0*t74*t224-2187.0/4096.0*t243-2187.0/4096.0*t244-2187.0/4096.0*t245-2187.0/4096.0*t78*t221-2187.0/4096.0*t254-2187.0/4096.0*t80*t224;
  values[20] = -6561.0/4096.0*t85*t221-6561.0/4096.0*t223-6561.0/4096.0*t15*t258-6561.0/4096.0*t260-6561.0/4096.0*t227-6561.0/4096.0*t261-6561.0/4096.0*t93*t221-6561.0/4096.0*t230-6561.0/4096.0*t30*t258;
  values[21] = 19683.0/4096.0*t97*t221+19683.0/4096.0*t234+19683.0/4096.0*t38*t258+19683.0/4096.0*t260+19683.0/4096.0*t227+19683.0/4096.0*t261+19683.0/4096.0*t100*t221+19683.0/4096.0*t237+19683.0/4096.0*t44*t258;
  values[22] = -19683.0/4096.0*t104*t221-19683.0/4096.0*t241-19683.0/4096.0*t53*t258-19683.0/4096.0*t272-19683.0/4096.0*t244-19683.0/4096.0*t273-19683.0/4096.0*t110*t221-19683.0/4096.0*t247-19683.0/4096.0*t66*t258;
  values[23] = 6561.0/4096.0*t114*t221+6561.0/4096.0*t251+6561.0/4096.0*t74*t258+6561.0/4096.0*t272+6561.0/4096.0*t244+6561.0/4096.0*t273+6561.0/4096.0*t117*t221+6561.0/4096.0*t254+6561.0/4096.0*t80*t258;
  values[24] = 6561.0/4096.0*t122*t282+6561.0/4096.0*t284+6561.0/4096.0*t128*t285+6561.0/4096.0*t287+6561.0/4096.0*t288+6561.0/4096.0*t289+6561.0/4096.0*t138*t282+6561.0/4096.0*t291+6561.0/4096.0*t140*t285;
  values[25] = -19683.0/4096.0*t144*t282-19683.0/4096.0*t295-19683.0/4096.0*t146*t285-19683.0/4096.0*t287-19683.0/4096.0*t288-19683.0/4096.0*t289-19683.0/4096.0*t149*t282-19683.0/4096.0*t298-19683.0/4096.0*t151*t285;
  values[26] = 19683.0/4096.0*t155*t282+19683.0/4096.0*t302+19683.0/4096.0*t157*t285+19683.0/4096.0*t304+19683.0/4096.0*t305+19683.0/4096.0*t306+19683.0/4096.0*t165*t282+19683.0/4096.0*t308+19683.0/4096.0*t167*t285;
  values[27] = -6561.0/4096.0*t171*t282-6561.0/4096.0*t312-6561.0/4096.0*t173*t285-6561.0/4096.0*t304-6561.0/4096.0*t305-6561.0/4096.0*t306-6561.0/4096.0*t176*t282-6561.0/4096.0*t315-6561.0/4096.0*t178*t285;
  values[28] = -2187.0/4096.0*t183*t282-2187.0/4096.0*t284-2187.0/4096.0*t128*t319-2187.0/4096.0*t321-2187.0/4096.0*t288-2187.0/4096.0*t322-2187.0/4096.0*t191*t282-2187.0/4096.0*t291-2187.0/4096.0*t140*t319;
  values[29] = 6561.0/4096.0*t195*t282+6561.0/4096.0*t295+6561.0/4096.0*t146*t319+6561.0/4096.0*t321+6561.0/4096.0*t288+6561.0/4096.0*t322+6561.0/4096.0*t198*t282+6561.0/4096.0*t298+6561.0/4096.0*t151*t319;
  values[30] = -6561.0/4096.0*t202*t282-6561.0/4096.0*t302-6561.0/4096.0*t157*t319-6561.0/4096.0*t333-6561.0/4096.0*t305-6561.0/4096.0*t334-6561.0/4096.0*t208*t282-6561.0/4096.0*t308-6561.0/4096.0*t167*t319;
  values[31] = 2187.0/4096.0*t212*t282+2187.0/4096.0*t312+2187.0/4096.0*t173*t319+2187.0/4096.0*t333+2187.0/4096.0*t305+2187.0/4096.0*t334+2187.0/4096.0*t215*t282+2187.0/4096.0*t315+2187.0/4096.0*t178*t319;
  values[32] = 2187.0/4096.0*t5*t348+2187.0/4096.0*t350+2187.0/4096.0*t15*t352+2187.0/4096.0*t354+2187.0/4096.0*t355+2187.0/4096.0*t356+2187.0/4096.0*t28*t348+2187.0/4096.0*t358+2187.0/4096.0*t30*t352;
  values[33] = -6561.0/4096.0*t36*t348-6561.0/4096.0*t362-6561.0/4096.0*t38*t352-6561.0/4096.0*t354-6561.0/4096.0*t355-6561.0/4096.0*t356-6561.0/4096.0*t42*t348-6561.0/4096.0*t365-6561.0/4096.0*t44*t352;
  values[34] = 6561.0/4096.0*t51*t348+6561.0/4096.0*t369+6561.0/4096.0*t53*t352+6561.0/4096.0*t371+6561.0/4096.0*t372+6561.0/4096.0*t373+6561.0/4096.0*t64*t348+6561.0/4096.0*t375+6561.0/4096.0*t66*t352;
  values[35] = -2187.0/4096.0*t72*t348-2187.0/4096.0*t379-2187.0/4096.0*t74*t352-2187.0/4096.0*t371-2187.0/4096.0*t372-2187.0/4096.0*t373-2187.0/4096.0*t78*t348-2187.0/4096.0*t382-2187.0/4096.0*t80*t352;
  values[36] = -6561.0/4096.0*t85*t348-6561.0/4096.0*t350-6561.0/4096.0*t15*t387-6561.0/4096.0*t389-6561.0/4096.0*t355-6561.0/4096.0*t390-6561.0/4096.0*t93*t348-6561.0/4096.0*t358-6561.0/4096.0*t30*t387;
  values[37] = 19683.0/4096.0*t97*t348+19683.0/4096.0*t362+19683.0/4096.0*t38*t387+19683.0/4096.0*t389+19683.0/4096.0*t355+19683.0/4096.0*t390+19683.0/4096.0*t100*t348+19683.0/4096.0*t365+19683.0/4096.0*t44*t387;
  values[38] = -19683.0/4096.0*t104*t348-19683.0/4096.0*t369-19683.0/4096.0*t53*t387-19683.0/4096.0*t401-19683.0/4096.0*t372-19683.0/4096.0*t402-19683.0/4096.0*t110*t348-19683.0/4096.0*t375-19683.0/4096.0*t66*t387;
  values[39] = 6561.0/4096.0*t114*t348+6561.0/4096.0*t379+6561.0/4096.0*t74*t387+6561.0/4096.0*t401+6561.0/4096.0*t372+6561.0/4096.0*t402+6561.0/4096.0*t117*t348+6561.0/4096.0*t382+6561.0/4096.0*t80*t387;
  values[40] = 6561.0/4096.0*t122*t412+6561.0/4096.0*t414+6561.0/4096.0*t128*t416+6561.0/4096.0*t418+6561.0/4096.0*t419+6561.0/4096.0*t420+6561.0/4096.0*t138*t412+6561.0/4096.0*t422+6561.0/4096.0*t140*t416;
  values[41] = -19683.0/4096.0*t144*t412-19683.0/4096.0*t426-19683.0/4096.0*t146*t416-19683.0/4096.0*t418-19683.0/4096.0*t419-19683.0/4096.0*t420-19683.0/4096.0*t149*t412-19683.0/4096.0*t429-19683.0/4096.0*t151*t416;
  values[42] = 19683.0/4096.0*t155*t412+19683.0/4096.0*t433+19683.0/4096.0*t157*t416+19683.0/4096.0*t435+19683.0/4096.0*t436+19683.0/4096.0*t437+19683.0/4096.0*t165*t412+19683.0/4096.0*t439+19683.0/4096.0*t167*t416;
  values[43] = -6561.0/4096.0*t171*t412-6561.0/4096.0*t443-6561.0/4096.0*t173*t416-6561.0/4096.0*t435-6561.0/4096.0*t436-6561.0/4096.0*t437-6561.0/4096.0*t176*t412-6561.0/4096.0*t446-6561.0/4096.0*t178*t416;
  values[44] = -2187.0/4096.0*t183*t412-2187.0/4096.0*t414-2187.0/4096.0*t128*t451-2187.0/4096.0*t453-2187.0/4096.0*t419-2187.0/4096.0*t454-2187.0/4096.0*t191*t412-2187.0/4096.0*t422-2187.0/4096.0*t140*t451;
  values[45] = 6561.0/4096.0*t195*t412+6561.0/4096.0*t426+6561.0/4096.0*t146*t451+6561.0/4096.0*t453+6561.0/4096.0*t419+6561.0/4096.0*t454+6561.0/4096.0*t198*t412+6561.0/4096.0*t429+6561.0/4096.0*t151*t451;
  values[46] = -6561.0/4096.0*t202*t412-6561.0/4096.0*t433-6561.0/4096.0*t157*t451-6561.0/4096.0*t465-6561.0/4096.0*t436-6561.0/4096.0*t466-6561.0/4096.0*t208*t412-6561.0/4096.0*t439-6561.0/4096.0*t167*t451;
  values[47] = 2187.0/4096.0*t212*t412+2187.0/4096.0*t443+2187.0/4096.0*t173*t451+2187.0/4096.0*t465+2187.0/4096.0*t436+2187.0/4096.0*t466+2187.0/4096.0*t215*t412+2187.0/4096.0*t446+2187.0/4096.0*t178*t451;
  values[48] = -729.0/4096.0*t5*t477-729.0/4096.0*t479-729.0/4096.0*t15*t480-729.0/4096.0*t482-729.0/4096.0*t483-729.0/4096.0*t484-729.0/4096.0*t28*t477-729.0/4096.0*t486-729.0/4096.0*t30*t480;
  values[49] = 2187.0/4096.0*t36*t477+2187.0/4096.0*t490+2187.0/4096.0*t38*t480+2187.0/4096.0*t482+2187.0/4096.0*t483+2187.0/4096.0*t484+2187.0/4096.0*t42*t477+2187.0/4096.0*t493+2187.0/4096.0*t44*t480;
  values[50] = -2187.0/4096.0*t51*t477-2187.0/4096.0*t497-2187.0/4096.0*t53*t480-2187.0/4096.0*t499-2187.0/4096.0*t500-2187.0/4096.0*t501-2187.0/4096.0*t64*t477-2187.0/4096.0*t503-2187.0/4096.0*t66*t480;
  values[51] = 729.0/4096.0*t72*t477+729.0/4096.0*t507+729.0/4096.0*t74*t480+729.0/4096.0*t499+729.0/4096.0*t500+729.0/4096.0*t501+729.0/4096.0*t78*t477+729.0/4096.0*t510+729.0/4096.0*t80*t480;
  values[52] = 2187.0/4096.0*t85*t477+2187.0/4096.0*t479+2187.0/4096.0*t15*t514+2187.0/4096.0*t516+2187.0/4096.0*t483+2187.0/4096.0*t517+2187.0/4096.0*t93*t477+2187.0/4096.0*t486+2187.0/4096.0*t30*t514;
  values[53] = -6561.0/4096.0*t97*t477-6561.0/4096.0*t490-6561.0/4096.0*t38*t514-6561.0/4096.0*t516-6561.0/4096.0*t483-6561.0/4096.0*t517-6561.0/4096.0*t100*t477-6561.0/4096.0*t493-6561.0/4096.0*t44*t514;
  values[54] = 6561.0/4096.0*t104*t477+6561.0/4096.0*t497+6561.0/4096.0*t53*t514+6561.0/4096.0*t528+6561.0/4096.0*t500+6561.0/4096.0*t529+6561.0/4096.0*t110*t477+6561.0/4096.0*t503+6561.0/4096.0*t66*t514;
  values[55] = -2187.0/4096.0*t114*t477-2187.0/4096.0*t507-2187.0/4096.0*t74*t514-2187.0/4096.0*t528-2187.0/4096.0*t500-2187.0/4096.0*t529-2187.0/4096.0*t117*t477-2187.0/4096.0*t510-2187.0/4096.0*t80*t514;
  values[56] = -2187.0/4096.0*t122*t538-2187.0/4096.0*t540-2187.0/4096.0*t128*t541-2187.0/4096.0*t543-2187.0/4096.0*t544-2187.0/4096.0*t545-2187.0/4096.0*t138*t538-2187.0/4096.0*t547-2187.0/4096.0*t140*t541;
  values[57] = 6561.0/4096.0*t144*t538+6561.0/4096.0*t551+6561.0/4096.0*t146*t541+6561.0/4096.0*t543+6561.0/4096.0*t544+6561.0/4096.0*t545+6561.0/4096.0*t149*t538+6561.0/4096.0*t554+6561.0/4096.0*t151*t541;
  values[58] = -6561.0/4096.0*t155*t538-6561.0/4096.0*t558-6561.0/4096.0*t157*t541-6561.0/4096.0*t560-6561.0/4096.0*t561-6561.0/4096.0*t562-6561.0/4096.0*t165*t538-6561.0/4096.0*t564-6561.0/4096.0*t167*t541;
  values[59] = 2187.0/4096.0*t171*t538+2187.0/4096.0*t568+2187.0/4096.0*t173*t541+2187.0/4096.0*t560+2187.0/4096.0*t561+2187.0/4096.0*t562+2187.0/4096.0*t176*t538+2187.0/4096.0*t571+2187.0/4096.0*t178*t541;
  values[60] = 729.0/4096.0*t183*t538+729.0/4096.0*t540+729.0/4096.0*t128*t575+729.0/4096.0*t577+729.0/4096.0*t544+729.0/4096.0*t578+729.0/4096.0*t191*t538+729.0/4096.0*t547+729.0/4096.0*t140*t575;
  values[61] = -2187.0/4096.0*t195*t538-2187.0/4096.0*t551-2187.0/4096.0*t146*t575-2187.0/4096.0*t577-2187.0/4096.0*t544-2187.0/4096.0*t578-2187.0/4096.0*t198*t538-2187.0/4096.0*t554-2187.0/4096.0*t151*t575;
  values[62] = 2187.0/4096.0*t202*t538+2187.0/4096.0*t558+2187.0/4096.0*t157*t575+2187.0/4096.0*t589+2187.0/4096.0*t561+2187.0/4096.0*t590+2187.0/4096.0*t208*t538+2187.0/4096.0*t564+2187.0/4096.0*t167*t575;
  values[63] = -729.0/4096.0*t212*t538-729.0/4096.0*t568-729.0/4096.0*t173*t575-729.0/4096.0*t589-729.0/4096.0*t561-729.0/4096.0*t590-729.0/4096.0*t215*t538-729.0/4096.0*t571-729.0/4096.0*t178*t575;
}

static void C_H_Q3_3D_DeriveXiZeta(double xi, double eta, double zeta,
                             double *values)
{
  double t1, t2, t3, t4, t5, t7, t8, t9, t10, t11, t12, t14, t15, t16, t17;
  double t18, t19, t21, t22, t23, t24, t25, t26, t27, t28, t30, t33, t34;
  double t35, t37, t39, t40, t42, t45, t46, t47, t48, t50, t52, t53, t54;
  double t55, t56, t57, t58, t59, t61, t64, t65, t66, t68, t70, t71, t73;
  double t77, t78, t80, t81, t82, t84, t85, t86, t88, t92, t95, t99, t101;
  double t102, t103, t105, t109, t112, t115, t116, t118, t119, t120, t122;
  double t123, t124, t126, t127, t128, t129, t130, t132, t135, t137, t139;
  double t141, t144, t146, t148, t149, t150, t151, t152, t154, t157, t159;
  double t161, t163, t167, t168, t170, t171, t172, t174, t175, t176, t178;
  double t182, t185, t189, t191, t192, t193, t195, t199, t202, t205, t206;
  double t207, t209, t210, t212, t213, t224, t225, t234, t236, t238, t239;
  double t250, t251, t260, t262, t264, t265, t276, t277, t286, t288, t290;
  double t291, t302, t303, t312, t313, t314, t315, t317, t318, t319, t320;
  double t321, t322, t324, t325, t326, t328, t332, t335, t339, t341, t342;
  double t343, t345, t349, t352, t355, t357, t358, t359, t361, t362, t363;
  double t365, t369, t372, t376, t378, t379, t380, t382, t386, t389, t392;
  double t394, t395, t396, t398, t399, t400, t402, t406, t409, t413, t415;
  double t416, t417, t419, t423, t426, t429, t431, t432, t433, t435, t436;
  double t437, t439, t443, t446, t450, t452, t453, t454, t456, t460, t463;
  double t466, t467, t468, t470, t471, t473, t474, t485, t486, t495, t497;
  double t499, t500, t511, t512, t521, t523, t525, t526, t537, t538, t547;
  double t549, t551, t552, t563, t564;

  t1 = xi+1.0/3.0;
  t2 = xi-1.0/3.0;
  t3 = t1*t2;
  t4 = eta-1.0;
  t5 = t3*t4;
  t7 = eta-1.0/3.0;
  t8 = (eta+1.0/3.0)*t7;
  t9 = zeta+1.0/3.0;
  t10 = zeta-1.0/3.0;
  t11 = t9*t10;
  t12 = t8*t11;
  t14 = zeta-1.0;
  t15 = t14*t10;
  t16 = t8*t15;
  t17 = t5*t16;
  t18 = t14*t9;
  t19 = t8*t18;
  t21 = xi-1.0;
  t22 = t21*t2;
  t23 = t22*t4;
  t24 = t23*t12;
  t25 = t23*t16;
  t26 = t23*t19;
  t27 = t21*t1;
  t28 = t27*t4;
  t30 = t28*t16;
  t33 = xi+1.0;
  t34 = t33*t2;
  t35 = t34*t4;
  t37 = t35*t16;
  t39 = t21*t33;
  t40 = t39*t4;
  t42 = t40*t16;
  t45 = -xi+1.0;
  t46 = -xi-1.0/3.0;
  t47 = t45*t46;
  t48 = t47*t4;
  t50 = t48*t16;
  t52 = -xi-1.0;
  t53 = t52*t46;
  t54 = t53*t4;
  t55 = t54*t12;
  t56 = t54*t16;
  t57 = t54*t19;
  t58 = t52*t45;
  t59 = t58*t4;
  t61 = t59*t16;
  t64 = -xi+1.0/3.0;
  t65 = t64*t46;
  t66 = t65*t4;
  t68 = t66*t16;
  t70 = t52*t64;
  t71 = t70*t4;
  t73 = t71*t16;
  t77 = (eta+1.0)*t7;
  t78 = t77*t11;
  t80 = t77*t15;
  t81 = t5*t80;
  t82 = t77*t18;
  t84 = t23*t78;
  t85 = t23*t80;
  t86 = t23*t82;
  t88 = t28*t80;
  t92 = t35*t80;
  t95 = t40*t80;
  t99 = t48*t80;
  t101 = t54*t78;
  t102 = t54*t80;
  t103 = t54*t82;
  t105 = t59*t80;
  t109 = t66*t80;
  t112 = t71*t80;
  t115 = -eta-1.0;
  t116 = t3*t115;
  t118 = -eta-1.0/3.0;
  t119 = (-eta+1.0)*t118;
  t120 = t119*t11;
  t122 = t119*t15;
  t123 = t116*t122;
  t124 = t119*t18;
  t126 = t22*t115;
  t127 = t126*t120;
  t128 = t126*t122;
  t129 = t126*t124;
  t130 = t27*t115;
  t132 = t130*t122;
  t135 = t34*t115;
  t137 = t135*t122;
  t139 = t39*t115;
  t141 = t139*t122;
  t144 = t47*t115;
  t146 = t144*t122;
  t148 = t53*t115;
  t149 = t148*t120;
  t150 = t148*t122;
  t151 = t148*t124;
  t152 = t58*t115;
  t154 = t152*t122;
  t157 = t65*t115;
  t159 = t157*t122;
  t161 = t70*t115;
  t163 = t161*t122;
  t167 = (-eta+1.0/3.0)*t118;
  t168 = t167*t11;
  t170 = t167*t15;
  t171 = t116*t170;
  t172 = t167*t18;
  t174 = t126*t168;
  t175 = t126*t170;
  t176 = t126*t172;
  t178 = t130*t170;
  t182 = t135*t170;
  t185 = t139*t170;
  t189 = t144*t170;
  t191 = t148*t168;
  t192 = t148*t170;
  t193 = t148*t172;
  t195 = t152*t170;
  t199 = t157*t170;
  t202 = t161*t170;
  t205 = zeta+1.0;
  t206 = t205*t10;
  t207 = t8*t206;
  t209 = t14*t205;
  t210 = t8*t209;
  t212 = t23*t207;
  t213 = t23*t210;
  t224 = t54*t207;
  t225 = t54*t210;
  t234 = t77*t206;
  t236 = t77*t209;
  t238 = t23*t234;
  t239 = t23*t236;
  t250 = t54*t234;
  t251 = t54*t236;
  t260 = t119*t206;
  t262 = t119*t209;
  t264 = t126*t260;
  t265 = t126*t262;
  t276 = t148*t260;
  t277 = t148*t262;
  t286 = t167*t206;
  t288 = t167*t209;
  t290 = t126*t286;
  t291 = t126*t288;
  t302 = t148*t286;
  t303 = t148*t288;
  t312 = -zeta+1.0;
  t313 = -zeta-1.0/3.0;
  t314 = t312*t313;
  t315 = t8*t314;
  t317 = -zeta-1.0;
  t318 = t317*t313;
  t319 = t8*t318;
  t320 = t5*t319;
  t321 = t317*t312;
  t322 = t8*t321;
  t324 = t23*t315;
  t325 = t23*t319;
  t326 = t23*t322;
  t328 = t28*t319;
  t332 = t35*t319;
  t335 = t40*t319;
  t339 = t48*t319;
  t341 = t54*t315;
  t342 = t54*t319;
  t343 = t54*t322;
  t345 = t59*t319;
  t349 = t66*t319;
  t352 = t71*t319;
  t355 = t77*t314;
  t357 = t77*t318;
  t358 = t5*t357;
  t359 = t77*t321;
  t361 = t23*t355;
  t362 = t23*t357;
  t363 = t23*t359;
  t365 = t28*t357;
  t369 = t35*t357;
  t372 = t40*t357;
  t376 = t48*t357;
  t378 = t54*t355;
  t379 = t54*t357;
  t380 = t54*t359;
  t382 = t59*t357;
  t386 = t66*t357;
  t389 = t71*t357;
  t392 = t119*t314;
  t394 = t119*t318;
  t395 = t116*t394;
  t396 = t119*t321;
  t398 = t126*t392;
  t399 = t126*t394;
  t400 = t126*t396;
  t402 = t130*t394;
  t406 = t135*t394;
  t409 = t139*t394;
  t413 = t144*t394;
  t415 = t148*t392;
  t416 = t148*t394;
  t417 = t148*t396;
  t419 = t152*t394;
  t423 = t157*t394;
  t426 = t161*t394;
  t429 = t167*t314;
  t431 = t167*t318;
  t432 = t116*t431;
  t433 = t167*t321;
  t435 = t126*t429;
  t436 = t126*t431;
  t437 = t126*t433;
  t439 = t130*t431;
  t443 = t135*t431;
  t446 = t139*t431;
  t450 = t144*t431;
  t452 = t148*t429;
  t453 = t148*t431;
  t454 = t148*t433;
  t456 = t152*t431;
  t460 = t157*t431;
  t463 = t161*t431;
  t466 = -zeta+1.0/3.0;
  t467 = t466*t313;
  t468 = t8*t467;
  t470 = t317*t466;
  t471 = t8*t470;
  t473 = t23*t468;
  t474 = t23*t471;
  t485 = t54*t468;
  t486 = t54*t471;
  t495 = t77*t467;
  t497 = t77*t470;
  t499 = t23*t495;
  t500 = t23*t497;
  t511 = t54*t495;
  t512 = t54*t497;
  t521 = t119*t467;
  t523 = t119*t470;
  t525 = t126*t521;
  t526 = t126*t523;
  t537 = t148*t521;
  t538 = t148*t523;
  t547 = t167*t467;
  t549 = t167*t470;
  t551 = t126*t547;
  t552 = t126*t549;
  t563 = t148*t547;
  t564 = t148*t549;

  values[0] = -729.0/4096.0*t5*t12-729.0/4096.0*t17-729.0/4096.0*t5*t19-729.0/4096.0*t24-729.0/4096.0*t25-729.0/4096.0*t26-729.0/4096.0*t28*t12-729.0/4096.0*t30-729.0/4096.0*t28*t19;
  values[1] = 2187.0/4096.0*t35*t12+2187.0/4096.0*t37+2187.0/4096.0*t35*t19+2187.0/4096.0*t24+2187.0/4096.0*t25+2187.0/4096.0*t26+2187.0/4096.0*t40*t12+2187.0/4096.0*t42+2187.0/4096.0*t40*t19;
  values[2] = -2187.0/4096.0*t48*t12-2187.0/4096.0*t50-2187.0/4096.0*t48*t19-2187.0/4096.0*t55-2187.0/4096.0*t56-2187.0/4096.0*t57-2187.0/4096.0*t59*t12-2187.0/4096.0*t61-2187.0/4096.0*t59*t19;
  values[3] = 729.0/4096.0*t66*t12+729.0/4096.0*t68+729.0/4096.0*t66*t19+729.0/4096.0*t55+729.0/4096.0*t56+729.0/4096.0*t57+729.0/4096.0*t71*t12+729.0/4096.0*t73+729.0/4096.0*t71*t19;
  values[4] = 2187.0/4096.0*t5*t78+2187.0/4096.0*t81+2187.0/4096.0*t5*t82+2187.0/4096.0*t84+2187.0/4096.0*t85+2187.0/4096.0*t86+2187.0/4096.0*t28*t78+2187.0/4096.0*t88+2187.0/4096.0*t28*t82;
  values[5] = -6561.0/4096.0*t35*t78-6561.0/4096.0*t92-6561.0/4096.0*t35*t82-6561.0/4096.0*t84-6561.0/4096.0*t85-6561.0/4096.0*t86-6561.0/4096.0*t40*t78-6561.0/4096.0*t95-6561.0/4096.0*t40*t82;
  values[6] = 6561.0/4096.0*t48*t78+6561.0/4096.0*t99+6561.0/4096.0*t48*t82+6561.0/4096.0*t101+6561.0/4096.0*t102+6561.0/4096.0*t103+6561.0/4096.0*t59*t78+6561.0/4096.0*t105+6561.0/4096.0*t59*t82;
  values[7] = -2187.0/4096.0*t66*t78-2187.0/4096.0*t109-2187.0/4096.0*t66*t82-2187.0/4096.0*t101-2187.0/4096.0*t102-2187.0/4096.0*t103-2187.0/4096.0*t71*t78-2187.0/4096.0*t112-2187.0/4096.0*t71*t82;
  values[8] = 2187.0/4096.0*t116*t120+2187.0/4096.0*t123+2187.0/4096.0*t116*t124+2187.0/4096.0*t127+2187.0/4096.0*t128+2187.0/4096.0*t129+2187.0/4096.0*t130*t120+2187.0/4096.0*t132+2187.0/4096.0*t130*t124;
  values[9] = -6561.0/4096.0*t135*t120-6561.0/4096.0*t137-6561.0/4096.0*t135*t124-6561.0/4096.0*t127-6561.0/4096.0*t128-6561.0/4096.0*t129-6561.0/4096.0*t139*t120-6561.0/4096.0*t141-6561.0/4096.0*t139*t124;
  values[10] = 6561.0/4096.0*t144*t120+6561.0/4096.0*t146+6561.0/4096.0*t144*t124+6561.0/4096.0*t149+6561.0/4096.0*t150+6561.0/4096.0*t151+6561.0/4096.0*t152*t120+6561.0/4096.0*t154+6561.0/4096.0*t152*t124;
  values[11] = -2187.0/4096.0*t157*t120-2187.0/4096.0*t159-2187.0/4096.0*t157*t124-2187.0/4096.0*t149-2187.0/4096.0*t150-2187.0/4096.0*t151-2187.0/4096.0*t161*t120-2187.0/4096.0*t163-2187.0/4096.0*t161*t124;
  values[12] = -729.0/4096.0*t116*t168-729.0/4096.0*t171-729.0/4096.0*t116*t172-729.0/4096.0*t174-729.0/4096.0*t175-729.0/4096.0*t176-729.0/4096.0*t130*t168-729.0/4096.0*t178-729.0/4096.0*t130*t172;
  values[13] = 2187.0/4096.0*t135*t168+2187.0/4096.0*t182+2187.0/4096.0*t135*t172+2187.0/4096.0*t174+2187.0/4096.0*t175+2187.0/4096.0*t176+2187.0/4096.0*t139*t168+2187.0/4096.0*t185+2187.0/4096.0*t139*t172;
  values[14] = -2187.0/4096.0*t144*t168-2187.0/4096.0*t189-2187.0/4096.0*t144*t172-2187.0/4096.0*t191-2187.0/4096.0*t192-2187.0/4096.0*t193-2187.0/4096.0*t152*t168-2187.0/4096.0*t195-2187.0/4096.0*t152*t172;
  values[15] = 729.0/4096.0*t157*t168+729.0/4096.0*t199+729.0/4096.0*t157*t172+729.0/4096.0*t191+729.0/4096.0*t192+729.0/4096.0*t193+729.0/4096.0*t161*t168+729.0/4096.0*t202+729.0/4096.0*t161*t172;
  values[16] = 2187.0/4096.0*t5*t207+2187.0/4096.0*t17+2187.0/4096.0*t5*t210+2187.0/4096.0*t212+2187.0/4096.0*t25+2187.0/4096.0*t213+2187.0/4096.0*t28*t207+2187.0/4096.0*t30+2187.0/4096.0*t28*t210;
  values[17] = -6561.0/4096.0*t35*t207-6561.0/4096.0*t37-6561.0/4096.0*t35*t210-6561.0/4096.0*t212-6561.0/4096.0*t25-6561.0/4096.0*t213-6561.0/4096.0*t40*t207-6561.0/4096.0*t42-6561.0/4096.0*t40*t210;
  values[18] = 6561.0/4096.0*t48*t207+6561.0/4096.0*t50+6561.0/4096.0*t48*t210+6561.0/4096.0*t224+6561.0/4096.0*t56+6561.0/4096.0*t225+6561.0/4096.0*t59*t207+6561.0/4096.0*t61+6561.0/4096.0*t59*t210;
  values[19] = -2187.0/4096.0*t66*t207-2187.0/4096.0*t68-2187.0/4096.0*t66*t210-2187.0/4096.0*t224-2187.0/4096.0*t56-2187.0/4096.0*t225-2187.0/4096.0*t71*t207-2187.0/4096.0*t73-2187.0/4096.0*t71*t210;
  values[20] = -6561.0/4096.0*t5*t234-6561.0/4096.0*t81-6561.0/4096.0*t5*t236-6561.0/4096.0*t238-6561.0/4096.0*t85-6561.0/4096.0*t239-6561.0/4096.0*t28*t234-6561.0/4096.0*t88-6561.0/4096.0*t28*t236;
  values[21] = 19683.0/4096.0*t35*t234+19683.0/4096.0*t92+19683.0/4096.0*t35*t236+19683.0/4096.0*t238+19683.0/4096.0*t85+19683.0/4096.0*t239+19683.0/4096.0*t40*t234+19683.0/4096.0*t95+19683.0/4096.0*t40*t236;
  values[22] = -19683.0/4096.0*t48*t234-19683.0/4096.0*t99-19683.0/4096.0*t48*t236-19683.0/4096.0*t250-19683.0/4096.0*t102-19683.0/4096.0*t251-19683.0/4096.0*t59*t234-19683.0/4096.0*t105-19683.0/4096.0*t59*t236;
  values[23] = 6561.0/4096.0*t66*t234+6561.0/4096.0*t109+6561.0/4096.0*t66*t236+6561.0/4096.0*t250+6561.0/4096.0*t102+6561.0/4096.0*t251+6561.0/4096.0*t71*t234+6561.0/4096.0*t112+6561.0/4096.0*t71*t236;
  values[24] = -6561.0/4096.0*t116*t260-6561.0/4096.0*t123-6561.0/4096.0*t116*t262-6561.0/4096.0*t264-6561.0/4096.0*t128-6561.0/4096.0*t265-6561.0/4096.0*t130*t260-6561.0/4096.0*t132-6561.0/4096.0*t130*t262;
  values[25] = 19683.0/4096.0*t135*t260+19683.0/4096.0*t137+19683.0/4096.0*t135*t262+19683.0/4096.0*t264+19683.0/4096.0*t128+19683.0/4096.0*t265+19683.0/4096.0*t139*t260+19683.0/4096.0*t141+19683.0/4096.0*t139*t262;
  values[26] = -19683.0/4096.0*t144*t260-19683.0/4096.0*t146-19683.0/4096.0*t144*t262-19683.0/4096.0*t276-19683.0/4096.0*t150-19683.0/4096.0*t277-19683.0/4096.0*t152*t260-19683.0/4096.0*t154-19683.0/4096.0*t152*t262;
  values[27] = 6561.0/4096.0*t157*t260+6561.0/4096.0*t159+6561.0/4096.0*t157*t262+6561.0/4096.0*t276+6561.0/4096.0*t150+6561.0/4096.0*t277+6561.0/4096.0*t161*t260+6561.0/4096.0*t163+6561.0/4096.0*t161*t262;
  values[28] = 2187.0/4096.0*t116*t286+2187.0/4096.0*t171+2187.0/4096.0*t116*t288+2187.0/4096.0*t290+2187.0/4096.0*t175+2187.0/4096.0*t291+2187.0/4096.0*t130*t286+2187.0/4096.0*t178+2187.0/4096.0*t130*t288;
  values[29] = -6561.0/4096.0*t135*t286-6561.0/4096.0*t182-6561.0/4096.0*t135*t288-6561.0/4096.0*t290-6561.0/4096.0*t175-6561.0/4096.0*t291-6561.0/4096.0*t139*t286-6561.0/4096.0*t185-6561.0/4096.0*t139*t288;
  values[30] = 6561.0/4096.0*t144*t286+6561.0/4096.0*t189+6561.0/4096.0*t144*t288+6561.0/4096.0*t302+6561.0/4096.0*t192+6561.0/4096.0*t303+6561.0/4096.0*t152*t286+6561.0/4096.0*t195+6561.0/4096.0*t152*t288;
  values[31] = -2187.0/4096.0*t157*t286-2187.0/4096.0*t199-2187.0/4096.0*t157*t288-2187.0/4096.0*t302-2187.0/4096.0*t192-2187.0/4096.0*t303-2187.0/4096.0*t161*t286-2187.0/4096.0*t202-2187.0/4096.0*t161*t288;
  values[32] = -2187.0/4096.0*t5*t315-2187.0/4096.0*t320-2187.0/4096.0*t5*t322-2187.0/4096.0*t324-2187.0/4096.0*t325-2187.0/4096.0*t326-2187.0/4096.0*t28*t315-2187.0/4096.0*t328-2187.0/4096.0*t28*t322;
  values[33] = 6561.0/4096.0*t35*t315+6561.0/4096.0*t332+6561.0/4096.0*t35*t322+6561.0/4096.0*t324+6561.0/4096.0*t325+6561.0/4096.0*t326+6561.0/4096.0*t40*t315+6561.0/4096.0*t335+6561.0/4096.0*t40*t322;
  values[34] = -6561.0/4096.0*t48*t315-6561.0/4096.0*t339-6561.0/4096.0*t48*t322-6561.0/4096.0*t341-6561.0/4096.0*t342-6561.0/4096.0*t343-6561.0/4096.0*t59*t315-6561.0/4096.0*t345-6561.0/4096.0*t59*t322;
  values[35] = 2187.0/4096.0*t66*t315+2187.0/4096.0*t349+2187.0/4096.0*t66*t322+2187.0/4096.0*t341+2187.0/4096.0*t342+2187.0/4096.0*t343+2187.0/4096.0*t71*t315+2187.0/4096.0*t352+2187.0/4096.0*t71*t322;
  values[36] = 6561.0/4096.0*t5*t355+6561.0/4096.0*t358+6561.0/4096.0*t5*t359+6561.0/4096.0*t361+6561.0/4096.0*t362+6561.0/4096.0*t363+6561.0/4096.0*t28*t355+6561.0/4096.0*t365+6561.0/4096.0*t28*t359;
  values[37] = -19683.0/4096.0*t35*t355-19683.0/4096.0*t369-19683.0/4096.0*t35*t359-19683.0/4096.0*t361-19683.0/4096.0*t362-19683.0/4096.0*t363-19683.0/4096.0*t40*t355-19683.0/4096.0*t372-19683.0/4096.0*t40*t359;
  values[38] = 19683.0/4096.0*t48*t355+19683.0/4096.0*t376+19683.0/4096.0*t48*t359+19683.0/4096.0*t378+19683.0/4096.0*t379+19683.0/4096.0*t380+19683.0/4096.0*t59*t355+19683.0/4096.0*t382+19683.0/4096.0*t59*t359;
  values[39] = -6561.0/4096.0*t66*t355-6561.0/4096.0*t386-6561.0/4096.0*t66*t359-6561.0/4096.0*t378-6561.0/4096.0*t379-6561.0/4096.0*t380-6561.0/4096.0*t71*t355-6561.0/4096.0*t389-6561.0/4096.0*t71*t359;
  values[40] = 6561.0/4096.0*t116*t392+6561.0/4096.0*t395+6561.0/4096.0*t116*t396+6561.0/4096.0*t398+6561.0/4096.0*t399+6561.0/4096.0*t400+6561.0/4096.0*t130*t392+6561.0/4096.0*t402+6561.0/4096.0*t130*t396;
  values[41] = -19683.0/4096.0*t135*t392-19683.0/4096.0*t406-19683.0/4096.0*t135*t396-19683.0/4096.0*t398-19683.0/4096.0*t399-19683.0/4096.0*t400-19683.0/4096.0*t139*t392-19683.0/4096.0*t409-19683.0/4096.0*t139*t396;
  values[42] = 19683.0/4096.0*t144*t392+19683.0/4096.0*t413+19683.0/4096.0*t144*t396+19683.0/4096.0*t415+19683.0/4096.0*t416+19683.0/4096.0*t417+19683.0/4096.0*t152*t392+19683.0/4096.0*t419+19683.0/4096.0*t152*t396;
  values[43] = -6561.0/4096.0*t157*t392-6561.0/4096.0*t423-6561.0/4096.0*t157*t396-6561.0/4096.0*t415-6561.0/4096.0*t416-6561.0/4096.0*t417-6561.0/4096.0*t161*t392-6561.0/4096.0*t426-6561.0/4096.0*t161*t396;
  values[44] = -2187.0/4096.0*t116*t429-2187.0/4096.0*t432-2187.0/4096.0*t116*t433-2187.0/4096.0*t435-2187.0/4096.0*t436-2187.0/4096.0*t437-2187.0/4096.0*t130*t429-2187.0/4096.0*t439-2187.0/4096.0*t130*t433;
  values[45] = 6561.0/4096.0*t135*t429+6561.0/4096.0*t443+6561.0/4096.0*t135*t433+6561.0/4096.0*t435+6561.0/4096.0*t436+6561.0/4096.0*t437+6561.0/4096.0*t139*t429+6561.0/4096.0*t446+6561.0/4096.0*t139*t433;
  values[46] = -6561.0/4096.0*t144*t429-6561.0/4096.0*t450-6561.0/4096.0*t144*t433-6561.0/4096.0*t452-6561.0/4096.0*t453-6561.0/4096.0*t454-6561.0/4096.0*t152*t429-6561.0/4096.0*t456-6561.0/4096.0*t152*t433;
  values[47] = 2187.0/4096.0*t157*t429+2187.0/4096.0*t460+2187.0/4096.0*t157*t433+2187.0/4096.0*t452+2187.0/4096.0*t453+2187.0/4096.0*t454+2187.0/4096.0*t161*t429+2187.0/4096.0*t463+2187.0/4096.0*t161*t433;
  values[48] = 729.0/4096.0*t5*t468+729.0/4096.0*t320+729.0/4096.0*t5*t471+729.0/4096.0*t473+729.0/4096.0*t325+729.0/4096.0*t474+729.0/4096.0*t28*t468+729.0/4096.0*t328+729.0/4096.0*t28*t471;
  values[49] = -2187.0/4096.0*t35*t468-2187.0/4096.0*t332-2187.0/4096.0*t35*t471-2187.0/4096.0*t473-2187.0/4096.0*t325-2187.0/4096.0*t474-2187.0/4096.0*t40*t468-2187.0/4096.0*t335-2187.0/4096.0*t40*t471;
  values[50] = 2187.0/4096.0*t48*t468+2187.0/4096.0*t339+2187.0/4096.0*t48*t471+2187.0/4096.0*t485+2187.0/4096.0*t342+2187.0/4096.0*t486+2187.0/4096.0*t59*t468+2187.0/4096.0*t345+2187.0/4096.0*t59*t471;
  values[51] = -729.0/4096.0*t66*t468-729.0/4096.0*t349-729.0/4096.0*t66*t471-729.0/4096.0*t485-729.0/4096.0*t342-729.0/4096.0*t486-729.0/4096.0*t71*t468-729.0/4096.0*t352-729.0/4096.0*t71*t471;
  values[52] = -2187.0/4096.0*t5*t495-2187.0/4096.0*t358-2187.0/4096.0*t5*t497-2187.0/4096.0*t499-2187.0/4096.0*t362-2187.0/4096.0*t500-2187.0/4096.0*t28*t495-2187.0/4096.0*t365-2187.0/4096.0*t28*t497;
  values[53] = 6561.0/4096.0*t35*t495+6561.0/4096.0*t369+6561.0/4096.0*t35*t497+6561.0/4096.0*t499+6561.0/4096.0*t362+6561.0/4096.0*t500+6561.0/4096.0*t40*t495+6561.0/4096.0*t372+6561.0/4096.0*t40*t497;
  values[54] = -6561.0/4096.0*t48*t495-6561.0/4096.0*t376-6561.0/4096.0*t48*t497-6561.0/4096.0*t511-6561.0/4096.0*t379-6561.0/4096.0*t512-6561.0/4096.0*t59*t495-6561.0/4096.0*t382-6561.0/4096.0*t59*t497;
  values[55] = 2187.0/4096.0*t66*t495+2187.0/4096.0*t386+2187.0/4096.0*t66*t497+2187.0/4096.0*t511+2187.0/4096.0*t379+2187.0/4096.0*t512+2187.0/4096.0*t71*t495+2187.0/4096.0*t389+2187.0/4096.0*t71*t497;
  values[56] = -2187.0/4096.0*t116*t521-2187.0/4096.0*t395-2187.0/4096.0*t116*t523-2187.0/4096.0*t525-2187.0/4096.0*t399-2187.0/4096.0*t526-2187.0/4096.0*t130*t521-2187.0/4096.0*t402-2187.0/4096.0*t130*t523;
  values[57] = 6561.0/4096.0*t135*t521+6561.0/4096.0*t406+6561.0/4096.0*t135*t523+6561.0/4096.0*t525+6561.0/4096.0*t399+6561.0/4096.0*t526+6561.0/4096.0*t139*t521+6561.0/4096.0*t409+6561.0/4096.0*t139*t523;
  values[58] = -6561.0/4096.0*t144*t521-6561.0/4096.0*t413-6561.0/4096.0*t144*t523-6561.0/4096.0*t537-6561.0/4096.0*t416-6561.0/4096.0*t538-6561.0/4096.0*t152*t521-6561.0/4096.0*t419-6561.0/4096.0*t152*t523;
  values[59] = 2187.0/4096.0*t157*t521+2187.0/4096.0*t423+2187.0/4096.0*t157*t523+2187.0/4096.0*t537+2187.0/4096.0*t416+2187.0/4096.0*t538+2187.0/4096.0*t161*t521+2187.0/4096.0*t426+2187.0/4096.0*t161*t523;
  values[60] = 729.0/4096.0*t116*t547+729.0/4096.0*t432+729.0/4096.0*t116*t549+729.0/4096.0*t551+729.0/4096.0*t436+729.0/4096.0*t552+729.0/4096.0*t130*t547+729.0/4096.0*t439+729.0/4096.0*t130*t549;
  values[61] = -2187.0/4096.0*t135*t547-2187.0/4096.0*t443-2187.0/4096.0*t135*t549-2187.0/4096.0*t551-2187.0/4096.0*t436-2187.0/4096.0*t552-2187.0/4096.0*t139*t547-2187.0/4096.0*t446-2187.0/4096.0*t139*t549;
  values[62] = 2187.0/4096.0*t144*t547+2187.0/4096.0*t450+2187.0/4096.0*t144*t549+2187.0/4096.0*t563+2187.0/4096.0*t453+2187.0/4096.0*t564+2187.0/4096.0*t152*t547+2187.0/4096.0*t456+2187.0/4096.0*t152*t549;
  values[63] = -729.0/4096.0*t157*t547-729.0/4096.0*t460-729.0/4096.0*t157*t549-729.0/4096.0*t563-729.0/4096.0*t453-729.0/4096.0*t564-729.0/4096.0*t161*t547-729.0/4096.0*t463-729.0/4096.0*t161*t549;
}

static void C_H_Q3_3D_DeriveEtaEta(double xi, double eta, double zeta,
                             double *values)
{
  double t1, t4, t5, t6, t7, t8, t10, t11, t12, t13, t14, t15, t16, t18;
  double t19, t20, t21, t25, t26, t28, t30, t33, t34, t35, t37, t41, t42;
  double t44, t46, t47, t48, t57, t58, t59, t60, t61, t62, t63, t65, t66;
  double t67, t68, t70, t72, t74, t76, t78, t80, t82, t83, t84, t94, t95;
  double t96, t97, t99, t100, t102, t104, t106, t108, t110, t112, t114;
  double t123, t124, t125, t127, t128, t130, t132, t134, t136, t138, t140;
  double t142, t151, t152, t154, t155, t156, t157, t158, t159, t161, t162;
  double t163, t165, t167, t169, t171, t173, t175, t177, t178, t187, t188;
  double t189, t190, t191, t193, t194, t195, t197, t199, t201, t203, t205;
  double t207, t209, t210, t220, t221, t222, t223, t225, t226, t228, t230;
  double t232, t234, t236, t238, t240, t249, t250, t251, t253, t254, t256;
  double t258, t260, t262, t264, t266, t268;

  t1 = xi-1.0;
  t4 = xi-1.0/3.0;
  t5 = t1*(xi+1.0/3.0)*t4;
  t6 = eta-1.0/3.0;
  t7 = zeta-1.0;
  t8 = t6*t7;
  t10 = zeta-1.0/3.0;
  t11 = (zeta+1.0/3.0)*t10;
  t12 = t8*t11;
  t13 = t5*t12;
  t14 = eta+1.0/3.0;
  t15 = t14*t7;
  t16 = t15*t11;
  t18 = eta-1.0;
  t19 = t18*t7;
  t20 = t19*t11;
  t21 = t5*t20;
  t25 = t1*(xi+1.0)*t4;
  t26 = t25*t12;
  t28 = t25*t20;
  t30 = -xi-1.0;
  t33 = -xi-1.0/3.0;
  t34 = t30*(-xi+1.0)*t33;
  t35 = t34*t12;
  t37 = t34*t20;
  t41 = t30*(-xi+1.0/3.0)*t33;
  t42 = t41*t12;
  t44 = t41*t20;
  t46 = eta+1.0;
  t47 = t46*t7;
  t48 = t47*t11;
  t57 = -eta-1.0/3.0;
  t58 = t57*t7;
  t59 = t58*t11;
  t60 = t5*t59;
  t61 = -eta+1.0;
  t62 = t61*t7;
  t63 = t62*t11;
  t65 = -eta-1.0;
  t66 = t65*t7;
  t67 = t66*t11;
  t68 = t5*t67;
  t70 = t25*t59;
  t72 = t25*t67;
  t74 = t34*t59;
  t76 = t34*t67;
  t78 = t41*t59;
  t80 = t41*t67;
  t82 = -eta+1.0/3.0;
  t83 = t82*t7;
  t84 = t83*t11;
  t94 = (zeta+1.0)*t10;
  t95 = t8*t94;
  t96 = t5*t95;
  t97 = t15*t94;
  t99 = t19*t94;
  t100 = t5*t99;
  t102 = t25*t95;
  t104 = t25*t99;
  t106 = t34*t95;
  t108 = t34*t99;
  t110 = t41*t95;
  t112 = t41*t99;
  t114 = t47*t94;
  t123 = t58*t94;
  t124 = t5*t123;
  t125 = t62*t94;
  t127 = t66*t94;
  t128 = t5*t127;
  t130 = t25*t123;
  t132 = t25*t127;
  t134 = t34*t123;
  t136 = t34*t127;
  t138 = t41*t123;
  t140 = t41*t127;
  t142 = t83*t94;
  t151 = -zeta-1.0;
  t152 = t6*t151;
  t154 = -zeta-1.0/3.0;
  t155 = (-zeta+1.0)*t154;
  t156 = t152*t155;
  t157 = t5*t156;
  t158 = t14*t151;
  t159 = t158*t155;
  t161 = t18*t151;
  t162 = t161*t155;
  t163 = t5*t162;
  t165 = t25*t156;
  t167 = t25*t162;
  t169 = t34*t156;
  t171 = t34*t162;
  t173 = t41*t156;
  t175 = t41*t162;
  t177 = t46*t151;
  t178 = t177*t155;
  t187 = t57*t151;
  t188 = t187*t155;
  t189 = t5*t188;
  t190 = t61*t151;
  t191 = t190*t155;
  t193 = t65*t151;
  t194 = t193*t155;
  t195 = t5*t194;
  t197 = t25*t188;
  t199 = t25*t194;
  t201 = t34*t188;
  t203 = t34*t194;
  t205 = t41*t188;
  t207 = t41*t194;
  t209 = t82*t151;
  t210 = t209*t155;
  t220 = (-zeta+1.0/3.0)*t154;
  t221 = t152*t220;
  t222 = t5*t221;
  t223 = t158*t220;
  t225 = t161*t220;
  t226 = t5*t225;
  t228 = t25*t221;
  t230 = t25*t225;
  t232 = t34*t221;
  t234 = t34*t225;
  t236 = t41*t221;
  t238 = t41*t225;
  t240 = t177*t220;
  t249 = t187*t220;
  t250 = t5*t249;
  t251 = t190*t220;
  t253 = t193*t220;
  t254 = t5*t253;
  t256 = t25*t249;
  t258 = t25*t253;
  t260 = t34*t249;
  t262 = t34*t253;
  t264 = t41*t249;
  t266 = t41*t253;
  t268 = t209*t220;

  values[0] = -729.0/2048.0*t13-729.0/2048.0*t5*t16-729.0/2048.0*t21;
  values[1] = 2187.0/2048.0*t26+2187.0/2048.0*t25*t16+2187.0/2048.0*t28;
  values[2] = 2187.0/2048.0*t35+2187.0/2048.0*t34*t16+2187.0/2048.0*t37;
  values[3] = -729.0/2048.0*t42-729.0/2048.0*t41*t16-729.0/2048.0*t44;
  values[4] = 2187.0/2048.0*t13+2187.0/2048.0*t5*t48+2187.0/2048.0*t21;
  values[5] = -6561.0/2048.0*t26-6561.0/2048.0*t25*t48-6561.0/2048.0*t28;
  values[6] = -6561.0/2048.0*t35-6561.0/2048.0*t34*t48-6561.0/2048.0*t37;
  values[7] = 2187.0/2048.0*t42+2187.0/2048.0*t41*t48+2187.0/2048.0*t44;
  values[8] = 2187.0/2048.0*t60+2187.0/2048.0*t5*t63+2187.0/2048.0*t68;
  values[9] = -6561.0/2048.0*t70-6561.0/2048.0*t25*t63-6561.0/2048.0*t72;
  values[10] = -6561.0/2048.0*t74-6561.0/2048.0*t34*t63-6561.0/2048.0*t76;
  values[11] = 2187.0/2048.0*t78+2187.0/2048.0*t41*t63+2187.0/2048.0*t80;
  values[12] = -729.0/2048.0*t60-729.0/2048.0*t5*t84-729.0/2048.0*t68;
  values[13] = 2187.0/2048.0*t70+2187.0/2048.0*t25*t84+2187.0/2048.0*t72;
  values[14] = 2187.0/2048.0*t74+2187.0/2048.0*t34*t84+2187.0/2048.0*t76;
  values[15] = -729.0/2048.0*t78-729.0/2048.0*t41*t84-729.0/2048.0*t80;
  values[16] = 2187.0/2048.0*t96+2187.0/2048.0*t5*t97+2187.0/2048.0*t100;
  values[17] = -6561.0/2048.0*t102-6561.0/2048.0*t25*t97-6561.0/2048.0*t104;
  values[18] = -6561.0/2048.0*t106-6561.0/2048.0*t34*t97-6561.0/2048.0*t108;
  values[19] = 2187.0/2048.0*t110+2187.0/2048.0*t41*t97+2187.0/2048.0*t112;
  values[20] = -6561.0/2048.0*t96-6561.0/2048.0*t5*t114-6561.0/2048.0*t100;
  values[21] = 19683.0/2048.0*t102+19683.0/2048.0*t25*t114+19683.0/2048.0*t104;
  values[22] = 19683.0/2048.0*t106+19683.0/2048.0*t34*t114+19683.0/2048.0*t108;
  values[23] = -6561.0/2048.0*t110-6561.0/2048.0*t41*t114-6561.0/2048.0*t112;
  values[24] = -6561.0/2048.0*t124-6561.0/2048.0*t5*t125-6561.0/2048.0*t128;
  values[25] = 19683.0/2048.0*t130+19683.0/2048.0*t25*t125+19683.0/2048.0*t132;
  values[26] = 19683.0/2048.0*t134+19683.0/2048.0*t34*t125+19683.0/2048.0*t136;
  values[27] = -6561.0/2048.0*t138-6561.0/2048.0*t41*t125-6561.0/2048.0*t140;
  values[28] = 2187.0/2048.0*t124+2187.0/2048.0*t5*t142+2187.0/2048.0*t128;
  values[29] = -6561.0/2048.0*t130-6561.0/2048.0*t25*t142-6561.0/2048.0*t132;
  values[30] = -6561.0/2048.0*t134-6561.0/2048.0*t34*t142-6561.0/2048.0*t136;
  values[31] = 2187.0/2048.0*t138+2187.0/2048.0*t41*t142+2187.0/2048.0*t140;
  values[32] = 2187.0/2048.0*t157+2187.0/2048.0*t5*t159+2187.0/2048.0*t163;
  values[33] = -6561.0/2048.0*t165-6561.0/2048.0*t25*t159-6561.0/2048.0*t167;
  values[34] = -6561.0/2048.0*t169-6561.0/2048.0*t34*t159-6561.0/2048.0*t171;
  values[35] = 2187.0/2048.0*t173+2187.0/2048.0*t41*t159+2187.0/2048.0*t175;
  values[36] = -6561.0/2048.0*t157-6561.0/2048.0*t5*t178-6561.0/2048.0*t163;
  values[37] = 19683.0/2048.0*t165+19683.0/2048.0*t25*t178+19683.0/2048.0*t167;
  values[38] = 19683.0/2048.0*t169+19683.0/2048.0*t34*t178+19683.0/2048.0*t171;
  values[39] = -6561.0/2048.0*t173-6561.0/2048.0*t41*t178-6561.0/2048.0*t175;
  values[40] = -6561.0/2048.0*t189-6561.0/2048.0*t5*t191-6561.0/2048.0*t195;
  values[41] = 19683.0/2048.0*t197+19683.0/2048.0*t25*t191+19683.0/2048.0*t199;
  values[42] = 19683.0/2048.0*t201+19683.0/2048.0*t34*t191+19683.0/2048.0*t203;
  values[43] = -6561.0/2048.0*t205-6561.0/2048.0*t41*t191-6561.0/2048.0*t207;
  values[44] = 2187.0/2048.0*t189+2187.0/2048.0*t5*t210+2187.0/2048.0*t195;
  values[45] = -6561.0/2048.0*t197-6561.0/2048.0*t25*t210-6561.0/2048.0*t199;
  values[46] = -6561.0/2048.0*t201-6561.0/2048.0*t34*t210-6561.0/2048.0*t203;
  values[47] = 2187.0/2048.0*t205+2187.0/2048.0*t41*t210+2187.0/2048.0*t207;
  values[48] = -729.0/2048.0*t222-729.0/2048.0*t5*t223-729.0/2048.0*t226;
  values[49] = 2187.0/2048.0*t228+2187.0/2048.0*t25*t223+2187.0/2048.0*t230;
  values[50] = 2187.0/2048.0*t232+2187.0/2048.0*t34*t223+2187.0/2048.0*t234;
  values[51] = -729.0/2048.0*t236-729.0/2048.0*t41*t223-729.0/2048.0*t238;
  values[52] = 2187.0/2048.0*t222+2187.0/2048.0*t5*t240+2187.0/2048.0*t226;
  values[53] = -6561.0/2048.0*t228-6561.0/2048.0*t25*t240-6561.0/2048.0*t230;
  values[54] = -6561.0/2048.0*t232-6561.0/2048.0*t34*t240-6561.0/2048.0*t234;
  values[55] = 2187.0/2048.0*t236+2187.0/2048.0*t41*t240+2187.0/2048.0*t238;
  values[56] = 2187.0/2048.0*t250+2187.0/2048.0*t5*t251+2187.0/2048.0*t254;
  values[57] = -6561.0/2048.0*t256-6561.0/2048.0*t25*t251-6561.0/2048.0*t258;
  values[58] = -6561.0/2048.0*t260-6561.0/2048.0*t34*t251-6561.0/2048.0*t262;
  values[59] = 2187.0/2048.0*t264+2187.0/2048.0*t41*t251+2187.0/2048.0*t266;
  values[60] = -729.0/2048.0*t250-729.0/2048.0*t5*t268-729.0/2048.0*t254;
  values[61] = 2187.0/2048.0*t256+2187.0/2048.0*t25*t268+2187.0/2048.0*t258;
  values[62] = 2187.0/2048.0*t260+2187.0/2048.0*t34*t268+2187.0/2048.0*t262;
  values[63] = -729.0/2048.0*t264-729.0/2048.0*t41*t268-729.0/2048.0*t266;
}

static void C_H_Q3_3D_DeriveEtaZeta(double xi, double eta, double zeta,
                             double *values)
{
  double t1, t4, t5, t6, t7, t8, t9, t10, t11, t12, t14, t15, t16, t17;
  double t18, t19, t21, t22, t23, t24, t25, t26, t27, t28, t29, t30, t32;
  double t33, t34, t39, t41, t43, t44, t45, t47, t50, t53, t54, t56, t58;
  double t59, t60, t62, t67, t69, t71, t72, t73, t75, t78, t79, t80, t82;
  double t83, t84, t86, t87, t89, t90, t91, t95, t98, t102, t105, t109;
  double t112, t115, t116, t117, t118, t120, t121, t122, t124, t125, t126;
  double t127, t128, t129, t130, t131, t132, t133, t135, t136, t137, t141;
  double t143, t144, t145, t147, t151, t153, t154, t155, t157, t161, t163;
  double t164, t165, t167, t170, t171, t172, t174, t175, t176, t178, t179;
  double t181, t182, t183, t187, t190, t194, t197, t201, t204, t207, t208;
  double t209, t211, t212, t214, t215, t216, t217, t218, t220, t225, t226;
  double t232, t233, t239, t240, t244, t246, t248, t250, t268, t270, t272;
  double t273, t274, t275, t276, t278, t283, t284, t290, t291, t297, t298;
  double t302, t304, t306, t308, t326, t327, t328, t329, t331, t332, t333;
  double t334, t335, t336, t338, t339, t340, t341, t342, t343, t344, t346;
  double t347, t348, t352, t354, t355, t356, t358, t362, t364, t365, t366;
  double t368, t372, t374, t375, t376, t378, t381, t383, t384, t385, t387;
  double t389, t390, t391, t395, t398, t402, t405, t409, t412, t415, t417;
  double t418, t419, t421, t422, t423, t424, t425, t426, t427, t429, t430;
  double t431, t435, t437, t438, t439, t441, t445, t447, t448, t449, t451;
  double t455, t457, t458, t459, t461, t464, t466, t467, t468, t470, t472;
  double t473, t474, t478, t481, t485, t488, t492, t495, t498, t499, t500;
  double t502, t503, t505, t506, t507, t508, t509, t511, t516, t517, t523;
  double t524, t530, t531, t535, t537, t539, t541, t559, t561, t563, t564;
  double t565, t566, t567, t569, t574, t575, t581, t582, t588, t589, t593;
  double t595, t597, t599;

  t1 = xi-1.0;
  t4 = xi-1.0/3.0;
  t5 = t1*(xi+1.0/3.0)*t4;
  t6 = eta+1.0/3.0;
  t7 = eta-1.0/3.0;
  t8 = t6*t7;
  t9 = zeta+1.0/3.0;
  t10 = zeta-1.0/3.0;
  t11 = t9*t10;
  t12 = t8*t11;
  t14 = zeta-1.0;
  t15 = t14*t10;
  t16 = t8*t15;
  t17 = t5*t16;
  t18 = t14*t9;
  t19 = t8*t18;
  t21 = eta-1.0;
  t22 = t21*t7;
  t23 = t22*t11;
  t24 = t5*t23;
  t25 = t22*t15;
  t26 = t5*t25;
  t27 = t22*t18;
  t28 = t5*t27;
  t29 = t21*t6;
  t30 = t29*t11;
  t32 = t29*t15;
  t33 = t5*t32;
  t34 = t29*t18;
  t39 = t1*(xi+1.0)*t4;
  t41 = t39*t16;
  t43 = t39*t23;
  t44 = t39*t25;
  t45 = t39*t27;
  t47 = t39*t32;
  t50 = -xi-1.0;
  t53 = -xi-1.0/3.0;
  t54 = t50*(-xi+1.0)*t53;
  t56 = t54*t16;
  t58 = t54*t23;
  t59 = t54*t25;
  t60 = t54*t27;
  t62 = t54*t32;
  t67 = t50*(-xi+1.0/3.0)*t53;
  t69 = t67*t16;
  t71 = t67*t23;
  t72 = t67*t25;
  t73 = t67*t27;
  t75 = t67*t32;
  t78 = eta+1.0;
  t79 = t78*t7;
  t80 = t79*t11;
  t82 = t79*t15;
  t83 = t5*t82;
  t84 = t79*t18;
  t86 = t21*t78;
  t87 = t86*t11;
  t89 = t86*t15;
  t90 = t5*t89;
  t91 = t86*t18;
  t95 = t39*t82;
  t98 = t39*t89;
  t102 = t54*t82;
  t105 = t54*t89;
  t109 = t67*t82;
  t112 = t67*t89;
  t115 = -eta+1.0;
  t116 = -eta-1.0/3.0;
  t117 = t115*t116;
  t118 = t117*t11;
  t120 = t117*t15;
  t121 = t5*t120;
  t122 = t117*t18;
  t124 = -eta-1.0;
  t125 = t124*t116;
  t126 = t125*t11;
  t127 = t5*t126;
  t128 = t125*t15;
  t129 = t5*t128;
  t130 = t125*t18;
  t131 = t5*t130;
  t132 = t124*t115;
  t133 = t132*t11;
  t135 = t132*t15;
  t136 = t5*t135;
  t137 = t132*t18;
  t141 = t39*t120;
  t143 = t39*t126;
  t144 = t39*t128;
  t145 = t39*t130;
  t147 = t39*t135;
  t151 = t54*t120;
  t153 = t54*t126;
  t154 = t54*t128;
  t155 = t54*t130;
  t157 = t54*t135;
  t161 = t67*t120;
  t163 = t67*t126;
  t164 = t67*t128;
  t165 = t67*t130;
  t167 = t67*t135;
  t170 = -eta+1.0/3.0;
  t171 = t170*t116;
  t172 = t171*t11;
  t174 = t171*t15;
  t175 = t5*t174;
  t176 = t171*t18;
  t178 = t124*t170;
  t179 = t178*t11;
  t181 = t178*t15;
  t182 = t5*t181;
  t183 = t178*t18;
  t187 = t39*t174;
  t190 = t39*t181;
  t194 = t54*t174;
  t197 = t54*t181;
  t201 = t67*t174;
  t204 = t67*t181;
  t207 = zeta+1.0;
  t208 = t207*t10;
  t209 = t8*t208;
  t211 = t14*t207;
  t212 = t8*t211;
  t214 = t22*t208;
  t215 = t5*t214;
  t216 = t22*t211;
  t217 = t5*t216;
  t218 = t29*t208;
  t220 = t29*t211;
  t225 = t39*t214;
  t226 = t39*t216;
  t232 = t54*t214;
  t233 = t54*t216;
  t239 = t67*t214;
  t240 = t67*t216;
  t244 = t79*t208;
  t246 = t79*t211;
  t248 = t86*t208;
  t250 = t86*t211;
  t268 = t117*t208;
  t270 = t117*t211;
  t272 = t125*t208;
  t273 = t5*t272;
  t274 = t125*t211;
  t275 = t5*t274;
  t276 = t132*t208;
  t278 = t132*t211;
  t283 = t39*t272;
  t284 = t39*t274;
  t290 = t54*t272;
  t291 = t54*t274;
  t297 = t67*t272;
  t298 = t67*t274;
  t302 = t171*t208;
  t304 = t171*t211;
  t306 = t178*t208;
  t308 = t178*t211;
  t326 = -zeta+1.0;
  t327 = -zeta-1.0/3.0;
  t328 = t326*t327;
  t329 = t8*t328;
  t331 = -zeta-1.0;
  t332 = t331*t327;
  t333 = t8*t332;
  t334 = t5*t333;
  t335 = t331*t326;
  t336 = t8*t335;
  t338 = t22*t328;
  t339 = t5*t338;
  t340 = t22*t332;
  t341 = t5*t340;
  t342 = t22*t335;
  t343 = t5*t342;
  t344 = t29*t328;
  t346 = t29*t332;
  t347 = t5*t346;
  t348 = t29*t335;
  t352 = t39*t333;
  t354 = t39*t338;
  t355 = t39*t340;
  t356 = t39*t342;
  t358 = t39*t346;
  t362 = t54*t333;
  t364 = t54*t338;
  t365 = t54*t340;
  t366 = t54*t342;
  t368 = t54*t346;
  t372 = t67*t333;
  t374 = t67*t338;
  t375 = t67*t340;
  t376 = t67*t342;
  t378 = t67*t346;
  t381 = t79*t328;
  t383 = t79*t332;
  t384 = t5*t383;
  t385 = t79*t335;
  t387 = t86*t328;
  t389 = t86*t332;
  t390 = t5*t389;
  t391 = t86*t335;
  t395 = t39*t383;
  t398 = t39*t389;
  t402 = t54*t383;
  t405 = t54*t389;
  t409 = t67*t383;
  t412 = t67*t389;
  t415 = t117*t328;
  t417 = t117*t332;
  t418 = t5*t417;
  t419 = t117*t335;
  t421 = t125*t328;
  t422 = t5*t421;
  t423 = t125*t332;
  t424 = t5*t423;
  t425 = t125*t335;
  t426 = t5*t425;
  t427 = t132*t328;
  t429 = t132*t332;
  t430 = t5*t429;
  t431 = t132*t335;
  t435 = t39*t417;
  t437 = t39*t421;
  t438 = t39*t423;
  t439 = t39*t425;
  t441 = t39*t429;
  t445 = t54*t417;
  t447 = t54*t421;
  t448 = t54*t423;
  t449 = t54*t425;
  t451 = t54*t429;
  t455 = t67*t417;
  t457 = t67*t421;
  t458 = t67*t423;
  t459 = t67*t425;
  t461 = t67*t429;
  t464 = t171*t328;
  t466 = t171*t332;
  t467 = t5*t466;
  t468 = t171*t335;
  t470 = t178*t328;
  t472 = t178*t332;
  t473 = t5*t472;
  t474 = t178*t335;
  t478 = t39*t466;
  t481 = t39*t472;
  t485 = t54*t466;
  t488 = t54*t472;
  t492 = t67*t466;
  t495 = t67*t472;
  t498 = -zeta+1.0/3.0;
  t499 = t498*t327;
  t500 = t8*t499;
  t502 = t331*t498;
  t503 = t8*t502;
  t505 = t22*t499;
  t506 = t5*t505;
  t507 = t22*t502;
  t508 = t5*t507;
  t509 = t29*t499;
  t511 = t29*t502;
  t516 = t39*t505;
  t517 = t39*t507;
  t523 = t54*t505;
  t524 = t54*t507;
  t530 = t67*t505;
  t531 = t67*t507;
  t535 = t79*t499;
  t537 = t79*t502;
  t539 = t86*t499;
  t541 = t86*t502;
  t559 = t117*t499;
  t561 = t117*t502;
  t563 = t125*t499;
  t564 = t5*t563;
  t565 = t125*t502;
  t566 = t5*t565;
  t567 = t132*t499;
  t569 = t132*t502;
  t574 = t39*t563;
  t575 = t39*t565;
  t581 = t54*t563;
  t582 = t54*t565;
  t588 = t67*t563;
  t589 = t67*t565;
  t593 = t171*t499;
  t595 = t171*t502;
  t597 = t178*t499;
  t599 = t178*t502;

  values[0] = -729.0/4096.0*t5*t12-729.0/4096.0*t17-729.0/4096.0*t5*t19-729.0/4096.0*t24-729.0/4096.0*t26-729.0/4096.0*t28-729.0/4096.0*t5*t30-729.0/4096.0*t33-729.0/4096.0*t5*t34;
  values[1] = 2187.0/4096.0*t39*t12+2187.0/4096.0*t41+2187.0/4096.0*t39*t19+2187.0/4096.0*t43+2187.0/4096.0*t44+2187.0/4096.0*t45+2187.0/4096.0*t39*t30+2187.0/4096.0*t47+2187.0/4096.0*t39*t34;
  values[2] = 2187.0/4096.0*t54*t12+2187.0/4096.0*t56+2187.0/4096.0*t54*t19+2187.0/4096.0*t58+2187.0/4096.0*t59+2187.0/4096.0*t60+2187.0/4096.0*t54*t30+2187.0/4096.0*t62+2187.0/4096.0*t54*t34;
  values[3] = -729.0/4096.0*t67*t12-729.0/4096.0*t69-729.0/4096.0*t67*t19-729.0/4096.0*t71-729.0/4096.0*t72-729.0/4096.0*t73-729.0/4096.0*t67*t30-729.0/4096.0*t75-729.0/4096.0*t67*t34;
  values[4] = 2187.0/4096.0*t5*t80+2187.0/4096.0*t83+2187.0/4096.0*t5*t84+2187.0/4096.0*t24+2187.0/4096.0*t26+2187.0/4096.0*t28+2187.0/4096.0*t5*t87+2187.0/4096.0*t90+2187.0/4096.0*t5*t91;
  values[5] = -6561.0/4096.0*t39*t80-6561.0/4096.0*t95-6561.0/4096.0*t39*t84-6561.0/4096.0*t43-6561.0/4096.0*t44-6561.0/4096.0*t45-6561.0/4096.0*t39*t87-6561.0/4096.0*t98-6561.0/4096.0*t39*t91;
  values[6] = -6561.0/4096.0*t54*t80-6561.0/4096.0*t102-6561.0/4096.0*t54*t84-6561.0/4096.0*t58-6561.0/4096.0*t59-6561.0/4096.0*t60-6561.0/4096.0*t54*t87-6561.0/4096.0*t105-6561.0/4096.0*t54*t91;
  values[7] = 2187.0/4096.0*t67*t80+2187.0/4096.0*t109+2187.0/4096.0*t67*t84+2187.0/4096.0*t71+2187.0/4096.0*t72+2187.0/4096.0*t73+2187.0/4096.0*t67*t87+2187.0/4096.0*t112+2187.0/4096.0*t67*t91;
  values[8] = -2187.0/4096.0*t5*t118-2187.0/4096.0*t121-2187.0/4096.0*t5*t122-2187.0/4096.0*t127-2187.0/4096.0*t129-2187.0/4096.0*t131-2187.0/4096.0*t5*t133-2187.0/4096.0*t136-2187.0/4096.0*t5*t137;
  values[9] = 6561.0/4096.0*t39*t118+6561.0/4096.0*t141+6561.0/4096.0*t39*t122+6561.0/4096.0*t143+6561.0/4096.0*t144+6561.0/4096.0*t145+6561.0/4096.0*t39*t133+6561.0/4096.0*t147+6561.0/4096.0*t39*t137;
  values[10] = 6561.0/4096.0*t54*t118+6561.0/4096.0*t151+6561.0/4096.0*t54*t122+6561.0/4096.0*t153+6561.0/4096.0*t154+6561.0/4096.0*t155+6561.0/4096.0*t54*t133+6561.0/4096.0*t157+6561.0/4096.0*t54*t137;
  values[11] = -2187.0/4096.0*t67*t118-2187.0/4096.0*t161-2187.0/4096.0*t67*t122-2187.0/4096.0*t163-2187.0/4096.0*t164-2187.0/4096.0*t165-2187.0/4096.0*t67*t133-2187.0/4096.0*t167-2187.0/4096.0*t67*t137;
  values[12] = 729.0/4096.0*t5*t172+729.0/4096.0*t175+729.0/4096.0*t5*t176+729.0/4096.0*t127+729.0/4096.0*t129+729.0/4096.0*t131+729.0/4096.0*t5*t179+729.0/4096.0*t182+729.0/4096.0*t5*t183;
  values[13] = -2187.0/4096.0*t39*t172-2187.0/4096.0*t187-2187.0/4096.0*t39*t176-2187.0/4096.0*t143-2187.0/4096.0*t144-2187.0/4096.0*t145-2187.0/4096.0*t39*t179-2187.0/4096.0*t190-2187.0/4096.0*t39*t183;
  values[14] = -2187.0/4096.0*t54*t172-2187.0/4096.0*t194-2187.0/4096.0*t54*t176-2187.0/4096.0*t153-2187.0/4096.0*t154-2187.0/4096.0*t155-2187.0/4096.0*t54*t179-2187.0/4096.0*t197-2187.0/4096.0*t54*t183;
  values[15] = 729.0/4096.0*t67*t172+729.0/4096.0*t201+729.0/4096.0*t67*t176+729.0/4096.0*t163+729.0/4096.0*t164+729.0/4096.0*t165+729.0/4096.0*t67*t179+729.0/4096.0*t204+729.0/4096.0*t67*t183;
  values[16] = 2187.0/4096.0*t5*t209+2187.0/4096.0*t17+2187.0/4096.0*t5*t212+2187.0/4096.0*t215+2187.0/4096.0*t26+2187.0/4096.0*t217+2187.0/4096.0*t5*t218+2187.0/4096.0*t33+2187.0/4096.0*t5*t220;
  values[17] = -6561.0/4096.0*t39*t209-6561.0/4096.0*t41-6561.0/4096.0*t39*t212-6561.0/4096.0*t225-6561.0/4096.0*t44-6561.0/4096.0*t226-6561.0/4096.0*t39*t218-6561.0/4096.0*t47-6561.0/4096.0*t39*t220;
  values[18] = -6561.0/4096.0*t54*t209-6561.0/4096.0*t56-6561.0/4096.0*t54*t212-6561.0/4096.0*t232-6561.0/4096.0*t59-6561.0/4096.0*t233-6561.0/4096.0*t54*t218-6561.0/4096.0*t62-6561.0/4096.0*t54*t220;
  values[19] = 2187.0/4096.0*t67*t209+2187.0/4096.0*t69+2187.0/4096.0*t67*t212+2187.0/4096.0*t239+2187.0/4096.0*t72+2187.0/4096.0*t240+2187.0/4096.0*t67*t218+2187.0/4096.0*t75+2187.0/4096.0*t67*t220;
  values[20] = -6561.0/4096.0*t5*t244-6561.0/4096.0*t83-6561.0/4096.0*t5*t246-6561.0/4096.0*t215-6561.0/4096.0*t26-6561.0/4096.0*t217-6561.0/4096.0*t5*t248-6561.0/4096.0*t90-6561.0/4096.0*t5*t250;
  values[21] = 19683.0/4096.0*t39*t244+19683.0/4096.0*t95+19683.0/4096.0*t39*t246+19683.0/4096.0*t225+19683.0/4096.0*t44+19683.0/4096.0*t226+19683.0/4096.0*t39*t248+19683.0/4096.0*t98+19683.0/4096.0*t39*t250;
  values[22] = 19683.0/4096.0*t54*t244+19683.0/4096.0*t102+19683.0/4096.0*t54*t246+19683.0/4096.0*t232+19683.0/4096.0*t59+19683.0/4096.0*t233+19683.0/4096.0*t54*t248+19683.0/4096.0*t105+19683.0/4096.0*t54*t250;
  values[23] = -6561.0/4096.0*t67*t244-6561.0/4096.0*t109-6561.0/4096.0*t67*t246-6561.0/4096.0*t239-6561.0/4096.0*t72-6561.0/4096.0*t240-6561.0/4096.0*t67*t248-6561.0/4096.0*t112-6561.0/4096.0*t67*t250;
  values[24] = 6561.0/4096.0*t5*t268+6561.0/4096.0*t121+6561.0/4096.0*t5*t270+6561.0/4096.0*t273+6561.0/4096.0*t129+6561.0/4096.0*t275+6561.0/4096.0*t5*t276+6561.0/4096.0*t136+6561.0/4096.0*t5*t278;
  values[25] = -19683.0/4096.0*t39*t268-19683.0/4096.0*t141-19683.0/4096.0*t39*t270-19683.0/4096.0*t283-19683.0/4096.0*t144-19683.0/4096.0*t284-19683.0/4096.0*t39*t276-19683.0/4096.0*t147-19683.0/4096.0*t39*t278;
  values[26] = -19683.0/4096.0*t54*t268-19683.0/4096.0*t151-19683.0/4096.0*t54*t270-19683.0/4096.0*t290-19683.0/4096.0*t154-19683.0/4096.0*t291-19683.0/4096.0*t54*t276-19683.0/4096.0*t157-19683.0/4096.0*t54*t278;
  values[27] = 6561.0/4096.0*t67*t268+6561.0/4096.0*t161+6561.0/4096.0*t67*t270+6561.0/4096.0*t297+6561.0/4096.0*t164+6561.0/4096.0*t298+6561.0/4096.0*t67*t276+6561.0/4096.0*t167+6561.0/4096.0*t67*t278;
  values[28] = -2187.0/4096.0*t5*t302-2187.0/4096.0*t175-2187.0/4096.0*t5*t304-2187.0/4096.0*t273-2187.0/4096.0*t129-2187.0/4096.0*t275-2187.0/4096.0*t5*t306-2187.0/4096.0*t182-2187.0/4096.0*t5*t308;
  values[29] = 6561.0/4096.0*t39*t302+6561.0/4096.0*t187+6561.0/4096.0*t39*t304+6561.0/4096.0*t283+6561.0/4096.0*t144+6561.0/4096.0*t284+6561.0/4096.0*t39*t306+6561.0/4096.0*t190+6561.0/4096.0*t39*t308;
  values[30] = 6561.0/4096.0*t54*t302+6561.0/4096.0*t194+6561.0/4096.0*t54*t304+6561.0/4096.0*t290+6561.0/4096.0*t154+6561.0/4096.0*t291+6561.0/4096.0*t54*t306+6561.0/4096.0*t197+6561.0/4096.0*t54*t308;
  values[31] = -2187.0/4096.0*t67*t302-2187.0/4096.0*t201-2187.0/4096.0*t67*t304-2187.0/4096.0*t297-2187.0/4096.0*t164-2187.0/4096.0*t298-2187.0/4096.0*t67*t306-2187.0/4096.0*t204-2187.0/4096.0*t67*t308;
  values[32] = -2187.0/4096.0*t5*t329-2187.0/4096.0*t334-2187.0/4096.0*t5*t336-2187.0/4096.0*t339-2187.0/4096.0*t341-2187.0/4096.0*t343-2187.0/4096.0*t5*t344-2187.0/4096.0*t347-2187.0/4096.0*t5*t348;
  values[33] = 6561.0/4096.0*t39*t329+6561.0/4096.0*t352+6561.0/4096.0*t39*t336+6561.0/4096.0*t354+6561.0/4096.0*t355+6561.0/4096.0*t356+6561.0/4096.0*t39*t344+6561.0/4096.0*t358+6561.0/4096.0*t39*t348;
  values[34] = 6561.0/4096.0*t54*t329+6561.0/4096.0*t362+6561.0/4096.0*t54*t336+6561.0/4096.0*t364+6561.0/4096.0*t365+6561.0/4096.0*t366+6561.0/4096.0*t54*t344+6561.0/4096.0*t368+6561.0/4096.0*t54*t348;
  values[35] = -2187.0/4096.0*t67*t329-2187.0/4096.0*t372-2187.0/4096.0*t67*t336-2187.0/4096.0*t374-2187.0/4096.0*t375-2187.0/4096.0*t376-2187.0/4096.0*t67*t344-2187.0/4096.0*t378-2187.0/4096.0*t67*t348;
  values[36] = 6561.0/4096.0*t5*t381+6561.0/4096.0*t384+6561.0/4096.0*t5*t385+6561.0/4096.0*t339+6561.0/4096.0*t341+6561.0/4096.0*t343+6561.0/4096.0*t5*t387+6561.0/4096.0*t390+6561.0/4096.0*t5*t391;
  values[37] = -19683.0/4096.0*t39*t381-19683.0/4096.0*t395-19683.0/4096.0*t39*t385-19683.0/4096.0*t354-19683.0/4096.0*t355-19683.0/4096.0*t356-19683.0/4096.0*t39*t387-19683.0/4096.0*t398-19683.0/4096.0*t39*t391;
  values[38] = -19683.0/4096.0*t54*t381-19683.0/4096.0*t402-19683.0/4096.0*t54*t385-19683.0/4096.0*t364-19683.0/4096.0*t365-19683.0/4096.0*t366-19683.0/4096.0*t54*t387-19683.0/4096.0*t405-19683.0/4096.0*t54*t391;
  values[39] = 6561.0/4096.0*t67*t381+6561.0/4096.0*t409+6561.0/4096.0*t67*t385+6561.0/4096.0*t374+6561.0/4096.0*t375+6561.0/4096.0*t376+6561.0/4096.0*t67*t387+6561.0/4096.0*t412+6561.0/4096.0*t67*t391;
  values[40] = -6561.0/4096.0*t5*t415-6561.0/4096.0*t418-6561.0/4096.0*t5*t419-6561.0/4096.0*t422-6561.0/4096.0*t424-6561.0/4096.0*t426-6561.0/4096.0*t5*t427-6561.0/4096.0*t430-6561.0/4096.0*t5*t431;
  values[41] = 19683.0/4096.0*t39*t415+19683.0/4096.0*t435+19683.0/4096.0*t39*t419+19683.0/4096.0*t437+19683.0/4096.0*t438+19683.0/4096.0*t439+19683.0/4096.0*t39*t427+19683.0/4096.0*t441+19683.0/4096.0*t39*t431;
  values[42] = 19683.0/4096.0*t54*t415+19683.0/4096.0*t445+19683.0/4096.0*t54*t419+19683.0/4096.0*t447+19683.0/4096.0*t448+19683.0/4096.0*t449+19683.0/4096.0*t54*t427+19683.0/4096.0*t451+19683.0/4096.0*t54*t431;
  values[43] = -6561.0/4096.0*t67*t415-6561.0/4096.0*t455-6561.0/4096.0*t67*t419-6561.0/4096.0*t457-6561.0/4096.0*t458-6561.0/4096.0*t459-6561.0/4096.0*t67*t427-6561.0/4096.0*t461-6561.0/4096.0*t67*t431;
  values[44] = 2187.0/4096.0*t5*t464+2187.0/4096.0*t467+2187.0/4096.0*t5*t468+2187.0/4096.0*t422+2187.0/4096.0*t424+2187.0/4096.0*t426+2187.0/4096.0*t5*t470+2187.0/4096.0*t473+2187.0/4096.0*t5*t474;
  values[45] = -6561.0/4096.0*t39*t464-6561.0/4096.0*t478-6561.0/4096.0*t39*t468-6561.0/4096.0*t437-6561.0/4096.0*t438-6561.0/4096.0*t439-6561.0/4096.0*t39*t470-6561.0/4096.0*t481-6561.0/4096.0*t39*t474;
  values[46] = -6561.0/4096.0*t54*t464-6561.0/4096.0*t485-6561.0/4096.0*t54*t468-6561.0/4096.0*t447-6561.0/4096.0*t448-6561.0/4096.0*t449-6561.0/4096.0*t54*t470-6561.0/4096.0*t488-6561.0/4096.0*t54*t474;
  values[47] = 2187.0/4096.0*t67*t464+2187.0/4096.0*t492+2187.0/4096.0*t67*t468+2187.0/4096.0*t457+2187.0/4096.0*t458+2187.0/4096.0*t459+2187.0/4096.0*t67*t470+2187.0/4096.0*t495+2187.0/4096.0*t67*t474;
  values[48] = 729.0/4096.0*t5*t500+729.0/4096.0*t334+729.0/4096.0*t5*t503+729.0/4096.0*t506+729.0/4096.0*t341+729.0/4096.0*t508+729.0/4096.0*t5*t509+729.0/4096.0*t347+729.0/4096.0*t5*t511;
  values[49] = -2187.0/4096.0*t39*t500-2187.0/4096.0*t352-2187.0/4096.0*t39*t503-2187.0/4096.0*t516-2187.0/4096.0*t355-2187.0/4096.0*t517-2187.0/4096.0*t39*t509-2187.0/4096.0*t358-2187.0/4096.0*t39*t511;
  values[50] = -2187.0/4096.0*t54*t500-2187.0/4096.0*t362-2187.0/4096.0*t54*t503-2187.0/4096.0*t523-2187.0/4096.0*t365-2187.0/4096.0*t524-2187.0/4096.0*t54*t509-2187.0/4096.0*t368-2187.0/4096.0*t54*t511;
  values[51] = 729.0/4096.0*t67*t500+729.0/4096.0*t372+729.0/4096.0*t67*t503+729.0/4096.0*t530+729.0/4096.0*t375+729.0/4096.0*t531+729.0/4096.0*t67*t509+729.0/4096.0*t378+729.0/4096.0*t67*t511;
  values[52] = -2187.0/4096.0*t5*t535-2187.0/4096.0*t384-2187.0/4096.0*t5*t537-2187.0/4096.0*t506-2187.0/4096.0*t341-2187.0/4096.0*t508-2187.0/4096.0*t5*t539-2187.0/4096.0*t390-2187.0/4096.0*t5*t541;
  values[53] = 6561.0/4096.0*t39*t535+6561.0/4096.0*t395+6561.0/4096.0*t39*t537+6561.0/4096.0*t516+6561.0/4096.0*t355+6561.0/4096.0*t517+6561.0/4096.0*t39*t539+6561.0/4096.0*t398+6561.0/4096.0*t39*t541;
  values[54] = 6561.0/4096.0*t54*t535+6561.0/4096.0*t402+6561.0/4096.0*t54*t537+6561.0/4096.0*t523+6561.0/4096.0*t365+6561.0/4096.0*t524+6561.0/4096.0*t54*t539+6561.0/4096.0*t405+6561.0/4096.0*t54*t541;
  values[55] = -2187.0/4096.0*t67*t535-2187.0/4096.0*t409-2187.0/4096.0*t67*t537-2187.0/4096.0*t530-2187.0/4096.0*t375-2187.0/4096.0*t531-2187.0/4096.0*t67*t539-2187.0/4096.0*t412-2187.0/4096.0*t67*t541;
  values[56] = 2187.0/4096.0*t5*t559+2187.0/4096.0*t418+2187.0/4096.0*t5*t561+2187.0/4096.0*t564+2187.0/4096.0*t424+2187.0/4096.0*t566+2187.0/4096.0*t5*t567+2187.0/4096.0*t430+2187.0/4096.0*t5*t569;
  values[57] = -6561.0/4096.0*t39*t559-6561.0/4096.0*t435-6561.0/4096.0*t39*t561-6561.0/4096.0*t574-6561.0/4096.0*t438-6561.0/4096.0*t575-6561.0/4096.0*t39*t567-6561.0/4096.0*t441-6561.0/4096.0*t39*t569;
  values[58] = -6561.0/4096.0*t54*t559-6561.0/4096.0*t445-6561.0/4096.0*t54*t561-6561.0/4096.0*t581-6561.0/4096.0*t448-6561.0/4096.0*t582-6561.0/4096.0*t54*t567-6561.0/4096.0*t451-6561.0/4096.0*t54*t569;
  values[59] = 2187.0/4096.0*t67*t559+2187.0/4096.0*t455+2187.0/4096.0*t67*t561+2187.0/4096.0*t588+2187.0/4096.0*t458+2187.0/4096.0*t589+2187.0/4096.0*t67*t567+2187.0/4096.0*t461+2187.0/4096.0*t67*t569;
  values[60] = -729.0/4096.0*t5*t593-729.0/4096.0*t467-729.0/4096.0*t5*t595-729.0/4096.0*t564-729.0/4096.0*t424-729.0/4096.0*t566-729.0/4096.0*t5*t597-729.0/4096.0*t473-729.0/4096.0*t5*t599;
  values[61] = 2187.0/4096.0*t39*t593+2187.0/4096.0*t478+2187.0/4096.0*t39*t595+2187.0/4096.0*t574+2187.0/4096.0*t438+2187.0/4096.0*t575+2187.0/4096.0*t39*t597+2187.0/4096.0*t481+2187.0/4096.0*t39*t599;
  values[62] = 2187.0/4096.0*t54*t593+2187.0/4096.0*t485+2187.0/4096.0*t54*t595+2187.0/4096.0*t581+2187.0/4096.0*t448+2187.0/4096.0*t582+2187.0/4096.0*t54*t597+2187.0/4096.0*t488+2187.0/4096.0*t54*t599;
  values[63] = -729.0/4096.0*t67*t593-729.0/4096.0*t492-729.0/4096.0*t67*t595-729.0/4096.0*t588-729.0/4096.0*t458-729.0/4096.0*t589-729.0/4096.0*t67*t597-729.0/4096.0*t495-729.0/4096.0*t67*t599;
}

static void C_H_Q3_3D_DeriveZetaZeta(double xi, double eta, double zeta,
                             double *values)
{
  double t1, t4, t5, t6, t8, t9, t10, t11, t12, t13, t14, t15, t16, t18;
  double t19, t20, t21, t25, t26, t28, t30, t33, t34, t35, t37, t41, t42;
  double t44, t47, t48, t49, t50, t52, t53, t55, t57, t59, t61, t63, t65;
  double t67, t69, t70, t71, t72, t73, t74, t75, t77, t78, t79, t81, t83;
  double t85, t87, t89, t91, t94, t95, t96, t97, t99, t100, t102, t104;
  double t106, t108, t110, t112, t114, t115, t116, t125, t134, t135, t144;
  double t153, t154, t155, t156, t157, t158, t159, t161, t162, t163, t164;
  double t166, t168, t170, t172, t174, t176, t178, t179, t180, t182, t183;
  double t185, t187, t189, t191, t193, t195, t197, t198, t199, t200, t201;
  double t203, t204, t205, t207, t209, t211, t213, t215, t217, t219, t220;
  double t221, t223, t224, t226, t228, t230, t232, t234, t236, t238, t239;
  double t240, t249, t258, t259, t268;

  t1 = xi-1.0;
  t4 = xi-1.0/3.0;
  t5 = t1*(xi+1.0/3.0)*t4;
  t6 = eta-1.0;
  t8 = t6*(eta+1.0/3.0);
  t9 = eta-1.0/3.0;
  t10 = zeta-1.0/3.0;
  t11 = t9*t10;
  t12 = t8*t11;
  t13 = t5*t12;
  t14 = zeta+1.0/3.0;
  t15 = t9*t14;
  t16 = t8*t15;
  t18 = zeta-1.0;
  t19 = t9*t18;
  t20 = t8*t19;
  t21 = t5*t20;
  t25 = t1*(xi+1.0)*t4;
  t26 = t25*t12;
  t28 = t25*t20;
  t30 = -xi-1.0;
  t33 = -xi-1.0/3.0;
  t34 = t30*(-xi+1.0)*t33;
  t35 = t34*t12;
  t37 = t34*t20;
  t41 = t30*(-xi+1.0/3.0)*t33;
  t42 = t41*t12;
  t44 = t41*t20;
  t47 = t6*(eta+1.0);
  t48 = t47*t11;
  t49 = t5*t48;
  t50 = t47*t15;
  t52 = t47*t19;
  t53 = t5*t52;
  t55 = t25*t48;
  t57 = t25*t52;
  t59 = t34*t48;
  t61 = t34*t52;
  t63 = t41*t48;
  t65 = t41*t52;
  t67 = -eta-1.0;
  t69 = t67*(-eta+1.0);
  t70 = -eta-1.0/3.0;
  t71 = t70*t10;
  t72 = t69*t71;
  t73 = t5*t72;
  t74 = t70*t14;
  t75 = t69*t74;
  t77 = t70*t18;
  t78 = t69*t77;
  t79 = t5*t78;
  t81 = t25*t72;
  t83 = t25*t78;
  t85 = t34*t72;
  t87 = t34*t78;
  t89 = t41*t72;
  t91 = t41*t78;
  t94 = t67*(-eta+1.0/3.0);
  t95 = t94*t71;
  t96 = t5*t95;
  t97 = t94*t74;
  t99 = t94*t77;
  t100 = t5*t99;
  t102 = t25*t95;
  t104 = t25*t99;
  t106 = t34*t95;
  t108 = t34*t99;
  t110 = t41*t95;
  t112 = t41*t99;
  t114 = zeta+1.0;
  t115 = t9*t114;
  t116 = t8*t115;
  t125 = t47*t115;
  t134 = t70*t114;
  t135 = t69*t134;
  t144 = t94*t134;
  t153 = -zeta-1.0/3.0;
  t154 = t9*t153;
  t155 = t8*t154;
  t156 = t5*t155;
  t157 = -zeta+1.0;
  t158 = t9*t157;
  t159 = t8*t158;
  t161 = -zeta-1.0;
  t162 = t9*t161;
  t163 = t8*t162;
  t164 = t5*t163;
  t166 = t25*t155;
  t168 = t25*t163;
  t170 = t34*t155;
  t172 = t34*t163;
  t174 = t41*t155;
  t176 = t41*t163;
  t178 = t47*t154;
  t179 = t5*t178;
  t180 = t47*t158;
  t182 = t47*t162;
  t183 = t5*t182;
  t185 = t25*t178;
  t187 = t25*t182;
  t189 = t34*t178;
  t191 = t34*t182;
  t193 = t41*t178;
  t195 = t41*t182;
  t197 = t70*t153;
  t198 = t69*t197;
  t199 = t5*t198;
  t200 = t70*t157;
  t201 = t69*t200;
  t203 = t70*t161;
  t204 = t69*t203;
  t205 = t5*t204;
  t207 = t25*t198;
  t209 = t25*t204;
  t211 = t34*t198;
  t213 = t34*t204;
  t215 = t41*t198;
  t217 = t41*t204;
  t219 = t94*t197;
  t220 = t5*t219;
  t221 = t94*t200;
  t223 = t94*t203;
  t224 = t5*t223;
  t226 = t25*t219;
  t228 = t25*t223;
  t230 = t34*t219;
  t232 = t34*t223;
  t234 = t41*t219;
  t236 = t41*t223;
  t238 = -zeta+1.0/3.0;
  t239 = t9*t238;
  t240 = t8*t239;
  t249 = t47*t239;
  t258 = t70*t238;
  t259 = t69*t258;
  t268 = t94*t258;

  values[0] = -729.0/2048.0*t13-729.0/2048.0*t5*t16-729.0/2048.0*t21;
  values[1] = 2187.0/2048.0*t26+2187.0/2048.0*t25*t16+2187.0/2048.0*t28;
  values[2] = 2187.0/2048.0*t35+2187.0/2048.0*t34*t16+2187.0/2048.0*t37;
  values[3] = -729.0/2048.0*t42-729.0/2048.0*t41*t16-729.0/2048.0*t44;
  values[4] = 2187.0/2048.0*t49+2187.0/2048.0*t5*t50+2187.0/2048.0*t53;
  values[5] = -6561.0/2048.0*t55-6561.0/2048.0*t25*t50-6561.0/2048.0*t57;
  values[6] = -6561.0/2048.0*t59-6561.0/2048.0*t34*t50-6561.0/2048.0*t61;
  values[7] = 2187.0/2048.0*t63+2187.0/2048.0*t41*t50+2187.0/2048.0*t65;
  values[8] = 2187.0/2048.0*t73+2187.0/2048.0*t5*t75+2187.0/2048.0*t79;
  values[9] = -6561.0/2048.0*t81-6561.0/2048.0*t25*t75-6561.0/2048.0*t83;
  values[10] = -6561.0/2048.0*t85-6561.0/2048.0*t34*t75-6561.0/2048.0*t87;
  values[11] = 2187.0/2048.0*t89+2187.0/2048.0*t41*t75+2187.0/2048.0*t91;
  values[12] = -729.0/2048.0*t96-729.0/2048.0*t5*t97-729.0/2048.0*t100;
  values[13] = 2187.0/2048.0*t102+2187.0/2048.0*t25*t97+2187.0/2048.0*t104;
  values[14] = 2187.0/2048.0*t106+2187.0/2048.0*t34*t97+2187.0/2048.0*t108;
  values[15] = -729.0/2048.0*t110-729.0/2048.0*t41*t97-729.0/2048.0*t112;
  values[16] = 2187.0/2048.0*t13+2187.0/2048.0*t5*t116+2187.0/2048.0*t21;
  values[17] = -6561.0/2048.0*t26-6561.0/2048.0*t25*t116-6561.0/2048.0*t28;
  values[18] = -6561.0/2048.0*t35-6561.0/2048.0*t34*t116-6561.0/2048.0*t37;
  values[19] = 2187.0/2048.0*t42+2187.0/2048.0*t41*t116+2187.0/2048.0*t44;
  values[20] = -6561.0/2048.0*t49-6561.0/2048.0*t5*t125-6561.0/2048.0*t53;
  values[21] = 19683.0/2048.0*t55+19683.0/2048.0*t25*t125+19683.0/2048.0*t57;
  values[22] = 19683.0/2048.0*t59+19683.0/2048.0*t34*t125+19683.0/2048.0*t61;
  values[23] = -6561.0/2048.0*t63-6561.0/2048.0*t41*t125-6561.0/2048.0*t65;
  values[24] = -6561.0/2048.0*t73-6561.0/2048.0*t5*t135-6561.0/2048.0*t79;
  values[25] = 19683.0/2048.0*t81+19683.0/2048.0*t25*t135+19683.0/2048.0*t83;
  values[26] = 19683.0/2048.0*t85+19683.0/2048.0*t34*t135+19683.0/2048.0*t87;
  values[27] = -6561.0/2048.0*t89-6561.0/2048.0*t41*t135-6561.0/2048.0*t91;
  values[28] = 2187.0/2048.0*t96+2187.0/2048.0*t5*t144+2187.0/2048.0*t100;
  values[29] = -6561.0/2048.0*t102-6561.0/2048.0*t25*t144-6561.0/2048.0*t104;
  values[30] = -6561.0/2048.0*t106-6561.0/2048.0*t34*t144-6561.0/2048.0*t108;
  values[31] = 2187.0/2048.0*t110+2187.0/2048.0*t41*t144+2187.0/2048.0*t112;
  values[32] = 2187.0/2048.0*t156+2187.0/2048.0*t5*t159+2187.0/2048.0*t164;
  values[33] = -6561.0/2048.0*t166-6561.0/2048.0*t25*t159-6561.0/2048.0*t168;
  values[34] = -6561.0/2048.0*t170-6561.0/2048.0*t34*t159-6561.0/2048.0*t172;
  values[35] = 2187.0/2048.0*t174+2187.0/2048.0*t41*t159+2187.0/2048.0*t176;
  values[36] = -6561.0/2048.0*t179-6561.0/2048.0*t5*t180-6561.0/2048.0*t183;
  values[37] = 19683.0/2048.0*t185+19683.0/2048.0*t25*t180+19683.0/2048.0*t187;
  values[38] = 19683.0/2048.0*t189+19683.0/2048.0*t34*t180+19683.0/2048.0*t191;
  values[39] = -6561.0/2048.0*t193-6561.0/2048.0*t41*t180-6561.0/2048.0*t195;
  values[40] = -6561.0/2048.0*t199-6561.0/2048.0*t5*t201-6561.0/2048.0*t205;
  values[41] = 19683.0/2048.0*t207+19683.0/2048.0*t25*t201+19683.0/2048.0*t209;
  values[42] = 19683.0/2048.0*t211+19683.0/2048.0*t34*t201+19683.0/2048.0*t213;
  values[43] = -6561.0/2048.0*t215-6561.0/2048.0*t41*t201-6561.0/2048.0*t217;
  values[44] = 2187.0/2048.0*t220+2187.0/2048.0*t5*t221+2187.0/2048.0*t224;
  values[45] = -6561.0/2048.0*t226-6561.0/2048.0*t25*t221-6561.0/2048.0*t228;
  values[46] = -6561.0/2048.0*t230-6561.0/2048.0*t34*t221-6561.0/2048.0*t232;
  values[47] = 2187.0/2048.0*t234+2187.0/2048.0*t41*t221+2187.0/2048.0*t236;
  values[48] = -729.0/2048.0*t156-729.0/2048.0*t5*t240-729.0/2048.0*t164;
  values[49] = 2187.0/2048.0*t166+2187.0/2048.0*t25*t240+2187.0/2048.0*t168;
  values[50] = 2187.0/2048.0*t170+2187.0/2048.0*t34*t240+2187.0/2048.0*t172;
  values[51] = -729.0/2048.0*t174-729.0/2048.0*t41*t240-729.0/2048.0*t176;
  values[52] = 2187.0/2048.0*t179+2187.0/2048.0*t5*t249+2187.0/2048.0*t183;
  values[53] = -6561.0/2048.0*t185-6561.0/2048.0*t25*t249-6561.0/2048.0*t187;
  values[54] = -6561.0/2048.0*t189-6561.0/2048.0*t34*t249-6561.0/2048.0*t191;
  values[55] = 2187.0/2048.0*t193+2187.0/2048.0*t41*t249+2187.0/2048.0*t195;
  values[56] = 2187.0/2048.0*t199+2187.0/2048.0*t5*t259+2187.0/2048.0*t205;
  values[57] = -6561.0/2048.0*t207-6561.0/2048.0*t25*t259-6561.0/2048.0*t209;
  values[58] = -6561.0/2048.0*t211-6561.0/2048.0*t34*t259-6561.0/2048.0*t213;
  values[59] = 2187.0/2048.0*t215+2187.0/2048.0*t41*t259+2187.0/2048.0*t217;
  values[60] = -729.0/2048.0*t220-729.0/2048.0*t5*t268-729.0/2048.0*t224;
  values[61] = 2187.0/2048.0*t226+2187.0/2048.0*t25*t268+2187.0/2048.0*t228;
  values[62] = 2187.0/2048.0*t230+2187.0/2048.0*t34*t268+2187.0/2048.0*t232;
  values[63] = -729.0/2048.0*t234-729.0/2048.0*t41*t268-729.0/2048.0*t236;
}

TBaseFunct3D *BF_C_H_Q3_3D_Obj = 
new TBaseFunct3D(64, BF_C_H_Q3_3D, BFUnitHexahedron, 
                 C_H_Q3_3D_Funct, C_H_Q3_3D_DeriveXi,
                 C_H_Q3_3D_DeriveEta, C_H_Q3_3D_DeriveZeta,
                 C_H_Q3_3D_DeriveXiXi, C_H_Q3_3D_DeriveXiEta,
                 C_H_Q3_3D_DeriveXiZeta, C_H_Q3_3D_DeriveEtaEta,
                 C_H_Q3_3D_DeriveEtaZeta, C_H_Q3_3D_DeriveZetaZeta,
                 3, 3,
                 0, NULL);