BF_C_Q_UL5S_2D.h

// ***********************************************************************
// Q5 element with bubbles, conforming, 2D
// ***********************************************************************

// base function values
static void C_Q_UL5S_2D_Funct(double xi, double eta, double *values)
{
  double t1, t2, t3, t4, t5, t6, t7, t8, t9, t10, t12, t15, t17, t20, t21;
  double t22, t24, t28, t29, t30, t31, t32, t33, t34, t35, t38, t39, t40, t41;
  double t43, t44, t45, t46, t47, t48, t49, t50, t51, t52, t53, t54, t55, t56;
  double t57, t58, t59, t61, t62, t63, t64, t65, t67, t69, t71, t72, t73, t74;
  double t75, t76, t77, t78, t79, t80, t81, t82, t83, t84, t85, t88, t90, t92;
  double t94, t95, t97, t98, t99, t100, t102, t103, t106, t107, t108, t109;
  double t110, t111, t112, t113, t114, t115, t116, t117, t119, t120, t122;
  double t123, t125, t126, t130, t131, t132, t133, t134, t135, t137, t139;
  double t140, t142, t143, t144, t147, t150, t151, t152, t157, t158, t159;
  double t162, t163, t164, t165, t166, t168, t169, t170, t171, t173, t174;
  double t175, t176, t177, t178, t179, t180, t181, t182, t183, t184, t185;
  double t189, t191, t192, t195, t197, t198, t199, t202, t203, t204, t206;
  double t207, t208, t209, t210, t211, t212, t213, t214, t215, t216, t217;
  double t218, t220, t221, t223, t228, t229, t230, t231, t232, t233, t234;
  double t235, t236, t237, t238, t240, t241, t242, t245, t248, t249, t252;
  double t253, t256, t259, t260, t262, t263, t264, t266, t269, t271, t272;
  double t273, t274, t275, t280, t281, t282, t285, t288, t289, t290, t291;
  double t292, t293, t294, t299, t301, t303, t310, t316, t318, t320, t321;
  double t322, t323, t324, t325, t330, t331, t332, t335, t336, t338, t339;
  double t340, t342, t344, t346, t350, t353, t355, t361, t363, t364, t366;
  double t367, t368, t369, t370, t371, t372, t373, t381, t382, t383, t384;
  double t385, t388, t390, t392, t393, t395, t396, t401, t404, t405, t406;
  double t407, t411, t413, t415, t416, t417, t418, t419, t420, t425, t426;
  double t428, t429, t430, t431, t432, t433, t434, t435, t436, t438, t439;
  double t442, t444, t445, t446, t449, t450, t451, t452, t453, t454, t456;
  double t458, t459, t467, t468, t469, t472, t473, t474, t475, t477, t479;
  double t486, t490, t492, t496, t500;

  t1 = 857.0/3072.0*xi;
  t2 = 857.0/3072.0*eta;
  t3 = xi*eta;
  t4 = 1049.0/2304.0*t3;
  t5 = xi*xi;
  t6 = eta*eta;
  t7 = t6*eta;
  t8 = t5*t7;
  t9 = 4375.0/9216.0*t8;
  t10 = t5*t5;
  t12 = t6*t6;
  t15 = 1.0-t5;
  t17 = t5*xi;
  t20 = 5.0/2.0*t17-3.0/2.0*xi;
  t21 = (1.0-eta)*t15*t20;
  t22 = 125.0/768.0*t21;
  t24 = 1.0-t6;
  t28 = 5.0/2.0*t7-3.0/2.0*eta;
  t29 = (1.0-xi)*t24*t28;
  t30 = 125.0/768.0*t29;
  t31 = t15*t24;
  t32 = t31*t7;
  t33 = 3125.0/3072.0*t32;
  t34 = t31*t17;
  t35 = 3125.0/3072.0*t34;
  t38 = 875.0/3072.0*t7;
  t39 = 31.0/1024.0+t1+t2-t4+t9+625.0/1536.0*t10+625.0/1536.0*t12+t22+t30-t33-t35-1325.0/3072.0*t5-1325.0/3072.0*t6-t38;
  t40 = 875.0/3072.0*t17;
  t41 = t5*t6;
  t43 = t17*eta;
  t44 = 125.0/576.0*t43;
  t45 = t17*t6;
  t46 = 4375.0/9216.0*t45;
  t47 = xi*t7;
  t48 = 125.0/576.0*t47;
  t49 = t17*t7;
  t50 = 625.0/2304.0*t49;
  t51 = xi*t6;
  t52 = 2875.0/9216.0*t51;
  t53 = t5*eta;
  t54 = 2875.0/9216.0*t53;
  t55 = t10*t7;
  t56 = 3125.0/4608.0*t55;
  t57 = xi*t12;
  t58 = 625.0/2304.0*t57;
  t59 = t5*t12;
  t61 = t17*t12;
  t62 = 3125.0/4608.0*t61;
  t63 = t10*eta;
  t64 = 625.0/2304.0*t63;
  t65 = t10*t6;
  t67 = t10*t12;
  t69 = -t40+4525.0/2304.0*t41+t44+t46+t48+t50-t52-t54-t56+t58-15625.0/9216.0*t59-t62+t64-15625.0/9216.0*t65+15625.0/9216.0*t67;
  t71 = 17525.0/12288.0*t53;
  t72 = 625.0/3072.0*t51;
  t73 = 50375.0/6144.0*t41;
  t74 = 875.0/768.0*t43;
  t75 = 625.0/3072.0*t45;
  t76 = 33125.0/12288.0*t8;
  t77 = 1375.0/1024.0*t17;
  t78 = 625.0/4096.0*t7;
  t79 = 1375.0/1024.0*xi;
  t80 = 425.0/4096.0*eta;
  t81 = 15625.0/3072.0*t34;
  t82 = -225.0/2048.0+t71+t72-t73-t74-t75-t76+t77+t78-t79-t80+t81;
  t83 = 625.0/768.0*t21;
  t84 = 15625.0/12288.0*t32;
  t85 = 875.0/768.0*t3;
  t88 = 8125.0/6144.0*t63;
  t90 = 15625.0/6144.0*t55;
  t92 = 78125.0/12288.0*t67;
  t94 = 3125.0/3072.0*t12;
  t95 = -t83+t84+t85+8575.0/4096.0*t5+6625.0/6144.0*t6-t88+21875.0/3072.0*t65+t90+90625.0/12288.0*t59-t92-8125.0/4096.0*t10-t94;
  t97 = 125.0/48.0*xi;
  t98 = 175.0/1536.0*eta;
  t99 = 5375.0/2304.0*t3;
  t100 = 4375.0/1152.0*t8;
  t102 = 625.0/384.0*t21;
  t103 = 15625.0/1536.0*t34;
  t106 = 625.0/1536.0*t7;
  t107 = 125.0/48.0*t17;
  t108 = 175.0/1024.0+t97+t98-t99+t100+11875.0/3072.0*t10+t94+t102-t103-775.0/192.0*t5-1375.0/1536.0*t6-t106-t107;
  t109 = 28625.0/2304.0*t41;
  t110 = 5375.0/2304.0*t43;
  t111 = 8125.0/4608.0*t45;
  t112 = 3125.0/2304.0*t47;
  t113 = 3125.0/2304.0*t49;
  t114 = 8125.0/4608.0*t51;
  t115 = 775.0/288.0*t53;
  t116 = 15625.0/4608.0*t55;
  t117 = 15625.0/4608.0*t57;
  t119 = 15625.0/4608.0*t61;
  t120 = 11875.0/4608.0*t63;
  t122 = 78125.0/9216.0*t67;
  t123 = t109+t110-t111+t112-t113+t114-t115-t116-t117-21875.0/2304.0*t59+t119+t120-53125.0/4608.0*t65+t122;
  t125 = 275.0/3072.0*eta;
  t126 = 6875.0/4608.0*t8;
  t130 = 125.0/1536.0*t6;
  t131 = 625.0/3072.0*t7;
  t132 = -275.0/2048.0-t97-t125+t99-t126-23125.0/6144.0*t10+3125.0/6144.0*t12-t102+t103+11975.0/3072.0*t5-t130-t131+t107;
  t133 = 20125.0/2304.0*t41;
  t134 = 11975.0/4608.0*t53;
  t135 = 15625.0/9216.0*t55;
  t137 = 23125.0/9216.0*t63;
  t139 = 78125.0/18432.0*t67;
  t140 = -t133-t110+t111-t112+t113-t114+t134+t135+t117+34375.0/9216.0*t59-t119-t137+40625.0/4608.0*t65-t139;
  t142 = 325.0/256.0*t53;
  t143 = 1625.0/512.0*t41;
  t144 = 25.0/512.0*eta;
  t147 = 625.0/512.0*t63;
  t150 = 75.0/1024.0-t142-t72+t143+t74+t75-t77+t79+t144-t81+t83-t85-975.0/512.0*t5-125.0/1024.0*t6+t147-3125.0/1024.0*t65+1875.0/1024.0*t10;
  t151 = 2803.0/12288.0*eta;
  t152 = 3125.0/36864.0*t8;
  t157 = (1.0+xi)*t24*t28;
  t158 = 125.0/768.0*t157;
  t159 = 3125.0/12288.0*t32;
  t162 = 2875.0/12288.0*t7;
  t163 = -19.0/1024.0-t1+t151+t4-t152-4375.0/12288.0*t10-625.0/6144.0*t12-t22+t158-t159+t35+4675.0/12288.0*t5+775.0/6144.0*t6-t162;
  t164 = 11575.0/18432.0*t41;
  t165 = 9125.0/36864.0*t53;
  t166 = 3125.0/18432.0*t55;
  t168 = 4375.0/18432.0*t63;
  t169 = 3125.0/9216.0*t65;
  t170 = 15625.0/36864.0*t67;
  t171 = t40-t164-t44-t46-t48-t50+t52+t165-t166-t58+3125.0/36864.0*t59+t62-t168+t169+t170;
  t173 = 625.0/12288.0*t53;
  t174 = 17525.0/12288.0*t51;
  t175 = 11375.0/12288.0*t41;
  t176 = 33125.0/12288.0*t45;
  t177 = 875.0/768.0*t47;
  t178 = 625.0/12288.0*t8;
  t179 = 625.0/4096.0*t17;
  t180 = 4875.0/4096.0*t7;
  t181 = 425.0/4096.0*xi;
  t182 = 4875.0/4096.0*eta;
  t183 = 15625.0/12288.0*t34;
  t184 = 625.0/768.0*t157;
  t185 = t173-t174+t175+t176+t177-t178-t179+t180+t181-t182-t183-t184;
  t189 = 8125.0/6144.0*t57;
  t191 = 15625.0/6144.0*t61;
  t192 = 78125.0/49152.0*t67;
  t195 = 825.0/16384.0+t84-t85-19625.0/24576.0*t5-4475.0/8192.0*t6+21875.0/24576.0*t65+t189-3125.0/24576.0*t59-t191-t192+34375.0/49152.0*t10+8125.0/16384.0*t12;
  t197 = 175.0/1536.0*xi;
  t198 = 7375.0/3072.0*eta;
  t199 = 10625.0/18432.0*t8;
  t202 = 625.0/384.0*t157;
  t203 = 15625.0/6144.0*t32;
  t204 = 125.0/1536.0*t5;
  t206 = 7375.0/3072.0*t7;
  t207 = 625.0/1536.0*t17;
  t208 = -175.0/4096.0-t197+t198+t99+t199+3125.0/12288.0*t10-11875.0/12288.0*t12+t202-t203+t204+775.0/768.0*t6-t206+t207;
  t209 = 2375.0/9216.0*t41;
  t210 = 3125.0/2304.0*t43;
  t211 = 4375.0/1152.0*t45;
  t212 = 5375.0/2304.0*t47;
  t213 = 775.0/288.0*t51;
  t214 = 10625.0/18432.0*t53;
  t215 = 15625.0/18432.0*t55;
  t216 = 11875.0/4608.0*t57;
  t217 = 3125.0/9216.0*t59;
  t218 = 15625.0/18432.0*t63;
  t220 = 78125.0/36864.0*t67;
  t221 = t209-t210-t211-t212+t113+t213-t214+t215-t216-t217+t119-t218-21875.0/9216.0*t65+t220;
  t223 = 275.0/3072.0*xi;
  t228 = 625.0/3072.0*t17;
  t229 = 275.0/8192.0+t223-t198-t99-t199+3125.0/24576.0*t10+23125.0/24576.0*t12-t202+t203+1625.0/12288.0*t5-11975.0/12288.0*t6+t206+t228;
  t230 = 19625.0/18432.0*t41;
  t231 = 6875.0/4608.0*t45;
  t232 = 11975.0/4608.0*t51;
  t233 = 23125.0/9216.0*t57;
  t234 = 34375.0/36864.0*t59;
  t235 = 15625.0/9216.0*t61;
  t236 = 34375.0/36864.0*t65;
  t237 = 78125.0/73728.0*t67;
  t238 = -t230+t210+t231+t212-t113-t232+t214-t215+t233+t234-t235+t218+t236-t237;
  t240 = 325.0/256.0*t51;
  t241 = 1625.0/2048.0*t41;
  t242 = 25.0/512.0*xi;
  t245 = 625.0/512.0*t57;
  t248 = -t173+t240+t241-t177+t178-t180-t242+t182+t184-t84-75.0/4096.0+t85-125.0/4096.0*t5+975.0/2048.0*t6-t245-3125.0/4096.0*t59-1875.0/4096.0*t12;
  t249 = 2803.0/12288.0*xi;
  t252 = (1.0+eta)*t15*t20;
  t253 = 125.0/768.0*t252;
  t256 = 3125.0/12288.0*t34;
  t259 = 121.0/16384.0-t249-t151-t4-t253+t152+4375.0/49152.0*t10+4375.0/49152.0*t12-t158+t159+t256-2225.0/24576.0*t5-2225.0/24576.0*t6+t162;
  t260 = 2875.0/12288.0*t17;
  t262 = 3125.0/36864.0*t45;
  t263 = 9125.0/36864.0*t51;
  t264 = 4375.0/18432.0*t57;
  t266 = 3125.0/18432.0*t61;
  t269 = t260-10475.0/36864.0*t41+t44+t262+t48+t50-t263-t165+t166+t264+15625.0/73728.0*t59+t266+t168+15625.0/73728.0*t65+15625.0/147456.0*t67;
  t271 = 625.0/12288.0*t51;
  t272 = 625.0/12288.0*t45;
  t273 = 4875.0/4096.0*t17;
  t274 = 4875.0/4096.0*xi;
  t275 = 625.0/768.0*t252;
  t280 = t142-t271+t241-t74+t272-t273+t274-t144-t183+t275-75.0/4096.0+t85+975.0/2048.0*t5-125.0/4096.0*t6-t147-3125.0/4096.0*t65-1875.0/4096.0*t10;
  t281 = 7375.0/3072.0*xi;
  t282 = 625.0/384.0*t252;
  t285 = 15625.0/6144.0*t34;
  t288 = 7375.0/3072.0*t17;
  t289 = 275.0/8192.0-t281+t125-t99-t282+t126+23125.0/24576.0*t10+3125.0/24576.0*t12+t285-11975.0/12288.0*t5+1625.0/12288.0*t6+t131+t288;
  t290 = 10625.0/18432.0*t45;
  t291 = 10625.0/18432.0*t51;
  t292 = 15625.0/18432.0*t57;
  t293 = 15625.0/18432.0*t61;
  t294 = -t230+t110-t290+t112-t113+t291-t134-t135+t292+t234-t293+t137+t236-t237;
  t299 = -175.0/4096.0+t281-t98+t99+t282-t100-11875.0/12288.0*t10+3125.0/12288.0*t12-t285+775.0/768.0*t5+t130+t106-t288;
  t301 = t209-t110+t290-t112+t113-t291+t115+t116-t292-21875.0/9216.0*t59+t293-t120-t169+t220;
  t303 = -t71+t271+t175+t74-t272+t76+t273-t78-t274+t80+t183-t275;
  t310 = 825.0/16384.0-t84-t85-4475.0/8192.0*t5-19625.0/24576.0*t6+t88-3125.0/24576.0*t65-t90+21875.0/24576.0*t59-t192+8125.0/16384.0*t10+34375.0/49152.0*t12;
  t316 = -19.0/1024.0+t249-t2+t4+t253-t9-625.0/6144.0*t10-4375.0/12288.0*t12-t30+t33-t256+775.0/6144.0*t5+4675.0/12288.0*t6+t38;
  t318 = -t260-t164-t44-t262-t48-t50+t263+t54+t56-t264+t217-t266-t64+3125.0/36864.0*t65+t170;
  t320 = 625.0/3072.0*t53;
  t321 = 625.0/3072.0*t8;
  t322 = 1375.0/1024.0*t7;
  t323 = 1375.0/1024.0*eta;
  t324 = 625.0/768.0*t29;
  t325 = 15625.0/3072.0*t32;
  t330 = 75.0/1024.0-t320-t240+t143+t177+t321-t322+t242+t323+t324-t325-t85-125.0/1024.0*t5-975.0/512.0*t6+t245-3125.0/1024.0*t59+1875.0/1024.0*t12;
  t331 = 125.0/48.0*eta;
  t332 = 8125.0/4608.0*t8;
  t335 = 625.0/384.0*t29;
  t336 = 15625.0/1536.0*t32;
  t338 = 125.0/48.0*t7;
  t339 = -275.0/2048.0-t223-t331+t99+t332+3125.0/6144.0*t10-23125.0/6144.0*t12-t335+t336-t204+11975.0/3072.0*t6+t338-t228;
  t340 = 8125.0/4608.0*t53;
  t342 = 15625.0/4608.0*t63;
  t344 = -t133-t210-t231-t212+t113+t232-t340-t116-t233+40625.0/4608.0*t59+t235+t342+34375.0/9216.0*t65-t139;
  t346 = 3125.0/3072.0*t10;
  t350 = 175.0/1024.0+t197+t331-t99-t332+t346+11875.0/3072.0*t12+t335-t336-1375.0/1536.0*t5-775.0/192.0*t6-t338-t207;
  t353 = t109+t210+t211+t212-t113-t213+t340+t116+t216-53125.0/4608.0*t59-t119-t342-21875.0/2304.0*t65+t122;
  t355 = -225.0/2048.0+t320+t174-t73-t176-t177-t321+t179+t322-t181-t323+t183;
  t361 = -t324+t325+t85+6625.0/6144.0*t5+8575.0/4096.0*t6+90625.0/12288.0*t65-t189+21875.0/3072.0*t59+t191-t92-t346-8125.0/4096.0*t12;
  t363 = 3125.0/12288.0*t53;
  t364 = 3125.0/12288.0*t51;
  t366 = 3125.0/12288.0*t45;
  t367 = 3125.0/12288.0*t8;
  t368 = 3125.0/12288.0*t17;
  t369 = 3125.0/12288.0*t7;
  t370 = 3125.0/12288.0*xi;
  t371 = 3125.0/12288.0*eta;
  t372 = 78125.0/12288.0*t34;
  t373 = 78125.0/12288.0*t32;
  t381 = -t363-t364+138125.0/4096.0*t41+t366+t367-t368-t369+t370+t371+5625.0/16384.0-t372-t373-41875.0/8192.0*t5-41875.0/8192.0*t6-234375.0/8192.0*t65-234375.0/8192.0*t59+390625.0/16384.0*t67+78125.0/16384.0*t10+78125.0/16384.0*t12;
  t382 = 78125.0/6144.0*t53;
  t383 = 138125.0/3072.0*t41;
  t384 = 78125.0/6144.0*t8;
  t385 = 78125.0/6144.0*t32;
  t388 = 78125.0/6144.0*t63;
  t390 = 78125.0/6144.0*t55;
  t392 = 390625.0/12288.0*t67;
  t393 = 15625.0/4096.0*t10;
  t395 = -t382-t383+t384-1875.0/4096.0+t385+4375.0/1024.0*t5+14375.0/1536.0*t6+t388+109375.0/3072.0*t65-t390+15625.0/384.0*t59-t392-t393-109375.0/12288.0*t12;
  t396 = 138125.0/6144.0*t41;
  t401 = 390625.0/24576.0*t67;
  t404 = 1875.0/8192.0+t382+t396-t384-t385+6875.0/4096.0*t5-104375.0/12288.0*t6-t388-171875.0/12288.0*t65+t390-296875.0/12288.0*t59+t401-15625.0/8192.0*t10+203125.0/24576.0*t12;
  t405 = 78125.0/6144.0*t51;
  t406 = 78125.0/6144.0*t45;
  t407 = 78125.0/6144.0*t34;
  t411 = 78125.0/6144.0*t57;
  t413 = 78125.0/6144.0*t61;
  t415 = 15625.0/4096.0*t12;
  t416 = -t405-t383+t406-1875.0/4096.0+t407+14375.0/1536.0*t5+4375.0/1024.0*t6+15625.0/384.0*t65+t411+109375.0/3072.0*t59-t413-t392-109375.0/12288.0*t10-t415;
  t417 = 3125.0/1536.0*xi;
  t418 = 3125.0/1536.0*eta;
  t419 = 15625.0/2304.0*t3;
  t420 = 21875.0/1152.0*t8;
  t425 = 3125.0/1536.0*t7;
  t426 = 3125.0/1536.0*t17;
  t428 = 625.0/1024.0-t417-t418+t419-t420+15625.0/3072.0*t10+15625.0/3072.0*t12-4375.0/768.0*t5-4375.0/768.0*t6+t425+t426+30625.0/576.0*t41;
  t429 = 15625.0/2304.0*t43;
  t430 = 21875.0/1152.0*t45;
  t431 = 15625.0/2304.0*t47;
  t432 = 15625.0/2304.0*t49;
  t433 = 21875.0/1152.0*t51;
  t434 = 21875.0/1152.0*t53;
  t435 = 78125.0/4608.0*t55;
  t436 = 78125.0/4608.0*t57;
  t438 = 78125.0/4608.0*t61;
  t439 = 78125.0/4608.0*t63;
  t442 = -t429-t430-t431+t432+t433+t434+t435-t436-109375.0/2304.0*t59+t438-t439-109375.0/2304.0*t65+390625.0/9216.0*t67;
  t444 = 3125.0/3072.0*xi;
  t445 = 15625.0/6144.0*t10;
  t446 = 15625.0/6144.0*t12;
  t449 = 3125.0/3072.0*t17;
  t450 = 48125.0/2304.0*t41;
  t451 = 625.0/2048.0-t444+t418-t419+t420+t445-t446-4375.0/1536.0*t5+6875.0/3072.0*t6-t425+t449-t450;
  t452 = 34375.0/4608.0*t45;
  t453 = 34375.0/4608.0*t51;
  t454 = 78125.0/9216.0*t57;
  t456 = 78125.0/9216.0*t61;
  t458 = 390625.0/18432.0*t67;
  t459 = t429+t452+t431-t432-t453-t434-t435+t454+109375.0/4608.0*t59-t456+t439+171875.0/9216.0*t65-t458;
  t467 = 1875.0/8192.0+t405+t396-t406-t407-104375.0/12288.0*t5+6875.0/4096.0*t6-296875.0/12288.0*t65-t411-171875.0/12288.0*t59+t413+t401+203125.0/24576.0*t10-15625.0/8192.0*t12;
  t468 = 3125.0/3072.0*eta;
  t469 = 34375.0/4608.0*t8;
  t472 = 3125.0/3072.0*t7;
  t473 = 625.0/2048.0+t417-t468-t419+t469-t445+t446+6875.0/3072.0*t5-4375.0/1536.0*t6+t472-t426-t450;
  t474 = 34375.0/4608.0*t53;
  t475 = 78125.0/9216.0*t55;
  t477 = 78125.0/9216.0*t63;
  t479 = t429+t430+t431-t432-t433-t474-t475+t436+171875.0/9216.0*t59-t438+t477+109375.0/4608.0*t65-t458;
  t486 = 625.0/4096.0+t444+t468+t419-t469-15625.0/12288.0*t10-15625.0/12288.0*t12+6875.0/6144.0*t5+6875.0/6144.0*t6-t472-t449+75625.0/9216.0*t41;
  t490 = -t429-t452-t431+t432+t453+t474+t475-t454-171875.0/18432.0*t59+t456-t477-171875.0/18432.0*t65+390625.0/36864.0*t67;
  t492 = 8125.0/2048.0*t41;
  t496 = -625.0/4096.0+t364-t492-t366+t368-t370+t372+8125.0/2048.0*t5+625.0/4096.0*t6+15625.0/4096.0*t65-t393;
  t500 = -625.0/4096.0+t363-t492-t367+t369-t371+t373+625.0/4096.0*t5+8125.0/2048.0*t6+15625.0/4096.0*t59-t415;

  values[0] = t39+t69;
  values[1] = t82+t95;
  values[2] = t108+t123;
  values[3] = t132+t140;
  values[4] = t150;
  values[5] = t163+t171;
  values[6] = t185+t195;
  values[7] = t208+t221;
  values[8] = t229+t238;
  values[9] = t248;
  values[10] = t259+t269;
  values[11] = t280;
  values[12] = t289+t294;
  values[13] = t299+t301;
  values[14] = t303+t310;
  values[15] = t316+t318;
  values[16] = t330;
  values[17] = t339+t344;
  values[18] = t350+t353;
  values[19] = t355+t361;
  values[20] = t381;
  values[21] = t395;
  values[22] = t404;
  values[23] = t416;
  values[24] = t428+t442;
  values[25] = t451+t459;
  values[26] = t467;
  values[27] = t473+t479;
  values[28] = t486+t490;
  values[29] = t496;
  values[30] = t500;
}

// values of the derivatives in xi direction
static void C_Q_UL5S_2D_DeriveXi(double xi, double eta, double *values)
{
  double t1, t2, t3, t4, t5, t6, t7, t8, t12, t13, t14, t15, t16, t17, t19;
  double t20, t21, t22, t24, t25, t26, t28, t32, t33, t34, t35, t37, t38, t39;
  double t41, t42, t45, t46, t47, t48, t49, t50, t51, t52, t53, t54, t56, t57;
  double t58, t59, t60, t61, t62, t65, t66, t68, t69, t70, t71, t72, t73, t74;
  double t75, t76, t78, t79, t80, t81, t83, t85, t87, t89, t91, t92, t93, t94;
  double t95, t96, t97, t98, t99, t100, t101, t102, t103, t105, t107, t108;
  double t109, t111, t112, t114, t115, t117, t118, t119, t122, t123, t125;
  double t127, t128, t130, t131, t133, t136, t138, t140, t141, t142, t143;
  double t144, t145, t146, t147, t149, t151, t152, t153, t154, t155, t156;
  double t157, t158, t159, t160, t163, t165, t166, t168, t169, t170, t171;
  double t172, t173, t174, t175, t177, t178, t179, t180, t181, t182, t183;
  double t184, t185, t187, t188, t190, t191, t192, t194, t195, t197, t198;
  double t199, t200, t201, t202, t204, t205, t207, t209, t210, t211, t213;
  double t215, t217, t218, t219, t221, t222, t224, t225, t226, t230, t231;
  double t233, t234, t235, t236, t237, t241, t242, t243, t244, t245, t247;
  double t248, t249, t250, t252, t253, t254, t257, t260, t267, t270, t273;
  double t275, t276, t277, t280, t281, t282, t284, t285, t287, t289, t291;
  double t292, t296, t301, t302, t303, t305, t306, t307, t308, t309, t310;
  double t316, t317, t318, t319, t320, t322, t324, t326, t327, t328, t329;
  double t333, t335, t336, t337, t338, t339, t342, t344, t346, t347, t349;
  double t350, t351, t352, t353, t354, t356, t358, t359, t360, t362, t363;
  double t366, t367, t368, t369, t370, t371, t373, t375, t377, t378, t383;
  double t384, t386, t387, t389, t391, t398, t399;

  t1 = xi*xi;
  t2 = 875.0/1024.0*t1;
  t3 = eta*eta;
  t4 = 2875.0/9216.0*t3;
  t5 = t3*eta;
  t6 = 125.0/576.0*t5;
  t7 = 1049.0/2304.0*eta;
  t8 = 1.0-t3;
  t12 = t8*(5.0/2.0*t5-3.0/2.0*eta);
  t13 = 125.0/768.0*t12;
  t14 = t1*eta;
  t15 = 125.0/192.0*t14;
  t16 = t3*t3;
  t17 = xi*t16;
  t19 = t1*t16;
  t20 = 3125.0/1536.0*t19;
  t21 = t1*xi;
  t22 = t21*t16;
  t24 = t1*t1;
  t25 = t24*t8;
  t26 = 3125.0/1536.0*t25;
  t28 = 1.0-eta;
  t32 = 5.0/2.0*t21-3.0/2.0*xi;
  t33 = t28*xi*t32;
  t34 = 125.0/384.0*t33;
  t35 = 1.0-t1;
  t37 = t35*t8*t1;
  t38 = 3125.0/1024.0*t37;
  t39 = -t2-t4+t6-t7-t13+t15-15625.0/4608.0*t17-t20+15625.0/2304.0*t22+t26-1325.0/1536.0*xi-t34-t38;
  t41 = xi*t8*t5;
  t42 = 3125.0/1536.0*t41;
  t45 = 15.0/2.0*t1-3.0/2.0;
  t46 = t28*t35*t45;
  t47 = 125.0/768.0*t46;
  t48 = xi*eta;
  t49 = 2875.0/4608.0*t48;
  t50 = t1*t3;
  t51 = 4375.0/3072.0*t50;
  t52 = t21*eta;
  t53 = 625.0/576.0*t52;
  t54 = t21*t3;
  t56 = t5*xi;
  t57 = 4375.0/4608.0*t56;
  t58 = t1*t5;
  t59 = 625.0/768.0*t58;
  t60 = t21*t5;
  t61 = 3125.0/1152.0*t60;
  t62 = xi*t3;
  t65 = 625.0/2304.0*t16;
  t66 = 857.0/3072.0+t42+t47-t49+t51+t53-15625.0/2304.0*t54+t57+t59-t61+4525.0/1152.0*t62+625.0/384.0*t21+t65;
  t68 = 17525.0/6144.0*t48;
  t69 = 625.0/3072.0*t3;
  t70 = 50375.0/3072.0*t62;
  t71 = 875.0/256.0*t14;
  t72 = 625.0/1024.0*t50;
  t73 = 33125.0/6144.0*t56;
  t74 = 4125.0/1024.0*t1;
  t75 = 15625.0/1536.0*t25;
  t76 = 15625.0/1024.0*t37;
  t78 = 625.0/384.0*t33;
  t79 = 625.0/768.0*t46;
  t80 = 15625.0/6144.0*t41;
  t81 = 875.0/768.0*eta;
  t83 = 8125.0/1536.0*t52;
  t85 = 15625.0/1536.0*t60;
  t87 = 78125.0/3072.0*t22;
  t89 = t78-t79-t80+t81+8575.0/2048.0*xi-t83+21875.0/768.0*t54+t85+90625.0/6144.0*t17-t87-8125.0/1024.0*t21;
  t91 = 5375.0/2304.0*eta;
  t92 = 15625.0/768.0*t25;
  t93 = 15625.0/512.0*t37;
  t94 = 28625.0/1152.0*t62;
  t95 = 5375.0/768.0*t14;
  t96 = 8125.0/1536.0*t50;
  t97 = 3125.0/2304.0*t5;
  t98 = 4375.0/576.0*t56;
  t99 = 3125.0/768.0*t58;
  t100 = 8125.0/4608.0*t3;
  t101 = 625.0/192.0*t33;
  t102 = 125.0/48.0-t91+t92-t93+t94+t95-t96+t97+t98-t99+t100-t101;
  t103 = 625.0/384.0*t46;
  t105 = 125.0/16.0*t1;
  t107 = 775.0/144.0*t48;
  t108 = 15625.0/1152.0*t60;
  t109 = 15625.0/4608.0*t16;
  t111 = 15625.0/1536.0*t19;
  t112 = 11875.0/1152.0*t52;
  t114 = 78125.0/2304.0*t22;
  t115 = t103+11875.0/768.0*t21-t105-775.0/96.0*xi-t107-t108-t109-21875.0/1152.0*t17+t111+t112-53125.0/1152.0*t54+t114;
  t117 = 20125.0/1152.0*t62;
  t118 = 6875.0/2304.0*t56;
  t119 = -125.0/48.0+t91-t92+t93-t117-t95+t96-t97-t118+t99-t100+t101;
  t122 = 11975.0/2304.0*t48;
  t123 = 15625.0/2304.0*t60;
  t125 = 23125.0/2304.0*t52;
  t127 = 78125.0/4608.0*t22;
  t128 = -t103-23125.0/1536.0*t21+t105+11975.0/1536.0*xi+t122+t123+t109+34375.0/4608.0*t17-t111-t125+40625.0/1152.0*t54-t127;
  t130 = 325.0/128.0*t48;
  t131 = 1625.0/256.0*t62;
  t133 = 625.0/128.0*t52;
  t136 = -t130-t69+t131+t71+t72-t74+1375.0/1024.0+t75-t76-t78+t79-t81-975.0/256.0*xi+t133-3125.0/256.0*t54+1875.0/256.0*t21;
  t138 = 15625.0/9216.0*t22;
  t140 = t2+t4-t6+t7+t13-t15+3125.0/18432.0*t17+t20+t138-t26+4675.0/6144.0*xi+t34+t38;
  t141 = 3125.0/6144.0*t41;
  t142 = 9125.0/18432.0*t48;
  t143 = 4375.0/4608.0*t52;
  t144 = 3125.0/2304.0*t54;
  t145 = 3125.0/18432.0*t56;
  t146 = 3125.0/4608.0*t60;
  t147 = 11575.0/9216.0*t62;
  t149 = -857.0/3072.0+t141-t47+t142-t51-t143+t144-t145-t59-t146-t147-4375.0/3072.0*t21-t65;
  t151 = 625.0/6144.0*t48;
  t152 = 17525.0/12288.0*t3;
  t153 = 11375.0/6144.0*t62;
  t154 = 33125.0/4096.0*t50;
  t155 = 875.0/768.0*t5;
  t156 = 625.0/6144.0*t56;
  t157 = 1875.0/4096.0*t1;
  t158 = 15625.0/6144.0*t25;
  t159 = 15625.0/4096.0*t37;
  t160 = 625.0/768.0*t12;
  t163 = 8125.0/6144.0*t16;
  t165 = 15625.0/2048.0*t19;
  t166 = 78125.0/12288.0*t22;
  t168 = t151-t152+t153+t154+t155-t156-t157+425.0/4096.0+t158-t159-t160-t80-t81-19625.0/12288.0*xi+21875.0/6144.0*t54+t163-3125.0/12288.0*t17-t165-t166+34375.0/12288.0*t21;
  t169 = 15625.0/3072.0*t41;
  t170 = 2375.0/4608.0*t62;
  t171 = 3125.0/768.0*t14;
  t172 = 4375.0/384.0*t50;
  t173 = 5375.0/2304.0*t5;
  t174 = 10625.0/9216.0*t56;
  t175 = 775.0/288.0*t3;
  t177 = -175.0/1536.0+t91+t169+t170-t171-t172-t173+t174+t99+t175+3125.0/3072.0*t21;
  t178 = 625.0/512.0*t1;
  t179 = 625.0/384.0*t12;
  t180 = 125.0/768.0*xi;
  t181 = 10625.0/9216.0*t48;
  t182 = 15625.0/4608.0*t60;
  t183 = 11875.0/4608.0*t16;
  t184 = 3125.0/4608.0*t17;
  t185 = 15625.0/4608.0*t52;
  t187 = 78125.0/9216.0*t22;
  t188 = t178+t179+t180-t181+t182-t183-t184+t111-t185-21875.0/2304.0*t54+t187;
  t190 = 19625.0/9216.0*t62;
  t191 = 6875.0/1536.0*t50;
  t192 = 11975.0/4608.0*t3;
  t194 = 275.0/3072.0-t91-t169-t190+t171+t191+t173-t174-t99-t192+3125.0/6144.0*t21;
  t195 = 625.0/1024.0*t1;
  t197 = 23125.0/9216.0*t16;
  t198 = 34375.0/18432.0*t17;
  t199 = 15625.0/3072.0*t19;
  t200 = 34375.0/9216.0*t54;
  t201 = 78125.0/18432.0*t22;
  t202 = t195-t179+1625.0/6144.0*xi+t181-t182+t197+t198-t199+t185+t200-t201;
  t204 = 325.0/256.0*t3;
  t205 = 1625.0/1024.0*t62;
  t207 = 625.0/512.0*t16;
  t209 = -t151+t204+t205-t155+t156-25.0/512.0+t160+t80+t81-125.0/2048.0*xi-t207-3125.0/2048.0*t17;
  t210 = 2875.0/4096.0*t1;
  t211 = 9125.0/36864.0*t3;
  t213 = 3125.0/6144.0*t19;
  t215 = 3125.0/6144.0*t25;
  t217 = 3125.0/4096.0*t37;
  t218 = -2803.0/12288.0+t210-t211+t6-t7-t13+t15+15625.0/36864.0*t17+t213+15625.0/36864.0*t22-t215-2225.0/12288.0*xi+t217;
  t219 = 1.0+eta;
  t221 = t219*xi*t32;
  t222 = 125.0/384.0*t221;
  t224 = t219*t35*t45;
  t225 = 125.0/768.0*t224;
  t226 = 3125.0/12288.0*t50;
  t230 = 4375.0/18432.0*t16;
  t231 = -t141+t222-t225-t142+t226+t143+15625.0/18432.0*t54+t145+t59+t146-10475.0/18432.0*t62+4375.0/12288.0*t21+t230;
  t233 = 625.0/12288.0*t3;
  t234 = 625.0/4096.0*t50;
  t235 = 14625.0/4096.0*t1;
  t236 = 625.0/384.0*t221;
  t237 = 625.0/768.0*t224;
  t241 = t130-t233+t205-t71+t234-t235+4875.0/4096.0+t158-t159-t236+t237+t81+975.0/1024.0*xi-t133-3125.0/1024.0*t54-1875.0/1024.0*t21;
  t242 = 15625.0/3072.0*t25;
  t243 = 15625.0/2048.0*t37;
  t244 = 10625.0/6144.0*t50;
  t245 = 10625.0/18432.0*t3;
  t247 = -7375.0/3072.0-t91-t242+t243-t190+t95-t244+t97+t118-t99+t245+23125.0/6144.0*t21;
  t248 = 7375.0/1024.0*t1;
  t249 = 625.0/192.0*t221;
  t250 = 625.0/384.0*t224;
  t252 = 15625.0/18432.0*t16;
  t253 = 15625.0/6144.0*t19;
  t254 = t248+t249-t250-11975.0/6144.0*xi-t122-t123+t252+t198-t253+t125+t200-t201;
  t257 = 7375.0/3072.0+t91+t242-t243+t170-t95+t244-t97-t98+t99-t245-11875.0/3072.0*t21;
  t260 = -t248-t249+t250+775.0/384.0*xi+t107+t108-t252-21875.0/4608.0*t17+t253-t112-t144+t187;
  t267 = t236-t237+t80-t81-4475.0/4096.0*xi+t83-3125.0/6144.0*t54-t85+21875.0/12288.0*t17-t166+8125.0/4096.0*t21;
  t270 = -t210+t211-t6+t7+t13-t15+t184-t213+t138+t215+775.0/3072.0*xi-t217-t42;
  t273 = 2803.0/12288.0-t222+t225+t49-t226-t53+3125.0/9216.0*t54-t57-t59+t61-t147-625.0/1536.0*t21-t230;
  t275 = 625.0/1536.0*t48;
  t276 = 625.0/1536.0*t56;
  t277 = 15625.0/1536.0*t41;
  t280 = -t275-t204+t131+t155+t276+25.0/512.0-t160+t277-t81-125.0/512.0*xi+t207-3125.0/512.0*t17;
  t281 = 15625.0/768.0*t41;
  t282 = 8125.0/2304.0*t56;
  t284 = -275.0/3072.0+t91-t281-t117-t171-t191-t173+t282+t99+t192+3125.0/1536.0*t21;
  t285 = 8125.0/2304.0*t48;
  t287 = 15625.0/1152.0*t52;
  t289 = -t195+t179-t180-t285-t108-t197+40625.0/2304.0*t17+t199+t287+34375.0/2304.0*t54-t127;
  t291 = 3125.0/768.0*t21;
  t292 = 175.0/1536.0-t91+t281+t94+t171+t172+t173-t282-t99-t175+t291;
  t296 = -t178-t179-1375.0/768.0*xi+t285+t108+t183-53125.0/2304.0*t17-t111-t287-21875.0/576.0*t54+t114;
  t301 = t275+t152-t70-t154-t155-t276+t157-425.0/4096.0-t158+t159+t160-t277+t81+6625.0/3072.0*xi+90625.0/3072.0*t54-t163+21875.0/1536.0*t17+t165-t87-t291;
  t302 = 3125.0/6144.0*t48;
  t303 = 3125.0/12288.0*t3;
  t305 = 3125.0/4096.0*t50;
  t306 = 3125.0/6144.0*t56;
  t307 = 3125.0/4096.0*t1;
  t308 = 78125.0/6144.0*t25;
  t309 = 78125.0/4096.0*t37;
  t310 = 78125.0/6144.0*t41;
  t316 = -t302-t303+138125.0/2048.0*t62+t305+t306-t307+3125.0/12288.0+t308-t309+t310-41875.0/4096.0*xi-234375.0/2048.0*t54-234375.0/4096.0*t17+390625.0/4096.0*t22+78125.0/4096.0*t21;
  t317 = 78125.0/3072.0*t48;
  t318 = 138125.0/1536.0*t62;
  t319 = 78125.0/3072.0*t56;
  t320 = 78125.0/3072.0*t41;
  t322 = 78125.0/1536.0*t52;
  t324 = 78125.0/1536.0*t60;
  t326 = 390625.0/3072.0*t22;
  t327 = 15625.0/1024.0*t21;
  t328 = -t317-t318+t319-t320+4375.0/512.0*xi+t322+109375.0/768.0*t54-t324+15625.0/192.0*t17-t326-t327;
  t329 = 138125.0/3072.0*t62;
  t333 = 390625.0/6144.0*t22;
  t335 = t317+t329-t319+t320+6875.0/2048.0*xi-t322-171875.0/3072.0*t54+t324-296875.0/6144.0*t17+t333-15625.0/2048.0*t21;
  t336 = 78125.0/6144.0*t3;
  t337 = 78125.0/2048.0*t50;
  t338 = 78125.0/3072.0*t25;
  t339 = 78125.0/2048.0*t37;
  t342 = 78125.0/6144.0*t16;
  t344 = 78125.0/2048.0*t19;
  t346 = -t336-t318+t337-t338+t339+14375.0/768.0*xi+15625.0/96.0*t54+t342+109375.0/1536.0*t17-t344-t326-109375.0/3072.0*t21;
  t347 = 15625.0/2304.0*eta;
  t349 = 15625.0/768.0*t14;
  t350 = 21875.0/384.0*t50;
  t351 = 15625.0/2304.0*t5;
  t352 = 21875.0/576.0*t56;
  t353 = 15625.0/768.0*t58;
  t354 = 21875.0/1152.0*t3;
  t356 = 3125.0/512.0*t1;
  t358 = 21875.0/576.0*t48;
  t359 = 78125.0/1152.0*t60;
  t360 = 78125.0/4608.0*t16;
  t362 = 78125.0/1536.0*t19;
  t363 = 78125.0/1152.0*t52;
  t366 = t347-3125.0/1536.0+30625.0/288.0*t62-t349-t350-t351-t352+t353+t354+15625.0/768.0*t21+t356-4375.0/384.0*xi+t358+t359-t360-109375.0/1152.0*t17+t362-t363-109375.0/576.0*t54+390625.0/2304.0*t22;
  t367 = 48125.0/1152.0*t62;
  t368 = 34375.0/1536.0*t50;
  t369 = 34375.0/4608.0*t3;
  t370 = 15625.0/1536.0*t21;
  t371 = 3125.0/1024.0*t1;
  t373 = 78125.0/9216.0*t16;
  t375 = 78125.0/3072.0*t19;
  t377 = 390625.0/4608.0*t22;
  t378 = -t347-3125.0/3072.0-t367+t349+t368+t351+t352-t353-t369+t370+t371-4375.0/768.0*xi-t358-t359+t373+109375.0/2304.0*t17-t375+t363+171875.0/2304.0*t54-t377;
  t383 = t336+t329-t337+t338-t339-104375.0/6144.0*xi-296875.0/3072.0*t54-t342-171875.0/6144.0*t17+t344+t333+203125.0/6144.0*t21;
  t384 = 34375.0/2304.0*t56;
  t386 = 34375.0/2304.0*t48;
  t387 = 78125.0/2304.0*t60;
  t389 = 78125.0/2304.0*t52;
  t391 = -t347+3125.0/1536.0-t367+t349+t350+t351+t384-t353-t354-t370-t356+6875.0/1536.0*xi-t386-t387+t360+171875.0/4608.0*t17-t362+t389+109375.0/1152.0*t54-t377;
  t398 = t347+3125.0/3072.0+75625.0/4608.0*t62-t349-t368-t351-t384+t353+t369-15625.0/3072.0*t21-t371+6875.0/3072.0*xi+t386+t387-t373-171875.0/9216.0*t17+t375-t389-171875.0/4608.0*t54+390625.0/9216.0*t22;
  t399 = 8125.0/1024.0*t62;

  values[0] = t39+t66;
  values[1] = -1375.0/1024.0+t68+t69-t70-t71-t72-t73+t74-t75+t76+t89;
  values[2] = t102+t115;
  values[3] = t119+t128;
  values[4] = t136;
  values[5] = t140+t149;
  values[6] = t168;
  values[7] = t177+t188;
  values[8] = t194+t202;
  values[9] = t209;
  values[10] = t218+t231;
  values[11] = t241;
  values[12] = t247+t254;
  values[13] = t257+t260;
  values[14] = -4875.0/4096.0-t68+t233+t153+t71-t234+t73+t235-t158+t159+t267;
  values[15] = t270+t273;
  values[16] = t280;
  values[17] = t284+t289;
  values[18] = t292+t296;
  values[19] = t301;
  values[20] = t316;
  values[21] = t328;
  values[22] = t335;
  values[23] = t346;
  values[24] = t366;
  values[25] = t378;
  values[26] = t383;
  values[27] = t391;
  values[28] = t398;
  values[29] = t303-t399-t305+t307-3125.0/12288.0-t308+t309+8125.0/1024.0*xi+15625.0/1024.0*t54-t327;
  values[30] = t302-t399-t306-t310+625.0/2048.0*xi+15625.0/2048.0*t17;
}

// values of the derivatives in eta direction
static void C_Q_UL5S_2D_DeriveEta(double xi, double eta, double *values)
{
  double t1, t2, t3, t4, t5, t6, t7, t8, t9, t11, t12, t13, t14, t15, t16;
  double t17, t18, t19, t20, t22, t23, t24, t27, t28, t29, t30, t34, t35, t36;
  double t38, t39, t40, t41, t42, t43, t44, t45, t49, t50, t52, t54, t55, t57;
  double t59, t60, t61, t62, t63, t64, t66, t68, t69, t70, t71, t72, t73, t74;
  double t75, t76, t77, t78, t79, t81, t83, t85, t86, t87, t88, t89, t90, t91;
  double t92, t93, t94, t95, t96, t97, t98, t99, t101, t102, t103, t105, t106;
  double t108, t109, t111, t112, t113, t115, t116, t117, t118, t120, t122;
  double t123, t125, t126, t128, t130, t131, t133, t134, t135, t136, t138;
  double t139, t140, t142, t143, t144, t146, t147, t148, t149, t150, t152;
  double t154, t155, t156, t157, t158, t159, t160, t161, t162, t164, t167;
  double t169, t170, t172, t174, t175, t176, t177, t178, t179, t180, t181;
  double t183, t184, t185, t186, t188, t189, t190, t191, t192, t194, t195;
  double t197, t198, t199, t201, t203, t204, t205, t206, t207, t208, t210;
  double t211, t213, t216, t217, t219, t220, t221, t222, t227, t229, t231;
  double t232, t235, t236, t237, t238, t240, t242, t243, t244, t247, t249;
  double t255, t256, t260, t262, t263, t264, t265, t266, t267, t268, t272;
  double t273, t274, t275, t277, t278, t279, t280, t282, t284, t286, t289;
  double t293, t300, t302, t303, t305, t306, t307, t308, t309, t310, t316;
  double t317, t318, t319, t320, t321, t323, t325, t327, t329, t330, t334;
  double t336, t337, t338, t339, t342, t344, t345, t346, t347, t349, t350;
  double t351, t352, t353, t354, t356, t358, t359, t360, t362, t363, t366;
  double t367, t368, t369, t370, t372, t374, t376, t377, t382, t383, t384;
  double t386, t387, t389, t391, t398, t399;

  t1 = xi*xi;
  t2 = eta*eta;
  t3 = t1*t2;
  t4 = 4375.0/3072.0*t3;
  t5 = t1*xi;
  t6 = t2*eta;
  t7 = t5*t6;
  t8 = 3125.0/1152.0*t7;
  t9 = t1*t6;
  t11 = xi*t6;
  t12 = 625.0/576.0*t11;
  t13 = t5*t2;
  t14 = 625.0/768.0*t13;
  t15 = t5*eta;
  t16 = 4375.0/4608.0*t15;
  t17 = xi*t2;
  t18 = 125.0/192.0*t17;
  t19 = 1.0-t1;
  t20 = 1.0-t2;
  t22 = t19*t20*t2;
  t23 = 3125.0/1024.0*t22;
  t24 = 1.0-xi;
  t27 = 15.0/2.0*t2-3.0/2.0;
  t28 = t24*t20*t27;
  t29 = 125.0/768.0*t28;
  t30 = 2875.0/9216.0*t1;
  t34 = 5.0/2.0*t6-3.0/2.0*eta;
  t35 = t24*eta*t34;
  t36 = 125.0/384.0*t35;
  t38 = t19*eta*t5;
  t39 = 3125.0/1536.0*t38;
  t40 = 875.0/1024.0*t2;
  t41 = t4-t8-15625.0/2304.0*t9+t12+t14+t16+t18-t23+t29-t30-t36+t39-t40;
  t42 = 125.0/576.0*t5;
  t43 = t2*t2;
  t44 = t19*t43;
  t45 = 3125.0/1536.0*t44;
  t49 = t19*(5.0/2.0*t5-3.0/2.0*xi);
  t50 = 125.0/768.0*t49;
  t52 = t1*eta;
  t54 = t1*t1;
  t55 = t54*t6;
  t57 = t54*eta;
  t59 = t54*t2;
  t60 = 3125.0/1536.0*t59;
  t61 = 625.0/2304.0*t54;
  t62 = xi*eta;
  t63 = 2875.0/4608.0*t62;
  t64 = 1049.0/2304.0*xi;
  t66 = 857.0/3072.0+t42+t45-t50+625.0/384.0*t6+4525.0/1152.0*t52+15625.0/2304.0*t55-15625.0/4608.0*t57-t60+t61-t63-t64-1325.0/1536.0*eta;
  t68 = 17525.0/12288.0*t1;
  t69 = 625.0/1536.0*t62;
  t70 = 50375.0/3072.0*t52;
  t71 = 875.0/768.0*t5;
  t72 = 625.0/1536.0*t15;
  t73 = 33125.0/4096.0*t3;
  t74 = 1875.0/4096.0*t2;
  t75 = 15625.0/1536.0*t38;
  t76 = 625.0/768.0*t49;
  t77 = 15625.0/6144.0*t44;
  t78 = 15625.0/4096.0*t22;
  t79 = 875.0/768.0*xi;
  t81 = 8125.0/6144.0*t54;
  t83 = 15625.0/2048.0*t59;
  t85 = 78125.0/3072.0*t55;
  t86 = 3125.0/768.0*t6;
  t87 = t68+t69-t70-t71-t72-t73+t74-425.0/4096.0-t75+t76-t77+t78+t79+6625.0/3072.0*eta-t81+21875.0/1536.0*t57+t83+90625.0/3072.0*t9-t85-t86;
  t88 = 5375.0/2304.0*xi;
  t89 = 15625.0/768.0*t38;
  t90 = 28625.0/1152.0*t52;
  t91 = 5375.0/2304.0*t5;
  t92 = 8125.0/2304.0*t15;
  t93 = 3125.0/768.0*t17;
  t94 = 4375.0/384.0*t3;
  t95 = 3125.0/768.0*t13;
  t96 = 8125.0/2304.0*t62;
  t97 = 625.0/384.0*t49;
  t98 = 175.0/1536.0-t88+t89+t90+t91-t92+t93+t94-t95+t96-t97;
  t99 = 625.0/512.0*t2;
  t101 = 775.0/288.0*t1;
  t102 = 15625.0/1536.0*t59;
  t103 = 15625.0/1152.0*t11;
  t105 = 15625.0/1152.0*t7;
  t106 = 11875.0/4608.0*t54;
  t108 = 78125.0/2304.0*t55;
  t109 = t86-t99-1375.0/768.0*eta-t101-t102-t103-21875.0/576.0*t9+t105+t106-53125.0/2304.0*t57+t108;
  t111 = 20125.0/1152.0*t52;
  t112 = 6875.0/1536.0*t3;
  t113 = -275.0/3072.0+t88-t89-t111-t91+t92-t93-t112+t95-t96+t97;
  t115 = 625.0/1024.0*t2;
  t116 = 125.0/768.0*eta;
  t117 = 11975.0/4608.0*t1;
  t118 = 15625.0/3072.0*t59;
  t120 = 23125.0/9216.0*t54;
  t122 = 78125.0/4608.0*t55;
  t123 = 3125.0/1536.0*t6-t115-t116+t117+t118+t103+34375.0/2304.0*t9-t105-t120+40625.0/2304.0*t57-t122;
  t125 = 325.0/256.0*t1;
  t126 = 1625.0/256.0*t52;
  t128 = 625.0/512.0*t54;
  t130 = -t125-t69+t126+t71+t72+25.0/512.0+t75-t76-t79-125.0/512.0*eta+t128-3125.0/512.0*t57;
  t131 = 3125.0/12288.0*t3;
  t133 = 3125.0/4096.0*t22;
  t134 = 9125.0/36864.0*t1;
  t135 = 2875.0/4096.0*t2;
  t136 = 1.0+xi;
  t138 = t136*eta*t34;
  t139 = 125.0/384.0*t138;
  t140 = -t131+t8+3125.0/9216.0*t9-t12-t14-t16-t18-t133+t134-t39-t135-t42-t139;
  t142 = t136*t20*t27;
  t143 = 125.0/768.0*t142;
  t144 = 3125.0/6144.0*t44;
  t146 = 11575.0/9216.0*t52;
  t147 = 15625.0/9216.0*t55;
  t148 = 3125.0/4608.0*t57;
  t149 = 3125.0/6144.0*t59;
  t150 = 4375.0/18432.0*t54;
  t152 = 2803.0/12288.0+t143+t144+t50-625.0/1536.0*t6-t146+t147+t148-t149-t150+t63+t64+775.0/3072.0*eta;
  t154 = 625.0/12288.0*t1;
  t155 = 17525.0/6144.0*t62;
  t156 = 11375.0/6144.0*t52;
  t157 = 33125.0/6144.0*t15;
  t158 = 875.0/256.0*t17;
  t159 = 625.0/4096.0*t3;
  t160 = 14625.0/4096.0*t2;
  t161 = 15625.0/6144.0*t38;
  t162 = 625.0/384.0*t138;
  t164 = 625.0/768.0*t142;
  t167 = 8125.0/1536.0*t11;
  t169 = 15625.0/1536.0*t7;
  t170 = 78125.0/12288.0*t55;
  t172 = -t164-t77+t78-t79-4475.0/4096.0*eta+21875.0/12288.0*t57+t167-3125.0/6144.0*t9-t169-t170+8125.0/4096.0*t6;
  t174 = 15625.0/3072.0*t44;
  t175 = 15625.0/2048.0*t22;
  t176 = 2375.0/4608.0*t52;
  t177 = 3125.0/2304.0*t5;
  t178 = 4375.0/576.0*t15;
  t179 = 5375.0/768.0*t17;
  t180 = 10625.0/6144.0*t3;
  t181 = 775.0/144.0*t62;
  t183 = 7375.0/3072.0+t88+t174-t175+t176-t177-t178-t179+t180+t95+t181-11875.0/3072.0*t6;
  t184 = 7375.0/1024.0*t2;
  t185 = 625.0/192.0*t138;
  t186 = 625.0/384.0*t142;
  t188 = 10625.0/18432.0*t1;
  t189 = 15625.0/6144.0*t59;
  t190 = 11875.0/1152.0*t11;
  t191 = 3125.0/2304.0*t9;
  t192 = 15625.0/18432.0*t54;
  t194 = 78125.0/9216.0*t55;
  t195 = -t184-t185+t186+775.0/384.0*eta-t188+t189-t190-t191+t105-t192-21875.0/4608.0*t57+t194;
  t197 = 19625.0/9216.0*t52;
  t198 = 6875.0/2304.0*t15;
  t199 = 11975.0/2304.0*t62;
  t201 = -7375.0/3072.0-t88-t174+t175-t197+t177+t198+t179-t180-t95-t199+23125.0/6144.0*t6;
  t203 = 23125.0/2304.0*t11;
  t204 = 34375.0/9216.0*t9;
  t205 = 15625.0/2304.0*t7;
  t206 = 34375.0/18432.0*t57;
  t207 = 78125.0/18432.0*t55;
  t208 = t184+t185-t186-11975.0/6144.0*eta+t188-t189+t203+t204-t205+t192+t206-t207;
  t210 = 325.0/128.0*t62;
  t211 = 1625.0/1024.0*t52;
  t213 = 625.0/128.0*t11;
  t216 = -t154+t210+t211-t158+t159-t160+4875.0/4096.0-t162+t164+t77-t78+t79+975.0/1024.0*eta-t213-3125.0/1024.0*t9-1875.0/1024.0*t6;
  t217 = 3125.0/4608.0*t7;
  t219 = 4375.0/4608.0*t11;
  t220 = 3125.0/18432.0*t15;
  t221 = 3125.0/6144.0*t38;
  t222 = t131+t217+15625.0/18432.0*t9+t219+t14+t220+t18+t133-t134-t221+t135+t42+t139;
  t227 = 9125.0/18432.0*t62;
  t229 = -2803.0/12288.0-t143-t144-t50+4375.0/12288.0*t6-10475.0/18432.0*t52+15625.0/36864.0*t55+15625.0/36864.0*t57+t149+t150-t227-t64-2225.0/12288.0*eta;
  t231 = 625.0/6144.0*t62;
  t232 = 625.0/6144.0*t15;
  t235 = t125-t231+t211-t71+t232-25.0/512.0+t161+t76+t79-125.0/2048.0*eta-t128-3125.0/2048.0*t57;
  t236 = 15625.0/3072.0*t38;
  t237 = 10625.0/9216.0*t15;
  t238 = 10625.0/9216.0*t62;
  t240 = 275.0/3072.0-t88-t236-t197+t91-t237+t93+t112-t95+t238+3125.0/6144.0*t6;
  t242 = 15625.0/4608.0*t11;
  t243 = 15625.0/4608.0*t7;
  t244 = t115-t97+1625.0/6144.0*eta-t117-t118+t242+t204-t243+t120+t206-t207;
  t247 = -175.0/1536.0+t88+t236+t176-t91+t237-t93-t94+t95-t238+3125.0/3072.0*t6;
  t249 = t99+t97+t116+t101+t102-t242-21875.0/2304.0*t9+t243-t106-t148+t194;
  t255 = -t68+t231+t156+t71-t232+t73-t74+425.0/4096.0-t161-t76+t77-t78-t79-19625.0/12288.0*eta+t81-3125.0/12288.0*t57-t83+21875.0/6144.0*t9-t170+34375.0/12288.0*t6;
  t256 = -t4-t217+t191-t219-t14-t220-t18+t23-t29+t30+t36+t221+t40;
  t260 = -857.0/3072.0-t42-t45+t50-4375.0/3072.0*t6-t146+t147+3125.0/18432.0*t57+t60-t61+t227+t64+4675.0/6144.0*eta;
  t262 = 625.0/3072.0*t1;
  t263 = 625.0/1024.0*t3;
  t264 = 4125.0/1024.0*t2;
  t265 = 625.0/384.0*t35;
  t266 = 625.0/768.0*t28;
  t267 = 15625.0/1536.0*t44;
  t268 = 15625.0/1024.0*t22;
  t272 = -t262-t210+t126+t158+t263-t264+1375.0/1024.0-t265+t266+t267-t268-t79-975.0/256.0*eta+t213-3125.0/256.0*t9+1875.0/256.0*t6;
  t273 = 15625.0/768.0*t44;
  t274 = 15625.0/512.0*t22;
  t275 = 8125.0/1536.0*t3;
  t277 = -125.0/48.0+t88-t273+t274-t111-t177-t198-t179+t275+t95+t199-23125.0/1536.0*t6;
  t278 = 125.0/16.0*t2;
  t279 = 625.0/192.0*t35;
  t280 = 625.0/384.0*t28;
  t282 = 8125.0/4608.0*t1;
  t284 = 15625.0/4608.0*t54;
  t286 = t278+t279-t280+11975.0/1536.0*eta-t282-t102-t203+40625.0/1152.0*t9+t205+t284+34375.0/4608.0*t57-t122;
  t289 = 125.0/48.0-t88+t273-t274+t90+t177+t178+t179-t275-t95-t181+11875.0/768.0*t6;
  t293 = -t278-t279+t280-775.0/96.0*eta+t282+t102+t190-53125.0/1152.0*t9-t105-t284-21875.0/1152.0*t57+t108;
  t300 = -t266-t267+t268+t79+8575.0/2048.0*eta+90625.0/6144.0*t57-t167+21875.0/768.0*t9+t169-t85-8125.0/1024.0*t6;
  t302 = 3125.0/12288.0*t1;
  t303 = 3125.0/6144.0*t62;
  t305 = 3125.0/6144.0*t15;
  t306 = 3125.0/4096.0*t3;
  t307 = 3125.0/4096.0*t2;
  t308 = 78125.0/6144.0*t38;
  t309 = 78125.0/6144.0*t44;
  t310 = 78125.0/4096.0*t22;
  t316 = -t302-t303+138125.0/2048.0*t52+t305+t306-t307+3125.0/12288.0+t308+t309-t310-41875.0/4096.0*eta-234375.0/4096.0*t57-234375.0/2048.0*t9+390625.0/4096.0*t55+78125.0/4096.0*t6;
  t317 = 78125.0/6144.0*t1;
  t318 = 138125.0/1536.0*t52;
  t319 = 78125.0/2048.0*t3;
  t320 = 78125.0/3072.0*t44;
  t321 = 78125.0/2048.0*t22;
  t323 = 78125.0/6144.0*t54;
  t325 = 78125.0/2048.0*t59;
  t327 = 390625.0/3072.0*t55;
  t329 = -t317-t318+t319-t320+t321+14375.0/768.0*eta+t323+109375.0/1536.0*t57-t325+15625.0/96.0*t9-t327-109375.0/3072.0*t6;
  t330 = 138125.0/3072.0*t52;
  t334 = 390625.0/6144.0*t55;
  t336 = t317+t330-t319+t320-t321-104375.0/6144.0*eta-t323-171875.0/6144.0*t57+t325-296875.0/3072.0*t9+t334+203125.0/6144.0*t6;
  t337 = 78125.0/3072.0*t62;
  t338 = 78125.0/3072.0*t15;
  t339 = 78125.0/3072.0*t38;
  t342 = 78125.0/1536.0*t11;
  t344 = 78125.0/1536.0*t7;
  t345 = 15625.0/1024.0*t6;
  t346 = -t337-t318+t338-t339+4375.0/512.0*eta+15625.0/192.0*t57+t342+109375.0/768.0*t9-t344-t327-t345;
  t347 = 15625.0/2304.0*xi;
  t349 = 15625.0/2304.0*t5;
  t350 = 21875.0/576.0*t15;
  t351 = 15625.0/768.0*t17;
  t352 = 21875.0/384.0*t3;
  t353 = 15625.0/768.0*t13;
  t354 = 21875.0/576.0*t62;
  t356 = 3125.0/512.0*t2;
  t358 = 21875.0/1152.0*t1;
  t359 = 78125.0/1536.0*t59;
  t360 = 78125.0/1152.0*t11;
  t362 = 78125.0/1152.0*t7;
  t363 = 78125.0/4608.0*t54;
  t366 = t347-3125.0/1536.0+30625.0/288.0*t52-t349-t350-t351-t352+t353+t354+15625.0/768.0*t6+t356-4375.0/384.0*eta+t358+t359-t360-109375.0/576.0*t9+t362-t363-109375.0/1152.0*t57+390625.0/2304.0*t55;
  t367 = 48125.0/1152.0*t52;
  t368 = 34375.0/2304.0*t15;
  t369 = 34375.0/2304.0*t62;
  t370 = 15625.0/1536.0*t6;
  t372 = 78125.0/2304.0*t11;
  t374 = 78125.0/2304.0*t7;
  t376 = 390625.0/4608.0*t55;
  t377 = -t347+3125.0/1536.0-t367+t349+t368+t351+t352-t353-t369-t370-t356+6875.0/1536.0*eta-t358-t359+t372+109375.0/1152.0*t9-t374+t363+171875.0/4608.0*t57-t376;
  t382 = t337+t330-t338+t339+6875.0/2048.0*eta-296875.0/6144.0*t57-t342-171875.0/3072.0*t9+t344+t334-15625.0/2048.0*t6;
  t383 = 34375.0/1536.0*t3;
  t384 = 3125.0/1024.0*t2;
  t386 = 34375.0/4608.0*t1;
  t387 = 78125.0/3072.0*t59;
  t389 = 78125.0/9216.0*t54;
  t391 = -t347-3125.0/3072.0-t367+t349+t350+t351+t383-t353-t354+t370+t384-4375.0/768.0*eta-t386-t387+t360+171875.0/2304.0*t9-t362+t389+109375.0/2304.0*t57-t376;
  t398 = t347+3125.0/3072.0+75625.0/4608.0*t52-t349-t368-t351-t383+t353+t369-15625.0/3072.0*t6-t384+6875.0/3072.0*eta+t386+t387-t372-171875.0/4608.0*t9+t374-t389-171875.0/9216.0*t57+390625.0/9216.0*t55;
  t399 = 8125.0/1024.0*t52;

  values[0] = t41+t66;
  values[1] = t87;
  values[2] = t98+t109;
  values[3] = t113+t123;
  values[4] = t130;
  values[5] = t140+t152;
  values[6] = -4875.0/4096.0+t154-t155+t156+t157+t158-t159+t160+t161+t162+t172;
  values[7] = t183+t195;
  values[8] = t201+t208;
  values[9] = t216;
  values[10] = t222+t229;
  values[11] = t235;
  values[12] = t240+t244;
  values[13] = t247+t249;
  values[14] = t255;
  values[15] = t256+t260;
  values[16] = t272;
  values[17] = t277+t286;
  values[18] = t289+t293;
  values[19] = -1375.0/1024.0+t262+t155-t70-t157-t158-t263+t264-t161+t265+t300;
  values[20] = t316;
  values[21] = t329;
  values[22] = t336;
  values[23] = t346;
  values[24] = t366;
  values[25] = t377;
  values[26] = t382;
  values[27] = t391;
  values[28] = t398;
  values[29] = t303-t399-t305-t308+625.0/2048.0*eta+15625.0/2048.0*t57;
  values[30] = t302-t399-t306+t307-3125.0/12288.0-t309+t310+8125.0/1024.0*eta+15625.0/1024.0*t9-t345;
}

// values of the derivatives in xi-xi  direction
static void C_Q_UL5S_2D_DeriveXiXi(double xi, double eta, double *values)
{
  double t1, t2, t4, t6, t7, t8, t9, t11, t12, t13, t14, t15, t16, t17, t18;
  double t19, t20, t21, t22, t24, t25, t26, t28, t29, t30, t31, t32, t34, t35;
  double t38, t39, t40, t43, t44, t45, t46, t48, t49, t50, t52, t54, t55, t56;
  double t57, t58, t59, t60, t61, t62, t63, t64, t65, t66, t68, t70, t72, t73;
  double t74, t75, t76, t77, t78, t79, t80, t81, t82, t84, t85, t86, t88, t89;
  double t91, t92, t93, t94, t96, t97, t99, t101, t102, t103, t104, t105, t108;
  double t109, t110, t111, t112, t113, t115, t116, t119, t120, t122, t123;
  double t124, t125, t126, t127, t128, t131, t132, t134, t135, t136, t137;
  double t138, t139, t141, t142, t143, t144, t145, t147, t148, t149, t150;
  double t152, t153, t154, t155, t156, t157, t158, t161, t163, t165, t168;
  double t170, t171, t173, t174, t175, t176, t177, t179, t181, t183, t184;
  double t185, t186, t187, t190, t191, t192, t193, t195, t196, t197, t198;
  double t199, t200, t203, t207, t211, t213, t214, t215, t218, t219, t221;
  double t223, t225, t226, t229, t232, t233, t235, t236, t237, t238, t239;
  double t240, t245, t246, t247, t248, t249, t250, t252, t254, t255, t256;
  double t257, t260, t262, t263, t264, t265, t268, t272, t273, t274, t275;
  double t277, t278, t279, t281, t282, t285, t286, t287, t288, t289, t291;
  double t293, t294, t299, t300, t301, t303, t305, t311, t312;

  t1 = 2875.0/4608.0*eta;
  t2 = 1.0-eta;
  t4 = xi*xi;
  t6 = 15.0/2.0*t4-3.0/2.0;
  t7 = t2*xi*t6;
  t8 = 125.0/192.0*t7;
  t9 = 1.0-t4;
  t11 = t2*t9*xi;
  t12 = 625.0/256.0*t11;
  t13 = eta*eta;
  t14 = 1.0-t13;
  t15 = t13*eta;
  t16 = t14*t15;
  t17 = 3125.0/1536.0*t16;
  t18 = xi*t13;
  t19 = 4375.0/1536.0*t18;
  t20 = t4*eta;
  t21 = 625.0/192.0*t20;
  t22 = t4*t13;
  t24 = 4375.0/4608.0*t15;
  t25 = xi*t15;
  t26 = 625.0/384.0*t25;
  t28 = t4*t15;
  t29 = 3125.0/384.0*t28;
  t30 = t4*xi;
  t31 = t30*t14;
  t32 = 21875.0/1536.0*t31;
  t34 = t9*t14*xi;
  t35 = 3125.0/512.0*t34;
  t38 = 5.0/2.0*t30-3.0/2.0*xi;
  t39 = t2*t38;
  t40 = 125.0/384.0*t39;
  t43 = 875.0/512.0*xi;
  t44 = xi*eta;
  t45 = 125.0/96.0*t44;
  t46 = t13*t13;
  t48 = xi*t46;
  t49 = 3125.0/768.0*t48;
  t50 = t4*t46;
  t52 = -t29+t32-t35-t40+4525.0/1152.0*t13+625.0/128.0*t4-t43+t45-15625.0/4608.0*t46-t49+15625.0/768.0*t50;
  t54 = 17525.0/6144.0*eta;
  t55 = 50375.0/3072.0*t13;
  t56 = 875.0/128.0*t44;
  t57 = 625.0/512.0*t18;
  t58 = 33125.0/6144.0*t15;
  t59 = 4125.0/512.0*xi;
  t60 = 109375.0/1536.0*t31;
  t61 = 15625.0/512.0*t34;
  t62 = 625.0/384.0*t39;
  t63 = 625.0/192.0*t7;
  t64 = 3125.0/256.0*t11;
  t65 = 15625.0/6144.0*t16;
  t66 = 8125.0/512.0*t20;
  t68 = 15625.0/512.0*t28;
  t70 = 78125.0/1024.0*t50;
  t72 = t54-t55-t56-t57-t58+t59-t60+t61+t62+t63-t64-t65+8575.0/2048.0-t66+21875.0/256.0*t22+t68+90625.0/6144.0*t46-t70-24375.0/1024.0*t4;
  t73 = 109375.0/768.0*t31;
  t74 = 15625.0/256.0*t34;
  t75 = 28625.0/1152.0*t13;
  t76 = 5375.0/384.0*t44;
  t77 = 8125.0/768.0*t18;
  t78 = 4375.0/576.0*t15;
  t79 = 3125.0/384.0*t25;
  t80 = 625.0/192.0*t39;
  t81 = 625.0/96.0*t7;
  t82 = 3125.0/128.0*t11;
  t84 = 125.0/8.0*xi;
  t85 = 775.0/144.0*eta;
  t86 = 15625.0/384.0*t28;
  t88 = 15625.0/768.0*t48;
  t89 = 11875.0/384.0*t20;
  t91 = 78125.0/768.0*t50;
  t92 = t73-t74+t75+t76-t77+t78-t79-t80-t81+t82+11875.0/256.0*t4-t84-775.0/96.0-t85-t86-21875.0/1152.0*t46+t88+t89-53125.0/384.0*t22+t91;
  t93 = 20125.0/1152.0*t13;
  t94 = 6875.0/2304.0*t15;
  t96 = 11975.0/2304.0*eta;
  t97 = 15625.0/768.0*t28;
  t99 = 23125.0/768.0*t20;
  t101 = 78125.0/1536.0*t50;
  t102 = -t73+t74-t93-t76+t77-t94+t79+t80+t81-t82-23125.0/512.0*t4+t84+11975.0/1536.0+t96+t97+34375.0/4608.0*t46-t88-t99+40625.0/384.0*t22-t101;
  t103 = 325.0/128.0*eta;
  t104 = 1625.0/256.0*t13;
  t105 = 1875.0/128.0*t20;
  t108 = -t103+t104+t56+t57-t59+t60-t61-t62-t63+t64-975.0/256.0+t105-9375.0/256.0*t22+5625.0/256.0*t4;
  t109 = 9125.0/18432.0*eta;
  t110 = 3125.0/6144.0*t16;
  t111 = 4375.0/1536.0*t20;
  t112 = 3125.0/768.0*t22;
  t113 = 3125.0/18432.0*t15;
  t115 = 3125.0/1536.0*t28;
  t116 = 11575.0/9216.0*t13;
  t119 = 15625.0/3072.0*t50;
  t120 = -t115-t32+t35+t40-t116-4375.0/1024.0*t4+t43-t45+3125.0/18432.0*t46+t49+t119;
  t122 = 625.0/6144.0*eta;
  t123 = 11375.0/6144.0*t13;
  t124 = 33125.0/2048.0*t18;
  t125 = 625.0/6144.0*t15;
  t126 = 1875.0/2048.0*xi;
  t127 = 109375.0/6144.0*t31;
  t128 = 15625.0/2048.0*t34;
  t131 = 15625.0/1024.0*t48;
  t132 = 78125.0/4096.0*t50;
  t134 = t122+t123+t124-t125-t126+t127-t128-t65-19625.0/12288.0+21875.0/2048.0*t22-3125.0/12288.0*t46-t131-t132+34375.0/4096.0*t4;
  t135 = 15625.0/3072.0*t16;
  t136 = 2375.0/4608.0*t13;
  t137 = 3125.0/384.0*t44;
  t138 = 4375.0/192.0*t18;
  t139 = 10625.0/9216.0*t15;
  t141 = 625.0/256.0*xi;
  t142 = 10625.0/9216.0*eta;
  t143 = 15625.0/1536.0*t28;
  t144 = 3125.0/4608.0*t46;
  t145 = 15625.0/1536.0*t20;
  t147 = 78125.0/3072.0*t50;
  t148 = t135+t136-t137-t138+t139+t79+3125.0/1024.0*t4+t141+125.0/768.0-t142+t143-t144+t88-t145-21875.0/768.0*t22+t147;
  t149 = 19625.0/9216.0*t13;
  t150 = 6875.0/768.0*t18;
  t152 = 625.0/512.0*xi;
  t153 = 34375.0/18432.0*t46;
  t154 = 15625.0/1536.0*t48;
  t155 = 34375.0/3072.0*t22;
  t156 = 78125.0/6144.0*t50;
  t157 = -t135-t149+t137+t150-t139-t79+3125.0/2048.0*t4+t152+1625.0/6144.0+t142-t143+t153-t154+t145+t155-t156;
  t158 = 1625.0/1024.0*t13;
  t161 = 3125.0/6144.0*t18;
  t163 = 21875.0/6144.0*t31;
  t165 = 3125.0/2048.0*t34;
  t168 = 1.0+eta;
  t170 = t168*xi*t6;
  t171 = 125.0/192.0*t170;
  t173 = t168*t9*xi;
  t174 = 625.0/256.0*t173;
  t175 = t168*t38;
  t176 = 125.0/384.0*t175;
  t177 = 2875.0/2048.0*xi;
  t179 = 3125.0/3072.0*t48;
  t181 = t165-10475.0/18432.0*t13+4375.0/4096.0*t4+t171-t174+t176+t177+t45+15625.0/36864.0*t46+t179+15625.0/12288.0*t50;
  t183 = 625.0/2048.0*t18;
  t184 = 14625.0/2048.0*xi;
  t185 = 625.0/384.0*t175;
  t186 = 625.0/192.0*t170;
  t187 = 3125.0/256.0*t173;
  t190 = t103+t158-t56+t183-t184+t127-t128-t185-t186+t187+975.0/1024.0-t105-9375.0/1024.0*t22-5625.0/1024.0*t4;
  t191 = 109375.0/3072.0*t31;
  t192 = 15625.0/1024.0*t34;
  t193 = 10625.0/3072.0*t18;
  t195 = 7375.0/512.0*xi;
  t196 = 625.0/192.0*t175;
  t197 = 625.0/96.0*t170;
  t198 = 3125.0/128.0*t173;
  t199 = 15625.0/3072.0*t48;
  t200 = -t191+t192-t149+t76-t193+t94-t79+23125.0/2048.0*t4+t195+t196+t197-t198-11975.0/6144.0-t96-t97+t153-t199+t99+t155-t156;
  t203 = t191-t192+t136-t76+t193-t78+t79-11875.0/1024.0*t4-t195-t196-t197+t198+775.0/384.0+t85+t86-21875.0/4608.0*t46+t199-t89-t112+t147;
  t207 = -t54+t123+t56-t183+t58+t184-t127+t128+t185+t186-t187+t65-4475.0/4096.0+t66-3125.0/2048.0*t22-t68+21875.0/12288.0*t46-t132+24375.0/4096.0*t4;
  t211 = -t165-t116-625.0/512.0*t4-t171+t174-t176-t177-t45+t144-t179+t119;
  t213 = 625.0/1536.0*eta;
  t214 = 625.0/1536.0*t15;
  t215 = 15625.0/1536.0*t16;
  t218 = 15625.0/768.0*t16;
  t219 = 8125.0/2304.0*t15;
  t221 = 8125.0/2304.0*eta;
  t223 = 15625.0/384.0*t20;
  t225 = -t218-t93-t137-t150+t219+t79+3125.0/512.0*t4-t152-125.0/768.0-t221-t86+40625.0/2304.0*t46+t154+t223+34375.0/768.0*t22-t101;
  t226 = 3125.0/256.0*t4;
  t229 = t218+t75+t137+t138-t219-t79+t226-t141-1375.0/768.0+t221+t86-53125.0/2304.0*t46-t88-t223-21875.0/192.0*t22+t91;
  t232 = t213-t55-t124-t214+t126-t127+t128-t215+6625.0/3072.0+90625.0/1024.0*t22+21875.0/1536.0*t46+t131-t70-t226;
  t233 = 3125.0/6144.0*eta;
  t235 = 3125.0/2048.0*t18;
  t236 = 3125.0/6144.0*t15;
  t237 = 3125.0/2048.0*xi;
  t238 = 546875.0/6144.0*t31;
  t239 = 78125.0/2048.0*t34;
  t240 = 78125.0/6144.0*t16;
  t245 = -t233+138125.0/2048.0*t13+t235+t236-t237+t238-t239+t240-41875.0/4096.0-703125.0/2048.0*t22-234375.0/4096.0*t46+1171875.0/4096.0*t50+234375.0/4096.0*t4;
  t246 = 78125.0/3072.0*eta;
  t247 = 138125.0/1536.0*t13;
  t248 = 78125.0/3072.0*t15;
  t249 = 78125.0/3072.0*t16;
  t250 = 78125.0/512.0*t20;
  t252 = 78125.0/512.0*t28;
  t254 = 390625.0/1024.0*t50;
  t255 = 46875.0/1024.0*t4;
  t256 = -t246-t247+t248-t249+4375.0/512.0+t250+109375.0/256.0*t22-t252+15625.0/192.0*t46-t254-t255;
  t257 = 138125.0/3072.0*t13;
  t260 = 390625.0/2048.0*t50;
  t262 = t246+t257-t248+t249+6875.0/2048.0-t250-171875.0/1024.0*t22+t252-296875.0/6144.0*t46+t260-46875.0/2048.0*t4;
  t263 = 78125.0/1024.0*t18;
  t264 = 546875.0/3072.0*t31;
  t265 = 78125.0/1024.0*t34;
  t268 = 78125.0/1024.0*t48;
  t272 = 15625.0/384.0*t44;
  t273 = 21875.0/192.0*t18;
  t274 = 21875.0/576.0*t15;
  t275 = 15625.0/384.0*t25;
  t277 = 3125.0/256.0*xi;
  t278 = 21875.0/576.0*eta;
  t279 = 78125.0/384.0*t28;
  t281 = 78125.0/768.0*t48;
  t282 = 78125.0/384.0*t20;
  t285 = 30625.0/288.0*t13-t272-t273-t274+t275+15625.0/256.0*t4+t277-4375.0/384.0+t278+t279-109375.0/1152.0*t46+t281-t282-109375.0/192.0*t22+390625.0/768.0*t50;
  t286 = 48125.0/1152.0*t13;
  t287 = 34375.0/768.0*t18;
  t288 = 15625.0/512.0*t4;
  t289 = 3125.0/512.0*xi;
  t291 = 78125.0/1536.0*t48;
  t293 = 390625.0/1536.0*t50;
  t294 = -t286+t272+t287+t274-t275+t288+t289-4375.0/768.0-t278-t279+109375.0/2304.0*t46-t291+t282+171875.0/768.0*t22-t293;
  t299 = 34375.0/2304.0*t15;
  t300 = 34375.0/2304.0*eta;
  t301 = 78125.0/768.0*t28;
  t303 = 78125.0/768.0*t20;
  t305 = -t286+t272+t273+t299-t275-t288-t277+6875.0/1536.0-t300-t301+171875.0/4608.0*t46-t281+t303+109375.0/384.0*t22-t293;
  t311 = 75625.0/4608.0*t13-t272-t287-t299+t275-15625.0/1024.0*t4-t289+6875.0/3072.0+t300+t301-171875.0/9216.0*t46+t291-t303-171875.0/1536.0*t22+390625.0/3072.0*t50;
  t312 = 8125.0/1024.0*t13;

  values[0] = -1325.0/1536.0-t1-t8+t12+t17+t19+t21-15625.0/768.0*t22+t24+t26+t52;
  values[1] = t72;
  values[2] = t92;
  values[3] = t102;
  values[4] = t108;
  values[5] = 4675.0/6144.0+t109+t8-t12+t110-t19-t111+t112-t113-t26+t120;
  values[6] = t134;
  values[7] = t148;
  values[8] = t157;
  values[9] = -t122+t158+t125+t65-125.0/2048.0-3125.0/2048.0*t46;
  values[10] = -2225.0/12288.0-t109-t110+t161+t111+15625.0/6144.0*t22+t113+t26+t115-t163+t181;
  values[11] = t190;
  values[12] = t200;
  values[13] = t203;
  values[14] = t207;
  values[15] = 775.0/3072.0+t1-t17-t161-t21+3125.0/3072.0*t22-t24-t26+t29+t163+t211;
  values[16] = -t213+t104+t214+t215-125.0/512.0-3125.0/512.0*t46;
  values[17] = t225;
  values[18] = t229;
  values[19] = t232;
  values[20] = t245;
  values[21] = t256;
  values[22] = t262;
  values[23] = -t247+t263-t264+t265+14375.0/768.0+15625.0/32.0*t22+109375.0/1536.0*t46-t268-t254-109375.0/1024.0*t4;
  values[24] = t285;
  values[25] = t294;
  values[26] = t257-t263+t264-t265-104375.0/6144.0-296875.0/1024.0*t22-171875.0/6144.0*t46+t268+t260+203125.0/2048.0*t4;
  values[27] = t305;
  values[28] = t311;
  values[29] = -t312-t235+t237-t238+t239+8125.0/1024.0+46875.0/1024.0*t22-t255;
  values[30] = t233-t312-t236-t240+625.0/2048.0+15625.0/2048.0*t46;
}

// values of the derivatives in xi-eta direction
static void C_Q_UL5S_2D_DeriveXiEta(double xi, double eta, double *values)
{
  double t1, t2, t3, t6, t7, t8, t9, t10, t11, t12, t14, t15, t16, t17, t18;
  double t19, t20, t22, t23, t24, t25, t26, t27, t29, t30, t31, t35, t36, t37;
  double t39, t40, t44, t45, t48, t49, t50, t51, t52, t53, t55, t56, t57, t59;
  double t60, t61, t62, t64, t65, t66, t67, t68, t69, t70, t71, t72, t73, t74;
  double t75, t76, t78, t80, t81, t82, t83, t84, t85, t86, t87, t88, t89, t90;
  double t91, t92, t93, t94, t95, t97, t98, t100, t101, t102, t103, t104, t105;
  double t107, t109, t110, t111, t112, t113, t115, t116, t117, t118, t119;
  double t120, t121, t122, t123, t124, t126, t127, t129, t130, t131, t132;
  double t133, t134, t135, t136, t137, t138, t140, t142, t143, t144, t145;
  double t146, t147, t148, t149, t150, t151, t152, t153, t154, t155, t156;
  double t157, t158, t159, t161, t162, t163, t164, t165, t166, t167, t168;
  double t169, t170, t171, t172, t173, t174, t176, t177, t179, t181, t182;
  double t183, t185, t187, t188, t190, t191, t193, t194, t195, t196, t197;
  double t198, t199, t200, t202, t205, t207, t208, t210, t211, t212, t213;
  double t215, t216, t217, t218, t219, t221, t223, t226, t229, t230, t231;
  double t233, t234, t235, t236, t237, t238, t242, t243, t244, t245, t246;
  double t247, t248, t250, t252, t254, t257, t259, t260, t261, t262, t264;
  double t266, t269, t270, t271, t272, t273, t274, t275, t276, t277, t279;
  double t280, t283, t284, t285, t286, t287, t289, t291, t292, t296, t297;
  double t298, t300, t302, t307, t308;

  t1 = 2875.0/4608.0*xi;
  t2 = xi*xi;
  t3 = 1.0-t2;
  t6 = t3*(15.0/2.0*t2-3.0/2.0);
  t7 = 125.0/768.0*t6;
  t8 = eta*eta;
  t9 = t8*t8;
  t10 = xi*t9;
  t11 = 3125.0/768.0*t10;
  t12 = 1.0-t8;
  t14 = xi*t12*t8;
  t15 = 3125.0/512.0*t14;
  t16 = t2*eta;
  t17 = 4375.0/1536.0*t16;
  t18 = t2*xi;
  t19 = 625.0/576.0*t18;
  t20 = t18*eta;
  t22 = xi*t8;
  t23 = 4375.0/1536.0*t22;
  t24 = t2*t8;
  t25 = 625.0/256.0*t24;
  t26 = t18*t8;
  t27 = 3125.0/384.0*t26;
  t29 = t3*eta*t2;
  t30 = 3125.0/512.0*t29;
  t31 = -1049.0/2304.0-t1-t7-t11+t15+t17+t19-15625.0/1152.0*t20+t23+t25-t27+t30;
  t35 = xi*(5.0/2.0*t18-3.0/2.0*xi);
  t36 = 125.0/384.0*t35;
  t37 = xi*eta;
  t39 = t8*eta;
  t40 = 625.0/576.0*t39;
  t44 = eta*(5.0/2.0*t39-3.0/2.0*eta);
  t45 = 125.0/384.0*t44;
  t48 = t12*(15.0/2.0*t8-3.0/2.0);
  t49 = 125.0/768.0*t48;
  t50 = 125.0/192.0*t8;
  t51 = 2875.0/4608.0*eta;
  t52 = 125.0/192.0*t2;
  t53 = xi*t39;
  t55 = t2*t39;
  t56 = 3125.0/384.0*t55;
  t57 = t18*t39;
  t59 = t2*t2;
  t60 = t59*eta;
  t61 = 3125.0/768.0*t60;
  t62 = t36+4525.0/576.0*t37+t40+t45-t49+t50-t51+t52-15625.0/1152.0*t53-t56+15625.0/576.0*t57-t61;
  t64 = 17525.0/6144.0*xi;
  t65 = 625.0/1536.0*eta;
  t66 = 50375.0/1536.0*t37;
  t67 = 875.0/256.0*t2;
  t68 = 625.0/512.0*t16;
  t69 = 33125.0/2048.0*t22;
  t70 = 15625.0/768.0*t60;
  t71 = 15625.0/512.0*t29;
  t72 = 625.0/384.0*t35;
  t73 = 625.0/768.0*t6;
  t74 = 15625.0/3072.0*t10;
  t75 = 15625.0/2048.0*t14;
  t76 = 8125.0/1536.0*t18;
  t78 = 15625.0/512.0*t26;
  t80 = 78125.0/768.0*t57;
  t81 = t64+t65-t66-t67-t68-t69+t70-t71-t72+t73+t74-t75+875.0/768.0-t76+21875.0/384.0*t20+t78+90625.0/1536.0*t53-t80;
  t82 = 15625.0/384.0*t60;
  t83 = 15625.0/256.0*t29;
  t84 = 28625.0/576.0*t37;
  t85 = 5375.0/768.0*t2;
  t86 = 8125.0/768.0*t16;
  t87 = 3125.0/768.0*t8;
  t88 = 4375.0/192.0*t22;
  t89 = 3125.0/256.0*t24;
  t90 = 8125.0/2304.0*eta;
  t91 = 625.0/192.0*t35;
  t92 = 625.0/384.0*t6;
  t93 = 775.0/144.0*xi;
  t94 = 15625.0/384.0*t26;
  t95 = 15625.0/1152.0*t39;
  t97 = 15625.0/384.0*t55;
  t98 = 11875.0/1152.0*t18;
  t100 = 78125.0/576.0*t57;
  t101 = -5375.0/2304.0-t82+t83+t84+t85-t86+t87+t88-t89+t90+t91-t92-t93-t94-t95-21875.0/288.0*t53+t97+t98-53125.0/576.0*t20+t100;
  t102 = 20125.0/576.0*t37;
  t103 = 6875.0/768.0*t22;
  t104 = 11975.0/2304.0*xi;
  t105 = 15625.0/768.0*t26;
  t107 = 23125.0/2304.0*t18;
  t109 = 78125.0/1152.0*t57;
  t110 = 5375.0/2304.0+t82-t83-t102-t85+t86-t87-t103+t89-t90-t91+t92+t104+t105+t95+34375.0/1152.0*t53-t97-t107+40625.0/576.0*t20-t109;
  t111 = 325.0/128.0*xi;
  t112 = 1625.0/128.0*t37;
  t113 = 625.0/128.0*t18;
  t115 = -t111-t65+t112+t67+t68-t70+t71+t72-t73-875.0/768.0+t113-3125.0/128.0*t20;
  t116 = 9125.0/18432.0*xi;
  t117 = 3125.0/3072.0*t10;
  t118 = 3125.0/2048.0*t14;
  t119 = 4375.0/4608.0*t18;
  t120 = 3125.0/1152.0*t20;
  t121 = 3125.0/6144.0*t22;
  t122 = 3125.0/1536.0*t26;
  t123 = 1049.0/2304.0+t116+t7-t117+t118-t17-t119+t120-t121-t25-t122-t30;
  t124 = 11575.0/4608.0*t37;
  t126 = 15625.0/2304.0*t57;
  t127 = -t36-t124-t40-t45+t49-t50+t51-t52+3125.0/4608.0*t53+t56+t126+t61;
  t129 = 625.0/6144.0*xi;
  t130 = 17525.0/6144.0*eta;
  t131 = 11375.0/3072.0*t37;
  t132 = 33125.0/2048.0*t16;
  t133 = 875.0/256.0*t8;
  t134 = 625.0/2048.0*t22;
  t135 = 15625.0/3072.0*t60;
  t136 = 15625.0/2048.0*t29;
  t137 = 625.0/384.0*t44;
  t138 = 625.0/768.0*t48;
  t140 = 8125.0/1536.0*t39;
  t142 = 15625.0/512.0*t55;
  t143 = 78125.0/3072.0*t57;
  t144 = t129-t130+t131+t132+t133-t134-t135+t136+t137-t138+t74-t75-875.0/768.0+21875.0/3072.0*t20+t140-3125.0/3072.0*t53-t142-t143;
  t145 = 15625.0/1536.0*t10;
  t146 = 15625.0/1024.0*t14;
  t147 = 2375.0/2304.0*t37;
  t148 = 3125.0/768.0*t2;
  t149 = 4375.0/192.0*t16;
  t150 = 5375.0/768.0*t8;
  t151 = 10625.0/3072.0*t22;
  t152 = 775.0/144.0*eta;
  t153 = 625.0/192.0*t44;
  t154 = 625.0/384.0*t48;
  t155 = 10625.0/9216.0*xi;
  t156 = 15625.0/1536.0*t26;
  t157 = 11875.0/1152.0*t39;
  t158 = 3125.0/1152.0*t53;
  t159 = 15625.0/4608.0*t18;
  t161 = 78125.0/2304.0*t57;
  t162 = 5375.0/2304.0-t145+t146+t147-t148-t149-t150+t151+t89+t152-t153+t154-t155+t156-t157-t158+t97-t159-21875.0/1152.0*t20+t161;
  t163 = 19625.0/4608.0*t37;
  t164 = 6875.0/768.0*t16;
  t165 = 11975.0/2304.0*eta;
  t166 = 23125.0/2304.0*t39;
  t167 = 34375.0/4608.0*t53;
  t168 = 15625.0/768.0*t55;
  t169 = 34375.0/4608.0*t20;
  t170 = 78125.0/4608.0*t57;
  t171 = -5375.0/2304.0+t145-t146-t163+t148+t164+t150-t151-t89-t165+t153-t154+t155-t156+t166+t167-t168+t159+t169-t170;
  t172 = 325.0/128.0*eta;
  t173 = 1625.0/512.0*t37;
  t174 = 625.0/128.0*t39;
  t176 = -t129+t172+t173-t133+t134-t137+t138-t74+t75+875.0/768.0-t174-3125.0/512.0*t53;
  t177 = 3125.0/6144.0*t16;
  t179 = 3125.0/2048.0*t29;
  t181 = -1049.0/2304.0-t116+t117-t118+t177+t119+15625.0/9216.0*t20+t121+t25+t122-t179-10475.0/9216.0*t37;
  t182 = 4375.0/4608.0*t39;
  t183 = 9125.0/18432.0*eta;
  t185 = 3125.0/1536.0*t55;
  t187 = 3125.0/3072.0*t60;
  t188 = t182+t45-t49+t50-t7+t36-t183+t52+15625.0/9216.0*t53+t185+15625.0/9216.0*t57+t187;
  t190 = 625.0/6144.0*eta;
  t191 = 625.0/2048.0*t16;
  t193 = t111-t190+t173-t67+t191-t135+t136-t72+t73+875.0/768.0-t113-3125.0/512.0*t20;
  t194 = 15625.0/1536.0*t60;
  t195 = 15625.0/1024.0*t29;
  t196 = 10625.0/3072.0*t16;
  t197 = 10625.0/9216.0*eta;
  t198 = 15625.0/4608.0*t39;
  t199 = 15625.0/1536.0*t55;
  t200 = -5375.0/2304.0+t194-t195-t163+t85-t196+t87+t103-t89+t197+t91-t92-t104-t105+t198+t167-t199+t107+t169-t170;
  t202 = 5375.0/2304.0-t194+t195+t147-t85+t196-t87-t88+t89-t197-t91+t92+t93+t94-t198-21875.0/1152.0*t53+t199-t98-t120+t161;
  t205 = -t64+t190+t131+t67-t191+t69+t135-t136+t72-t73-t74+t75-875.0/768.0+t76-3125.0/3072.0*t20-t78+21875.0/3072.0*t53-t143;
  t207 = 1049.0/2304.0+t1+t11-t15-t177-t19+3125.0/4608.0*t20-t23-t25+t27+t179-t124;
  t208 = -t182-t45+t49-t50+t7-t36+t183-t52+t158-t185+t126-t187;
  t210 = 625.0/1536.0*xi;
  t211 = 625.0/512.0*t22;
  t212 = 15625.0/768.0*t10;
  t213 = 15625.0/512.0*t14;
  t215 = -t210-t172+t112+t133+t211+t137-t138-t212+t213-875.0/768.0+t174-3125.0/128.0*t53;
  t216 = 15625.0/384.0*t10;
  t217 = 15625.0/256.0*t14;
  t218 = 8125.0/768.0*t22;
  t219 = 8125.0/2304.0*xi;
  t221 = 15625.0/1152.0*t18;
  t223 = 5375.0/2304.0+t216-t217-t102-t148-t164-t150+t218+t89+t165-t153+t154-t219-t94-t166+40625.0/576.0*t53+t168+t221+34375.0/1152.0*t20-t109;
  t226 = -5375.0/2304.0-t216+t217+t84+t148+t149+t150-t218-t89-t152+t153-t154+t219+t94+t157-53125.0/576.0*t53-t97-t221-21875.0/288.0*t20+t100;
  t229 = t210+t130-t66-t132-t133-t211+t135-t136-t137+t138+t212-t213+875.0/768.0+90625.0/1536.0*t20-t140+21875.0/384.0*t53+t142-t80;
  t230 = 3125.0/6144.0*xi;
  t231 = 3125.0/6144.0*eta;
  t233 = 3125.0/2048.0*t16;
  t234 = 3125.0/2048.0*t22;
  t235 = 78125.0/3072.0*t60;
  t236 = 78125.0/2048.0*t29;
  t237 = 78125.0/3072.0*t10;
  t238 = 78125.0/2048.0*t14;
  t242 = -t230-t231+138125.0/1024.0*t37+t233+t234-t235+t236-t237+t238-234375.0/1024.0*t20-234375.0/1024.0*t53+390625.0/1024.0*t57;
  t243 = 78125.0/3072.0*xi;
  t244 = 138125.0/768.0*t37;
  t245 = 78125.0/1024.0*t22;
  t246 = 78125.0/1536.0*t10;
  t247 = 78125.0/1024.0*t14;
  t248 = 78125.0/1536.0*t18;
  t250 = 78125.0/512.0*t26;
  t252 = 390625.0/768.0*t57;
  t254 = 138125.0/1536.0*t37;
  t257 = 390625.0/1536.0*t57;
  t259 = 78125.0/3072.0*eta;
  t260 = 78125.0/1024.0*t16;
  t261 = 78125.0/1536.0*t60;
  t262 = 78125.0/1024.0*t29;
  t264 = 78125.0/1536.0*t39;
  t266 = 78125.0/512.0*t55;
  t269 = 15625.0/768.0*t2;
  t270 = 21875.0/192.0*t16;
  t271 = 15625.0/768.0*t8;
  t272 = 21875.0/192.0*t22;
  t273 = 15625.0/256.0*t24;
  t274 = 21875.0/576.0*eta;
  t275 = 21875.0/576.0*xi;
  t276 = 78125.0/384.0*t26;
  t277 = 78125.0/1152.0*t39;
  t279 = 78125.0/384.0*t55;
  t280 = 78125.0/1152.0*t18;
  t283 = 15625.0/2304.0+30625.0/144.0*t37-t269-t270-t271-t272+t273+t274+t275+t276-t277-109375.0/288.0*t53+t279-t280-109375.0/288.0*t20+390625.0/576.0*t57;
  t284 = 48125.0/576.0*t37;
  t285 = 34375.0/768.0*t16;
  t286 = 34375.0/2304.0*eta;
  t287 = 78125.0/2304.0*t39;
  t289 = 78125.0/768.0*t55;
  t291 = 390625.0/1152.0*t57;
  t292 = -15625.0/2304.0-t284+t269+t285+t271+t272-t273-t286-t275-t276+t287+109375.0/576.0*t53-t289+t280+171875.0/1152.0*t20-t291;
  t296 = 34375.0/768.0*t22;
  t297 = 34375.0/2304.0*xi;
  t298 = 78125.0/768.0*t26;
  t300 = 78125.0/2304.0*t18;
  t302 = -15625.0/2304.0-t284+t269+t270+t271+t296-t273-t274-t297-t298+t277+171875.0/1152.0*t53-t279+t300+109375.0/576.0*t20-t291;
  t307 = 15625.0/2304.0+75625.0/2304.0*t37-t269-t285-t271-t296+t273+t286+t297+t298-t287-171875.0/2304.0*t53+t289-t300-171875.0/2304.0*t20+390625.0/2304.0*t57;
  t308 = 8125.0/512.0*t37;

  values[0] = t31+t62;
  values[1] = t81;
  values[2] = t101;
  values[3] = t110;
  values[4] = t115;
  values[5] = t123+t127;
  values[6] = t144;
  values[7] = t162;
  values[8] = t171;
  values[9] = t176;
  values[10] = t181+t188;
  values[11] = t193;
  values[12] = t200;
  values[13] = t202;
  values[14] = t205;
  values[15] = t207+t208;
  values[16] = t215;
  values[17] = t223;
  values[18] = t226;
  values[19] = t229;
  values[20] = t242;
  values[21] = -t243-t244+t245+t246-t247+t248+109375.0/384.0*t20-t250+15625.0/48.0*t53-t252;
  values[22] = t243+t254-t245-t246+t247-t248-171875.0/1536.0*t20+t250-296875.0/1536.0*t53+t257;
  values[23] = -t259-t244+t260+t261-t262+15625.0/48.0*t20+t264+109375.0/384.0*t53-t266-t252;
  values[24] = t283;
  values[25] = t292;
  values[26] = t259+t254-t260-t261+t262-296875.0/1536.0*t20-t264-171875.0/1536.0*t53+t266+t257;
  values[27] = t302;
  values[28] = t307;
  values[29] = t231-t308-t233+t235-t236+15625.0/512.0*t20;
  values[30] = t230-t308-t234+t237-t238+15625.0/512.0*t53;
}

// values of the derivatives in eta-eta direction
static void C_Q_UL5S_2D_DeriveEtaEta(double xi, double eta, double *values)
{
  double t1, t2, t3, t4, t5, t6, t7, t8, t9, t10, t11, t12, t14, t15, t16;
  double t17, t18, t19, t20, t22, t23, t25, t26, t27, t30, t31, t32, t34, t35;
  double t38, t39, t40, t41, t42, t44, t46, t47, t50, t51, t52, t54, t55, t56;
  double t57, t58, t59, t60, t61, t63, t65, t66, t67, t68, t69, t70, t71, t72;
  double t73, t74, t75, t76, t77, t79, t81, t82, t83, t84, t86, t87, t90, t91;
  double t92, t95, t97, t98, t101, t103, t104, t106, t107, t108, t109, t110;
  double t111, t112, t113, t114, t115, t117, t118, t119, t120, t121, t122;
  double t123, t124, t125, t126, t128, t130, t131, t133, t134, t135, t136;
  double t137, t138, t139, t140, t142, t143, t144, t145, t146, t147, t148;
  double t150, t151, t152, t153, t154, t156, t157, t158, t159, t160, t161;
  double t162, t163, t164, t167, t168, t169, t170, t172, t174, t179, t181;
  double t182, t185, t186, t187, t189, t190, t191, t194, t198, t202, t204;
  double t205, t206, t207, t208, t209, t210, t213, t214, t215, t216, t218;
  double t219, t220, t221, t224, t228, t232, t233, t235, t236, t237, t238;
  double t239, t240, t245, t246, t247, t248, t249, t251, t253, t256, t259;
  double t262, t263, t264, t266, t268, t269, t270, t272, t273, t274, t275;
  double t276, t278, t279, t280, t282, t285, t286, t287, t288, t289, t290;
  double t292, t294, t295, t299, t300, t301, t302, t305, t311, t312;

  t1 = 2875.0/4608.0*xi;
  t2 = xi*xi;
  t3 = t2*eta;
  t4 = 4375.0/1536.0*t3;
  t5 = t2*xi;
  t6 = 4375.0/4608.0*t5;
  t7 = t5*eta;
  t8 = 625.0/384.0*t7;
  t9 = eta*eta;
  t10 = xi*t9;
  t11 = 625.0/192.0*t10;
  t12 = t2*t9;
  t14 = t5*t9;
  t15 = 3125.0/384.0*t14;
  t16 = 1.0-t2;
  t17 = t9*eta;
  t18 = t16*t17;
  t19 = 21875.0/1536.0*t18;
  t20 = 1.0-t9;
  t22 = t16*t20*eta;
  t23 = 3125.0/512.0*t22;
  t25 = xi*eta;
  t26 = 125.0/96.0*t25;
  t27 = 1.0-xi;
  t30 = 15.0/2.0*t9-3.0/2.0;
  t31 = t27*eta*t30;
  t32 = 125.0/192.0*t31;
  t34 = t27*t20*eta;
  t35 = 625.0/256.0*t34;
  t38 = 5.0/2.0*t17-3.0/2.0*eta;
  t39 = t27*t38;
  t40 = 125.0/384.0*t39;
  t41 = t16*t5;
  t42 = 3125.0/1536.0*t41;
  t44 = 875.0/512.0*eta;
  t46 = t2*t2;
  t47 = t46*t9;
  t50 = t46*eta;
  t51 = 3125.0/768.0*t50;
  t52 = t26-t32+t35-t40+t42+625.0/128.0*t9-t44+4525.0/1152.0*t2+15625.0/768.0*t47-15625.0/4608.0*t46-t51;
  t54 = 625.0/1536.0*xi;
  t55 = 50375.0/3072.0*t2;
  t56 = 625.0/1536.0*t5;
  t57 = 33125.0/2048.0*t3;
  t58 = 1875.0/2048.0*eta;
  t59 = 15625.0/1536.0*t41;
  t60 = 109375.0/6144.0*t18;
  t61 = 15625.0/2048.0*t22;
  t63 = 15625.0/1024.0*t50;
  t65 = 78125.0/1024.0*t47;
  t66 = 3125.0/256.0*t9;
  t67 = t54-t55-t56-t57+t58-t59-t60+t61+6625.0/3072.0+21875.0/1536.0*t46+t63+90625.0/1024.0*t12-t65-t66;
  t68 = 15625.0/768.0*t41;
  t69 = 28625.0/1152.0*t2;
  t70 = 8125.0/2304.0*t5;
  t71 = 3125.0/384.0*t25;
  t72 = 4375.0/192.0*t3;
  t73 = 3125.0/384.0*t7;
  t74 = 8125.0/2304.0*xi;
  t75 = 625.0/256.0*eta;
  t76 = 15625.0/768.0*t50;
  t77 = 15625.0/384.0*t10;
  t79 = 15625.0/384.0*t14;
  t81 = 78125.0/768.0*t47;
  t82 = t68+t69-t70+t71+t72-t73+t74+t66-t75-1375.0/768.0-t76-t77-21875.0/192.0*t12+t79-53125.0/2304.0*t46+t81;
  t83 = 20125.0/1152.0*t2;
  t84 = 6875.0/768.0*t3;
  t86 = 625.0/512.0*eta;
  t87 = 15625.0/1536.0*t50;
  t90 = 78125.0/1536.0*t47;
  t91 = -t68-t83+t70-t71-t84+t73-t74+3125.0/512.0*t9-t86-125.0/768.0+t87+t77+34375.0/768.0*t12-t79+40625.0/2304.0*t46-t90;
  t92 = 1625.0/256.0*t2;
  t95 = 3125.0/6144.0*t3;
  t97 = 21875.0/6144.0*t18;
  t98 = 3125.0/2048.0*t22;
  t101 = 1.0+xi;
  t103 = t101*eta*t30;
  t104 = 125.0/192.0*t103;
  t106 = t101*t20*eta;
  t107 = 625.0/256.0*t106;
  t108 = t101*t38;
  t109 = 125.0/384.0*t108;
  t110 = 2875.0/2048.0*eta;
  t111 = 11575.0/9216.0*t2;
  t112 = 15625.0/3072.0*t47;
  t113 = 3125.0/4608.0*t46;
  t114 = 3125.0/3072.0*t50;
  t115 = -t26-t42-625.0/512.0*t9-t104+t107-t109-t110-t111+t112+t113-t114;
  t117 = 17525.0/6144.0*xi;
  t118 = 11375.0/6144.0*t2;
  t119 = 33125.0/6144.0*t5;
  t120 = 875.0/128.0*t25;
  t121 = 625.0/2048.0*t3;
  t122 = 14625.0/2048.0*eta;
  t123 = 15625.0/6144.0*t41;
  t124 = 625.0/384.0*t108;
  t125 = 625.0/192.0*t103;
  t126 = 3125.0/256.0*t106;
  t128 = 8125.0/512.0*t10;
  t130 = 15625.0/512.0*t14;
  t131 = 78125.0/4096.0*t47;
  t133 = -t117+t118+t119+t120-t121+t122+t123+t124+t125-t126-t60+t61-4475.0/4096.0+21875.0/12288.0*t46+t128-3125.0/2048.0*t12-t130-t131+24375.0/4096.0*t9;
  t134 = 109375.0/3072.0*t18;
  t135 = 15625.0/1024.0*t22;
  t136 = 2375.0/4608.0*t2;
  t137 = 4375.0/576.0*t5;
  t138 = 5375.0/384.0*t25;
  t139 = 10625.0/3072.0*t3;
  t140 = 775.0/144.0*xi;
  t142 = 7375.0/512.0*eta;
  t143 = 625.0/192.0*t108;
  t144 = 625.0/96.0*t103;
  t145 = 3125.0/128.0*t106;
  t146 = 15625.0/3072.0*t50;
  t147 = 11875.0/384.0*t10;
  t148 = 3125.0/768.0*t12;
  t150 = 78125.0/3072.0*t47;
  t151 = t134-t135+t136-t137-t138+t139+t73+t140-11875.0/1024.0*t9-t142-t143-t144+t145+775.0/384.0+t146-t147-t148+t79-21875.0/4608.0*t46+t150;
  t152 = 19625.0/9216.0*t2;
  t153 = 6875.0/2304.0*t5;
  t154 = 11975.0/2304.0*xi;
  t156 = 23125.0/768.0*t10;
  t157 = 34375.0/3072.0*t12;
  t158 = 15625.0/768.0*t14;
  t159 = 34375.0/18432.0*t46;
  t160 = 78125.0/6144.0*t47;
  t161 = -t134+t135-t152+t153+t138-t139-t73-t154+23125.0/2048.0*t9+t142+t143+t144-t145-11975.0/6144.0-t146+t156+t157-t158+t159-t160;
  t162 = 325.0/128.0*xi;
  t163 = 1625.0/1024.0*t2;
  t164 = 1875.0/128.0*t10;
  t167 = t162+t163-t120+t121-t122-t124-t125+t126+t60-t61+975.0/1024.0-t164-9375.0/1024.0*t12-5625.0/1024.0*t9;
  t168 = 9125.0/18432.0*xi;
  t169 = 3125.0/18432.0*t5;
  t170 = 4375.0/1536.0*t10;
  t172 = 3125.0/1536.0*t14;
  t174 = 3125.0/6144.0*t41;
  t179 = t26-t174+4375.0/4096.0*t9+t104-t107+t109+t110-10475.0/18432.0*t2+15625.0/12288.0*t47+15625.0/36864.0*t46+t114;
  t181 = 625.0/6144.0*xi;
  t182 = 625.0/6144.0*t5;
  t185 = 15625.0/3072.0*t41;
  t186 = 10625.0/9216.0*t5;
  t187 = 10625.0/9216.0*xi;
  t189 = 15625.0/1536.0*t10;
  t190 = 15625.0/1536.0*t14;
  t191 = -t185-t152-t186+t71+t84-t73+t187+3125.0/2048.0*t9+t86+1625.0/6144.0-t87+t189+t157-t190+t159-t160;
  t194 = t185+t136+t186-t71-t72+t73-t187+3125.0/1024.0*t9+t75+125.0/768.0+t76-t189-21875.0/768.0*t12+t190-t113+t150;
  t198 = t181+t118-t182+t57-t58-t123+t60-t61-19625.0/12288.0-3125.0/12288.0*t46-t63+21875.0/2048.0*t12-t131+34375.0/4096.0*t9;
  t202 = -t26+t32-t35+t40+t174-4375.0/1024.0*t9+t44-t111+t112+3125.0/18432.0*t46+t51;
  t204 = 625.0/512.0*t3;
  t205 = 4125.0/512.0*eta;
  t206 = 625.0/384.0*t39;
  t207 = 625.0/192.0*t31;
  t208 = 3125.0/256.0*t34;
  t209 = 109375.0/1536.0*t18;
  t210 = 15625.0/512.0*t22;
  t213 = -t162+t92+t120+t204-t205-t206-t207+t208+t209-t210-975.0/256.0+t164-9375.0/256.0*t12+5625.0/256.0*t9;
  t214 = 109375.0/768.0*t18;
  t215 = 15625.0/256.0*t22;
  t216 = 8125.0/768.0*t3;
  t218 = 125.0/8.0*eta;
  t219 = 625.0/192.0*t39;
  t220 = 625.0/96.0*t31;
  t221 = 3125.0/128.0*t34;
  t224 = -t214+t215-t83-t153-t138+t216+t73+t154-23125.0/512.0*t9+t218+t219+t220-t221+11975.0/1536.0-t76-t156+40625.0/384.0*t12+t158+34375.0/4608.0*t46-t90;
  t228 = t214-t215+t69+t137+t138-t216-t73-t140+11875.0/256.0*t9-t218-t219-t220+t221-775.0/96.0+t76+t147-53125.0/384.0*t12-t79-21875.0/1152.0*t46+t81;
  t232 = t117-t55-t119-t120-t204+t205-t123+t206+t207-t208-t209+t210+8575.0/2048.0+90625.0/6144.0*t46-t128+21875.0/256.0*t12+t130-t65-24375.0/1024.0*t9;
  t233 = 3125.0/6144.0*xi;
  t235 = 3125.0/6144.0*t5;
  t236 = 3125.0/2048.0*t3;
  t237 = 3125.0/2048.0*eta;
  t238 = 78125.0/6144.0*t41;
  t239 = 546875.0/6144.0*t18;
  t240 = 78125.0/2048.0*t22;
  t245 = -t233+138125.0/2048.0*t2+t235+t236-t237+t238+t239-t240-41875.0/4096.0-234375.0/4096.0*t46-703125.0/2048.0*t12+1171875.0/4096.0*t47+234375.0/4096.0*t9;
  t246 = 138125.0/1536.0*t2;
  t247 = 78125.0/1024.0*t3;
  t248 = 546875.0/3072.0*t18;
  t249 = 78125.0/1024.0*t22;
  t251 = 78125.0/1024.0*t50;
  t253 = 390625.0/1024.0*t47;
  t256 = 138125.0/3072.0*t2;
  t259 = 390625.0/2048.0*t47;
  t262 = 78125.0/3072.0*xi;
  t263 = 78125.0/3072.0*t5;
  t264 = 78125.0/3072.0*t41;
  t266 = 78125.0/512.0*t10;
  t268 = 78125.0/512.0*t14;
  t269 = 46875.0/1024.0*t9;
  t270 = -t262-t246+t263-t264+4375.0/512.0+15625.0/192.0*t46+t266+109375.0/256.0*t12-t268-t253-t269;
  t272 = 21875.0/576.0*t5;
  t273 = 15625.0/384.0*t25;
  t274 = 21875.0/192.0*t3;
  t275 = 15625.0/384.0*t7;
  t276 = 21875.0/576.0*xi;
  t278 = 3125.0/256.0*eta;
  t279 = 78125.0/768.0*t50;
  t280 = 78125.0/384.0*t10;
  t282 = 78125.0/384.0*t14;
  t285 = 30625.0/288.0*t2-t272-t273-t274+t275+t276+15625.0/256.0*t9+t278-4375.0/384.0+t279-t280-109375.0/192.0*t12+t282-109375.0/1152.0*t46+390625.0/768.0*t47;
  t286 = 48125.0/1152.0*t2;
  t287 = 34375.0/2304.0*t5;
  t288 = 34375.0/2304.0*xi;
  t289 = 15625.0/512.0*t9;
  t290 = 78125.0/768.0*t10;
  t292 = 78125.0/768.0*t14;
  t294 = 390625.0/1536.0*t47;
  t295 = -t286+t287+t273+t274-t275-t288-t289-t278+6875.0/1536.0-t279+t290+109375.0/384.0*t12-t292+171875.0/4608.0*t46-t294;
  t299 = t262+t256-t263+t264+6875.0/2048.0-296875.0/6144.0*t46-t266-171875.0/1024.0*t12+t268+t259-46875.0/2048.0*t9;
  t300 = 34375.0/768.0*t3;
  t301 = 3125.0/512.0*eta;
  t302 = 78125.0/1536.0*t50;
  t305 = -t286+t272+t273+t300-t275-t276+t289+t301-4375.0/768.0-t302+t280+171875.0/768.0*t12-t282+109375.0/2304.0*t46-t294;
  t311 = 75625.0/4608.0*t2-t287-t273-t300+t275+t288-15625.0/1024.0*t9-t301+6875.0/3072.0+t302-t290-171875.0/1536.0*t12+t292-171875.0/9216.0*t46+390625.0/3072.0*t47;
  t312 = 8125.0/1024.0*t2;

  values[0] = -1325.0/1536.0-t1+t4+t6+t8+t11-15625.0/768.0*t12-t15+t19-t23+t52;
  values[1] = t67;
  values[2] = t82;
  values[3] = t91;
  values[4] = -t54+t92+t56+t59-125.0/512.0-3125.0/512.0*t46;
  values[5] = 775.0/3072.0+t1-t95-t6-t8-t11+3125.0/3072.0*t12+t15+t97-t98+t115;
  values[6] = t133;
  values[7] = t151;
  values[8] = t161;
  values[9] = t167;
  values[10] = -2225.0/12288.0-t168+t95+t169+t8+t170+15625.0/6144.0*t12+t172-t97+t98+t179;
  values[11] = -t181+t163+t182+t123-125.0/2048.0-3125.0/2048.0*t46;
  values[12] = t191;
  values[13] = t194;
  values[14] = t198;
  values[15] = 4675.0/6144.0+t168-t4-t169-t8-t170+t148-t172-t19+t23+t202;
  values[16] = t213;
  values[17] = t224;
  values[18] = t228;
  values[19] = t232;
  values[20] = t245;
  values[21] = -t246+t247-t248+t249+14375.0/768.0+109375.0/1536.0*t46-t251+15625.0/32.0*t12-t253-109375.0/1024.0*t9;
  values[22] = t256-t247+t248-t249-104375.0/6144.0-171875.0/6144.0*t46+t251-296875.0/1024.0*t12+t259+203125.0/2048.0*t9;
  values[23] = t270;
  values[24] = t285;
  values[25] = t295;
  values[26] = t299;
  values[27] = t305;
  values[28] = t311;
  values[29] = t233-t312-t235-t238+625.0/2048.0+15625.0/2048.0*t46;
  values[30] = -t312-t236+t237-t239+t240+8125.0/1024.0+46875.0/1024.0*t12-t269;
}

// ***********************************************************************

TBaseFunct2D *BF_C_Q_UL5S_2D_Obj = new TBaseFunct2D
        (31, BF_C_Q_UL5S_2D, BFUnitSquare, 
         C_Q_UL5S_2D_Funct, C_Q_UL5S_2D_DeriveXi,
         C_Q_UL5S_2D_DeriveEta, C_Q_UL5S_2D_DeriveXiXi,
         C_Q_UL5S_2D_DeriveXiEta, C_Q_UL5S_2D_DeriveEtaEta, 5, 5,
         0, NULL);