parmoon/html/_b_f___c___q___u_l6_s__2_d_8h_source.html

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

 // UL6S element with bubbles, conforming, 2D

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


 // base function values

 static void C_Q_UL6S_2D_Funct(double xi, double eta, double *values)

 {

   double t1 = eta*eta;

   double t2 = t1*t1;

   double t3 = 153.0/280.0*t2;

   double t4 = xi*xi;

   double t5 = 3303.0/5600.0*t4;

   double t7 = 3303.0/5600.0*t1;

   double t8 = t4*xi;

   double t9 = 3.0/4.0*t8;

   double t10 = t1*eta;

   double t11 = 3.0/4.0*t10;

   double t12 = xi*t10;

   double t13 = 27.0/128.0*t12;

   double t14 = t4*t4;

   double t15 = t14*t1;

   double t16 = 261.0/1280.0*t15;

   double t17 = t14*xi;

   double t18 = t17*eta;

   double t20 = t14*eta;

   double t21 = 2007.0/4480.0*t20;

   double t22 = t8*eta;

   double t23 = 27.0/128.0*t22;

   double t24 = t17*t1;

   double t27 = 1.0-t4;

   double t28 = 1.0-t1;

   double t29 = t27*t28;

   double t30 = t29*t2;

   double t31 = 243.0/160.0*t30;

   double t32 = t2*eta;

   double t33 = 27.0/40.0*t32;

   double t34 = t8*t32;

   double t36 = t14*t32;

   double t37 = 243.0/256.0*t36;

   double t38 = xi*t32;

   double t40 = t17*t2;

   double t41 = 243.0/256.0*t40;

   double t42 = t4*t32;

   double t44 = t3-t5-663.0/5600.0*xi-t7+t9+t11+t13-t16-81.0/640.0*t18-t21+t23+1431.0/1280.0*t24-663.0/5600.0*eta+t31-t33-567.0/256.0*t34-t37-81.0/640.0*t38-t41+1431.0/1280.0*t42;

   double t45 = t4*t2;

   double t46 = 261.0/1280.0*t45;

   double t47 = t14*t2;

   double t48 = 81.0/256.0*t47;

   double t49 = t8*t2;

   double t50 = 189.0/256.0*t49;

   double t51 = t14*t10;

   double t52 = 189.0/256.0*t51;

   double t53 = t17*t10;

   double t55 = xi*t2;

   double t56 = 2007.0/4480.0*t55;

   double t57 = t8*t10;

   double t59 = t4*t10;

   double t61 = t17*t32;

   double t63 = xi*eta;

   double t65 = t4*eta;

   double t67 = xi*t1;

   double t69 = t8*t1;

   double t71 = t4*t1;

   double t73 = 153.0/280.0*t14;

   double t74 = 27.0/40.0*t17;

   double t79 = 3.0/8.0+35.0/8.0*t14-15.0/4.0*t4;

   double t80 = (1.0-eta)*t27*t79;

   double t81 = 81.0/700.0*t80;

   double t82 = t29*t14;

   double t83 = 243.0/160.0*t82;

   double t88 = 3.0/8.0+35.0/8.0*t2-15.0/4.0*t1;

   double t89 = (1.0-xi)*t28*t88;

   double t90 = 81.0/700.0*t89;

   double t91 = 243.0/2800.0-t46+t48+t50+t52-567.0/256.0*t53-t56+441.0/256.0*t57-309.0/256.0*t59+729.0/256.0*t61-19.0/320.0*t63+12127.0/22400.0*t65+12127.0/22400.0*t67-309.0/256.0*t69+437.0/1280.0*t71+t73-t74-t81+t83-t90;

   double t93 = 81.0/400.0*t2;

   double t94 = 9909.0/2800.0*t4;

   double t95 = 9.0/20.0*xi;

   double t96 = 81.0/400.0*t1;

   double t97 = 9.0/2.0*t8;

   double t98 = 27.0/20.0*t10;

   double t100 = 81.0/80.0*t15;

   double t102 = 1269.0/560.0*t20;

   double t106 = 729.0/400.0*t30;

   double t107 = 243.0/200.0*t32;

   double t110 = -729.0/2800.0-t93+t94+t95+t96-t97-t98+189.0/160.0*t12+t100+81.0/80.0*t18+t102-9.0/8.0*t22-1053.0/160.0*t24+1107.0/2800.0*eta-t106+t107+243.0/16.0*t34-243.0/160.0*t38;

   double t111 = 729.0/160.0*t40;

   double t112 = 243.0/200.0*t42;

   double t113 = 81.0/400.0*t45;

   double t114 = 81.0/16.0*t49;

   double t116 = 81.0/160.0*t55;

   double t117 = 189.0/16.0*t57;

   double t118 = 27.0/20.0*t59;

   double t119 = 2187.0/160.0*t61;

   double t120 = 9.0/80.0*t63;

   double t124 = 243.0/200.0*t71;

   double t125 = 459.0/140.0*t14;

   double t126 = 81.0/20.0*t17;

   double t127 = 243.0/350.0*t80;

   double t128 = 729.0/80.0*t82;

   double t129 = t111-t112+t113-t114+1701.0/160.0*t53+t116-t117+t118-t119+t120-1863.0/700.0*t65-117.0/160.0*t67+117.0/16.0*t69-t124-t125+t126+t127-t128;

   double t131 = 9909.0/1120.0*t4;

   double t132 = 9.0/8.0*xi;

   double t133 = 45.0/4.0*t8;

   double t134 = 567.0/128.0*t12;

   double t135 = 567.0/256.0*t15;

   double t137 = 4509.0/896.0*t20;

   double t140 = 729.0/1120.0*eta;

   double t142 = 2187.0/256.0*t36;

   double t143 = 729.0/128.0*t38;

   double t144 = 2187.0/256.0*t40;

   double t145 = 2187.0/256.0*t42;

   double t146 = 729.0/256.0*t45;

   double t147 = 729.0/1120.0-t131-t132+t133-t134-t135-405.0/128.0*t18-t137+423.0/128.0*t22+4131.0/256.0*t24-t140-8019.0/256.0*t34+t142+t143-t144-t145+t146;

   double t148 = 729.0/256.0*t47;

   double t149 = 2673.0/256.0*t49;

   double t150 = 1701.0/256.0*t51;

   double t151 = 5103.0/256.0*t53;

   double t152 = 243.0/128.0*t55;

   double t153 = 6237.0/256.0*t57;

   double t154 = 1701.0/256.0*t59;

   double t155 = 6561.0/256.0*t61;

   double t156 = 9.0/64.0*t63;

   double t157 = 25461.0/4480.0*t65;

   double t160 = 567.0/256.0*t71;

   double t161 = 459.0/56.0*t14;

   double t162 = 81.0/8.0*t17;

   double t163 = 243.0/140.0*t80;

   double t164 = 729.0/32.0*t82;

   double t165 = -t148+t149-t150-t151-t152+t153+t154+t155-t156+t157+243.0/128.0*t67-4617.0/256.0*t69+t160+t161-t162-t163+t164;

   double t167 = 9.0/16.0*t2;

   double t169 = 3.0/2.0*xi;

   double t170 = t1/16.0;

   double t171 = 15.0*t8;

   double t184 = 243.0/32.0*t40;

   double t186 = t167+3303.0/280.0*t4+t169-t170-t171+21.0/16.0*t10+63.0/32.0*t12+45.0/16.0*t15+81.0/16.0*t18+45.0/7.0*t20-45.0/8.0*t22-675.0/32.0*t24+26.0/35.0*eta-27.0/16.0*t32+405.0/16.0*t34-243.0/16.0*t36-81.0/32.0*t38+t184+135.0/8.0*t42;

   double t187 = 45.0/8.0*t45;

   double t188 = 81.0/16.0*t47;

   double t189 = 135.0/16.0*t49;

   double t192 = 27.0/32.0*t55;

   double t193 = 315.0/16.0*t57;

   double t195 = 729.0/32.0*t61;

   double t196 = 9.0/16.0*t63;

   double t200 = 11.0/4.0*t71;

   double t202 = 27.0/2.0*t17;

   double t205 = -243.0/280.0-t187+t188-t189+189.0/16.0*t51+567.0/32.0*t53+t192-t193-105.0/8.0*t59-t195+t196-251.0/35.0*t65-75.0/32.0*t67+375.0/16.0*t69-t200-153.0/14.0*t14+t202+81.0/35.0*t80-243.0/8.0*t82;

   double t211 = 1701.0/256.0*t34;

   double t213 = 729.0/256.0*t40;

   double t214 = 729.0/1120.0-t131-t132+t133+189.0/128.0*t12-t135-567.0/128.0*t18-t137+657.0/128.0*t22+3969.0/256.0*t24-t140-t211+t142-243.0/128.0*t38-t213-t145+t146;

   double t215 = 567.0/256.0*t49;

   double t216 = 1701.0/256.0*t53;

   double t217 = 81.0/128.0*t55;

   double t218 = 1323.0/256.0*t57;

   double t219 = 2187.0/256.0*t61;

   double t220 = 45.0/64.0*t63;

   double t223 = -t148+t215-t150-t216+t217+t218+t154+t219-t220+t157+207.0/128.0*t67-4383.0/256.0*t69+t160+t161-t162-t163+t164;

   double t225 = 729.0/2800.0*eta;

   double t226 = 9.0/40.0*t63;

   double t227 = 3537.0/1400.0*t65;

   double t228 = 9.0/4.0*t22;

   double t229 = 81.0/40.0*t18;

   double t231 = 81.0/80.0*t71;

   double t234 = t95+t225+t226-t227-t228+t102+t229-27.0/40.0*t67-t231+27.0/4.0*t69+t100-243.0/40.0*t24+t94-t97-t125+t126+t127-t128-729.0/2800.0;

   double t235 = 72.0/175.0*t2;

   double t237 = 2547.0/5600.0*t1;

   double t238 = 3.0/20.0*t10;

   double t240 = 243.0/640.0*t18;

   double t244 = 243.0/800.0*t30;

   double t245 = 27.0/200.0*t32;

   double t247 = 243.0/640.0*t38;

   double t248 = 243.0/1280.0*t40;

   double t250 = t235-t5-177.0/5600.0*xi-t237+t9-t238-261.0/640.0*t12-t16-t240-t21+45.0/128.0*t22+1269.0/1280.0*t24-159.0/5600.0*eta+t244+t245-81.0/256.0*t34-t37+t247+t248+1971.0/6400.0*t42;

   double t251 = 441.0/6400.0*t45;

   double t252 = 27.0/256.0*t49;

   double t254 = 1629.0/4480.0*t55;

   double t255 = 63.0/256.0*t57;

   double t257 = 729.0/1280.0*t61;

   double t258 = t63/320.0;

   double t262 = 1321.0/6400.0*t71;

   double t265 = (1.0+xi)*t28*t88;

   double t266 = 81.0/700.0*t265;

   double t267 = 243.0/2800.0-t251+t48+t252+t52+567.0/1280.0*t53+t254+t255-393.0/1280.0*t59-t257+t258+10111.0/22400.0*t65-6877.0/22400.0*t67-291.0/256.0*t69+t262+t73-t74-t81+t83-t266;

   double t269 = 432.0/175.0*t2;

   double t270 = 81.0/400.0*t4;

   double t272 = 7641.0/2800.0*t1;

   double t273 = 27.0/20.0*t8;

   double t274 = 9.0/10.0*t10;

   double t276 = 81.0/400.0*t15;

   double t278 = 81.0/160.0*t20;

   double t280 = 243.0/200.0*t24;

   double t281 = 9.0/100.0*eta;

   double t282 = 81.0/100.0*t32;

   double t284 = 729.0/160.0*t36;

   double t285 = -729.0/2800.0-t269+t270-351.0/2800.0*xi+t272-t273+t274+99.0/40.0*t12+t276+243.0/800.0*t18+t278+27.0/160.0*t22-t280-t281-t106-t282+243.0/160.0*t34+t284;

   double t288 = 81.0/16.0*t51;

   double t290 = 1269.0/560.0*t55;

   double t291 = 27.0/16.0*t57;

   double t293 = 2187.0/800.0*t61;

   double t294 = 99.0/400.0*t63;

   double t297 = 27.0/20.0*t69;

   double t298 = 81.0/200.0*t71;

   double t299 = 81.0/400.0*t14;

   double t300 = 243.0/200.0*t17;

   double t301 = 729.0/400.0*t82;

   double t302 = 243.0/350.0*t265;

   double t303 = -891.0/400.0*t38-1377.0/800.0*t42+t113-t288-243.0/80.0*t53-t290-t291+153.0/80.0*t59+t293-t294-153.0/800.0*t65+837.0/350.0*t67+t297-t298-t299+t300-t301+t302;

   double t305 = 216.0/35.0*t2;

   double t306 = 729.0/1120.0*xi;

   double t307 = 7641.0/1120.0*t1;

   double t308 = 9.0/4.0*t10;

   double t310 = 729.0/256.0*t15;

   double t312 = 243.0/128.0*t20;

   double t314 = 2187.0/1280.0*t24;

   double t315 = 9.0/40.0*eta;

   double t316 = 729.0/160.0*t30;

   double t317 = 81.0/40.0*t32;

   double t320 = 729.0/1120.0+t305+t306-t307-t308-3897.0/640.0*t12+t310-729.0/640.0*t18-t312-81.0/128.0*t22+t314+t315+t316+t317-729.0/256.0*t34-t142+3483.0/640.0*t38;

   double t321 = 2187.0/1280.0*t40;

   double t323 = 243.0/1280.0*t45;

   double t324 = 243.0/256.0*t49;

   double t325 = 2673.0/256.0*t51;

   double t327 = 25947.0/4480.0*t55;

   double t328 = 891.0/256.0*t57;

   double t330 = 6561.0/1280.0*t61;

   double t331 = 207.0/320.0*t63;

   double t333 = 28863.0/4480.0*t67;

   double t334 = 243.0/256.0*t69;

   double t335 = 243.0/1280.0*t71;

   double t336 = 243.0/140.0*t265;

   double t337 = -t321+5103.0/1280.0*t42-t323-t148-t324+t325+8019.0/1280.0*t53+t327+t328-1161.0/256.0*t59-t330+t331+351.0/640.0*t65-t333+t334+t335-t336;

   double t340 = t4/16.0;

   double t346 = 45.0/8.0*t15;

   double t348 = 27.0/32.0*t20;

   double t355 = 243.0/32.0*t36;

   double t358 = -288.0/35.0*t2-t340-25.0/28.0*xi+2547.0/280.0*t1+3.0/16.0*t8+3.0*t10+63.0/8.0*t12-t346+81.0/160.0*t18+t348+9.0/32.0*t22-27.0/8.0*t24-3.0/10.0*eta-243.0/40.0*t30-27.0/10.0*t32+81.0/32.0*t34+t355-567.0/80.0*t38+243.0/80.0*t40;

   double t362 = 135.0/16.0*t51;

   double t365 = 45.0/16.0*t57;

   double t367 = 729.0/160.0*t61;

   double t368 = 63.0/80.0*t63;

   double t372 = t71/20.0;

   double t373 = 9.0/16.0*t14;

   double t376 = -243.0/280.0-783.0/160.0*t42+9.0/80.0*t45+t188+27.0/16.0*t49-t362-81.0/16.0*t53-1089.0/140.0*t55-t365+87.0/16.0*t59+t367-t368-87.0/160.0*t65+607.0/70.0*t67-15.0/8.0*t69-t372+t373+27.0/80.0*t17+81.0/35.0*t265;

   double t379 = 81.0/128.0*t20;

   double t381 = 729.0/256.0*t36;

   double t383 = 729.0/1120.0+t305+t306-t307-t308-3663.0/640.0*t12+t310+t240+t379+t23+t314+t315+t316+t317-243.0/256.0*t34-t381+3321.0/640.0*t38;

   double t385 = 567.0/256.0*t51;

   double t387 = 189.0/256.0*t57;

   double t389 = 2187.0/1280.0*t61;

   double t390 = 171.0/320.0*t63;

   double t392 = -t321+4293.0/1280.0*t42-t323-t148-t324+t385+1701.0/1280.0*t53+t327+t387-927.0/256.0*t59-t389+t390+171.0/640.0*t65-t333+t334+t335-t336;

   double t394 = 729.0/2800.0*xi;

   double t396 = 3537.0/1400.0*t67;

   double t397 = 81.0/400.0*t71;

   double t398 = 9.0/4.0*t12;

   double t399 = 81.0/40.0*t38;

   double t400 = -t394-t281-t226-27.0/200.0*t65+t396-t397+t398+t118-t290+t113-t399-t112+t272+t274-t269-t282+t302-t106-729.0/2800.0;

   double t401 = 2547.0/5600.0*t4;

   double t403 = 3.0/20.0*t8;

   double t405 = 441.0/6400.0*t15;

   double t407 = 1629.0/4480.0*t20;

   double t413 = (1.0+eta)*t27*t79;

   double t414 = 81.0/700.0*t413;

   double t416 = 243.0/1280.0*t36;

   double t418 = t235-t401+327.0/5600.0*xi-t237-t403-t238-243.0/640.0*t12-t405+1053.0/3200.0*t18+t407-243.0/640.0*t22+1161.0/6400.0*t24+327.0/5600.0*eta-t414+t244+t245+81.0/1280.0*t34+t416+1053.0/3200.0*t38+t248;

   double t420 = 27.0/256.0*t51;

   double t430 = 72.0/175.0*t14;

   double t431 = 27.0/200.0*t17;

   double t432 = 243.0/800.0*t82;

   double t433 = 243.0/2800.0+1161.0/6400.0*t42-t251+t48+t252+t420+81.0/1280.0*t53+t254+9.0/256.0*t57-303.0/1280.0*t59+729.0/6400.0*t61+121.0/1600.0*t63-8893.0/22400.0*t65-8893.0/22400.0*t67-303.0/1280.0*t69+457.0/6400.0*t71+t430+t431+t432-t266;

   double t435 = 9.0/100.0*xi;

   double t437 = 7641.0/2800.0*t4;

   double t438 = 9.0/10.0*t8;

   double t439 = 432.0/175.0*t14;

   double t440 = 81.0/100.0*t17;

   double t441 = 243.0/350.0*t413;

   double t442 = -t435-t225-t226+t227+t228-t102-t229-27.0/200.0*t67-t397+t297+t276-t280+t437+t438-t439-t440+t441-t301-729.0/2800.0;

   double t443 = 7641.0/1120.0*t4;

   double t444 = 9.0/40.0*xi;

   double t445 = 9.0/4.0*t8;

   double t446 = 243.0/1280.0*t15;

   double t448 = 25947.0/4480.0*t20;

   double t451 = 243.0/140.0*t413;

   double t453 = 2187.0/1280.0*t36;

   double t454 = 2187.0/1280.0*t42;

   double t455 = 729.0/1120.0-t443+t444-t445+t13-t446+3321.0/640.0*t18+t448-3663.0/640.0*t22+4293.0/1280.0*t24+t140-t451+1701.0/1280.0*t34-t453+t247-t213+t454;

   double t456 = 243.0/256.0*t51;

   double t458 = 243.0/256.0*t59;

   double t459 = 28863.0/4480.0*t65;

   double t462 = 216.0/35.0*t14;

   double t463 = 81.0/40.0*t17;

   double t464 = 729.0/160.0*t82;

   double t465 = t146-t148+t215-t456-243.0/256.0*t53+t217+t387+t458-t389+t390-t459+171.0/640.0*t67-927.0/256.0*t69+t335+t462+t463+t464;

   double t483 = t167+2547.0/280.0*t4-3.0/10.0*xi-t170+3.0*t8+3.0/16.0*t10+9.0/32.0*t12+9.0/80.0*t15-567.0/80.0*t18-1089.0/140.0*t20+63.0/8.0*t22-783.0/160.0*t24-25.0/28.0*eta+81.0/35.0*t413+27.0/80.0*t32-81.0/16.0*t34+243.0/80.0*t36+81.0/160.0*t38+t184;

   double t494 = -243.0/280.0-27.0/8.0*t42-t187+t188-t189+27.0/16.0*t51+81.0/32.0*t53+t192-t365-15.0/8.0*t59+t367-t368+607.0/70.0*t65-87.0/160.0*t67+87.0/16.0*t69-t372-288.0/35.0*t14-27.0/10.0*t17-243.0/40.0*t82;

   double t502 = 729.0/1120.0-t443+t444-t445-81.0/128.0*t12-t446+3483.0/640.0*t18+t448-3897.0/640.0*t22+5103.0/1280.0*t24+t140-t451+8019.0/1280.0*t34-t453-729.0/640.0*t38-t144+t454;

   double t506 = t146-t148+t149-t456-729.0/256.0*t53-t152+t328+t458-t330+t331-t459+351.0/640.0*t67-1161.0/256.0*t69+t335+t462+t463+t464;

   double t514 = -729.0/2800.0-t93+t437-t435+t96+t438-t98+27.0/160.0*t12+t276-891.0/400.0*t18-t102+99.0/40.0*t22-1377.0/800.0*t24-351.0/2800.0*eta+t441-t106+t107-243.0/80.0*t34;

   double t520 = 243.0/800.0*t38+t111-t112+t113-t114+243.0/160.0*t53+t116-t291+t118+t293-t294+837.0/350.0*t65-153.0/800.0*t67+153.0/80.0*t69-t298-t439-t440-t301;

   double t528 = t3-t401-159.0/5600.0*xi-t7-t403+t11+45.0/128.0*t12-t405+t240+t407-261.0/640.0*t22+1971.0/6400.0*t24-177.0/5600.0*eta-t414+t31-t33+567.0/1280.0*t34+t416-t247-t41;

   double t535 = 243.0/2800.0+1269.0/1280.0*t42-t46+t48+t50+t420-81.0/256.0*t53-t56+t255-291.0/256.0*t59-t257+t258-6877.0/22400.0*t65+10111.0/22400.0*t67-393.0/1280.0*t69+t262+t430+t431+t432-t90;

   double t537 = 9.0/20.0*eta;

   double t540 = 81.0/80.0*t45;

   double t542 = 9909.0/2800.0*t1;

   double t543 = 9.0/2.0*t10;

   double t544 = 459.0/140.0*t2;

   double t545 = 81.0/20.0*t32;

   double t546 = 243.0/350.0*t89;

   double t547 = 729.0/80.0*t30;

   double t548 = t394+t537+t226-27.0/40.0*t65-t396-t231-t398+27.0/4.0*t59+t290+t540+t399-243.0/40.0*t42+t542-t543-t544+t545+t546-t547-729.0/2800.0;

   double t549 = 459.0/56.0*t2;

   double t550 = 9909.0/1120.0*t1;

   double t551 = 45.0/4.0*t10;

   double t555 = 2187.0/256.0*t24;

   double t556 = 9.0/8.0*eta;

   double t557 = 729.0/32.0*t30;

   double t558 = 81.0/8.0*t32;

   double t560 = 729.0/1120.0+t549-t306-t550+t551+657.0/128.0*t12+t310-243.0/128.0*t18+t379+189.0/128.0*t22-t555-t556+t557-t558-t211-t381-567.0/128.0*t38;

   double t562 = 567.0/256.0*t45;

   double t563 = 1701.0/256.0*t49;

   double t564 = 4509.0/896.0*t55;

   double t567 = 25461.0/4480.0*t67;

   double t568 = 1701.0/256.0*t69;

   double t569 = 243.0/140.0*t89;

   double t570 = t144+3969.0/256.0*t42-t562-t148-t563+t385-t216-t564+t218-4383.0/256.0*t59+t219-t220+207.0/128.0*t65+t567+t568+t160-t569;

   double t576 = 15.0*t10;

   double t581 = 3.0/2.0*eta;

   double t583 = 27.0/2.0*t32;

   double t587 = -153.0/14.0*t2-t340+26.0/35.0*xi+3303.0/280.0*t1+21.0/16.0*t8-t576-45.0/8.0*t12-t346-81.0/32.0*t18+t348+63.0/32.0*t22+135.0/8.0*t24+t581-243.0/8.0*t30+t583+567.0/32.0*t34+t355+81.0/16.0*t38-243.0/16.0*t40;

   double t599 = -243.0/280.0-675.0/32.0*t42+45.0/16.0*t45+t188+189.0/16.0*t49-t362+405.0/16.0*t53+45.0/7.0*t55-t193+375.0/16.0*t59-t195+t196-75.0/32.0*t65-251.0/35.0*t67-105.0/8.0*t69-t200+t373-27.0/16.0*t17+81.0/35.0*t89;

   double t602 = 729.0/128.0*t18;

   double t603 = 567.0/128.0*t22;

   double t604 = 5103.0/256.0*t34;

   double t606 = 729.0/1120.0+t549-t306-t550+t551+423.0/128.0*t12+t310+t602-t312-t603-t555-t556+t557-t558-t604-t142-405.0/128.0*t38;

   double t611 = t144+4131.0/256.0*t42-t562-t148-t563+t325-8019.0/256.0*t53-t564+t153-4617.0/256.0*t59+t155-t156+243.0/128.0*t65+t567+t568+t160-t569;

   double t618 = -729.0/2800.0-t544+t270+1107.0/2800.0*xi+t542-t273-t543-9.0/8.0*t12+t276-243.0/160.0*t18+t278+189.0/160.0*t22-t280+t537-t547+t545+1701.0/160.0*t34+t284;

   double t625 = 81.0/80.0*t38-1053.0/160.0*t42+t540-t288+243.0/16.0*t53+t290-t117+117.0/16.0*t59-t119+t120-117.0/160.0*t65-1863.0/700.0*t67+t297-t124-t299+t300-t301+t546;

   double t627 = 243.0/200.0*t2;

   double t628 = 243.0/200.0*t4;

   double t629 = 81.0/100.0*xi;

   double t630 = 243.0/200.0*t1;

   double t631 = 81.0/10.0*t8;

   double t632 = 81.0/10.0*t10;

   double t634 = 243.0/200.0*t15;

   double t637 = 729.0/100.0*t24;

   double t638 = 81.0/100.0*eta;

   double t639 = 2187.0/200.0*t30;

   double t640 = 729.0/100.0*t32;

   double t642 = t627-t628-t629-t630+t631+t632-81.0/10.0*t12-t634+729.0/100.0*t18-81.0/10.0*t22+t637-t638+t639-t640-729.0/10.0*t34;

   double t644 = 729.0/100.0*t42;

   double t645 = 243.0/200.0*t45;

   double t648 = 81.0/10.0*t59;

   double t651 = 81.0/100.0*t65;

   double t652 = 81.0/100.0*t67;

   double t653 = 81.0/10.0*t69;

   double t655 = 243.0/200.0*t14;

   double t656 = 729.0/100.0*t17;

   double t657 = 2187.0/200.0*t82;

   double t658 = 729.0/100.0*t38+t644-t645-729.0/10.0*t53+81.0*t57-t648+6561.0/100.0*t61+81.0/100.0*t63+t651+t652-t653+243.0/100.0*t71+t655-t656+t657;

   double t660 = 243.0/80.0*t2;

   double t661 = 243.0/80.0*t1;

   double t662 = 81.0/4.0*t10;

   double t666 = 6561.0/160.0*t24;

   double t667 = 81.0/40.0*eta;

   double t668 = 2187.0/80.0*t30;

   double t669 = 729.0/40.0*t32;

   double t672 = 6561.0/160.0*t40;

   double t673 = -t660+t661-t662+2673.0/160.0*t12-2187.0/80.0*t18+243.0/8.0*t22+t666+t667-t668+t669+2187.0/16.0*t34-2187.0/160.0*t38-t672;

   double t674 = 729.0/40.0*t42;

   double t675 = 243.0/80.0*t45;

   double t676 = 729.0/16.0*t49;

   double t678 = 729.0/160.0*t55;

   double t679 = 2673.0/16.0*t57;

   double t680 = 81.0/4.0*t59;

   double t681 = 19683.0/160.0*t61;

   double t682 = 243.0/80.0*t63;

   double t683 = 81.0/40.0*t65;

   double t684 = 729.0/160.0*t67;

   double t685 = 729.0/16.0*t69;

   double t686 = 243.0/80.0*t71;

   double t687 = -t674+t675+t676+24057.0/160.0*t53-t678-t679+t680-t681-t682-t683+t684-t685-t686;

   double t694 = 27.0/2.0*t12;

   double t695 = 243.0/20.0*t18;

   double t696 = 27.0/2.0*t22;

   double t701 = 243.0/2.0*t34;

   double t702 = 243.0/20.0*t38;

   double t703 = 81.0/20.0*t2+9.0/10.0*xi-81.0/20.0*t1-9.0*t8+27.0*t10-t694+t695-t696-81.0*t24-27.0/10.0*eta+729.0/20.0*t30-243.0/10.0*t32-t701+t702;

   double t708 = 243.0/2.0*t53;

   double t710 = 135.0*t57;

   double t712 = 2187.0/20.0*t61;

   double t713 = 27.0/20.0*t63;

   double t717 = 81.0/20.0*t71;

   double t719 = 729.0/10.0*t40+243.0/10.0*t42-81.0/20.0*t45-81.0*t49-t708+81.0/10.0*t55+t710-27.0*t59+t712+t713+27.0/10.0*t65-9.0*t67+90.0*t69+t717+81.0/10.0*t17;

   double t726 = -t660+t661-t662+567.0/160.0*t12+729.0/80.0*t18-81.0/8.0*t22+t666+t667-t668+t669+729.0/16.0*t34-729.0/160.0*t38-t672;

   double t728 = 567.0/16.0*t57;

   double t729 = 6561.0/160.0*t61;

   double t730 = 81.0/80.0*t63;

   double t731 = -t674+t675+t676+5103.0/160.0*t53-t678-t728+t680-t729+t730-t683+t684-t685-t686;

   double t733 = 243.0/80.0*t4;

   double t734 = 81.0/40.0*xi;

   double t735 = 81.0/4.0*t8;

   double t737 = 243.0/80.0*t15;

   double t739 = 729.0/160.0*t20;

   double t741 = 729.0/40.0*t24;

   double t743 = 6561.0/160.0*t36;

   double t745 = 6561.0/160.0*t42;

   double t746 = t733+t734-t735+243.0/8.0*t12+t737-2187.0/160.0*t18-t739+2673.0/160.0*t22-t741+24057.0/160.0*t34-t743-2187.0/80.0*t38+t745;

   double t747 = 729.0/16.0*t51;

   double t749 = 729.0/16.0*t59;

   double t750 = 729.0/160.0*t65;

   double t751 = 81.0/40.0*t67;

   double t752 = 81.0/4.0*t69;

   double t753 = 243.0/80.0*t14;

   double t754 = 729.0/40.0*t17;

   double t755 = 2187.0/80.0*t82;

   double t756 = t747+2187.0/16.0*t53-t679-t749-t681-t682+t750-t751+t752-t686-t753+t754-t755;

   double t759 = 6561.0/256.0*t15;

   double t761 = 2187.0/128.0*t20;

   double t763 = 19683.0/256.0*t24;

   double t765 = 19683.0/256.0*t36;

   double t767 = 19683.0/256.0*t40;

   double t768 = 19683.0/256.0*t42;

   double t769 = 6561.0/256.0*t45;

   double t770 = -8019.0/128.0*t12-t759+6561.0/128.0*t18+t761-8019.0/128.0*t22-t763-72171.0/256.0*t34+t765+6561.0/128.0*t38+t767-t768-t769;

   double t771 = 6561.0/256.0*t47;

   double t772 = 24057.0/256.0*t49;

   double t773 = 24057.0/256.0*t51;

   double t775 = 2187.0/128.0*t55;

   double t777 = 24057.0/256.0*t59;

   double t780 = 2187.0/128.0*t65;

   double t781 = 2187.0/128.0*t67;

   double t782 = 24057.0/256.0*t69;

   double t783 = 6561.0/256.0*t71;

   double t784 = t771-t772-t773-72171.0/256.0*t53+t775+88209.0/256.0*t57+t777+59049.0/256.0*t61+729.0/64.0*t63-t780-t781+t782+t783;

   double t786 = 81.0/16.0*t4;

   double t790 = 405.0/8.0*t15;

   double t792 = 243.0/32.0*t20;

   double t796 = 2187.0/32.0*t36;

   double t799 = 2187.0/32.0*t42;

   double t800 = 729.0/16.0*t45;

   double t801 = t786-27.0/8.0*xi+297.0/16.0*t8+405.0/8.0*t12+t790-729.0/32.0*t18-t792+891.0/32.0*t22+1215.0/8.0*t24+8019.0/32.0*t34-t796-729.0/16.0*t38-2187.0/16.0*t40+t799+t800;

   double t802 = 729.0/16.0*t47;

   double t804 = 1215.0/16.0*t51;

   double t807 = 4455.0/16.0*t57;

   double t808 = 1215.0/16.0*t59;

   double t809 = 6561.0/32.0*t61;

   double t810 = 81.0/16.0*t63;

   double t811 = 243.0/32.0*t65;

   double t814 = 405.0/8.0*t71;

   double t815 = 81.0/16.0*t14;

   double t817 = -t802+2673.0/16.0*t49+t804+3645.0/16.0*t53-243.0/8.0*t55-t807-t808-t809-t810+t811+135.0/4.0*t67-1485.0/8.0*t69-t814-t815-243.0/16.0*t17;

   double t820 = 2187.0/128.0*t18;

   double t821 = 729.0/128.0*t20;

   double t824 = 6561.0/256.0*t36;

   double t825 = 2187.0/128.0*t38;

   double t826 = 6561.0/256.0*t42;

   double t827 = -1701.0/128.0*t12-t759-t820-t821+2673.0/128.0*t22-t763-24057.0/256.0*t34+t824+t825+t767-t826-t769;

   double t828 = 5103.0/256.0*t51;

   double t830 = 18711.0/256.0*t57;

   double t831 = 5103.0/256.0*t59;

   double t832 = 19683.0/256.0*t61;

   double t833 = 243.0/64.0*t63;

   double t834 = 729.0/128.0*t65;

   double t835 = t771-t772-t828-15309.0/256.0*t53+t775+t830+t831+t832-t833+t834-t781+t782+t783;

   double t847 = -81.0/20.0*t4-27.0/10.0*xi+27.0*t8-9.0*t10-t694-81.0/20.0*t15+t695+81.0/10.0*t20-t696+243.0/10.0*t24+9.0/10.0*eta+81.0/10.0*t32-t701+729.0/10.0*t36;

   double t857 = t702-81.0*t42-81.0*t51-t708+t710+90.0*t59+t712+t713-9.0*t65+27.0/10.0*t67-27.0*t69+t717+81.0/20.0*t14-243.0/10.0*t17+729.0/20.0*t82;

   double t859 = 81.0/16.0*t2;

   double t860 = 81.0/16.0*t1;

   double t863 = 729.0/16.0*t15;

   double t867 = 2187.0/32.0*t24;

   double t873 = 2187.0/32.0*t40;

   double t874 = -t859+t860+297.0/16.0*t10+891.0/32.0*t12+t863-729.0/16.0*t18-243.0/8.0*t20+405.0/8.0*t22+t867-27.0/8.0*eta-243.0/16.0*t32+3645.0/16.0*t34-2187.0/16.0*t36-729.0/32.0*t38-t873;

   double t876 = 405.0/8.0*t45;

   double t877 = 1215.0/16.0*t49;

   double t880 = 243.0/32.0*t55;

   double t883 = 243.0/32.0*t67;

   double t884 = 1215.0/16.0*t69;

   double t885 = 1215.0/8.0*t42+t876-t802+t877+2673.0/16.0*t51+8019.0/32.0*t53-t880-t807-1485.0/8.0*t59-t809-t810+135.0/4.0*t65+t883-t884-t814;

   double t899 = 1.0+9.0*t2-10.0*t4+t169-10.0*t1-t171-t576-45.0/2.0*t12-90.0*t15+81.0/4.0*t18+27.0/2.0*t20-45.0/2.0*t22-135.0*t24+t581+t583-405.0/2.0*t34+243.0/2.0*t36+81.0/4.0*t38;

   double t917 = 243.0/2.0*t40-135.0*t42-90.0*t45+81.0*t47-135.0*t49-135.0*t51-405.0/2.0*t53+27.0/2.0*t55+225.0*t57+150.0*t59+729.0/4.0*t61+9.0/4.0*t63-15.0*t65-15.0*t67+150.0*t69+100.0*t71+9.0*t14+t202;

   double t928 = -t859+t860+63.0/16.0*t10+189.0/32.0*t12+t863+243.0/16.0*t18+81.0/8.0*t20-135.0/8.0*t22+t867+t556-81.0/16.0*t32+1215.0/16.0*t34-729.0/16.0*t36-243.0/32.0*t38-t873;

   double t932 = 945.0/16.0*t57;

   double t934 = 2187.0/32.0*t61;

   double t935 = 27.0/16.0*t63;

   double t937 = 405.0/8.0*t42+t876-t802+t877+567.0/16.0*t51+1701.0/32.0*t53-t880-t932-315.0/8.0*t59-t934+t935-45.0/4.0*t65+t883-t884-t814;

   double t944 = t733+t734-t735-81.0/8.0*t12+t737-729.0/160.0*t18-t739+567.0/160.0*t22-t741+5103.0/160.0*t34-t743+729.0/80.0*t38+t745;

   double t946 = t747+729.0/16.0*t53-t728-t749-t729+t730+t750-t751+t752-t686-t753+t754-t755;

   double t950 = 6561.0/256.0*t24;

   double t952 = 6561.0/256.0*t40;

   double t953 = 2673.0/128.0*t12-t759+t820+t761-1701.0/128.0*t22-t950-15309.0/256.0*t34+t765-t825+t952-t768-t769;

   double t954 = 5103.0/256.0*t49;

   double t956 = 729.0/128.0*t55;

   double t957 = 729.0/128.0*t67;

   double t958 = 5103.0/256.0*t69;

   double t959 = t771-t954-t773-24057.0/256.0*t53-t956+t830+t777+t832-t833-t780+t957+t958+t783;

   double t969 = t786+t132+63.0/16.0*t8-135.0/8.0*t12+t790-243.0/32.0*t18-t792+189.0/32.0*t22+405.0/8.0*t24+1701.0/32.0*t34-t796+243.0/16.0*t38-729.0/16.0*t40+t799+t800;

   double t976 = -t802+567.0/16.0*t49+t804+1215.0/16.0*t53+81.0/8.0*t55-t932-t808-t934+t935+t811-45.0/4.0*t67-315.0/8.0*t69-t814-t815-81.0/16.0*t17;

   double t978 = t134-t759-t602-t821+t603-t950-t604+t824-t143+t952-t826-t769;

   double t981 = t771-t954-t828-t151-t956+3969.0/256.0*t57+t831+t155+81.0/64.0*t63+t834+t957+t958+t783;

   double t983 = -t629+t652+t124-t653-t634+t637-t628+t631+t655-t656+t657;

   double t984 = -t638+t651+t124-t648-t645+t644-t630+t632+t627-t640+t639;


   values[0] = t44+t91;

   values[1] = t110+t129;

   values[2] = t147+t165;

   values[3] = t186+t205;

   values[4] = t214+t223;

   values[5] = t234;

   values[6] = t250+t267;

   values[7] = t285+t303;

   values[8] = t320+t337;

   values[9] = t358+t376;

   values[10] = t383+t392;

   values[11] = t400;

   values[12] = t418+t433;

   values[13] = t442;

   values[14] = t455+t465;

   values[15] = t483+t494;

   values[16] = t502+t506;

   values[17] = t514+t520;

   values[18] = t528+t535;

   values[19] = t548;

   values[20] = t560+t570;

   values[21] = t587+t599;

   values[22] = t606+t611;

   values[23] = t618+t625;

   values[24] = t642+t658;

   values[25] = t673+t687;

   values[26] = t703+t719;

   values[27] = t726+t731;

   values[28] = t746+t756;

   values[29] = t770+t784;

   values[30] = t801+t817;

   values[31] = t827+t835;

   values[32] = t847+t857;

   values[33] = t874+t885;

   values[34] = t899+t917;

   values[35] = t928+t937;

   values[36] = t944+t946;

   values[37] = t953+t959;

   values[38] = t969+t976;

   values[39] = t978+t981;

   values[40] = t983;

   values[41] = t984;

 }


 // values of the derivatives in xi direction

 static void C_Q_UL6S_2D_DeriveXi(double xi, double eta, double *values)

 {

   double t1 = eta*eta;

   double t2 = t1*eta;

   double t3 = 27.0/128.0*t2;

   double t4 = t1*t1;

   double t5 = t4*eta;

   double t7 = 2007.0/4480.0*t4;

   double t8 = 1.0-eta;

   double t9 = xi*xi;

   double t10 = 1.0-t9;

   double t12 = t9*xi;

   double t15 = 35.0/2.0*t12-15.0/2.0*xi;

   double t16 = t8*t10*t15;

   double t17 = 81.0/700.0*t16;

   double t18 = 1.0-t1;

   double t20 = xi*t18*t4;

   double t21 = 243.0/80.0*t20;

   double t22 = t9*t9;

   double t23 = t22*eta;

   double t25 = t12*eta;

   double t26 = 2007.0/1120.0*t25;

   double t27 = xi*t2;

   double t29 = t22*t1;

   double t33 = t9*t5;

   double t35 = xi*t5;

   double t37 = t22*t5;

   double t39 = t12*t5;

   double t40 = 243.0/64.0*t39;

   double t42 = t22*xi*t18;

   double t43 = 243.0/80.0*t42;

   double t47 = t18*(3.0/8.0+35.0/8.0*t4-15.0/4.0*t1);

   double t48 = 81.0/700.0*t47;

   double t49 = t4*xi;

   double t50 = 261.0/640.0*t49;

   double t51 = t3-81.0/640.0*t5-t7-t17-t21-81.0/128.0*t23-t26-309.0/128.0*t27+1431.0/256.0*t29+12127.0/22400.0*t1-19.0/320.0*eta-1701.0/256.0*t33+1431.0/640.0*t35+3645.0/256.0*t37-t40-t43+t48-t50;

   double t52 = t22*t2;

   double t54 = t12*t4;

   double t55 = 81.0/64.0*t54;

   double t56 = t22*t4;

   double t57 = 1215.0/256.0*t56;

   double t58 = t9*t4;

   double t59 = 567.0/256.0*t58;

   double t60 = xi*eta;

   double t62 = t9*t2;

   double t64 = t12*t2;

   double t65 = 189.0/64.0*t64;

   double t66 = t9*eta;

   double t67 = 81.0/128.0*t66;

   double t68 = 153.0/70.0*t12;

   double t70 = t10*t18*t12;

   double t71 = 243.0/40.0*t70;

   double t72 = xi*t1;

   double t74 = t9*t1;

   double t76 = t12*t1;

   double t77 = 261.0/320.0*t76;

   double t78 = 27.0/8.0*t22;

   double t79 = 9.0/4.0*t9;

   double t80 = 3303.0/2800.0*xi;

   double t84 = 3.0/8.0+35.0/8.0*t22-15.0/4.0*t9;

   double t85 = t8*xi*t84;

   double t86 = 81.0/350.0*t85;

   double t87 = -663.0/5600.0-2835.0/256.0*t52+t55-t57+t59+12127.0/11200.0*t60+1323.0/256.0*t62+t65+t67+t68+t71+437.0/640.0*t72-927.0/256.0*t74-t77-t78+t79-t80+t86;

   double t91 = 81.0/160.0*t4;

   double t92 = 243.0/350.0*t16;

   double t93 = 729.0/200.0*t20;

   double t95 = 1269.0/140.0*t25;

   double t96 = 27.0/10.0*t27;

   double t99 = 9.0/80.0*eta;

   double t101 = 243.0/100.0*t35;

   double t102 = 2187.0/32.0*t37;

   double t103 = 729.0/40.0*t42;

   double t104 = 81.0/200.0*t49;

   double t105 = 189.0/160.0*t2-243.0/160.0*t5+t91+t92+t93+81.0/16.0*t23+t95+t96-1053.0/32.0*t29-117.0/160.0*t1+t99+729.0/16.0*t33-t101-t102+t103+t104;

   double t107 = 729.0/32.0*t56;

   double t108 = 243.0/16.0*t58;

   double t110 = 567.0/16.0*t62;

   double t112 = 459.0/35.0*t12;

   double t113 = 729.0/20.0*t70;

   double t114 = 243.0/100.0*t72;

   double t116 = 81.0/20.0*t76;

   double t117 = 81.0/4.0*t22;

   double t118 = 27.0/2.0*t9;

   double t119 = 9909.0/1400.0*xi;

   double t120 = 243.0/175.0*t85;

   double t121 = 9.0/20.0+1701.0/32.0*t52+t107-t108-1863.0/350.0*t60-t110-27.0/8.0*t66-t112-t113-t114+351.0/16.0*t74+t116+t117-t118+t119-t120;

   double t123 = 567.0/128.0*t2;

   double t124 = 729.0/128.0*t5;

   double t125 = 243.0/128.0*t4;

   double t126 = 243.0/140.0*t16;

   double t128 = 4509.0/224.0*t25;

   double t129 = 1701.0/128.0*t27;

   double t132 = 9.0/64.0*eta;

   double t134 = 2187.0/128.0*t35;

   double t135 = 32805.0/256.0*t37;

   double t136 = 2187.0/64.0*t39;

   double t137 = 729.0/16.0*t42;

   double t138 = 729.0/128.0*t49;

   double t139 = 25515.0/256.0*t52;

   double t140 = -t123+t124-t125-t126-2025.0/128.0*t23-t128+t129+20655.0/256.0*t29+243.0/128.0*t1-t132-24057.0/256.0*t33-t134+t135+t136-t137+t138-t139;

   double t141 = 729.0/64.0*t54;

   double t142 = 10935.0/256.0*t56;

   double t143 = 8019.0/256.0*t58;

   double t144 = 25461.0/2240.0*t60;

   double t145 = 18711.0/256.0*t62;

   double t146 = 1701.0/64.0*t64;

   double t148 = 459.0/14.0*t12;

   double t149 = 729.0/8.0*t70;

   double t150 = 567.0/128.0*t72;

   double t152 = 567.0/64.0*t76;

   double t153 = 405.0/8.0*t22;

   double t154 = 135.0/4.0*t9;

   double t155 = 9909.0/560.0*xi;

   double t156 = 243.0/70.0*t85;

   double t157 = -9.0/8.0-t141-t142+t143+t144+t145-t146+1269.0/128.0*t66+t148+t149+t150-13851.0/256.0*t74-t152-t153+t154-t155+t156;

   double t161 = 27.0/32.0*t4;

   double t168 = 9.0/16.0*eta;

   double t171 = 3645.0/32.0*t37;

   double t174 = 45.0/4.0*t49;

   double t175 = 3.0/2.0+63.0/32.0*t2-81.0/32.0*t5+t161+81.0/35.0*t16+405.0/16.0*t23+180.0/7.0*t25-105.0/4.0*t27-3375.0/32.0*t29-75.0/32.0*t1+t168+1215.0/16.0*t33+135.0/4.0*t35-t171-243.0/4.0*t39+243.0/4.0*t42-t174;

   double t177 = 81.0/4.0*t54;

   double t178 = 1215.0/32.0*t56;

   double t179 = 405.0/16.0*t58;

   double t181 = 945.0/16.0*t62;

   double t186 = 11.0/2.0*t72;

   double t189 = 135.0/2.0*t22;

   double t190 = 45.0*t9;

   double t193 = 2835.0/32.0*t52+t177+t178-t179-502.0/35.0*t60-t181+189.0/4.0*t64-135.0/8.0*t66-306.0/7.0*t12-243.0/2.0*t70-t186+1125.0/16.0*t74+45.0/4.0*t76+t189-t190+3303.0/140.0*xi-162.0/35.0*t85;

   double t197 = 81.0/128.0*t4;

   double t201 = 45.0/64.0*eta;

   double t202 = 5103.0/256.0*t33;

   double t203 = 10935.0/256.0*t37;

   double t204 = 8505.0/256.0*t52;

   double t205 = 189.0/128.0*t2-243.0/128.0*t5+t197-t126-2835.0/128.0*t23-t128+t129+19845.0/256.0*t29+207.0/128.0*t1-t201-t202-t134+t203+t136-t137+t138-t204;

   double t206 = 3645.0/256.0*t56;

   double t207 = 1701.0/256.0*t58;

   double t208 = 3969.0/256.0*t62;

   double t211 = -9.0/8.0-t141-t206+t207+t144+t208-t146+1971.0/128.0*t66+t148+t149+t150-13149.0/256.0*t74-t152-t153+t154-t155+t156;

   double t213 = 9.0/40.0*eta;

   double t214 = 3537.0/700.0*t60;

   double t215 = 27.0/4.0*t66;

   double t216 = 81.0/8.0*t23;

   double t218 = 81.0/40.0*t72;

   double t221 = 9.0/20.0+t213-t214-t215+t95+t216-27.0/40.0*t1-t218+81.0/4.0*t74+t116-243.0/8.0*t29+t119-t118-t112+t117-t120+t92+t103-t113;

   double t223 = 243.0/640.0*t5;

   double t224 = 1629.0/4480.0*t4;

   double t225 = 243.0/400.0*t20;

   double t226 = 243.0/128.0*t23;

   double t230 = eta/320.0;

   double t233 = 729.0/256.0*t37;

   double t234 = 441.0/3200.0*t49;

   double t235 = -261.0/640.0*t2+t223+t224-t17-t225-t226-t26-393.0/640.0*t27+1269.0/256.0*t29-6877.0/22400.0*t1+t230-243.0/256.0*t33+1971.0/3200.0*t35-t233-t40-t43-t48-t234;

   double t237 = 243.0/256.0*t56;

   double t238 = 81.0/256.0*t58;

   double t240 = 189.0/256.0*t62;

   double t242 = 1321.0/3200.0*t72;

   double t244 = -177.0/5600.0+567.0/256.0*t52+t55+t237+t238+10111.0/11200.0*t60+t240+t65+135.0/128.0*t66+t68+t71+t242-873.0/256.0*t74-t77-t78+t79-t80+t86;

   double t248 = 1269.0/560.0*t4;

   double t250 = 81.0/40.0*t25;

   double t252 = 243.0/40.0*t29;

   double t254 = 99.0/400.0*eta;

   double t257 = 2187.0/160.0*t37;

   double t258 = 729.0/40.0*t39;

   double t259 = 729.0/200.0*t42;

   double t260 = 99.0/40.0*t2-891.0/400.0*t5-t248+t93+243.0/160.0*t23+t250+153.0/40.0*t27-t252+837.0/350.0*t1-t254+729.0/160.0*t33-1377.0/400.0*t35+t257+t258+t259;

   double t261 = 243.0/350.0*t47;

   double t264 = 81.0/16.0*t62;

   double t265 = 81.0/4.0*t64;

   double t267 = 81.0/100.0*t12;

   double t268 = 729.0/100.0*t70;

   double t269 = 81.0/100.0*t72;

   double t270 = 81.0/20.0*t74;

   double t271 = 81.0/100.0*t76;

   double t272 = 243.0/40.0*t22;

   double t273 = 81.0/20.0*t9;

   double t274 = 81.0/200.0*xi;

   double t275 = -351.0/2800.0+t261+t104-243.0/16.0*t52-153.0/400.0*t60-t264-t265+81.0/160.0*t66-t267-t268-t269+t270+t271+t272-t273+t274;

   double t279 = 25947.0/4480.0*t4;

   double t280 = 729.0/80.0*t20;

   double t282 = 243.0/32.0*t25;

   double t284 = 2187.0/256.0*t29;

   double t285 = 28863.0/4480.0*t1;

   double t286 = 207.0/320.0*eta;

   double t289 = 6561.0/256.0*t37;

   double t290 = -3897.0/640.0*t2+3483.0/640.0*t5+t279-t280-729.0/128.0*t23-t282-1161.0/128.0*t27+t284-t285+t286-2187.0/256.0*t33+5103.0/640.0*t35-t289-t136;

   double t291 = 243.0/140.0*t47;

   double t292 = 243.0/640.0*t49;

   double t294 = 2187.0/256.0*t56;

   double t295 = 729.0/256.0*t58;

   double t297 = 2673.0/256.0*t62;

   double t298 = 2673.0/64.0*t64;

   double t300 = 243.0/640.0*t72;

   double t301 = 729.0/256.0*t74;

   double t302 = 729.0/64.0*t76;

   double t303 = 729.0/1120.0-t291-t292+8019.0/256.0*t52-t141-t294-t295+351.0/320.0*t60+t297+t298-243.0/128.0*t66+t300+t301+t302;

   double t310 = 27.0/8.0*t25;

   double t314 = 63.0/80.0*eta;

   double t317 = 729.0/32.0*t37;

   double t318 = 243.0/8.0*t39;

   double t319 = 81.0/35.0*t47;

   double t321 = 63.0/8.0*t2-567.0/80.0*t5-1089.0/140.0*t4+243.0/20.0*t20+81.0/32.0*t23+t310+87.0/8.0*t27-135.0/8.0*t29+607.0/70.0*t1-t314+243.0/32.0*t33-783.0/80.0*t35+t317+t318+t319+9.0/40.0*t49;

   double t326 = 135.0/16.0*t62;

   double t327 = 135.0/4.0*t64;

   double t329 = 9.0/4.0*t12;

   double t330 = t72/10.0;

   double t332 = 45.0/2.0*t76;

   double t335 = xi/8.0;

   double t336 = -25.0/28.0-405.0/16.0*t52+t177+243.0/16.0*t56+81.0/16.0*t58-87.0/80.0*t60-t326-t327+27.0/32.0*t66+t329-t330-45.0/8.0*t74-t332+27.0/16.0*t22+9.0/16.0*t9-t335;

   double t340 = 81.0/32.0*t25;

   double t342 = 171.0/320.0*eta;

   double t345 = 2187.0/256.0*t37;

   double t346 = 729.0/64.0*t39;

   double t347 = -3663.0/640.0*t2+3321.0/640.0*t5+t279-t280+t226+t340-927.0/128.0*t27+t284-t285+t342-729.0/256.0*t33+4293.0/640.0*t35-t345-t346;

   double t350 = 567.0/256.0*t62;

   double t351 = 567.0/64.0*t64;

   double t352 = 729.0/1120.0-t291-t292+1701.0/256.0*t52-t141-t294-t295+171.0/320.0*t60+t350+t351+t67+t300+t301+t302;

   double t355 = 3537.0/1400.0*t1;

   double t356 = 81.0/200.0*t72;

   double t357 = 9.0/4.0*t2;

   double t358 = 81.0/40.0*t5;

   double t359 = -729.0/2800.0-t213-27.0/100.0*t60+t355-t356+t357+t96-t248+t104-t358-t101+t261+t93;

   double t363 = 1629.0/1120.0*t25;

   double t367 = 1.0+eta;

   double t369 = t367*xi*t84;

   double t370 = 81.0/350.0*t369;

   double t372 = t367*t10*t15;

   double t373 = 81.0/700.0*t372;

   double t378 = 243.0/320.0*t39;

   double t379 = 243.0/400.0*t42;

   double t380 = 327.0/5600.0-243.0/640.0*t2+1053.0/3200.0*t5+t224-t225+1053.0/640.0*t23+t363-303.0/640.0*t27+1161.0/1280.0*t29-8893.0/22400.0*t1+t370-t373+121.0/1600.0*eta+243.0/1280.0*t33+1161.0/3200.0*t35+729.0/1280.0*t37+t378-t379;

   double t384 = 27.0/64.0*t64;

   double t386 = 288.0/175.0*t12;

   double t387 = 243.0/200.0*t70;

   double t390 = 441.0/1600.0*t76;

   double t391 = 27.0/40.0*t22;

   double t392 = 9.0/20.0*t9;

   double t393 = 2547.0/2800.0*xi;

   double t394 = -t48-t234+81.0/256.0*t52+t55+t237+t238-8893.0/11200.0*t60+27.0/256.0*t62+t384-729.0/640.0*t66+t386+t387+457.0/3200.0*t72-909.0/1280.0*t74-t390+t391-t392-t393;

   double t397 = 7641.0/1400.0*xi;

   double t398 = 27.0/10.0*t9;

   double t399 = 1728.0/175.0*t12;

   double t400 = 81.0/20.0*t22;

   double t401 = 243.0/175.0*t369;

   double t402 = 243.0/350.0*t372;

   double t403 = -9.0/100.0-t213+t214+t215-t95-t216-27.0/200.0*t1-t356+t270+t271-t252+t397+t398-t399-t400-t401+t402+t259-t268;

   double t405 = 25947.0/1120.0*t25;

   double t406 = 243.0/128.0*t27;

   double t409 = 243.0/70.0*t369;

   double t410 = 243.0/140.0*t372;

   double t412 = 2187.0/640.0*t35;

   double t413 = 2187.0/320.0*t39;

   double t414 = 729.0/80.0*t42;

   double t415 = t3+t223+t197+3321.0/128.0*t23+t405+t406+4293.0/256.0*t29+171.0/640.0*t1+t409-t410+t342+5103.0/1280.0*t33+t412-t345-t413-t414+t138;

   double t417 = 28863.0/2240.0*t60;

   double t418 = 243.0/64.0*t64;

   double t420 = 864.0/35.0*t12;

   double t421 = 729.0/40.0*t70;

   double t423 = 243.0/320.0*t76;

   double t424 = 81.0/8.0*t22;

   double t425 = 27.0/4.0*t9;

   double t426 = 7641.0/560.0*xi;

   double t427 = 9.0/40.0-1215.0/256.0*t52-t141-t206+t207-t417+t350-t418-10989.0/640.0*t66+t420+t421+t300-2781.0/256.0*t74-t423+t424-t425-t426;

   double t442 = -3.0/10.0+9.0/32.0*t2+81.0/160.0*t5+t161-567.0/16.0*t23-1089.0/35.0*t25-15.0/4.0*t27-783.0/32.0*t29-87.0/160.0*t1-162.0/35.0*t369+81.0/35.0*t372-t314-243.0/16.0*t33-27.0/4.0*t35+t317+243.0/20.0*t39+243.0/20.0*t42;

   double t454 = -t174+405.0/32.0*t52+t177+t178-t179+607.0/35.0*t60-t326+27.0/4.0*t64+189.0/8.0*t66-1152.0/35.0*t12-243.0/10.0*t70-t330+261.0/16.0*t74+9.0/20.0*t76-27.0/2.0*t22+9.0*t9+2547.0/140.0*xi;

   double t462 = -81.0/128.0*t2-729.0/640.0*t5-t125+3483.0/128.0*t23+t405+t406+5103.0/256.0*t29+351.0/640.0*t1+t409-t410+t286+24057.0/1280.0*t33+t412-t289-t413-t414+t138;

   double t466 = 9.0/40.0-3645.0/256.0*t52-t141-t142+t143-t417+t297-t418-11691.0/640.0*t66+t420+t421+t300-3483.0/256.0*t74-t423+t424-t425-t426;

   double t474 = -9.0/100.0+27.0/160.0*t2+243.0/800.0*t5+t91+t93-891.0/80.0*t23-t95+t96-1377.0/160.0*t29-153.0/800.0*t1-t401+t402-t254-729.0/80.0*t33-t101+t257;

   double t479 = t259+t104+243.0/32.0*t52+t107-t108+837.0/175.0*t60-t264+297.0/40.0*t66-t399-t268-t269+459.0/80.0*t74+t271-t400+t398+t397;

   double t487 = 45.0/128.0*t2-t223-t7-t21+t226+t363-291.0/128.0*t27+1971.0/1280.0*t29+10111.0/22400.0*t1+t370-t373+t230+1701.0/1280.0*t33+1269.0/640.0*t35-t233+t378-t379+t48;

   double t492 = -159.0/5600.0-t50-405.0/256.0*t52+t55-t57+t59-6877.0/11200.0*t60+t240+t384-783.0/640.0*t66+t386+t387+t242-1179.0/1280.0*t74-t390+t391-t392-t393;

   double t496 = 81.0/40.0*t49;

   double t498 = 729.0/40.0*t20;

   double t499 = 729.0/2800.0+t213-27.0/20.0*t60-t355-t218-t357+27.0/2.0*t27+t248+t496+t358-243.0/20.0*t35-t261+t498;

   double t502 = 4509.0/896.0*t4;

   double t503 = 729.0/16.0*t20;

   double t506 = 10935.0/256.0*t29;

   double t507 = 25461.0/4480.0*t1;

   double t509 = 657.0/128.0*t2-567.0/128.0*t5-t502-t503-1215.0/128.0*t23+t340-4383.0/128.0*t27-t506+t507-t201-t202+3969.0/128.0*t35+t203-t346;

   double t510 = 567.0/128.0*t49;

   double t511 = 5103.0/256.0*t58;

   double t514 = 5103.0/256.0*t74;

   double t515 = -729.0/1120.0+t291-t510-t204-t141+t142-t511+207.0/64.0*t60+t208+t351+567.0/128.0*t66+t150+t514+t302;

   double t528 = -45.0/8.0*t2+81.0/16.0*t5+45.0/7.0*t4+243.0/4.0*t20-405.0/32.0*t23+t310+375.0/8.0*t27+675.0/8.0*t29-251.0/35.0*t1+t168+1701.0/32.0*t33-675.0/16.0*t35-t171+t318-t319+45.0/8.0*t49;

   double t537 = 26.0/35.0+2025.0/16.0*t52+t177-1215.0/16.0*t56+567.0/16.0*t58-75.0/16.0*t60-t181-t327+189.0/32.0*t66+t329-t186-315.0/8.0*t74-t332-135.0/16.0*t22+63.0/16.0*t9-t335;

   double t541 = 3645.0/128.0*t23;

   double t543 = 15309.0/256.0*t33;

   double t545 = 423.0/128.0*t2-405.0/128.0*t5-t502-t503+t541-t282-4617.0/128.0*t27-t506+t507-t132-t543+4131.0/128.0*t35+t135-t136;

   double t548 = 1701.0/128.0*t66;

   double t549 = -729.0/1120.0+t291-t510-40095.0/256.0*t52-t141+t142-t511+243.0/64.0*t60+t145+t298-t548+t150+t514+t302;

   double t558 = -9.0/8.0*t2+81.0/80.0*t5+t248+t498-243.0/32.0*t23+t250+117.0/8.0*t27-t252-1863.0/700.0*t1+t99+5103.0/160.0*t33-1053.0/80.0*t35-t102+t258+t259;

   double t562 = 1107.0/2800.0-t261+t496+1215.0/16.0*t52-117.0/80.0*t60-t110-t265+567.0/160.0*t66-t267-t268-t114+t270+t271+t272-t273+t274;

   double t566 = 2187.0/100.0*t20;

   double t568 = 81.0/5.0*t27;

   double t569 = 729.0/20.0*t29;

   double t570 = 81.0/100.0*t1;

   double t573 = 729.0/50.0*t35;

   double t575 = 2187.0/100.0*t42;

   double t576 = 243.0/100.0*t49;

   double t577 = -81.0/10.0*t2+729.0/100.0*t5-t566+729.0/20.0*t23-t568+t569+t570+81.0/100.0*eta-2187.0/10.0*t33+t573+6561.0/20.0*t37-t575-t576;

   double t579 = 81.0/50.0*t60;

   double t582 = 243.0/50.0*t12;

   double t583 = 2187.0/50.0*t70;

   double t585 = 243.0/10.0*t74;

   double t586 = 243.0/50.0*t76;

   double t587 = 729.0/20.0*t22;

   double t588 = 243.0/10.0*t9;

   double t589 = 243.0/100.0*xi;

   double t590 = -81.0/100.0-729.0/2.0*t52+t579+243.0*t62-243.0/10.0*t66+t582+t583+243.0/50.0*t72-t585-t586-t587+t588-t589;

   double t592 = 6561.0/32.0*t29;

   double t596 = 2187.0/40.0*t20;

   double t599 = 6561.0/32.0*t56;

   double t600 = 729.0/20.0*t35;

   double t601 = 243.0/40.0*t49;

   double t603 = 2187.0/16.0*t58;

   double t605 = 729.0/160.0*t4;

   double t606 = 8019.0/16.0*t62;

   double t607 = 81.0/2.0*t27;

   double t608 = 19683.0/32.0*t37;

   double t609 = 243.0/80.0*eta;

   double t610 = 81.0/20.0*t60;

   double t611 = 729.0/160.0*t1;

   double t612 = 2187.0/16.0*t74;

   double t613 = 243.0/40.0*t72;

   double t614 = t603+24057.0/32.0*t52-t605-t606+t607-t608-t609-t610+t611-t612-t613;

   double t617 = 81.0/2.0*t66;

   double t618 = 243.0/4.0*t23;

   double t619 = 27.0/2.0*t2;

   double t623 = 729.0/2.0*t33;

   double t624 = 243.0/20.0*t5;

   double t627 = 9.0/10.0-405.0*t29-t617+t618-t619-27.0*t9+81.0/2.0*t22-729.0/10.0*t20-t623+t624+729.0/2.0*t56+243.0/5.0*t35;

   double t630 = 1215.0/2.0*t52;

   double t632 = 405.0*t62;

   double t634 = 2187.0/4.0*t37;

   double t635 = 27.0/20.0*eta;

   double t639 = 81.0/10.0*t72;

   double t640 = -81.0/10.0*t49-243.0*t58-t630+81.0/10.0*t4+t632-54.0*t27+t634+t635+27.0/5.0*t60-9.0*t1+270.0*t74+t639;

   double t649 = 1701.0/16.0*t62;

   double t650 = 6561.0/32.0*t37;

   double t651 = 81.0/80.0*eta;

   double t652 = t603+5103.0/32.0*t52-t605-t649+t607-t650+t651-t610+t611-t612-t613;

   double t657 = 729.0/40.0*t25;

   double t658 = 729.0/8.0*t27;

   double t659 = 729.0/8.0*t29;

   double t660 = 81.0/40.0*t1;

   double t662 = 6561.0/80.0*t35;

   double t663 = 6561.0/40.0*t39;

   double t664 = 2187.0/40.0*t42;

   double t665 = 243.0/8.0*t2-2187.0/80.0*t5-2187.0/32.0*t23-t657-t658-t659-t660-t609+72171.0/160.0*t33+t662-t608-t663+t664;

   double t667 = 729.0/80.0*t60;

   double t668 = 729.0/4.0*t64;

   double t670 = 243.0/20.0*t12;

   double t671 = 2187.0/20.0*t70;

   double t672 = 243.0/4.0*t74;

   double t673 = 243.0/20.0*t76;

   double t674 = 729.0/8.0*t22;

   double t675 = 243.0/4.0*t9;

   double t676 = 243.0/40.0*xi;

   double t677 = 81.0/40.0+10935.0/16.0*t52+t667-t606+t668+8019.0/160.0*t66-t670-t671-t613+t672+t673+t674-t675+t676;

   double t681 = 2187.0/128.0*t4;

   double t683 = 2187.0/32.0*t25;

   double t684 = 24057.0/128.0*t27;

   double t685 = 98415.0/256.0*t29;

   double t686 = 2187.0/128.0*t1;

   double t689 = 19683.0/128.0*t35;

   double t691 = -8019.0/128.0*t2+6561.0/128.0*t5+t681+32805.0/128.0*t23+t683+t684-t685-t686+729.0/64.0*eta-216513.0/256.0*t33-t689+295245.0/256.0*t37;

   double t692 = 19683.0/64.0*t39;

   double t693 = 6561.0/128.0*t49;

   double t695 = 6561.0/64.0*t54;

   double t696 = 98415.0/256.0*t56;

   double t697 = 72171.0/256.0*t58;

   double t698 = 2187.0/64.0*t60;

   double t700 = 24057.0/64.0*t64;

   double t702 = 6561.0/128.0*t72;

   double t703 = 72171.0/256.0*t74;

   double t704 = 6561.0/64.0*t76;

   double t705 = t692-t693-360855.0/256.0*t52+t695+t696-t697-t698+264627.0/256.0*t62-t700-24057.0/128.0*t66+t702+t703-t704;

   double t711 = 243.0/8.0*t25;

   double t712 = 1215.0/8.0*t27;

   double t715 = 81.0/16.0*eta;

   double t717 = 2187.0/16.0*t35;

   double t718 = 32805.0/32.0*t37;

   double t719 = 2187.0/8.0*t39;

   double t720 = 729.0/8.0*t49;

   double t721 = -27.0/8.0+405.0/8.0*t2-729.0/16.0*t5-243.0/8.0*t4-3645.0/32.0*t23-t711-t712+6075.0/8.0*t29+135.0/4.0*t1-t715+24057.0/32.0*t33+t717-t718-t719+t720;

   double t723 = 729.0/4.0*t54;

   double t726 = 243.0/16.0*t60;

   double t727 = 13365.0/16.0*t62;

   double t728 = 1215.0/4.0*t64;

   double t730 = 81.0/4.0*t12;

   double t731 = 405.0/4.0*t72;

   double t733 = 405.0/2.0*t76;

   double t736 = 81.0/8.0*xi;

   double t737 = 18225.0/16.0*t52-t723-10935.0/16.0*t56+8019.0/16.0*t58+t726-t727+t728+2673.0/32.0*t66-t730-t731-4455.0/8.0*t74+t733-1215.0/16.0*t22+891.0/16.0*t9+t736;

   double t740 = 2187.0/128.0*t5;

   double t741 = 10935.0/128.0*t23;

   double t742 = 729.0/32.0*t25;

   double t743 = 5103.0/128.0*t27;

   double t744 = 243.0/64.0*eta;

   double t746 = 6561.0/128.0*t35;

   double t747 = 98415.0/256.0*t37;

   double t748 = -1701.0/128.0*t2+t740+t681-t741-t742+t743-t685-t686-t744-72171.0/256.0*t33-t746+t747;

   double t749 = 6561.0/64.0*t39;

   double t751 = 729.0/64.0*t60;

   double t752 = 56133.0/256.0*t62;

   double t753 = 5103.0/64.0*t64;

   double t755 = t749-t693-76545.0/256.0*t52+t695+t696-t697+t751+t752-t753+8019.0/128.0*t66+t702+t703-t704;

   double t764 = -t619+t624+t618+162.0/5.0*t25+180.0*t27+243.0/2.0*t29+27.0/10.0*t1+t635-t623-162.0*t35+t634+1458.0/5.0*t39-729.0/10.0*t42;

   double t774 = -27.0/10.0-t630-18.0*t60+t632-324.0*t64-t617+81.0/5.0*t12+729.0/5.0*t70+t639-81.0*t74-81.0/5.0*t76-243.0/2.0*t22+81.0*t9-81.0/10.0*xi;

   double t778 = 243.0/32.0*t4;

   double t782 = 10935.0/32.0*t29;

   double t783 = 243.0/32.0*t1;

   double t786 = 891.0/32.0*t2-729.0/32.0*t5-t778-3645.0/16.0*t23-243.0/2.0*t25-1485.0/4.0*t27+t782+t783-t715+10935.0/16.0*t33+1215.0/4.0*t35-t718;

   double t788 = 405.0/4.0*t49;

   double t790 = 10935.0/32.0*t56;

   double t791 = 3645.0/16.0*t58;

   double t795 = 3645.0/16.0*t74;

   double t796 = 729.0/4.0*t76;

   double t797 = -2187.0/4.0*t39+t788+40095.0/32.0*t52-t723-t790+t791+135.0/2.0*t60-t727+2673.0/4.0*t64+1215.0/8.0*t66-t731-t795+t796;

   double t813 = 3.0/2.0-45.0/2.0*t2+81.0/4.0*t5+27.0/2.0*t4+405.0/4.0*t23+54.0*t25+300.0*t27-675.0*t29-15.0*t1+9.0/4.0*eta-1215.0/2.0*t33-270.0*t35+3645.0/4.0*t37+486.0*t39-180.0*t49;

   double t827 = -2025.0/2.0*t52+324.0*t54+1215.0/2.0*t56-405.0*t58-30.0*t60+675.0*t62-540.0*t64-135.0/2.0*t66+36.0*t12+200.0*t72+450.0*t74-360.0*t76+t189-t190-20.0*xi;

   double t834 = 27.0/16.0*eta;

   double t837 = 10935.0/32.0*t37;

   double t838 = 189.0/32.0*t2-243.0/32.0*t5-t778+1215.0/16.0*t23+81.0/2.0*t25-315.0/4.0*t27+t782+t783+t834+3645.0/16.0*t33+405.0/4.0*t35-t837;

   double t842 = 2835.0/16.0*t62;

   double t845 = -729.0/4.0*t39+t788+8505.0/32.0*t52-t723-t790+t791-45.0/2.0*t60-t842+567.0/4.0*t64-405.0/8.0*t66-t731-t795+t796;

   double t851 = -81.0/8.0*t2+729.0/80.0*t5-729.0/32.0*t23-t657-t658-t659-t660+t651+15309.0/160.0*t33+t662-t650-t663+t664;

   double t854 = 81.0/40.0+3645.0/16.0*t52+t667-t649+t668+1701.0/160.0*t66-t670-t671-t613+t672+t673+t674-t675+t676;

   double t857 = 729.0/128.0*t4;

   double t858 = 32805.0/256.0*t29;

   double t859 = 729.0/128.0*t1;

   double t861 = 2673.0/128.0*t2-t740-t857+t741+t683+t684-t858+t859-t744-45927.0/256.0*t33-t689+t747;

   double t863 = 32805.0/256.0*t56;

   double t864 = 15309.0/256.0*t58;

   double t866 = 15309.0/256.0*t74;

   double t867 = t692-t693-120285.0/256.0*t52+t695+t863-t864-t698+t752-t700-5103.0/128.0*t66+t702+t866-t704;

   double t876 = 9.0/8.0-135.0/8.0*t2+243.0/16.0*t5+81.0/8.0*t4-1215.0/32.0*t23-t711-t712+2025.0/8.0*t29-45.0/4.0*t1+t834+5103.0/32.0*t33+t717-t837-t719+t720;

   double t884 = 6075.0/16.0*t52-t723-3645.0/16.0*t56+1701.0/16.0*t58+t726-t842+t728+567.0/32.0*t66-t730-t731-945.0/8.0*t74+t733-405.0/16.0*t22+189.0/16.0*t9+t736;

   double t887 = t123-t124-t857-t541-t742+t743-t858+t859+81.0/64.0*eta-t543-t746+t135;

   double t889 = t749-t693-t139+t695+t863-t864+t751+11907.0/256.0*t62-t753+t548+t702+t866-t704;

   double t891 = -81.0/100.0+t570+t114-t585-t586+t569-t589+t588+t582-t587-t575+t583;


   values[0] = t51+t87;

   values[1] = t105+t121;

   values[2] = t140+t157;

   values[3] = t175+t193;

   values[4] = t205+t211;

   values[5] = t221;

   values[6] = t235+t244;

   values[7] = t260+t275;

   values[8] = t290+t303;

   values[9] = t321+t336;

   values[10] = t347+t352;

   values[11] = t359;

   values[12] = t380+t394;

   values[13] = t403;

   values[14] = t415+t427;

   values[15] = t442+t454;

   values[16] = t462+t466;

   values[17] = t474+t479;

   values[18] = t487+t492;

   values[19] = t499;

   values[20] = t509+t515;

   values[21] = t528+t537;

   values[22] = t545+t549;

   values[23] = t558+t562;

   values[24] = t577+t590;

   values[25] = t592+729.0/8.0*t66-2187.0/16.0*t23+2673.0/160.0*t2+t596+6561.0/16.0*t33-2187.0/160.0*t5-t599-t600+t601+t614;

   values[26] = t627+t640;

   values[27] = t592-243.0/8.0*t66+729.0/16.0*t23+567.0/160.0*t2+t596+2187.0/16.0*t33-729.0/160.0*t5-t599-t600+t601+t652;

   values[28] = t665+t677;

   values[29] = t691+t705;

   values[30] = t721+t737;

   values[31] = t748+t755;

   values[32] = t764+t774;

   values[33] = t786+t797;

   values[34] = t813+t827;

   values[35] = t838+t845;

   values[36] = t851+t854;

   values[37] = t861+t867;

   values[38] = t876+t884;

   values[39] = t887+t889;

   values[40] = t891;

   values[41] = t579+t114-t568-t576+t573-t566;

 }


 // values of the derivatives in eta direction

 static void C_Q_UL6S_2D_DeriveEta(double xi, double eta, double *values)

 {

   double t1 = xi*xi;

   double t2 = t1*t1;

   double t3 = t2*xi;

   double t4 = eta*eta;

   double t5 = t4*eta;

   double t6 = t3*t5;

   double t7 = 243.0/64.0*t6;

   double t8 = 3303.0/2800.0*eta;

   double t10 = 1.0-t1;

   double t12 = t10*eta*t2;

   double t13 = 243.0/80.0*t12;

   double t14 = 1.0-xi;

   double t16 = t4*t4;

   double t19 = 3.0/8.0+35.0/8.0*t16-15.0/4.0*t4;

   double t20 = t14*eta*t19;

   double t21 = 81.0/350.0*t20;

   double t22 = 1.0-t4;

   double t26 = 35.0/2.0*t5-15.0/2.0*eta;

   double t27 = t14*t22*t26;

   double t28 = 81.0/700.0*t27;

   double t29 = t1*t16;

   double t32 = 9.0/4.0*t4;

   double t33 = 153.0/70.0*t5;

   double t35 = t1*xi;

   double t36 = 27.0/128.0*t35;

   double t37 = 2007.0/4480.0*t2;

   double t38 = 27.0/8.0*t16;

   double t42 = t10*(3.0/8.0+35.0/8.0*t2-15.0/4.0*t1);

   double t43 = 81.0/700.0*t42;

   double t44 = xi*t4;

   double t45 = 81.0/128.0*t44;

   double t46 = t35*t4;

   double t48 = -663.0/5600.0-t7-t8-19.0/320.0*xi-t13+t21-t28+1431.0/256.0*t29+12127.0/22400.0*t1+t32+t33-81.0/640.0*t3+t36-t37-t38+t43+t45+1323.0/256.0*t46;

   double t49 = t1*t4;

   double t52 = t10*t22*t5;

   double t53 = 243.0/40.0*t52;

   double t55 = t10*t16*eta;

   double t56 = 243.0/80.0*t55;

   double t57 = t1*eta;

   double t59 = xi*eta;

   double t61 = t1*t5;

   double t62 = 261.0/320.0*t61;

   double t63 = t35*t5;

   double t64 = 189.0/64.0*t63;

   double t65 = xi*t16;

   double t67 = t2*t5;

   double t68 = 81.0/64.0*t67;

   double t69 = t2*t16;

   double t70 = 1215.0/256.0*t69;

   double t71 = t3*t16;

   double t73 = t3*t4;

   double t75 = t35*eta;

   double t77 = t2*eta;

   double t78 = 261.0/640.0*t77;

   double t79 = t3*eta;

   double t81 = t2*t4;

   double t82 = 567.0/256.0*t81;

   double t83 = xi*t5;

   double t84 = 2007.0/1120.0*t83;

   double t85 = t35*t16;

   double t87 = -927.0/256.0*t49+t53-t56+437.0/640.0*t57+12127.0/11200.0*t59-t62+t64-81.0/128.0*t65+t68-t70+3645.0/256.0*t71-1701.0/256.0*t73-309.0/128.0*t75-t78+1431.0/640.0*t79+t82-t84-2835.0/256.0*t85;

   double t89 = 729.0/40.0*t6;

   double t90 = 81.0/200.0*eta;

   double t91 = 9.0/80.0*xi;

   double t92 = 729.0/40.0*t12;

   double t93 = 243.0/40.0*t29;

   double t95 = 81.0/20.0*t4;

   double t96 = 81.0/100.0*t5;

   double t99 = 1269.0/560.0*t2;

   double t100 = 243.0/40.0*t16;

   double t101 = 243.0/350.0*t42;

   double t103 = 1107.0/2800.0+t89+t90+t91+t92-t93-1863.0/700.0*t1-t95-t96+81.0/80.0*t3-9.0/8.0*t35+t99+t100-t101+567.0/160.0*t44;

   double t104 = 567.0/16.0*t46;

   double t105 = 81.0/20.0*t49;

   double t106 = 729.0/100.0*t52;

   double t107 = 729.0/200.0*t55;

   double t108 = 243.0/100.0*t57;

   double t110 = 81.0/100.0*t61;

   double t111 = 81.0/4.0*t63;

   double t113 = 2187.0/32.0*t71;

   double t116 = 81.0/40.0*t77;

   double t118 = 81.0/40.0*t83;

   double t120 = -t104+t105-t106+t107-t108-117.0/80.0*t59+t110-t111-243.0/32.0*t65-t113+5103.0/160.0*t73+117.0/8.0*t75+t116-1053.0/80.0*t79+t118+1215.0/16.0*t85;

   double t122 = 2187.0/64.0*t6;

   double t123 = 9.0/64.0*xi;

   double t124 = 729.0/16.0*t12;

   double t125 = 10935.0/256.0*t29;

   double t126 = 25461.0/4480.0*t1;

   double t129 = 4509.0/896.0*t2;

   double t130 = 243.0/140.0*t42;

   double t131 = 1701.0/128.0*t44;

   double t132 = 18711.0/256.0*t46;

   double t133 = 5103.0/256.0*t49;

   double t134 = 567.0/128.0*t57;

   double t135 = -729.0/1120.0-t122-t123-t124-t125+t126-405.0/128.0*t3+423.0/128.0*t35-t129+t130-t131+t132+t133+t134;

   double t137 = 729.0/64.0*t61;

   double t138 = 2673.0/64.0*t63;

   double t139 = 3645.0/128.0*t65;

   double t140 = 729.0/64.0*t67;

   double t141 = 10935.0/256.0*t69;

   double t142 = 32805.0/256.0*t71;

   double t143 = 15309.0/256.0*t73;

   double t145 = 567.0/128.0*t77;

   double t147 = 5103.0/256.0*t81;

   double t148 = 243.0/32.0*t83;

   double t150 = 243.0/64.0*t59+t137+t138+t139-t140+t141+t142-t143-4617.0/128.0*t75-t145+4131.0/128.0*t79-t147-t148-40095.0/256.0*t85;

   double t152 = 243.0/8.0*t6;

   double t153 = eta/8.0;

   double t154 = 9.0/16.0*xi;

   double t159 = 9.0/4.0*t5;

   double t164 = 81.0/35.0*t42;

   double t166 = 945.0/16.0*t46;

   double t167 = 26.0/35.0+t152-t153+t154+243.0/4.0*t12+675.0/8.0*t29-251.0/35.0*t1+63.0/16.0*t4+t159+81.0/16.0*t3-45.0/8.0*t35+45.0/7.0*t2-135.0/16.0*t16-t164+189.0/32.0*t44-t166;

   double t169 = 11.0/2.0*t57;

   double t171 = 45.0/2.0*t61;

   double t172 = 135.0/4.0*t63;

   double t174 = 81.0/4.0*t67;

   double t176 = 3645.0/32.0*t71;

   double t182 = 27.0/8.0*t83;

   double t184 = -315.0/8.0*t49-t169-75.0/16.0*t59-t171-t172-405.0/32.0*t65+t174-1215.0/16.0*t69-t176+1701.0/32.0*t73+375.0/8.0*t75+45.0/8.0*t77-675.0/16.0*t79+567.0/16.0*t81+t182+2025.0/16.0*t85;

   double t186 = 729.0/64.0*t6;

   double t187 = 45.0/64.0*xi;

   double t191 = 3969.0/256.0*t46;

   double t192 = -729.0/1120.0-t186-t187-t124-t125+t126-567.0/128.0*t3+657.0/128.0*t35-t129+t130+567.0/128.0*t44+t191+t133+t134;

   double t194 = 567.0/64.0*t63;

   double t196 = 10935.0/256.0*t71;

   double t197 = 5103.0/256.0*t73;

   double t200 = 81.0/32.0*t83;

   double t201 = 8505.0/256.0*t85;

   double t202 = 207.0/64.0*t59+t137+t194-1215.0/128.0*t65-t140+t141+t196-t197-4383.0/128.0*t75-t145+3969.0/128.0*t79-t147+t200-t201;

   double t204 = 9.0/40.0*xi;

   double t205 = 3537.0/1400.0*t1;

   double t206 = 9.0/4.0*t35;

   double t207 = 81.0/40.0*t3;

   double t209 = 81.0/40.0*t57;

   double t212 = 729.0/2800.0+t204-t205-t206+t99+t207-27.0/20.0*t59-t209+27.0/2.0*t75+t116-243.0/20.0*t79-t101+t92;

   double t213 = 243.0/320.0*t6;

   double t214 = 2547.0/2800.0*eta;

   double t215 = xi/320.0;

   double t218 = 9.0/20.0*t4;

   double t219 = 288.0/175.0*t5;

   double t220 = 243.0/640.0*t3;

   double t222 = 27.0/40.0*t16;

   double t224 = 189.0/256.0*t46;

   double t226 = 1.0+xi;

   double t228 = t226*eta*t19;

   double t229 = 81.0/350.0*t228;

   double t230 = -159.0/5600.0+t213-t214+t215-t13+1971.0/1280.0*t29+10111.0/22400.0*t1-t218+t219-t220+45.0/128.0*t35-t37+t222+t43-783.0/640.0*t44+t224-1179.0/1280.0*t49+t229;

   double t232 = t226*t22*t26;

   double t233 = 81.0/700.0*t232;

   double t234 = 243.0/200.0*t52;

   double t235 = 243.0/400.0*t55;

   double t236 = 1321.0/3200.0*t57;

   double t238 = 441.0/1600.0*t61;

   double t239 = 27.0/64.0*t63;

   double t240 = 243.0/128.0*t65;

   double t241 = 729.0/256.0*t71;

   double t245 = 1629.0/1120.0*t83;

   double t247 = -t233+t234-t235+t236-6877.0/11200.0*t59-t238+t239+t240+t68-t70-t241+1701.0/1280.0*t73-291.0/128.0*t75-t78+1269.0/640.0*t79+t82+t245-405.0/256.0*t85;

   double t249 = 7641.0/1400.0*eta;

   double t250 = 99.0/400.0*xi;

   double t251 = 729.0/200.0*t12;

   double t254 = 27.0/10.0*t4;

   double t255 = 1728.0/175.0*t5;

   double t258 = 81.0/160.0*t2;

   double t259 = 81.0/20.0*t16;

   double t261 = 81.0/16.0*t46;

   double t263 = 243.0/175.0*t228;

   double t264 = -9.0/100.0+t249-t250+t251-1377.0/160.0*t29-153.0/800.0*t1+t254-t255+243.0/800.0*t3+27.0/160.0*t35+t258-t259+297.0/40.0*t44-t261+459.0/80.0*t49-t263;

   double t265 = 243.0/350.0*t232;

   double t266 = 81.0/100.0*t57;

   double t269 = 729.0/32.0*t69;

   double t270 = 2187.0/160.0*t71;

   double t272 = 27.0/10.0*t75;

   double t273 = 81.0/200.0*t77;

   double t274 = 243.0/100.0*t79;

   double t275 = 243.0/16.0*t81;

   double t276 = 1269.0/140.0*t83;

   double t278 = t265-t106+t107-t266+837.0/175.0*t59+t110-891.0/80.0*t65+t269+t270-729.0/80.0*t73+t272+t273-t274-t275-t276+243.0/32.0*t85;

   double t280 = 2187.0/320.0*t6;

   double t281 = 7641.0/560.0*eta;

   double t282 = 207.0/320.0*xi;

   double t285 = 27.0/4.0*t4;

   double t286 = 864.0/35.0*t5;

   double t289 = 243.0/128.0*t2;

   double t290 = 81.0/8.0*t16;

   double t292 = 2673.0/256.0*t46;

   double t294 = 243.0/70.0*t228;

   double t295 = 243.0/140.0*t232;

   double t296 = 9.0/40.0-t280-t281+t282+5103.0/256.0*t29+351.0/640.0*t1-t285+t286-729.0/640.0*t3-81.0/128.0*t35-t289+t290-11691.0/640.0*t44+t292-3483.0/256.0*t49+t294-t295;

   double t297 = 729.0/40.0*t52;

   double t298 = 729.0/80.0*t55;

   double t299 = 243.0/640.0*t57;

   double t300 = 28863.0/2240.0*t59;

   double t301 = 243.0/320.0*t61;

   double t302 = 243.0/64.0*t63;

   double t304 = 6561.0/256.0*t71;

   double t306 = 243.0/128.0*t75;

   double t307 = 729.0/128.0*t77;

   double t308 = 2187.0/640.0*t79;

   double t309 = 8019.0/256.0*t81;

   double t310 = 25947.0/1120.0*t83;

   double t312 = t297-t298+t299-t300-t301-t302+3483.0/128.0*t65-t140-t141-t304+24057.0/1280.0*t73+t306+t307+t308+t309+t310-3645.0/256.0*t85;

   double t316 = 63.0/80.0*xi;

   double t323 = 27.0/32.0*t2;

   double t326 = 135.0/16.0*t46;

   double t330 = -3.0/10.0+243.0/20.0*t6+2547.0/140.0*eta-t316-783.0/32.0*t29-87.0/160.0*t1+9.0*t4-1152.0/35.0*t5+81.0/160.0*t3+9.0/32.0*t35+t323-27.0/2.0*t16+189.0/8.0*t44-t326+261.0/16.0*t49-162.0/35.0*t228+81.0/35.0*t232;

   double t333 = t57/10.0;

   double t338 = 1215.0/32.0*t69;

   double t339 = 729.0/32.0*t71;

   double t342 = 45.0/4.0*t77;

   double t344 = 405.0/16.0*t81;

   double t347 = -243.0/10.0*t52+243.0/20.0*t55-t333+607.0/35.0*t59+9.0/20.0*t61+27.0/4.0*t63-567.0/16.0*t65+t174+t338+t339-243.0/16.0*t73-15.0/4.0*t75-t342-27.0/4.0*t79-t344-1089.0/35.0*t83+405.0/32.0*t85;

   double t349 = 171.0/320.0*xi;

   double t352 = 81.0/128.0*t2;

   double t354 = 567.0/256.0*t46;

   double t356 = 9.0/40.0-t280-t281+t349+4293.0/256.0*t29+171.0/640.0*t1-t285+t286+t220+t36+t352+t290-10989.0/640.0*t44+t354-2781.0/256.0*t49+t294-t295;

   double t358 = 3645.0/256.0*t69;

   double t359 = 2187.0/256.0*t71;

   double t361 = 1701.0/256.0*t81;

   double t363 = t297-t298+t299-t300-t301-t302+3321.0/128.0*t65-t140-t358-t359+5103.0/1280.0*t73+t306+t307+t308+t361+t310-1215.0/256.0*t85;

   double t366 = 3537.0/700.0*t59;

   double t367 = 81.0/200.0*t57;

   double t368 = 27.0/4.0*t44;

   double t369 = 81.0/8.0*t65;

   double t370 = -9.0/100.0-t204-27.0/200.0*t1+t366-t367+t368+t105-t276+t110-t369-t93+t249+t254-t255-t259-t263+t265+t107-t106;

   double t372 = 243.0/400.0*t12;

   double t377 = 1629.0/4480.0*t2;

   double t381 = 327.0/5600.0+t213-t214+121.0/1600.0*xi-t372+1161.0/1280.0*t29-8893.0/22400.0*t1-t218+t219+1053.0/3200.0*t3-243.0/640.0*t35+t377+t222-t43-729.0/640.0*t44+27.0/256.0*t46-909.0/1280.0*t49+t229;

   double t385 = 243.0/256.0*t69;

   double t389 = 441.0/3200.0*t77;

   double t391 = 81.0/256.0*t81;

   double t393 = -t233+t234-t235+457.0/3200.0*t57-8893.0/11200.0*t59-t238+t239+1053.0/640.0*t65+t68+t385+729.0/1280.0*t71+243.0/1280.0*t73-303.0/640.0*t75-t389+1161.0/3200.0*t79+t391+t245+81.0/256.0*t85;

   double t396 = -729.0/2800.0-t204+t205+t206-t99-t207-27.0/100.0*t59-t367+t272+t273-t274+t101+t251;

   double t397 = 729.0/80.0*t12;

   double t398 = 2187.0/256.0*t29;

   double t399 = 28863.0/4480.0*t1;

   double t402 = 25947.0/4480.0*t2;

   double t403 = 729.0/256.0*t49;

   double t404 = 729.0/1120.0-t186+t349-t397+t398-t399+3321.0/640.0*t3-3663.0/640.0*t35+t402-t130+t45+t354+t403+t299;

   double t406 = 2187.0/256.0*t69;

   double t409 = 243.0/640.0*t77;

   double t411 = 729.0/256.0*t81;

   double t413 = 171.0/320.0*t59+t137+t194+t240-t140-t406-t359-729.0/256.0*t73-927.0/128.0*t75-t409+4293.0/640.0*t79-t411+t200+1701.0/256.0*t85;

   double t424 = -25.0/28.0+t152-t153-t316+243.0/20.0*t12-135.0/8.0*t29+607.0/70.0*t1+9.0/16.0*t4+t159-567.0/80.0*t3+63.0/8.0*t35-1089.0/140.0*t2+27.0/16.0*t16+t164+27.0/32.0*t44-t326;

   double t435 = -45.0/8.0*t49-t333-87.0/80.0*t59-t171-t172+81.0/32.0*t65+t174+243.0/16.0*t69+t339+243.0/32.0*t73+87.0/8.0*t75+9.0/40.0*t77-783.0/80.0*t79+81.0/16.0*t81+t182-405.0/16.0*t85;

   double t440 = 729.0/1120.0-t122+t282-t397+t398-t399+3483.0/640.0*t3-3897.0/640.0*t35+t402-t130-243.0/128.0*t44+t292+t403+t299;

   double t447 = 351.0/320.0*t59+t137+t138-729.0/128.0*t65-t140-t406-t304-2187.0/256.0*t73-1161.0/128.0*t75-t409+5103.0/640.0*t79-t411-t148+8019.0/256.0*t85;

   double t453 = -351.0/2800.0+t89+t90-t250+t251-t93+837.0/350.0*t1-t95-t96-891.0/400.0*t3+99.0/40.0*t35-t99+t100+t101+81.0/160.0*t44;

   double t460 = -t261+t105-t106+t107-t266-153.0/400.0*t59+t110-t111+243.0/160.0*t65+t270+729.0/160.0*t73+153.0/40.0*t75+t273-1377.0/400.0*t79+t118-243.0/16.0*t85;

   double t467 = -t7-t8+t215-t372+t21-t28+1269.0/256.0*t29-6877.0/22400.0*t1+t32+t33+t220-261.0/640.0*t35+t377-t38-t43+135.0/128.0*t44+t224-873.0/256.0*t49;

   double t473 = -177.0/5600.0+t53-t56+t236+10111.0/11200.0*t59-t62+t64-t240+t68+t385-t241-243.0/256.0*t73-393.0/640.0*t75-t389+1971.0/3200.0*t79+t391-t84+567.0/256.0*t85;

   double t477 = 81.0/20.0*t61;

   double t479 = 9909.0/1400.0*eta;

   double t480 = 27.0/2.0*t4;

   double t481 = 459.0/35.0*t5;

   double t482 = 81.0/4.0*t16;

   double t483 = 243.0/175.0*t20;

   double t484 = 243.0/350.0*t27;

   double t485 = 729.0/40.0*t55;

   double t486 = 729.0/20.0*t52;

   double t487 = 9.0/20.0+t204-27.0/40.0*t1-t366-t209-t368+81.0/4.0*t49+t276+t477+t369-243.0/8.0*t29+t479-t480-t481+t482-t483+t484+t485-t486;

   double t488 = 9909.0/560.0*eta;

   double t489 = 243.0/70.0*t20;

   double t490 = 243.0/140.0*t27;

   double t493 = 135.0/4.0*t4;

   double t494 = 459.0/14.0*t5;

   double t497 = 405.0/8.0*t16;

   double t500 = -9.0/8.0+t122-t488-t187+t489-t490+19845.0/256.0*t29+207.0/128.0*t1+t493+t494-243.0/128.0*t3+189.0/128.0*t35+t352-t497+1971.0/128.0*t44+t191-13149.0/256.0*t49;

   double t501 = 729.0/8.0*t52;

   double t502 = 729.0/16.0*t55;

   double t503 = 25461.0/2240.0*t59;

   double t504 = 567.0/64.0*t61;

   double t505 = 1701.0/64.0*t63;

   double t507 = 1701.0/128.0*t75;

   double t508 = 2187.0/128.0*t79;

   double t509 = 4509.0/224.0*t83;

   double t510 = t501-t502+t134+t503-t504-t505-2835.0/128.0*t65-t140-t358+t196-t197+t507+t307-t508+t361-t509-t201;

   double t518 = 45.0*t4;

   double t522 = 135.0/2.0*t16;

   double t525 = 3.0/2.0-243.0/4.0*t6+3303.0/140.0*eta+t154-162.0/35.0*t20+81.0/35.0*t27-3375.0/32.0*t29-75.0/32.0*t1-t518-306.0/7.0*t5-81.0/32.0*t3+63.0/32.0*t35+t323+t522-135.0/8.0*t44-t166+1125.0/16.0*t49;

   double t537 = -243.0/2.0*t52+243.0/4.0*t55-t169-502.0/35.0*t59+45.0/4.0*t61+189.0/4.0*t63+405.0/16.0*t65+t174+t338-t176+1215.0/16.0*t73-105.0/4.0*t75-t342+135.0/4.0*t79-t344+180.0/7.0*t83+2835.0/32.0*t85;

   double t541 = 729.0/128.0*t3;

   double t542 = 567.0/128.0*t35;

   double t545 = -9.0/8.0+t122-t488-t123+t489-t490+20655.0/256.0*t29+243.0/128.0*t1+t493+t494+t541-t542-t289-t497+1269.0/128.0*t44+t132-13851.0/256.0*t49;

   double t548 = 25515.0/256.0*t85;

   double t549 = t501-t502+t134+t503-t504-t505-2025.0/128.0*t65-t140-t141+t142-24057.0/256.0*t73+t507+t307-t508+t309-t509-t548;

   double t556 = 9.0/20.0+t479+t91+t251-t483+t484-1053.0/32.0*t29-117.0/160.0*t1-t480-t481-243.0/160.0*t3+189.0/160.0*t35+t258+t482-27.0/8.0*t44-t104;

   double t562 = 351.0/16.0*t49-t486+t485-t108-1863.0/350.0*t59+t477+81.0/16.0*t65+t269-t113+729.0/16.0*t73+t272+t273-t274-t275+t276+1701.0/32.0*t85;

   double t564 = 243.0/100.0*eta;

   double t566 = 2187.0/100.0*t12;

   double t567 = 729.0/20.0*t29;

   double t568 = 81.0/100.0*t1;

   double t569 = 243.0/10.0*t4;

   double t570 = 243.0/50.0*t5;

   double t573 = 729.0/20.0*t16;

   double t576 = -81.0/100.0-t564+81.0/100.0*xi-t566+t567+t568+t569+t570+729.0/100.0*t3-81.0/10.0*t35-t573-243.0/10.0*t44+243.0*t46;

   double t577 = 243.0/10.0*t49;

   double t578 = 2187.0/50.0*t52;

   double t579 = 2187.0/100.0*t55;

   double t581 = 81.0/50.0*t59;

   double t582 = 243.0/50.0*t61;

   double t586 = 81.0/5.0*t75;

   double t587 = 243.0/100.0*t77;

   double t588 = 729.0/50.0*t79;

   double t590 = -t577+t578-t579+243.0/50.0*t57+t581-t582+729.0/20.0*t65+6561.0/20.0*t71-2187.0/10.0*t73-t586-t587+t588-729.0/2.0*t85;

   double t592 = 6561.0/40.0*t6;

   double t593 = 243.0/40.0*eta;

   double t594 = 243.0/80.0*xi;

   double t595 = 729.0/8.0*t29;

   double t596 = 81.0/40.0*t1;

   double t597 = 243.0/4.0*t4;

   double t598 = 243.0/20.0*t5;

   double t601 = 729.0/8.0*t16;

   double t603 = 8019.0/16.0*t46;

   double t604 = 81.0/40.0-t592+t593-t594-t595-t596-t597-t598-2187.0/80.0*t3+243.0/8.0*t35+t601+8019.0/160.0*t44-t603;

   double t605 = 243.0/4.0*t49;

   double t606 = 2187.0/20.0*t52;

   double t607 = 2187.0/40.0*t55;

   double t608 = 243.0/40.0*t57;

   double t609 = 729.0/80.0*t59;

   double t610 = 243.0/20.0*t61;

   double t611 = 729.0/4.0*t63;

   double t613 = 19683.0/32.0*t71;

   double t615 = 729.0/8.0*t75;

   double t616 = 6561.0/80.0*t79;

   double t617 = 729.0/40.0*t83;

   double t619 = t605-t606+t607-t608+t609+t610+t611-2187.0/32.0*t65-t613+72171.0/160.0*t73-t615+t616-t617+10935.0/16.0*t85;

   double t623 = 27.0/20.0*xi;

   double t628 = 243.0/20.0*t3;

   double t629 = 27.0/2.0*t35;

   double t631 = 81.0/2.0*t44;

   double t632 = 405.0*t46;

   double t633 = -27.0/10.0+1458.0/5.0*t6-81.0/10.0*eta+t623+243.0/2.0*t29+27.0/10.0*t1+81.0*t4+81.0/5.0*t5+t628-t629-243.0/2.0*t16-t631+t632;

   double t637 = 81.0/10.0*t57;

   double t641 = 243.0/4.0*t65;

   double t642 = 2187.0/4.0*t71;

   double t643 = 729.0/2.0*t73;

   double t647 = 1215.0/2.0*t85;

   double t648 = -81.0*t49+729.0/5.0*t52-729.0/10.0*t55+t637-18.0*t59-81.0/5.0*t61-324.0*t63+t641+t642-t643+180.0*t75-162.0*t79+162.0/5.0*t83-t647;

   double t650 = 81.0/80.0*xi;

   double t654 = 1701.0/16.0*t46;

   double t655 = 81.0/40.0-t592+t593+t650-t595-t596-t597-t598+729.0/80.0*t3-81.0/8.0*t35+t601+1701.0/160.0*t44-t654;

   double t657 = 6561.0/32.0*t71;

   double t660 = t605-t606+t607-t608+t609+t610+t611-729.0/32.0*t65-t657+15309.0/160.0*t73-t615+t616-t617+3645.0/16.0*t85;

   double t662 = 2187.0/40.0*t12;

   double t663 = 729.0/20.0*t79;

   double t665 = 729.0/160.0*t2;

   double t667 = 243.0/40.0*t77;

   double t670 = 6561.0/32.0*t69;

   double t673 = 6561.0/32.0*t29;

   double t674 = 2187.0/16.0*t81;

   double t676 = 2187.0/16.0*t49;

   double t677 = 729.0/160.0*t1;

   double t678 = 81.0/20.0*t59;

   double t679 = 81.0/2.0*t75;

   double t680 = t673+t674+6561.0/16.0*t73-t603-t676-t613-t594+t677-t678+t679-t608;

   double t682 = 19683.0/64.0*t6;

   double t684 = 98415.0/256.0*t29;

   double t685 = 2187.0/128.0*t1;

   double t688 = 2187.0/128.0*t2;

   double t691 = 72171.0/256.0*t49;

   double t692 = 6561.0/128.0*t57;

   double t693 = 2187.0/64.0*t59;

   double t694 = t682+729.0/64.0*xi-t684-t685+6561.0/128.0*t3-8019.0/128.0*t35+t688-24057.0/128.0*t44+264627.0/256.0*t46+t691+t692-t693;

   double t695 = 6561.0/64.0*t61;

   double t696 = 24057.0/64.0*t63;

   double t698 = 6561.0/64.0*t67;

   double t699 = 98415.0/256.0*t69;

   double t702 = 24057.0/128.0*t75;

   double t703 = 6561.0/128.0*t77;

   double t704 = 19683.0/128.0*t79;

   double t705 = 72171.0/256.0*t81;

   double t706 = 2187.0/32.0*t83;

   double t708 = -t695-t696+32805.0/128.0*t65+t698+t699+295245.0/256.0*t71-216513.0/256.0*t73+t702-t703-t704-t705+t706-360855.0/256.0*t85;

   double t711 = 81.0/16.0*xi;

   double t712 = 10935.0/32.0*t29;

   double t713 = 243.0/32.0*t1;

   double t716 = 243.0/32.0*t2;

   double t718 = 13365.0/16.0*t46;

   double t719 = 3645.0/16.0*t49;

   double t720 = 405.0/4.0*t57;

   double t722 = -2187.0/4.0*t6-t711+t712+t713-729.0/32.0*t3+891.0/32.0*t35-t716+1215.0/8.0*t44-t718-t719-t720+135.0/2.0*t59;

   double t723 = 729.0/4.0*t61;

   double t726 = 729.0/4.0*t67;

   double t727 = 10935.0/32.0*t69;

   double t728 = 32805.0/32.0*t71;

   double t731 = 405.0/4.0*t77;

   double t733 = 3645.0/16.0*t81;

   double t736 = t723+2673.0/4.0*t63-3645.0/16.0*t65-t726-t727-t728+10935.0/16.0*t73-1485.0/4.0*t75+t731+1215.0/4.0*t79+t733-243.0/2.0*t83+40095.0/32.0*t85;

   double t738 = 243.0/64.0*xi;

   double t739 = 32805.0/256.0*t29;

   double t740 = 729.0/128.0*t1;

   double t741 = 2187.0/128.0*t3;

   double t743 = 729.0/128.0*t2;

   double t745 = 56133.0/256.0*t46;

   double t746 = 15309.0/256.0*t49;

   double t747 = t682-t738-t739+t740-t741+2673.0/128.0*t35-t743-5103.0/128.0*t44+t745+t746+t692-t693;

   double t748 = 10935.0/128.0*t65;

   double t749 = 32805.0/256.0*t69;

   double t750 = 98415.0/256.0*t71;

   double t752 = 15309.0/256.0*t81;

   double t754 = -t695-t696+t748+t698+t749+t750-45927.0/256.0*t73+t702-t703-t704-t752+t706-120285.0/256.0*t85;

   double t763 = 9.0/10.0-729.0/10.0*t12+243.0/5.0*t79-t629+81.0/10.0*t2+t628-81.0/10.0*t77-t631+81.0/2.0*t16-27.0*t4-t647+729.0/2.0*t69;

   double t770 = t641-405.0*t29-243.0*t81-t643+t632+270.0*t49+t642+t623-9.0*t1+27.0/5.0*t59-54.0*t75+t637;

   double t772 = 2187.0/8.0*t6;

   double t773 = 81.0/8.0*eta;

   double t777 = 81.0/4.0*t5;

   double t784 = -27.0/8.0-t772+t773-t711+6075.0/8.0*t29+135.0/4.0*t1+891.0/16.0*t4-t777-729.0/16.0*t3+405.0/8.0*t35-243.0/8.0*t2-1215.0/16.0*t16+2673.0/32.0*t44-t718-4455.0/8.0*t49;

   double t785 = 243.0/16.0*t59;

   double t786 = 405.0/2.0*t61;

   double t787 = 1215.0/4.0*t63;

   double t791 = 1215.0/8.0*t75;

   double t792 = 729.0/8.0*t77;

   double t793 = 2187.0/16.0*t79;

   double t795 = 243.0/8.0*t83;

   double t797 = -t720+t785+t786+t787-3645.0/32.0*t65-t726-10935.0/16.0*t69-t728+24057.0/32.0*t73-t791+t792+t793+8019.0/16.0*t81-t795+18225.0/16.0*t85;

   double t811 = 3.0/2.0+486.0*t6-20.0*eta+9.0/4.0*xi-675.0*t29-15.0*t1-t518+36.0*t5+81.0/4.0*t3-45.0/2.0*t35+27.0/2.0*t2+t522-135.0/2.0*t44+675.0*t46+450.0*t49;

   double t827 = 200.0*t57-30.0*t59-360.0*t61-540.0*t63+405.0/4.0*t65+324.0*t67+1215.0/2.0*t69+3645.0/4.0*t71-1215.0/2.0*t73+300.0*t75-180.0*t77-270.0*t79-405.0*t81+54.0*t83-2025.0/2.0*t85;

   double t829 = 27.0/16.0*xi;

   double t838 = 2835.0/16.0*t46;

   double t840 = 9.0/8.0-t772+t773+t829+2025.0/8.0*t29-45.0/4.0*t1+189.0/16.0*t4-t777+243.0/16.0*t3-135.0/8.0*t35+81.0/8.0*t2-405.0/16.0*t16+567.0/32.0*t44-t838-945.0/8.0*t49;

   double t843 = 10935.0/32.0*t71;

   double t847 = -t720+t785+t786+t787-1215.0/32.0*t65-t726-3645.0/16.0*t69-t843+5103.0/32.0*t73-t791+t792+t793+1701.0/16.0*t81-t795+6075.0/16.0*t85;

   double t856 = t673+t674+2187.0/16.0*t73-t654-t676-t657+t650+t677-t678+t679-t608;

   double t858 = 6561.0/64.0*t6;

   double t861 = 729.0/64.0*t59;

   double t862 = t858-t738-t684-t685+t741-1701.0/128.0*t35+t688+8019.0/128.0*t44+t745+t691+t692+t861;

   double t863 = 5103.0/64.0*t63;

   double t865 = 5103.0/128.0*t75;

   double t866 = 6561.0/128.0*t79;

   double t867 = 729.0/32.0*t83;

   double t869 = -t695-t863-t748+t698+t699+t750-72171.0/256.0*t73+t865-t703-t866-t705-t867-76545.0/256.0*t85;

   double t876 = -729.0/4.0*t6+t829+t712+t713-243.0/32.0*t3+189.0/32.0*t35-t716-405.0/8.0*t44-t838-t719-t720-45.0/2.0*t59;

   double t884 = t723+567.0/4.0*t63+1215.0/16.0*t65-t726-t727-t843+3645.0/16.0*t73-315.0/4.0*t75+t731+405.0/4.0*t79+t733+81.0/2.0*t83+8505.0/32.0*t85;

   double t888 = t858+81.0/64.0*xi-t739+t740-t541+t542-t743+t131+11907.0/256.0*t46+t746+t692+t861;

   double t889 = -t695-t863-t139+t698+t749+t142-t143+t865-t703-t866-t752-t867-t548;

   double t892 = -81.0/100.0+t568+t108-t577-t582+t567-t564+t569+t570-t573-t579+t578;


   values[0] = t48+t87;

   values[1] = t103+t120;

   values[2] = t135+t150;

   values[3] = t167+t184;

   values[4] = t192+t202;

   values[5] = t212;

   values[6] = t230+t247;

   values[7] = t264+t278;

   values[8] = t296+t312;

   values[9] = t330+t347;

   values[10] = t356+t363;

   values[11] = t370;

   values[12] = t381+t393;

   values[13] = t396;

   values[14] = t404+t413;

   values[15] = t424+t435;

   values[16] = t440+t447;

   values[17] = t453+t460;

   values[18] = t467+t473;

   values[19] = t487;

   values[20] = t500+t510;

   values[21] = t525+t537;

   values[22] = t545+t549;

   values[23] = t556+t562;

   values[24] = t576+t590;

   values[25] = t604+t619;

   values[26] = t633+t648;

   values[27] = t655+t660;

   values[28] = t662-t663+2673.0/160.0*t35-t665-2187.0/160.0*t3+t667+729.0/8.0*t44+24057.0/32.0*t85-t670-2187.0/16.0*t65+t680;

   values[29] = t694+t708;

   values[30] = t722+t736;

   values[31] = t747+t754;

   values[32] = t763+t770;

   values[33] = t784+t797;

   values[34] = t811+t827;

   values[35] = t840+t847;

   values[36] = t662-t663+567.0/160.0*t35-t665-729.0/160.0*t3+t667-243.0/8.0*t44+5103.0/32.0*t85-t670+729.0/16.0*t65+t856;

   values[37] = t862+t869;

   values[38] = t876+t884;

   values[39] = t888+t889;

   values[40] = t581+t108-t586-t587+t588-t566;

   values[41] = t892;

 }


 // values of the derivatives in xi-xi  direction

 static void C_Q_UL6S_2D_DeriveXiXi(double xi, double eta, double *values)

 {

   double t1 = xi*xi;

   double t2 = 459.0/70.0*t1;

   double t3 = eta*eta;

   double t5 = 9.0/2.0*xi;

   double t6 = t3*t3;

   double t7 = t6*eta;

   double t10 = t1*xi;

   double t11 = 27.0/2.0*t10;

   double t12 = t3*eta;

   double t14 = 261.0/640.0*t6;

   double t15 = t10*t7;

   double t17 = xi*t7;

   double t19 = t1*t7;

   double t20 = 729.0/64.0*t19;

   double t21 = t1*t6;

   double t22 = 243.0/64.0*t21;

   double t23 = t10*t6;

   double t24 = 1215.0/64.0*t23;

   double t25 = xi*t6;

   double t26 = 567.0/128.0*t25;

   double t27 = -3303.0/2800.0+t2+437.0/640.0*t3+t5+1431.0/640.0*t7+12127.0/11200.0*eta-t11-309.0/128.0*t12-t14+3645.0/64.0*t15-1701.0/128.0*t17-t20+t22-t24+t26;

   double t28 = t10*t12;

   double t30 = t1*t12;

   double t31 = 567.0/64.0*t30;

   double t32 = xi*eta;

   double t33 = 81.0/64.0*t32;

   double t34 = t1*eta;

   double t35 = 6021.0/1120.0*t34;

   double t36 = 1.0-t3;

   double t37 = t36*t6;

   double t38 = 243.0/80.0*t37;

   double t39 = 1.0-eta;

   double t40 = t1*t1;

   double t43 = 3.0/8.0+35.0/8.0*t40-15.0/4.0*t1;

   double t44 = t39*t43;

   double t45 = 81.0/350.0*t44;

   double t46 = t40*t36;

   double t47 = 2187.0/80.0*t46;

   double t48 = xi*t3;

   double t50 = t10*t3;

   double t52 = t1*t3;

   double t53 = 783.0/320.0*t52;

   double t54 = xi*t12;

   double t56 = t10*eta;

   double t61 = 35.0/2.0*t10-15.0/2.0*xi;

   double t62 = t39*xi*t61;

   double t63 = 81.0/175.0*t62;

   double t64 = 1.0-t1;

   double t66 = t64*t36*t1;

   double t67 = 729.0/40.0*t66;

   double t70 = 105.0/2.0*t1-15.0/2.0;

   double t71 = t39*t64*t70;

   double t72 = 81.0/700.0*t71;

   double t73 = -2835.0/64.0*t28+t31+t33-t35-t38+t45-t47-927.0/128.0*t48+1431.0/64.0*t50-t53+1323.0/128.0*t54-81.0/32.0*t56+t63+t67-t72;

   double t75 = 1377.0/35.0*t1;

   double t76 = 243.0/100.0*t3;

   double t77 = 27.0*xi;

   double t78 = 243.0/100.0*t7;

   double t80 = 81.0*t10;

   double t81 = 27.0/10.0*t12;

   double t82 = 81.0/200.0*t6;

   double t83 = 2187.0/8.0*t15;

   double t85 = 729.0/8.0*t23;

   double t86 = 243.0/8.0*t25;

   double t87 = 9909.0/1400.0-t75-t76-t77-t78-1863.0/350.0*eta+t80+t81+t82-t83+729.0/8.0*t17+t85-t86;

   double t90 = 3807.0/140.0*t34;

   double t91 = 729.0/200.0*t37;

   double t92 = 243.0/175.0*t44;

   double t93 = 6561.0/40.0*t46;

   double t96 = 243.0/20.0*t52;

   double t97 = 567.0/8.0*t54;

   double t99 = 486.0/175.0*t62;

   double t100 = 2187.0/20.0*t66;

   double t101 = 243.0/350.0*t71;

   double t102 = 1701.0/8.0*t28-27.0/4.0*t32+t90+t91-t92+t93+351.0/8.0*t48-1053.0/8.0*t50+t96-t97+81.0/4.0*t56-t99-t100+t101;

   double t104 = 1377.0/14.0*t1;

   double t105 = 567.0/128.0*t3;

   double t106 = 135.0/2.0*xi;

   double t107 = 2187.0/128.0*t7;

   double t108 = 25461.0/2240.0*eta;

   double t109 = 405.0/2.0*t10;

   double t110 = 1701.0/128.0*t12;

   double t111 = 729.0/128.0*t6;

   double t112 = 32805.0/64.0*t15;

   double t114 = 6561.0/64.0*t19;

   double t115 = 2187.0/64.0*t21;

   double t116 = 10935.0/64.0*t23;

   double t117 = -9909.0/560.0+t104+t105+t106-t107+t108-t109+t110+t111+t112-24057.0/128.0*t17+t114-t115-t116;

   double t118 = 8019.0/128.0*t25;

   double t119 = 25515.0/64.0*t28;

   double t120 = 5103.0/64.0*t30;

   double t122 = 13527.0/224.0*t34;

   double t123 = 243.0/70.0*t44;

   double t124 = 6561.0/16.0*t46;

   double t127 = 1701.0/64.0*t52;

   double t128 = 18711.0/128.0*t54;

   double t130 = 243.0/35.0*t62;

   double t131 = 2187.0/8.0*t66;

   double t132 = 243.0/140.0*t71;

   double t133 = t118-t119-t120+1269.0/64.0*t32-t122+t123-t124-13851.0/128.0*t48+20655.0/64.0*t50-t127+t128-2025.0/32.0*t56+t130+t131-t132;

   double t136 = 11.0/2.0*t3;

   double t137 = 90.0*xi;

   double t140 = 270.0*t10;

   double t142 = 45.0/4.0*t6;

   double t143 = 3645.0/8.0*t15;

   double t146 = 243.0/4.0*t21;

   double t147 = 1215.0/8.0*t23;

   double t148 = 3303.0/140.0-918.0/7.0*t1-t136-t137+135.0/4.0*t7-502.0/35.0*eta+t140-105.0/4.0*t12-t142-t143+1215.0/8.0*t17-729.0/4.0*t19+t146+t147;

   double t149 = 405.0/8.0*t25;

   double t159 = 945.0/8.0*t54;

   double t164 = -t149+2835.0/8.0*t28+567.0/4.0*t30-135.0/4.0*t32+540.0/7.0*t34-162.0/35.0*t44+2187.0/4.0*t46+1125.0/8.0*t48-3375.0/8.0*t50+135.0/4.0*t52-t159+405.0/4.0*t56-324.0/35.0*t62-729.0/2.0*t66+81.0/35.0*t71;

   double t166 = 10935.0/64.0*t15;

   double t167 = 5103.0/128.0*t17;

   double t168 = 3645.0/64.0*t23;

   double t169 = -9909.0/560.0+t104+t105+t106-t107+t108-t109+t110+t111+t166-t167+t114-t115-t168;

   double t170 = 1701.0/128.0*t25;

   double t171 = 8505.0/64.0*t28;

   double t175 = 3969.0/128.0*t54;

   double t177 = t170-t171-t120+1971.0/64.0*t32-t122+t123-t124-13149.0/128.0*t48+19845.0/64.0*t50-t127+t175-2835.0/32.0*t56+t130+t131-t132;

   double t179 = 3537.0/700.0*eta;

   double t180 = 27.0/2.0*t32;

   double t181 = 81.0/2.0*t56;

   double t182 = 81.0/40.0*t3;

   double t185 = -t179-t180+t90+t181-t182+81.0/2.0*t48+t96-243.0/2.0*t50+9909.0/1400.0-t77-t75+t80-t92-t99+t101+t93-t100;

   double t186 = 1321.0/3200.0*t3;

   double t190 = 441.0/3200.0*t6;

   double t191 = 729.0/64.0*t15;

   double t193 = 243.0/64.0*t23;

   double t194 = 81.0/128.0*t25;

   double t195 = -3303.0/2800.0+t2+t186+t5+1971.0/3200.0*t7+10111.0/11200.0*eta-t11-393.0/640.0*t12-t190-t191-243.0/128.0*t17-t20+t22+t193+t194;

   double t198 = 243.0/400.0*t37;

   double t201 = 189.0/128.0*t54;

   double t202 = 243.0/32.0*t56;

   double t203 = 567.0/64.0*t28+t31+135.0/64.0*t32-t35-t198+t45-t47-873.0/128.0*t48+1269.0/64.0*t50-t53+t201-t202+t63+t67-t72;

   double t205 = 2187.0/100.0*t66;

   double t207 = 81.0/8.0*t54;

   double t209 = 243.0/100.0*t1;

   double t210 = 243.0/10.0*t10;

   double t211 = 81.0/100.0*t3;

   double t213 = 81.0/10.0*xi;

   double t215 = 2187.0/40.0*t15;

   double t216 = -t205+243.0/40.0*t56-t207+t82-1377.0/400.0*t7-t209+t210-t211+153.0/40.0*t12-t213-153.0/400.0*eta+t215;

   double t218 = 2187.0/40.0*t19;

   double t220 = 243.0/4.0*t30;

   double t222 = 243.0/40.0*t34;

   double t223 = 6561.0/200.0*t46;

   double t224 = 81.0/10.0*t48;

   double t225 = 243.0/10.0*t50;

   double t226 = 243.0/100.0*t52;

   double t227 = 81.0/200.0+729.0/80.0*t17+t218-243.0/4.0*t28-t220+81.0/80.0*t32+t222+t91+t223+t224-t225+t226;

   double t229 = 729.0/32.0*t34;

   double t231 = 2187.0/64.0*t50;

   double t233 = 6561.0/64.0*t15;

   double t236 = 729.0/128.0*t25;

   double t237 = 2187.0/64.0*t23;

   double t240 = 243.0/640.0*t6;

   double t241 = 8019.0/64.0*t30;

   double t242 = 2673.0/128.0*t54;

   double t245 = 729.0/80.0*t37;

   double t246 = 243.0/640.0*t3;

   double t247 = 2187.0/64.0*t52;

   double t248 = 729.0/128.0*t48;

   double t249 = -t115+8019.0/64.0*t28-t240+t241+t242+351.0/320.0*eta-243.0/64.0*t32-t245+t246+t247+t248;

   double t251 = 27.0/4.0*t1;

   double t252 = t3/10.0;

   double t259 = 729.0/8.0*t15;

   double t261 = 729.0/8.0*t19;

   double t262 = -1.0/8.0+t251-t252+9.0/8.0*xi-783.0/80.0*t7-87.0/80.0*eta+27.0/4.0*t10+87.0/8.0*t12+9.0/40.0*t6+t259+243.0/16.0*t17+t261;

   double t266 = 405.0/4.0*t30;

   double t268 = 81.0/8.0*t34;

   double t272 = 135.0/2.0*t52;

   double t273 = 135.0/8.0*t54;

   double t275 = t146+243.0/4.0*t23+81.0/8.0*t25-405.0/4.0*t28-t266+27.0/16.0*t32+t268+243.0/20.0*t37-45.0/4.0*t48-135.0/2.0*t50-t272-t273+81.0/8.0*t56;

   double t277 = 243.0/32.0*t34;

   double t279 = 2187.0/64.0*t19;

   double t280 = 2187.0/64.0*t15;

   double t285 = 1701.0/64.0*t30;

   double t286 = 567.0/128.0*t54;

   double t288 = -t115+1701.0/64.0*t28-t240+t285+t286+171.0/320.0*eta+t33-t245+t246+t247+t248;

   double t291 = 81.0/200.0*t3;

   double t293 = 864.0/175.0*t1;

   double t295 = 1.0+eta;

   double t296 = t295*t43;

   double t297 = 81.0/350.0*t296;

   double t298 = 9.0/10.0*xi;

   double t301 = 27.0/10.0*t10;

   double t305 = 729.0/320.0*t19;

   double t306 = t293+457.0/3200.0*t3+t297-t298+1161.0/3200.0*t7-8893.0/11200.0*eta+t301-303.0/640.0*t12-t190+729.0/320.0*t15+243.0/640.0*t17+t305+t22+t193+t194;

   double t308 = 81.0/64.0*t30;

   double t310 = 4887.0/1120.0*t34;

   double t311 = 2187.0/400.0*t46;

   double t314 = 1323.0/1600.0*t52;

   double t317 = 729.0/200.0*t66;

   double t319 = t295*xi*t61;

   double t320 = 81.0/175.0*t319;

   double t322 = t295*t64*t70;

   double t323 = 81.0/700.0*t322;

   double t324 = -2547.0/2800.0+81.0/64.0*t28+t308-729.0/320.0*t32+t310-t198-t311-909.0/640.0*t48+1161.0/320.0*t50-t314+27.0/128.0*t54+1053.0/160.0*t56+t317+t320-t323;

   double t326 = 27.0/5.0*xi;

   double t327 = 5184.0/175.0*t1;

   double t328 = 81.0/5.0*t10;

   double t329 = 243.0/175.0*t296;

   double t330 = 486.0/175.0*t319;

   double t331 = 243.0/350.0*t322;

   double t332 = t179+t180-t90-t181-t291+t224+t226-t225+7641.0/1400.0+t326-t327-t328-t329-t330+t331+t223-t205;

   double t333 = 2592.0/35.0*t1;

   double t334 = 243.0/70.0*t296;

   double t335 = 27.0/2.0*xi;

   double t336 = 2187.0/640.0*t7;

   double t337 = 28863.0/2240.0*eta;

   double t338 = 81.0/2.0*t10;

   double t339 = 243.0/128.0*t12;

   double t341 = 6561.0/320.0*t19;

   double t342 = -7641.0/560.0+t333+t246+t334-t335+t336-t337+t338+t339+t111-t280+5103.0/640.0*t17-t341-t115;

   double t344 = 729.0/64.0*t30;

   double t346 = 77841.0/1120.0*t34;

   double t347 = 6561.0/80.0*t46;

   double t350 = 729.0/320.0*t52;

   double t352 = 2187.0/40.0*t66;

   double t353 = 243.0/35.0*t319;

   double t354 = 243.0/140.0*t322;

   double t355 = -t168+t170-1215.0/64.0*t28-t344-10989.0/320.0*t32+t346-t347-2781.0/128.0*t48+4293.0/64.0*t50-t350+t286+3321.0/32.0*t56+t352+t353-t354;

   double t366 = 2547.0/140.0-3456.0/35.0*t1-t252-162.0/35.0*t296+18.0*xi-27.0/4.0*t7+607.0/35.0*eta-54.0*t10-15.0/4.0*t12-t142+t259-243.0/8.0*t17+729.0/20.0*t19+t146;

   double t379 = t147-t149+405.0/8.0*t28+81.0/4.0*t30+189.0/4.0*t32-3267.0/35.0*t34+2187.0/20.0*t46+261.0/8.0*t48-783.0/8.0*t50+27.0/20.0*t52-t273-567.0/4.0*t56-729.0/10.0*t66-324.0/35.0*t319+81.0/35.0*t322;

   double t382 = -7641.0/560.0+t333+t246+t334-t335+t336-t337+t338+t339+t111-t233+24057.0/640.0*t17-t341-t115;

   double t388 = -t116+t118-3645.0/64.0*t28-t344-11691.0/320.0*t32+t346-t347-3483.0/128.0*t48+5103.0/64.0*t50-t350+t242+3483.0/32.0*t56+t352+t353-t354;

   double t392 = 7641.0/1400.0-t327-t211-t329+t326-t78+837.0/175.0*eta-t328+t81+t82+t215-729.0/40.0*t17+t85;

   double t398 = -t86+243.0/8.0*t28+297.0/20.0*t32-t90+t91+t223+459.0/40.0*t48-1377.0/40.0*t50+t226-t207-891.0/20.0*t56-t205-t330+t331;

   double t404 = t293+t186+t297-t298+1269.0/640.0*t7-6877.0/11200.0*eta+t301-291.0/128.0*t12-t14-t191+1701.0/640.0*t17+t305+t22-t24+t26;

   double t409 = -2547.0/2800.0-405.0/64.0*t28+t308-783.0/320.0*t32+t310-t38-t311-1179.0/640.0*t48+1971.0/320.0*t50-t314+t201+t202+t317+t320-t323;

   double t413 = 81.0/40.0*t6;

   double t415 = 729.0/40.0*t37;

   double t418 = 10935.0/64.0*t50;

   double t421 = 5103.0/128.0*t25;

   double t423 = 567.0/128.0*t6;

   double t426 = 729.0/16.0*t37;

   double t427 = 5103.0/128.0*t48;

   double t428 = -t115-t171-t423+t285+t175+207.0/64.0*eta+567.0/64.0*t32-t426+t105+t247+t427;

   double t437 = -1.0/8.0+t251-t136+63.0/8.0*xi-675.0/16.0*t7-75.0/16.0*eta-135.0/4.0*t10+375.0/8.0*t12+45.0/8.0*t6-t143+1701.0/16.0*t17+t261;

   double t446 = t146-1215.0/4.0*t23+567.0/8.0*t25+2025.0/4.0*t28-t266+189.0/16.0*t32+t268+243.0/4.0*t37-315.0/4.0*t48+675.0/2.0*t50-t272-t159-405.0/8.0*t56;

   double t448 = 3645.0/32.0*t56;

   double t451 = 15309.0/128.0*t17;

   double t455 = 1701.0/64.0*t32;

   double t456 = -t115-40095.0/64.0*t28-t423+t241+t128+243.0/64.0*eta-t455-t426+t105+t247+t427;

   double t462 = -t205-243.0/8.0*t56-t97+t413-1053.0/80.0*t7-t209+t210-t76+117.0/8.0*t12-t213-117.0/80.0*eta-t83;

   double t466 = 81.0/200.0+5103.0/80.0*t17+t218+1215.0/4.0*t28-t220+567.0/80.0*t32+t222+t415+t223+t224-t225+t226;

   double t468 = 19683.0/100.0*t46;

   double t469 = 6561.0/50.0*t66;

   double t471 = 729.0/5.0*t50;

   double t472 = 81.0/5.0*t12;

   double t473 = 243.0/5.0*xi;

   double t474 = 729.0/50.0*t1;

   double t475 = 729.0/5.0*t10;

   double t478 = 729.0/50.0*t7;

   double t481 = 243.0/100.0*t6;

   double t483 = 81.0/50.0*eta;

   double t485 = 2187.0/100.0*t37;

   double t487 = 729.0/50.0*t52;

   double t488 = 243.0/5.0*t48;

   double t489 = t478-2187.0/5.0*t17-1458.0*t28-t481+486.0*t54+t483-243.0/5.0*t32-t485+243.0/50.0*t3-t487-t488;

   double t491 = 6561.0/8.0*t50;

   double t494 = 2187.0/40.0*t37;

   double t496 = 6561.0/8.0*t23;

   double t497 = 729.0/20.0*t7;

   double t498 = 243.0/40.0*t6;

   double t499 = 2187.0/8.0*t25;

   double t501 = 8019.0/8.0*t54;

   double t502 = 81.0/2.0*t12;

   double t503 = 19683.0/8.0*t15;

   double t504 = 81.0/20.0*eta;

   double t505 = 2187.0/8.0*t48;

   double t506 = 243.0/40.0*t3;

   double t507 = t491+729.0/4.0*t32-2187.0/4.0*t56+t494+6561.0/8.0*t17-t496-t497+t498+t499+24057.0/8.0*t28-t501+t502-t503-t504-t505-t506;

   double t509 = 81.0*t32;

   double t510 = 243.0*t56;

   double t514 = 729.0*t17;

   double t519 = 2430.0*t28;

   double t520 = 810.0*t54;

   double t522 = 2187.0*t15;

   double t525 = 81.0/10.0*t3;

   double t526 = -1620.0*t50-t509+t510-54.0*xi+162.0*t10-729.0/10.0*t37-t514+1458.0*t23+243.0/5.0*t7-81.0/10.0*t6-486.0*t25-t519+t520-54.0*t12+t522+27.0/5.0*eta+540.0*t48+t525;

   double t531 = 1701.0/8.0*t54;

   double t532 = 6561.0/8.0*t15;

   double t533 = t491-243.0/4.0*t32+729.0/4.0*t56+t494+2187.0/8.0*t17-t496-t497+t498+t499+5103.0/8.0*t28-t531+t502-t532-t504-t505-t506;

   double t534 = 19683.0/40.0*t46;

   double t535 = 6561.0/20.0*t66;

   double t536 = 2187.0/40.0*t34;

   double t538 = 729.0/2.0*t50;

   double t539 = 729.0/8.0*t12;

   double t540 = 243.0/2.0*xi;

   double t541 = 729.0/20.0*t1;

   double t542 = 729.0/2.0*t10;

   double t543 = 19683.0/40.0*t19;

   double t544 = 243.0/40.0+t534-t535-t536-2187.0/8.0*t56-t538-t539-t540-t541+t542-t543;

   double t545 = 6561.0/80.0*t7;

   double t548 = 2187.0/4.0*t30;

   double t549 = 729.0/80.0*eta;

   double t551 = 729.0/20.0*t52;

   double t552 = 243.0/2.0*t48;

   double t553 = -t503+t545+72171.0/80.0*t17+10935.0/4.0*t28+t548-t501+t549+8019.0/80.0*t32-t506+t551+t552;

   double t555 = 6561.0/32.0*t34;

   double t557 = 98415.0/64.0*t50;

   double t558 = 24057.0/128.0*t12;

   double t559 = 59049.0/64.0*t19;

   double t561 = 19683.0/128.0*t7;

   double t563 = 72171.0/128.0*t25;

   double t564 = 98415.0/64.0*t23;

   double t565 = 19683.0/64.0*t21;

   double t567 = 6561.0/128.0*t6;

   double t568 = 72171.0/64.0*t30;

   double t570 = 2187.0/64.0*eta;

   double t572 = 6561.0/128.0*t3;

   double t573 = 19683.0/64.0*t52;

   double t574 = 72171.0/128.0*t48;

   double t575 = t555+32805.0/32.0*t56-t557+t558+t559+295245.0/64.0*t15-t561-216513.0/128.0*t17-t563+t564+t565-360855.0/64.0*t28-t567-t568+264627.0/128.0*t54-t570-24057.0/64.0*t32+t572-t573+t574;

   double t576 = 729.0/8.0*t34;

   double t579 = 1215.0/8.0*t12;

   double t581 = 243.0/4.0*t1;

   double t583 = 6561.0/8.0*t19;

   double t584 = 32805.0/8.0*t15;

   double t585 = 2187.0/16.0*t7;

   double t587 = 81.0/8.0-t576-3645.0/8.0*t56+6075.0/2.0*t50-t579+891.0/8.0*xi-t581-1215.0/4.0*t10-t583-t584+t585+24057.0/16.0*t17;

   double t590 = 2187.0/4.0*t21;

   double t592 = 729.0/8.0*t6;

   double t593 = 3645.0/4.0*t30;

   double t594 = 13365.0/8.0*t54;

   double t595 = 243.0/16.0*eta;

   double t597 = 405.0/4.0*t3;

   double t598 = 1215.0/2.0*t52;

   double t600 = 8019.0/8.0*t25-10935.0/4.0*t23-t590+18225.0/4.0*t28+t592+t593-t594+t595+2673.0/16.0*t32-t597+t598-4455.0/4.0*t48;

   double t602 = 2187.0/32.0*t34;

   double t603 = 10935.0/32.0*t56;

   double t604 = 5103.0/128.0*t12;

   double t605 = 19683.0/64.0*t19;

   double t606 = 98415.0/64.0*t15;

   double t607 = 6561.0/128.0*t7;

   double t610 = 15309.0/64.0*t30;

   double t611 = 56133.0/128.0*t54;

   double t612 = 729.0/64.0*eta;

   double t614 = -t602-t603-t557+t604+t605+t606-t607-72171.0/128.0*t17-t563+t564+t565-76545.0/64.0*t28-t567-t610+t611+t612+8019.0/64.0*t32+t572-t573+t574;

   double t624 = -81.0/10.0-6561.0/10.0*t46+2187.0/5.0*t66+486.0/5.0*t34+t510+486.0*t50+180.0*t12+162.0*xi+243.0/5.0*t1-486.0*t10+4374.0/5.0*t19;

   double t630 = t522-162.0*t7-t514-t519-972.0*t30+t520-18.0*eta-t509+t525-243.0/5.0*t52-162.0*t48;

   double t634 = 10935.0/8.0*t50;

   double t639 = 3645.0/8.0*t25;

   double t640 = 10935.0/8.0*t23;

   double t642 = 405.0/4.0*t6;

   double t646 = 2187.0/4.0*t52;

   double t647 = 3645.0/8.0*t48;

   double t648 = -729.0/2.0*t34-3645.0/4.0*t56+t634-1485.0/4.0*t12-6561.0/4.0*t19-t584+1215.0/4.0*t7+10935.0/8.0*t17+t639-t640-t590+40095.0/8.0*t28+t642+8019.0/4.0*t30-t594+135.0/2.0*eta+1215.0/4.0*t32-t597+t646-t647;

   double t658 = -20.0+162.0*t34+405.0*t56-2700.0*t50+300.0*t12-t137+108.0*t1+t140+1458.0*t19+3645.0*t15-270.0*t7-1215.0*t17;

   double t671 = -810.0*t25+2430.0*t23+972.0*t21-4050.0*t28-180.0*t6-1620.0*t30+1350.0*t54-30.0*eta-135.0*t32+200.0*t3-1080.0*t52+900.0*t48;

   double t677 = 10935.0/8.0*t15;

   double t682 = 2835.0/8.0*t54;

   double t685 = 243.0/2.0*t34+1215.0/4.0*t56+t634-315.0/4.0*t12-2187.0/4.0*t19-t677+405.0/4.0*t7+3645.0/8.0*t17+t639-t640-t590+8505.0/8.0*t28+t642+1701.0/4.0*t30-t682-45.0/2.0*eta-405.0/4.0*t32-t597+t646-t647;

   double t687 = 243.0/40.0+t534-t535-t536-729.0/8.0*t56-t538-t539-t540-t541+t542-t543;

   double t691 = -t532+t545+15309.0/80.0*t17+3645.0/4.0*t28+t548-t531+t549+1701.0/80.0*t32-t506+t551+t552;

   double t693 = 32805.0/64.0*t50;

   double t695 = 15309.0/128.0*t25;

   double t696 = 32805.0/64.0*t23;

   double t699 = 15309.0/128.0*t48;

   double t700 = t555+t603-t693+t558+t559+t606-t561-45927.0/128.0*t17-t695+t696+t565-120285.0/64.0*t28-t567-t568+t611-t570-5103.0/64.0*t32+t572-t573+t699;

   double t706 = 81.0/8.0-t576-1215.0/8.0*t56+2025.0/2.0*t50-t579+189.0/8.0*xi-t581-405.0/4.0*t10-t583-t677+t585+5103.0/16.0*t17;

   double t712 = 1701.0/8.0*t25-3645.0/4.0*t23-t590+6075.0/4.0*t28+t592+t593-t682+t595+567.0/16.0*t32-t597+t598-945.0/4.0*t48;

   double t715 = -t602-t448-t693+t604+t605+t112-t607-t451-t695+t696+t565-t119-t567-t610+11907.0/128.0*t54+t612+t455+t572-t573+t699;


   values[0] = t27+t73;

   values[1] = t87+t102;

   values[2] = t117+t133;

   values[3] = t148+t164;

   values[4] = t169+t177;

   values[5] = t185;

   values[6] = t195+t203;

   values[7] = t216+t227;

   values[8] = -t229-729.0/32.0*t56+t231-1161.0/128.0*t12-t114-t233+5103.0/640.0*t7-2187.0/128.0*t17-t236-t237+t249;

   values[9] = t262+t275;

   values[10] = t277+t202+t231-927.0/128.0*t12-t279-t280+4293.0/640.0*t7-729.0/128.0*t17-t236-t237+t288;

   values[11] = -27.0/100.0*eta-t291+t81+t82-t78+t91;

   values[12] = t306+t324;

   values[13] = t332;

   values[14] = t342+t355;

   values[15] = t366+t379;

   values[16] = t382+t388;

   values[17] = t392+t398;

   values[18] = t404+t409;

   values[19] = -27.0/20.0*eta-t182+27.0/2.0*t12+t413-243.0/20.0*t7+t415;

   values[20] = t277-1215.0/32.0*t56-t418-4383.0/128.0*t12-t279+t166+3969.0/128.0*t7-t167-t421+t116+t428;

   values[21] = t437+t446;

   values[22] = -t229+t448-t418-4617.0/128.0*t12-t114+t112+4131.0/128.0*t7-t451-t421+t116+t456;

   values[23] = t462+t466;

   values[24] = -243.0/100.0-t468+t469+729.0/5.0*t56+t471-t472+t473+t474-t475+6561.0/5.0*t15+t489;

   values[25] = t507;

   values[26] = t526;

   values[27] = t533;

   values[28] = t544+t553;

   values[29] = t575;

   values[30] = t587+t600;

   values[31] = t614;

   values[32] = t624+t630;

   values[33] = t648;

   values[34] = t658+t671;

   values[35] = t685;

   values[36] = t687+t691;

   values[37] = t700;

   values[38] = t706+t712;

   values[39] = t715;

   values[40] = t76-t488-t487+t471-243.0/100.0+t473+t474-t475-t468+t469;

   values[41] = t483+t76-t472-t481+t478-t485;

 }


 // values of the derivatives in xi-eta direction

 static void C_Q_UL6S_2D_DeriveXiEta(double xi, double eta, double *values)

 {

   double t1 = eta*eta;

   double t2 = t1*eta;

   double t3 = xi*t2;

   double t4 = 261.0/160.0*t3;

   double t5 = xi*xi;

   double t6 = t1*t1;

   double t7 = t5*t6;

   double t9 = t5*t5;

   double t10 = t9*t6;

   double t12 = t5*xi;

   double t13 = t12*t6;

   double t14 = 1215.0/64.0*t13;

   double t15 = t9*t2;

   double t16 = 1215.0/64.0*t15;

   double t17 = xi*t6;

   double t19 = t12*t2;

   double t20 = 81.0/16.0*t19;

   double t21 = t5*t2;

   double t22 = 567.0/64.0*t21;

   double t23 = xi*eta;

   double t25 = t5*eta;

   double t27 = 1.0-t5;

   double t31 = t27*(35.0/2.0*t12-15.0/2.0*xi);

   double t32 = 81.0/700.0*t31;

   double t33 = 1.0-t1;

   double t37 = t33*(35.0/2.0*t2-15.0/2.0*eta);

   double t38 = 81.0/700.0*t37;

   double t42 = eta*(3.0/8.0+35.0/8.0*t6-15.0/4.0*t1);

   double t43 = 81.0/350.0*t42;

   double t47 = xi*(3.0/8.0+35.0/8.0*t9-15.0/4.0*t5);

   double t48 = 81.0/350.0*t47;

   double t49 = t1*xi;

   double t51 = -19.0/320.0-t4-8505.0/256.0*t7+18225.0/256.0*t10-t14-t16+1431.0/128.0*t17+t20+t22+437.0/320.0*t23-927.0/128.0*t25+t32+t38-t43-t48-927.0/128.0*t49;

   double t52 = t1*t12;

   double t53 = 567.0/64.0*t52;

   double t54 = t5*t1;

   double t57 = 81.0/128.0*t1;

   double t58 = 2007.0/1120.0*t2;

   double t59 = 2007.0/1120.0*t12;

   double t61 = 81.0/128.0*t5;

   double t62 = t9*eta;

   double t66 = t12*eta;

   double t67 = 261.0/160.0*t66;

   double t68 = t9*t1;

   double t71 = t9*xi*eta;

   double t72 = 243.0/40.0*t71;

   double t74 = xi*t6*eta;

   double t75 = 243.0/40.0*t74;

   double t77 = xi*t33*t2;

   double t78 = 243.0/20.0*t77;

   double t80 = t27*eta*t12;

   double t81 = 243.0/20.0*t80;

   double t82 = t53+3969.0/256.0*t54-81.0/128.0*t9+t57-t58-t59-81.0/128.0*t6+t61+1431.0/128.0*t62+12127.0/11200.0*xi+12127.0/11200.0*eta-t67-8505.0/256.0*t68+t72+t75-t78-t81;

   double t84 = 81.0/50.0*t3;

   double t86 = 10935.0/32.0*t10;

   double t87 = 729.0/8.0*t15;

   double t88 = 243.0/20.0*t17;

   double t89 = 243.0/4.0*t21;

   double t90 = 243.0/50.0*t23;

   double t92 = 243.0/350.0*t31;

   double t93 = 243.0/175.0*t47;

   double t94 = 81.0/10.0*t49;

   double t95 = 1701.0/16.0*t54;

   double t98 = t84+3645.0/16.0*t7-t86+t87-t88-t89-t90+351.0/8.0*t25-t92+t93+t94-t95+81.0/16.0*t9+567.0/160.0*t1;

   double t99 = 81.0/40.0*t2;

   double t100 = 1269.0/140.0*t12;

   double t106 = 81.0/10.0*t66;

   double t108 = 729.0/20.0*t71;

   double t109 = 729.0/100.0*t74;

   double t110 = 729.0/50.0*t77;

   double t111 = 729.0/10.0*t80;

   double t112 = 9.0/80.0+t99+t100-243.0/32.0*t6-27.0/8.0*t5-1053.0/16.0*t62-1863.0/350.0*xi-117.0/80.0*eta+t106+5103.0/32.0*t68-t108-t109+t110+t111;

   double t114 = 729.0/32.0*t3;

   double t116 = 164025.0/256.0*t10;

   double t117 = 10935.0/64.0*t13;

   double t118 = 10935.0/64.0*t15;

   double t119 = 10935.0/128.0*t17;

   double t120 = 729.0/16.0*t19;

   double t121 = 8019.0/64.0*t21;

   double t122 = 567.0/64.0*t23;

   double t124 = 243.0/140.0*t31;

   double t125 = 243.0/70.0*t47;

   double t126 = 5103.0/128.0*t49;

   double t127 = -9.0/64.0+t114-120285.0/256.0*t7+t116+t117-t118-t119-t120+t121+t122-13851.0/128.0*t25+t124-t125+t126;

   double t128 = 5103.0/64.0*t52;

   double t129 = 56133.0/256.0*t54;

   double t131 = 1701.0/128.0*t1;

   double t132 = 243.0/32.0*t2;

   double t133 = 4509.0/224.0*t12;

   double t134 = 3645.0/128.0*t6;

   double t137 = 25461.0/2240.0*xi;

   double t139 = 567.0/32.0*t66;

   double t140 = 76545.0/256.0*t68;

   double t141 = 729.0/8.0*t71;

   double t142 = 729.0/4.0*t80;

   double t143 = -t128+t129-2025.0/128.0*t9-t131-t132-t133+t134+1269.0/128.0*t5+20655.0/128.0*t62+t137+243.0/64.0*eta-t139-t140+t141-t142;

   double t145 = 45.0*t3;

   double t147 = 18225.0/32.0*t10;

   double t149 = 1215.0/8.0*t15;

   double t151 = 81.0*t19;

   double t152 = 405.0/4.0*t21;

   double t153 = 11.0*t23;

   double t155 = 81.0/35.0*t31;

   double t156 = 162.0/35.0*t47;

   double t158 = 9.0/16.0-t145+6075.0/16.0*t7-t147-1215.0/4.0*t13+t149+675.0/4.0*t17+t151-t152-t153+1125.0/8.0*t25-t155+t156-315.0/4.0*t49;

   double t160 = 2835.0/16.0*t54;

   double t163 = 27.0/8.0*t2;

   double t174 = 567.0/4.0*t52-t160+405.0/16.0*t9+189.0/32.0*t1+t163+180.0/7.0*t12-405.0/32.0*t6-135.0/8.0*t5-3375.0/16.0*t62-502.0/35.0*xi-75.0/16.0*eta+45.0/2.0*t66+8505.0/32.0*t68-243.0/2.0*t71+243.0*t80;

   double t176 = 25515.0/256.0*t7;

   double t177 = 54675.0/256.0*t10;

   double t178 = 3645.0/64.0*t15;

   double t179 = 1701.0/64.0*t21;

   double t181 = t114-t176+t177+t117-t178-t119-t120+t179+t122-13149.0/128.0*t25+t124-t125+t126-t128;

   double t182 = 11907.0/256.0*t54;

   double t185 = 81.0/32.0*t2;

   double t190 = 25515.0/256.0*t68;

   double t191 = -45.0/64.0+t182-2835.0/128.0*t9+567.0/128.0*t1+t185-t133-1215.0/128.0*t6+1971.0/128.0*t5+19845.0/128.0*t62+t137+207.0/64.0*eta-t139-t190+t141-t142;

   double t193 = 3537.0/700.0*xi;

   double t194 = 27.0/4.0*t5;

   double t195 = 81.0/8.0*t9;

   double t197 = 81.0/20.0*t23;

   double t200 = 9.0/40.0-t193-t194+t100+t195-27.0/20.0*eta-t197+81.0/2.0*t25+t106-243.0/4.0*t62+t93-t92-t108+t111;

   double t201 = 441.0/800.0*t3;

   double t203 = 3645.0/256.0*t10;

   double t204 = 243.0/64.0*t15;

   double t206 = 81.0/64.0*t21;

   double t207 = 1321.0/1600.0*t23;

   double t210 = -t201-1215.0/256.0*t7-t203-t14+t204+1971.0/640.0*t17+t20+t206+t207-873.0/128.0*t25+t32-t38+t43-t48-1179.0/640.0*t49+t53;

   double t211 = 567.0/256.0*t54;

   double t212 = 243.0/128.0*t9;

   double t214 = 1629.0/1120.0*t2;

   double t215 = 243.0/128.0*t6;

   double t221 = 243.0/200.0*t74;

   double t222 = 243.0/100.0*t77;

   double t223 = 1.0/320.0+t211-t212-783.0/640.0*t1+t214-t59+t215+135.0/128.0*t5+1269.0/128.0*t62+10111.0/11200.0*xi-6877.0/11200.0*eta-t67+1701.0/256.0*t68+t72+t221-t222-t81;

   double t226 = 2187.0/32.0*t10;

   double t227 = 729.0/8.0*t13;

   double t229 = 81.0/50.0*t23;

   double t230 = 81.0/10.0*t25;

   double t231 = 243.0/350.0*t37;

   double t232 = 243.0/175.0*t42;

   double t234 = 243.0/4.0*t52;

   double t235 = 243.0/16.0*t54;

   double t238 = t84+729.0/32.0*t7+t226+t227-1377.0/80.0*t17-t229+t230+t231-t232+459.0/40.0*t49-t234-t235+243.0/160.0*t9+297.0/40.0*t1;

   double t239 = 1269.0/140.0*t2;

   double t240 = 81.0/40.0*t12;

   double t243 = 243.0/20.0*t62;

   double t246 = 81.0/50.0*t66;

   double t248 = 729.0/100.0*t71;

   double t249 = 729.0/50.0*t80;

   double t250 = -99.0/400.0-t239+t240-891.0/80.0*t6+81.0/160.0*t5-t243-153.0/400.0*xi+837.0/175.0*eta+t246-729.0/16.0*t68-t248-t109+t110+t249;

   double t252 = 243.0/160.0*t3;

   double t254 = 32805.0/256.0*t10;

   double t255 = 2187.0/64.0*t15;

   double t257 = 729.0/64.0*t21;

   double t258 = 243.0/320.0*t23;

   double t259 = 729.0/128.0*t25;

   double t260 = 243.0/140.0*t37;

   double t261 = 243.0/70.0*t42;

   double t263 = 207.0/320.0-t252-10935.0/256.0*t7-t254-t117-t255+5103.0/128.0*t17-t120-t257+t258+t259-t260+t261-3483.0/128.0*t49;

   double t264 = 8019.0/64.0*t52;

   double t265 = 8019.0/256.0*t54;

   double t268 = 25947.0/1120.0*t2;

   double t269 = 243.0/32.0*t12;

   double t272 = 2187.0/128.0*t62;

   double t274 = 28863.0/2240.0*eta;

   double t275 = 729.0/32.0*t66;

   double t277 = 729.0/40.0*t74;

   double t278 = 729.0/20.0*t77;

   double t279 = t264+t265-729.0/128.0*t9-11691.0/640.0*t1+t268-t269+3483.0/128.0*t6-243.0/128.0*t5+t272+351.0/320.0*xi-t274+t275+24057.0/256.0*t68+t277-t278;

   double t283 = 3645.0/32.0*t10;

   double t284 = 1215.0/8.0*t13;

   double t288 = t23/5.0;

   double t290 = 81.0/35.0*t37;

   double t291 = 162.0/35.0*t42;

   double t293 = -63.0/80.0+9.0/10.0*t3+1215.0/32.0*t7+t283+t284+243.0/4.0*t15-783.0/16.0*t17+t151+81.0/4.0*t21-t288-45.0/4.0*t25+t290-t291+261.0/8.0*t49;

   double t294 = 405.0/4.0*t52;

   double t295 = 405.0/16.0*t54;

   double t299 = 27.0/8.0*t12;

   double t305 = 45.0*t66;

   double t309 = -t294-t295+81.0/32.0*t9+189.0/8.0*t1-1089.0/35.0*t2+t299-567.0/16.0*t6+27.0/32.0*t5-135.0/4.0*t62-87.0/80.0*xi+607.0/35.0*eta-t305-1215.0/16.0*t68-243.0/10.0*t74+243.0/5.0*t77;

   double t312 = 10935.0/256.0*t10;

   double t313 = 3645.0/64.0*t13;

   double t316 = 1701.0/64.0*t52;

   double t317 = -t252-3645.0/256.0*t7-t312-t313-t255+4293.0/128.0*t17-t120-t257+t258+t259-t260+t261-2781.0/128.0*t49+t316;

   double t318 = 1701.0/256.0*t54;

   double t320 = 81.0/32.0*t12;

   double t324 = 171.0/320.0+t318+t212-10989.0/640.0*t1+t268+t320+3321.0/128.0*t6+t61+t272+171.0/320.0*xi-t274+t275+5103.0/256.0*t68+t277-t278;

   double t327 = 3537.0/700.0*eta;

   double t328 = 81.0/100.0*t23;

   double t329 = 27.0/4.0*t1;

   double t330 = 81.0/8.0*t6;

   double t331 = -9.0/40.0-27.0/100.0*xi+t327-t328+t329+t94-t239+t84-t330-t88-t232+t231-t109+t110;

   double t334 = 243.0/64.0*t13;

   double t339 = 121.0/1600.0-t201+243.0/256.0*t7+729.0/256.0*t10+t334+t204+1161.0/640.0*t17+t20+t206+457.0/1600.0*t23-909.0/640.0*t25-t32-t38+t43+t48-909.0/640.0*t49;

   double t340 = 81.0/64.0*t52;

   double t344 = 1629.0/1120.0*t12;

   double t350 = 441.0/800.0*t66;

   double t352 = 243.0/200.0*t71;

   double t353 = 243.0/100.0*t80;

   double t354 = t340+81.0/256.0*t54+1053.0/640.0*t9-729.0/640.0*t1+t214+t344+1053.0/640.0*t6-729.0/640.0*t5+1161.0/640.0*t62-8893.0/11200.0*xi-8893.0/11200.0*eta-t350+243.0/256.0*t68+t352+t221-t222-t353;

   double t357 = -9.0/40.0+t193+t194-t100-t195-27.0/100.0*eta-t328+t230+t246-t243-t93+t92-t248+t249;

   double t359 = 2187.0/64.0*t13;

   double t360 = 2187.0/128.0*t17;

   double t362 = 729.0/128.0*t49;

   double t363 = 729.0/64.0*t52;

   double t364 = t114+5103.0/256.0*t7-t312-t359-t178+t360-t120+t179+t258-2781.0/128.0*t25-t124+t125+t362-t363;

   double t366 = 25947.0/1120.0*t12;

   double t369 = 28863.0/2240.0*xi;

   double t371 = 243.0/160.0*t66;

   double t373 = 729.0/40.0*t71;

   double t374 = 729.0/20.0*t80;

   double t375 = 171.0/320.0+t318+3321.0/128.0*t9+t57+t185+t366+t215-10989.0/640.0*t5+4293.0/128.0*t62-t369+171.0/320.0*eta-t371-3645.0/256.0*t68+t373-t374;

   double t382 = -63.0/80.0-t145-1215.0/16.0*t7+t283+243.0/4.0*t13+t149-135.0/4.0*t17+t151-t152-t288+261.0/8.0*t25+t155-t156-45.0/4.0*t49;

   double t396 = 81.0/4.0*t52-t295-567.0/16.0*t9+27.0/32.0*t1+t163-1089.0/35.0*t12+81.0/32.0*t6+189.0/8.0*t5-783.0/16.0*t62+607.0/35.0*xi-87.0/80.0*eta+9.0/10.0*t66+1215.0/32.0*t68-243.0/10.0*t71+243.0/5.0*t80;

   double t400 = 207.0/320.0+t114+24057.0/256.0*t7-t254-t359-t118+t360-t120+t121+t258-3483.0/128.0*t25-t124+t125+t362;

   double t408 = -t363+t265+3483.0/128.0*t9-243.0/128.0*t1-t132+t366-729.0/128.0*t6-11691.0/640.0*t5+5103.0/128.0*t62-t369+351.0/320.0*eta-t371-10935.0/256.0*t68+t373-t374;

   double t414 = t84-729.0/16.0*t7+t226+t87-t88-t89-t229+459.0/40.0*t25+t92-t93+t94-t235-891.0/80.0*t9+81.0/160.0*t1;

   double t421 = -99.0/400.0+t99-t100+243.0/160.0*t6+297.0/40.0*t5-1377.0/80.0*t62+837.0/175.0*xi-153.0/400.0*eta+t246+729.0/32.0*t68-t248-t109+t110+t249;

   double t427 = -t4+1701.0/256.0*t7-t203+t334-t16+1269.0/128.0*t17+t20+t22+t207-1179.0/640.0*t25-t32+t38-t43+t48-873.0/128.0*t49+t340;

   double t434 = 1.0/320.0+t211+t212+135.0/128.0*t1-t58+t344-t215-783.0/640.0*t5+1971.0/640.0*t62-6877.0/11200.0*xi+10111.0/11200.0*eta-t350-1215.0/256.0*t68+t352+t75-t78-t353;

   double t438 = 81.0/10.0*t3;

   double t440 = 729.0/20.0*t74;

   double t441 = 729.0/10.0*t77;

   double t442 = 9.0/40.0-27.0/20.0*xi-t327-t197-t329+81.0/2.0*t49+t239+t438+t330-243.0/4.0*t17+t232-t231-t440+t441;

   double t443 = 567.0/32.0*t3;

   double t445 = 5103.0/64.0*t21;

   double t446 = 5103.0/128.0*t25;

   double t448 = -t443-t176+t177-t313+t118+19845.0/128.0*t17-t120-t445+t122+t446+t260-t261-13149.0/128.0*t49+t316;

   double t451 = 4509.0/224.0*t2;

   double t454 = 10935.0/128.0*t62;

   double t456 = 25461.0/2240.0*eta;

   double t457 = 729.0/8.0*t74;

   double t458 = 729.0/4.0*t77;

   double t459 = -45.0/64.0+t182-1215.0/128.0*t9+1971.0/128.0*t1-t451+t320-2835.0/128.0*t6+567.0/128.0*t5-t454+207.0/64.0*xi+t456+t275-t190+t457-t458;

   double t468 = 9.0/16.0+45.0/2.0*t3+8505.0/32.0*t7-t147+t284-1215.0/4.0*t15-3375.0/16.0*t17+t151+567.0/4.0*t21-t153-315.0/4.0*t25-t290+t291+1125.0/8.0*t49;

   double t480 = -t294-t160-405.0/32.0*t9-135.0/8.0*t1+180.0/7.0*t2+t299+405.0/16.0*t6+189.0/32.0*t5+675.0/4.0*t62-75.0/16.0*xi-502.0/35.0*eta-t305+6075.0/16.0*t68-243.0/2.0*t74+243.0*t77;

   double t482 = 76545.0/256.0*t7;

   double t485 = -9.0/64.0-t443-t482+t116-t117+t118+20655.0/128.0*t17-t120-t445+t122+t446+t260-t261-13851.0/128.0*t49;

   double t486 = 3645.0/128.0*t9;

   double t489 = 1701.0/128.0*t5;

   double t492 = t264+t129+t486+1269.0/128.0*t1-t451-t269-2025.0/128.0*t6-t489-t454+243.0/64.0*xi+t456+t275-120285.0/256.0*t68+t457-t458;

   double t499 = t438+5103.0/32.0*t7-t86+t227-1053.0/16.0*t17-t90+t230-t231+t232+351.0/8.0*t49-t234-t95-243.0/32.0*t9-27.0/8.0*t1;

   double t505 = 9.0/80.0+t239+t240+81.0/16.0*t6+567.0/160.0*t5-t243-117.0/80.0*xi-1863.0/350.0*eta+t246+3645.0/16.0*t68-t248-t440+t441+t249;

   double t507 = 2187.0/25.0*t80;

   double t509 = 729.0/10.0*t62;

   double t510 = 243.0/5.0*t49;

   double t512 = 81.0/50.0*eta;

   double t514 = 2187.0/50.0*t71;

   double t516 = 729.0/10.0*t17;

   double t517 = 81.0/100.0-t507+729.0/20.0*t9+t509-t510+729.0/20.0*t6+t512-243.0/10.0*t1+t514+6561.0/4.0*t10+t516;

   double t520 = 243.0/25.0*t3;

   double t522 = 81.0/50.0*xi;

   double t524 = 2187.0/50.0*t74;

   double t525 = 2187.0/25.0*t77;

   double t527 = 243.0/25.0*t66;

   double t528 = 243.0/5.0*t25;

   double t529 = -2187.0/2.0*t7-2187.0/2.0*t68-t520+729.0*t54+t522-243.0/10.0*t5+t524-t525+243.0/25.0*t23-t527-t528;

   double t531 = 6561.0/16.0*t62;

   double t535 = 2187.0/20.0*t74;

   double t536 = 2187.0/10.0*t77;

   double t539 = 6561.0/8.0*t15;

   double t540 = 729.0/4.0*t17;

   double t541 = 243.0/10.0*t3;

   double t542 = t531+729.0/8.0*t5-2187.0/16.0*t9+8019.0/160.0*t1-t535+t536+32805.0/16.0*t7-2187.0/32.0*t6-t539-t540+t541;

   double t543 = 2187.0/4.0*t21;

   double t545 = 729.0/40.0*t2;

   double t546 = 24057.0/16.0*t54;

   double t547 = 243.0/2.0*t49;

   double t548 = 98415.0/32.0*t10;

   double t549 = 81.0/20.0*xi;

   double t550 = 729.0/80.0*eta;

   double t551 = 2187.0/8.0*t25;

   double t552 = 243.0/20.0*t23;

   double t553 = -243.0/80.0+t543+72171.0/32.0*t68-t545-t546+t547-t548-t549+t550-t551-t552;

   double t556 = 81.0/2.0*t5;

   double t557 = 243.0/4.0*t9;

   double t558 = 81.0/2.0*t1;

   double t561 = 3645.0/2.0*t7;

   double t562 = 243.0/4.0*t6;

   double t566 = -810.0*t62-t556+t557-t558+729.0/5.0*t74-1458.0/5.0*t77-t561+t562+1458.0*t15+243.0*t17-162.0/5.0*t3;

   double t568 = 3645.0/2.0*t68;

   double t570 = 1215.0*t54;

   double t572 = 10935.0/4.0*t10;

   double t576 = 81.0/5.0*t23;

   double t577 = 27.0/20.0-972.0*t21-t568+162.0/5.0*t2+t570-162.0*t49+t572+27.0/5.0*xi-18.0*eta+540.0*t25+t576;

   double t584 = t531-243.0/8.0*t5+729.0/16.0*t9+1701.0/160.0*t1-t535+t536+10935.0/16.0*t7-729.0/32.0*t6-t539-t540+t541;

   double t586 = 5103.0/16.0*t54;

   double t587 = 32805.0/32.0*t10;

   double t588 = 81.0/80.0+t543+15309.0/32.0*t68-t545-t586+t547-t587-t549+t550-t551-t552;

   double t590 = 2187.0/10.0*t80;

   double t591 = 729.0/40.0*t12;

   double t593 = 729.0/4.0*t62;

   double t594 = 2187.0/8.0*t49;

   double t596 = 81.0/20.0*eta;

   double t598 = 2187.0/20.0*t71;

   double t599 = 6561.0/8.0*t13;

   double t600 = -243.0/80.0+t590-t591-2187.0/32.0*t9-t593-t594-2187.0/16.0*t6-t596+729.0/8.0*t1-t598-t599;

   double t601 = 6561.0/16.0*t17;

   double t604 = 2187.0/4.0*t52;

   double t605 = 729.0/80.0*xi;

   double t607 = 243.0/10.0*t66;

   double t608 = 243.0/2.0*t25;

   double t609 = -t548+t601+72171.0/32.0*t7+32805.0/16.0*t68+t604-t546+t605+8019.0/160.0*t5-t552+t607+t608;

   double t611 = 6561.0/32.0*t3;

   double t614 = 98415.0/64.0*t13;

   double t615 = 98415.0/64.0*t15;

   double t616 = 98415.0/128.0*t17;

   double t617 = 6561.0/16.0*t19;

   double t618 = 72171.0/64.0*t21;

   double t619 = 6561.0/64.0*t23;

   double t620 = 72171.0/128.0*t25;

   double t621 = 72171.0/128.0*t49;

   double t622 = 729.0/64.0-t611-1082565.0/256.0*t7+1476225.0/256.0*t10+t614+t615-t616+t617-t618+t619+t620+t621;

   double t623 = 72171.0/64.0*t52;

   double t627 = 2187.0/32.0*t2;

   double t628 = 2187.0/32.0*t12;

   double t631 = 98415.0/128.0*t62;

   double t632 = 2187.0/64.0*xi;

   double t633 = 2187.0/64.0*eta;

   double t634 = 6561.0/32.0*t66;

   double t636 = -t623+793881.0/256.0*t54+32805.0/128.0*t9-24057.0/128.0*t1+t627+t628+32805.0/128.0*t6-24057.0/128.0*t5-t631-t632-t633-t634-1082565.0/256.0*t68;

   double t638 = 729.0/2.0*t3;

   double t640 = 164025.0/32.0*t10;

   double t641 = 10935.0/8.0*t13;

   double t643 = 10935.0/16.0*t17;

   double t644 = 729.0*t19;

   double t646 = 405.0/2.0*t23;

   double t648 = 3645.0/8.0*t49;

   double t649 = -81.0/16.0+t638+120285.0/32.0*t7-t640-t641-10935.0/4.0*t15+t643-t644+8019.0/4.0*t21-t646-4455.0/4.0*t25-t648;

   double t650 = 3645.0/4.0*t52;

   double t651 = 40095.0/16.0*t54;

   double t655 = 243.0/8.0*t12;

   double t659 = 243.0/16.0*xi;

   double t661 = 405.0*t66;

   double t663 = t650-t651-3645.0/32.0*t9+1215.0/8.0*t1-243.0/2.0*t2-t655-3645.0/16.0*t6+2673.0/32.0*t5+6075.0/4.0*t62+t659+135.0/2.0*eta+t661+54675.0/16.0*t68;

   double t666 = 492075.0/256.0*t10;

   double t667 = 32805.0/64.0*t13;

   double t668 = 32805.0/128.0*t17;

   double t669 = 15309.0/128.0*t49;

   double t670 = 15309.0/64.0*t52;

   double t671 = -t611-360855.0/256.0*t7+t666+t667+t615-t668+t617-t618+t619+t620+t669-t670;

   double t672 = 168399.0/256.0*t54;

   double t673 = 10935.0/128.0*t9;

   double t675 = 729.0/32.0*t12;

   double t676 = 10935.0/128.0*t6;

   double t678 = 729.0/64.0*xi;

   double t680 = -243.0/64.0+t672-t673-5103.0/128.0*t1+t627-t675+t676+8019.0/128.0*t5-t631+t678-t633-t634-229635.0/256.0*t68;

   double t689 = 27.0/20.0-1458.0/5.0*t80+162.0/5.0*t12+t557+243.0*t62+540.0*t49+t562+27.0/5.0*eta-t558+729.0/5.0*t71+1458.0*t13;

   double t695 = t572-810.0*t17-t561-t568-972.0*t52+t570-18.0*xi-t556+t576-162.0/5.0*t66-162.0*t25;

   double t697 = 405.0*t3;

   double t700 = 10935.0/8.0*t15;

   double t702 = 3645.0/4.0*t21;

   double t703 = 3645.0/8.0*t25;

   double t705 = -81.0/16.0+t697+54675.0/16.0*t7-t640-10935.0/4.0*t13-t700+6075.0/4.0*t17-t644+t702-t646-t703-4455.0/4.0*t49;

   double t709 = 243.0/8.0*t2;

   double t713 = 10935.0/16.0*t62;

   double t715 = 243.0/16.0*eta;

   double t716 = 729.0/2.0*t66;

   double t718 = 8019.0/4.0*t52-t651-3645.0/16.0*t9+2673.0/32.0*t1-t709-243.0/2.0*t12-3645.0/32.0*t6+1215.0/8.0*t5+t713+135.0/2.0*xi+t715+t716+120285.0/32.0*t68;

   double t731 = 9.0/4.0-720.0*t3-6075.0/2.0*t7+18225.0/4.0*t10+2430.0*t13+2430.0*t15-1350.0*t17+1296.0*t19-1620.0*t21+400.0*t23+900.0*t25+900.0*t49;

   double t745 = -1620.0*t52+2025.0*t54+405.0/4.0*t9-135.0/2.0*t1+54.0*t2+54.0*t12+405.0/4.0*t6-135.0/2.0*t5-1350.0*t62-30.0*xi-30.0*eta-720.0*t66-6075.0/2.0*t68;

   double t748 = 54675.0/32.0*t10;

   double t752 = 27.0/16.0+t697+18225.0/16.0*t7-t748-3645.0/4.0*t13-t700+2025.0/4.0*t17-t644+t702-t646-t703-945.0/4.0*t49;

   double t754 = 8505.0/16.0*t54;

   double t762 = 1701.0/4.0*t52-t754+1215.0/16.0*t9+567.0/32.0*t1-t709+81.0/2.0*t12-1215.0/32.0*t6-405.0/8.0*t5+t713-45.0/2.0*xi+t715+t716+25515.0/32.0*t68;

   double t767 = 81.0/80.0+t590-t591-729.0/32.0*t9-t593-t594+729.0/16.0*t6-t596-243.0/8.0*t1-t598-t599;

   double t771 = -t587+t601+15309.0/32.0*t7+10935.0/16.0*t68+t604-t586+t605+1701.0/160.0*t5-t552+t607+t608;

   double t774 = 32805.0/64.0*t15;

   double t775 = 15309.0/64.0*t21;

   double t776 = 15309.0/128.0*t25;

   double t777 = -t611-229635.0/256.0*t7+t666+t614+t774-t616+t617-t775+t619+t776+t621-t623;

   double t779 = 729.0/32.0*t2;

   double t781 = 32805.0/128.0*t62;

   double t782 = 729.0/64.0*eta;

   double t784 = -243.0/64.0+t672+t673+8019.0/128.0*t1-t779+t628-t676-5103.0/128.0*t5-t781-t632+t782-t634-360855.0/256.0*t68;

   double t790 = 27.0/16.0+t638+25515.0/32.0*t7-t748-t641-3645.0/4.0*t15+t643-t644+1701.0/4.0*t21-t646-945.0/4.0*t25-t648;

   double t799 = t650-t754-1215.0/32.0*t9-405.0/8.0*t1+81.0/2.0*t2-t655+1215.0/16.0*t6+567.0/32.0*t5+2025.0/4.0*t62+t659-45.0/2.0*eta+t661+18225.0/16.0*t68;

   double t801 = 81.0/64.0-t611-t482+t116+t667+t774-t668+t617-t775+t619+t776+t669;

   double t803 = -t670+35721.0/256.0*t54-t486+t131-t779-t675-t134+t489-t781+t678+t782-t634-t140;


   values[0] = t51+t82;

   values[1] = t98+t112;

   values[2] = t127+t143;

   values[3] = t158+t174;

   values[4] = t181+t191;

   values[5] = t200;

   values[6] = t210+t223;

   values[7] = t238+t250;

   values[8] = t263+t279;

   values[9] = t293+t309;

   values[10] = t317+t324;

   values[11] = t331;

   values[12] = t339+t354;

   values[13] = t357;

   values[14] = t364+t375;

   values[15] = t382+t396;

   values[16] = t400+t408;

   values[17] = t414+t421;

   values[18] = t427+t434;

   values[19] = t442;

   values[20] = t448+t459;

   values[21] = t468+t480;

   values[22] = t485+t492;

   values[23] = t499+t505;

   values[24] = t517+t529;

   values[25] = t542+t553;

   values[26] = t566+t577;

   values[27] = t584+t588;

   values[28] = t600+t609;

   values[29] = t622+t636;

   values[30] = t649+t663;

   values[31] = t671+t680;

   values[32] = t689+t695;

   values[33] = t705+t718;

   values[34] = t731+t745;

   values[35] = t752+t762;

   values[36] = t767+t771;

   values[37] = t777+t784;

   values[38] = t790+t799;

   values[39] = t801+t803;

   values[40] = t512+t90-t528-t527+t509+t514-t507;

   values[41] = t522+t90-t510-t520+t516+t524-t525;

 }


 // values of the derivatives in eta-eta direction

 static void C_Q_UL6S_2D_DeriveEtaEta(double xi, double eta, double *values)

 {

   double t1 = xi*xi;

   double t2 = eta*eta;

   double t3 = t1*t2;

   double t4 = 783.0/320.0*t3;

   double t5 = t1*xi;

   double t6 = t5*t2;

   double t7 = 567.0/64.0*t6;

   double t8 = 1.0-t1;

   double t9 = t2*t2;

   double t10 = t8*t9;

   double t11 = 2187.0/80.0*t10;

   double t12 = t1*t1;

   double t13 = t8*t12;

   double t14 = 243.0/80.0*t13;

   double t15 = 1.0-xi;

   double t18 = 3.0/8.0+35.0/8.0*t9-15.0/4.0*t2;

   double t19 = t15*t18;

   double t20 = 81.0/350.0*t19;

   double t21 = t12*xi;

   double t22 = t21*t2;

   double t23 = 729.0/64.0*t22;

   double t24 = t5*eta;

   double t26 = t12*eta;

   double t27 = 567.0/128.0*t26;

   double t28 = t21*eta;

   double t30 = t12*t2;

   double t31 = 243.0/64.0*t30;

   double t32 = t2*eta;

   double t33 = xi*t32;

   double t35 = t12*t32;

   double t36 = 1215.0/64.0*t35;

   double t37 = t21*t32;

   double t39 = t5*t32;

   double t41 = -3303.0/2800.0-t4+t7-t11-t14+t20-t23+1323.0/128.0*t24+t27-1701.0/128.0*t28+t31-81.0/32.0*t33-t36+3645.0/64.0*t37-2835.0/64.0*t39;

   double t42 = t1*t32;

   double t44 = xi*eta;

   double t45 = 81.0/64.0*t44;

   double t46 = t1*eta;

   double t48 = xi*t2;

   double t49 = 6021.0/1120.0*t48;

   double t51 = 9.0/2.0*eta;

   double t52 = 27.0/2.0*t32;

   double t54 = 1.0-t2;

   double t56 = t8*t54*t2;

   double t57 = 729.0/40.0*t56;

   double t60 = 105.0/2.0*t2-15.0/2.0;

   double t61 = t15*t54*t60;

   double t62 = 81.0/700.0*t61;

   double t66 = 35.0/2.0*t32-15.0/2.0*eta;

   double t67 = t15*eta*t66;

   double t68 = 81.0/175.0*t67;

   double t70 = 459.0/70.0*t2;

   double t72 = 261.0/640.0*t12;

   double t73 = 1431.0/64.0*t42+t45-927.0/128.0*t46-t49+12127.0/11200.0*xi+t51-t52+437.0/640.0*t1+t57-t62+t68+1431.0/640.0*t21+t70-309.0/128.0*t5-t72;

   double t75 = 2187.0/40.0*t22;

   double t76 = 567.0/8.0*t24;

   double t79 = 243.0/100.0*t1;

   double t81 = 81.0/40.0*t12;

   double t83 = 243.0/100.0*t2;

   double t84 = 243.0/10.0*t32;

   double t86 = 81.0/10.0*eta;

   double t87 = t75-t76+5103.0/80.0*t28-243.0/8.0*t33-t79+117.0/8.0*t5+t81-1053.0/80.0*t21-t83+t84-117.0/80.0*xi-t86;

   double t88 = 2187.0/8.0*t37;

   double t90 = 243.0/10.0*t42;

   double t91 = 2187.0/100.0*t56;

   double t93 = 81.0/10.0*t46;

   double t94 = 729.0/40.0*t13;

   double t95 = 6561.0/200.0*t10;

   double t96 = 243.0/40.0*t48;

   double t97 = 243.0/4.0*t6;

   double t98 = 243.0/100.0*t3;

   double t99 = 81.0/200.0-t88+1215.0/4.0*t39-t90-t91+567.0/80.0*t44+t93+t94+t95+t96-t97+t98;

   double t101 = 15309.0/128.0*t28;

   double t103 = 567.0/128.0*t12;

   double t105 = 5103.0/128.0*t26;

   double t106 = 729.0/32.0*t48;

   double t107 = 32805.0/64.0*t37;

   double t108 = 10935.0/64.0*t42;

   double t109 = 10935.0/64.0*t35;

   double t112 = 2187.0/64.0*t30;

   double t113 = 6561.0/64.0*t22;

   double t114 = 3645.0/32.0*t33;

   double t115 = 8019.0/64.0*t6;

   double t116 = 2187.0/64.0*t3;

   double t118 = 567.0/128.0*t1;

   double t119 = 729.0/16.0*t13;

   double t120 = 1701.0/64.0*t44;

   double t121 = 18711.0/128.0*t24;

   double t122 = 5103.0/128.0*t46;

   double t123 = -t112-t113+t114+t115+t116+243.0/64.0*xi+t118-t119-t120+t121+t122;

   double t125 = 135.0/2.0*t3;

   double t126 = 405.0/4.0*t6;

   double t128 = 729.0/8.0*t22;

   double t129 = 945.0/8.0*t24;

   double t132 = 243.0/4.0*t30;

   double t135 = 3645.0/8.0*t37;

   double t136 = -1.0/8.0-t125-t126+243.0/4.0*t13+t128-t129+567.0/8.0*t26+1701.0/16.0*t28+t132-405.0/8.0*t33-1215.0/4.0*t35-t135;

   double t141 = 81.0/8.0*t48;

   double t145 = 11.0/2.0*t1;

   double t147 = 27.0/4.0*t2;

   double t150 = 2025.0/4.0*t39+675.0/2.0*t42+189.0/16.0*t44-315.0/4.0*t46+t141-75.0/16.0*xi+63.0/8.0*eta-135.0/4.0*t32-t145-675.0/16.0*t21+t147+375.0/8.0*t5+45.0/8.0*t12;

   double t152 = 5103.0/128.0*t28;

   double t155 = 243.0/32.0*t48;

   double t156 = 10935.0/64.0*t37;

   double t157 = 8505.0/64.0*t39;

   double t159 = 2187.0/64.0*t22;

   double t161 = 1701.0/64.0*t6;

   double t164 = 3969.0/128.0*t24;

   double t165 = -t112-t159-1215.0/32.0*t33+t161+t116+207.0/64.0*xi+t118-t119+567.0/64.0*t44+t164+t122;

   double t168 = 81.0/40.0*t1;

   double t172 = 1.0+xi;

   double t173 = t172*t18;

   double t174 = 81.0/350.0*t173;

   double t176 = t172*eta*t66;

   double t177 = 81.0/175.0*t176;

   double t179 = t172*t54*t60;

   double t180 = 81.0/700.0*t179;

   double t181 = 1323.0/1600.0*t3;

   double t182 = 81.0/64.0*t6;

   double t183 = 2187.0/400.0*t10;

   double t184 = 729.0/320.0*t22;

   double t185 = 189.0/128.0*t24;

   double t187 = 243.0/32.0*t33;

   double t188 = 729.0/64.0*t37;

   double t189 = t174+t177-t180-t181+t182-t183-t14+t184+t185+t27+1701.0/640.0*t28+t31+t187-t36-t188;

   double t194 = 4887.0/1120.0*t48;

   double t196 = 9.0/10.0*eta;

   double t197 = 27.0/10.0*t32;

   double t198 = 1321.0/3200.0*t1;

   double t199 = 729.0/200.0*t56;

   double t201 = 864.0/175.0*t2;

   double t203 = -2547.0/2800.0-405.0/64.0*t39+1971.0/320.0*t42-783.0/320.0*t44-1179.0/640.0*t46+t194-6877.0/11200.0*xi-t196+t197+t198+t199+1269.0/640.0*t21+t201-291.0/128.0*t5-t72;

   double t205 = 243.0/175.0*t173;

   double t206 = 486.0/175.0*t176;

   double t207 = 243.0/350.0*t179;

   double t208 = 729.0/200.0*t13;

   double t209 = 81.0/8.0*t24;

   double t210 = 243.0/8.0*t26;

   double t213 = 729.0/8.0*t35;

   double t214 = 2187.0/40.0*t37;

   double t215 = 7641.0/1400.0-t205-t206+t207+t98+t95+t208-t209-t210-729.0/40.0*t28-891.0/20.0*t33+t213+t214;

   double t220 = 3807.0/140.0*t48;

   double t222 = 27.0/5.0*eta;

   double t223 = 81.0/5.0*t32;

   double t224 = 81.0/100.0*t1;

   double t225 = 243.0/100.0*t21;

   double t226 = 5184.0/175.0*t2;

   double t227 = 27.0/10.0*t5;

   double t228 = 81.0/200.0*t12;

   double t229 = 243.0/8.0*t39-1377.0/40.0*t42+297.0/20.0*t44+459.0/40.0*t46-t220+837.0/175.0*xi+t222-t223-t224-t91-t225-t226+t227+t228;

   double t231 = 243.0/70.0*t173;

   double t232 = 243.0/35.0*t176;

   double t233 = 243.0/140.0*t179;

   double t234 = 729.0/320.0*t3;

   double t235 = 729.0/64.0*t6;

   double t236 = 6561.0/80.0*t10;

   double t237 = 6561.0/320.0*t22;

   double t238 = 2673.0/128.0*t24;

   double t239 = 8019.0/128.0*t26;

   double t242 = -7641.0/560.0+t231+t232-t233-t234-t235-t236-t237+t238+t239+24057.0/640.0*t28-t112+3483.0/32.0*t33-t109;

   double t243 = 6561.0/64.0*t37;

   double t248 = 77841.0/1120.0*t48;

   double t249 = 28863.0/2240.0*xi;

   double t250 = 27.0/2.0*eta;

   double t251 = 81.0/2.0*t32;

   double t252 = 243.0/640.0*t1;

   double t253 = 2187.0/40.0*t56;

   double t254 = 2187.0/640.0*t21;

   double t255 = 2592.0/35.0*t2;

   double t256 = 243.0/128.0*t5;

   double t257 = 729.0/128.0*t12;

   double t258 = -t243-3645.0/64.0*t39+5103.0/64.0*t42-11691.0/320.0*t44-3483.0/128.0*t46+t248-t249-t250+t251+t252+t253+t254+t255+t256+t257;

   double t267 = 135.0/8.0*t24;

   double t268 = 405.0/8.0*t26;

   double t271 = 1215.0/8.0*t35;

   double t272 = 2547.0/140.0-162.0/35.0*t173-324.0/35.0*t176+81.0/35.0*t179+27.0/20.0*t3+81.0/4.0*t6+2187.0/20.0*t10+729.0/20.0*t22-t267-t268-243.0/8.0*t28+t132-567.0/4.0*t33+t271;

   double t273 = 729.0/8.0*t37;

   double t282 = t1/10.0;

   double t287 = 45.0/4.0*t12;

   double t288 = t273+405.0/8.0*t39-783.0/8.0*t42+189.0/4.0*t44+261.0/8.0*t46-3267.0/35.0*t48+607.0/35.0*xi+18.0*eta-54.0*t32-t282-729.0/10.0*t56-27.0/4.0*t21-3456.0/35.0*t2-15.0/4.0*t5-t287;

   double t290 = 567.0/128.0*t24;

   double t291 = 1701.0/128.0*t26;

   double t294 = 3645.0/64.0*t35;

   double t295 = -7641.0/560.0+t231+t232-t233-t234-t235-t236-t237+t290+t291+5103.0/640.0*t28-t112+3321.0/32.0*t33-t294;

   double t296 = 2187.0/64.0*t37;

   double t301 = -t296-1215.0/64.0*t39+4293.0/64.0*t42-10989.0/320.0*t44-2781.0/128.0*t46+t248-t249-t250+t251+t252+t253+t254+t255+t256+t257;

   double t303 = 3537.0/700.0*xi;

   double t304 = 81.0/200.0*t1;

   double t305 = 27.0/2.0*t44;

   double t306 = 81.0/2.0*t33;

   double t307 = t303-t304+t305+t93-t220+t98-t306-t90+7641.0/1400.0+t222-t226-t223-t205-t206+t207+t95-t91;

   double t308 = 243.0/400.0*t13;

   double t310 = 81.0/128.0*t26;

   double t313 = 243.0/64.0*t35;

   double t315 = t174+t177-t180-t181+t182-t183-t308+t184+27.0/128.0*t24+t310+243.0/640.0*t28+t31+1053.0/160.0*t33+t313+729.0/320.0*t37;

   double t324 = 441.0/3200.0*t12;

   double t325 = -2547.0/2800.0+81.0/64.0*t39+1161.0/320.0*t42-729.0/320.0*t44-909.0/640.0*t46+t194-8893.0/11200.0*xi-t196+t197+457.0/3200.0*t1+t199+1161.0/3200.0*t21+t201-303.0/640.0*t5-t324;

   double t331 = 243.0/640.0*t12;

   double t333 = 729.0/128.0*t26;

   double t334 = 2187.0/64.0*t42;

   double t335 = 2187.0/64.0*t35;

   double t339 = 729.0/80.0*t13;

   double t340 = 729.0/128.0*t46;

   double t341 = -t112-t159+t187+t161+t116+171.0/320.0*xi+t252-t339+t45+t290+t340;

   double t348 = -1.0/8.0-t125-t126+243.0/20.0*t13+t128-t267+81.0/8.0*t26+243.0/16.0*t28+t132+81.0/8.0*t33+243.0/4.0*t35+t273;

   double t359 = -405.0/4.0*t39-135.0/2.0*t42+27.0/16.0*t44-45.0/4.0*t46+t141-87.0/80.0*xi+9.0/8.0*eta+27.0/4.0*t32-t282-783.0/80.0*t21+t147+87.0/8.0*t5+9.0/40.0*t12;

   double t369 = -t112-t113-729.0/32.0*t33+t115+t116+351.0/320.0*xi+t252-t339-243.0/64.0*t44+t238+t340;

   double t376 = t75-t209+729.0/80.0*t28+243.0/40.0*t33-t224+153.0/40.0*t5+t228-1377.0/400.0*t21-t83+t84-153.0/400.0*xi-t86;

   double t379 = 81.0/200.0+t214-243.0/4.0*t39-t90-t91+81.0/80.0*t44+t93+t208+t95+t96-t97+t98;

   double t383 = -3303.0/2800.0-t4+t7-t11-t308+t20-t23+t185+t310-243.0/128.0*t28+t31-t187+t313-t188+567.0/64.0*t39;

   double t390 = 1269.0/64.0*t42+135.0/64.0*t44-873.0/128.0*t46-t49+10111.0/11200.0*xi+t51-t52+t198+t57-t62+t68+1971.0/3200.0*t21+t70-393.0/640.0*t5-t324;

   double t393 = 243.0/20.0*t3;

   double t395 = 27.0*eta;

   double t396 = 1377.0/35.0*t2;

   double t397 = 81.0*t32;

   double t398 = 243.0/175.0*t19;

   double t399 = 486.0/175.0*t67;

   double t400 = 243.0/350.0*t61;

   double t401 = 6561.0/40.0*t10;

   double t402 = 2187.0/20.0*t56;

   double t403 = -t303-t168-t305+81.0/2.0*t46+t220+t393+t306-243.0/2.0*t42+9909.0/1400.0-t395-t396+t397-t398-t399+t400+t401-t402;

   double t404 = 1701.0/64.0*t3;

   double t405 = 5103.0/64.0*t6;

   double t406 = 6561.0/16.0*t10;

   double t407 = 243.0/70.0*t19;

   double t409 = -9909.0/560.0-t404-t405-t406+t407+t113+t164+t291-t152-t112-2835.0/32.0*t33-t294+t156-t157;

   double t413 = 13527.0/224.0*t48;

   double t414 = 25461.0/2240.0*xi;

   double t415 = 135.0/2.0*eta;

   double t416 = 405.0/2.0*t32;

   double t417 = 2187.0/8.0*t56;

   double t418 = 243.0/140.0*t61;

   double t419 = 243.0/35.0*t67;

   double t420 = 2187.0/128.0*t21;

   double t421 = 1377.0/14.0*t2;

   double t422 = 1701.0/128.0*t5;

   double t423 = 19845.0/64.0*t42+1971.0/64.0*t44-13149.0/128.0*t46-t413+t414+t415-t416+t118+t417-t418+t419-t420+t421+t422+t257;

   double t433 = 3303.0/140.0+135.0/4.0*t3+567.0/4.0*t6+2187.0/4.0*t10-162.0/35.0*t19-729.0/4.0*t22-t129-t268+1215.0/8.0*t28+t132+405.0/4.0*t33+t271-t135+2835.0/8.0*t39;

   double t439 = 90.0*eta;

   double t440 = 270.0*t32;

   double t447 = -3375.0/8.0*t42-135.0/4.0*t44+1125.0/8.0*t46+540.0/7.0*t48-502.0/35.0*xi-t439+t440-t145-729.0/2.0*t56+81.0/35.0*t61-324.0/35.0*t67+135.0/4.0*t21-918.0/7.0*t2-105.0/4.0*t5-t287;

   double t451 = 25515.0/64.0*t39;

   double t452 = -9909.0/560.0-t404-t405-t406+t407+t113+t121+t239-24057.0/128.0*t28-t112-2025.0/32.0*t33-t109+t107-t451;

   double t456 = 20655.0/64.0*t42+1269.0/64.0*t44-13851.0/128.0*t46-t413+t414+t415-t416+t118+t417-t418+t419-t420+t421+t422+t257;

   double t462 = 9909.0/1400.0+t393+t401+t208-t398-t76-t210+729.0/8.0*t28+81.0/4.0*t33+t213-t88+1701.0/8.0*t39-1053.0/8.0*t42;

   double t466 = -27.0/4.0*t44+351.0/8.0*t46+t220-1863.0/350.0*xi-t395+t397-t79-t402+t400-t399-t225-t396+t227+t228;

   double t468 = 19683.0/100.0*t10;

   double t469 = 6561.0/50.0*t56;

   double t471 = 81.0/5.0*t5;

   double t472 = 243.0/100.0*t12;

   double t473 = 729.0/50.0*t21;

   double t474 = 729.0/5.0*t32;

   double t475 = 243.0/5.0*eta;

   double t476 = 729.0/50.0*t2;

   double t479 = 729.0/5.0*t42;

   double t482 = 729.0/50.0*t3;

   double t483 = 81.0/50.0*xi;

   double t485 = 2187.0/100.0*t13;

   double t488 = 243.0/5.0*t46;

   double t489 = 6561.0/5.0*t37+t479-1458.0*t39+729.0/5.0*t33-t482+t483+243.0/50.0*t1-t485-243.0/5.0*t44+486.0*t24-t488;

   double t491 = 19683.0/40.0*t10;

   double t492 = 6561.0/20.0*t56;

   double t494 = 729.0/8.0*t5;

   double t495 = 6561.0/80.0*t21;

   double t496 = 2187.0/40.0*t48;

   double t497 = 729.0/2.0*t32;

   double t498 = 243.0/2.0*eta;

   double t499 = 729.0/20.0*t2;

   double t500 = 19683.0/8.0*t37;

   double t501 = 243.0/40.0+t491-t492+72171.0/80.0*t28-t494+t495-t496+t497-t498-t499-t500;

   double t502 = 729.0/2.0*t42;

   double t504 = 19683.0/40.0*t22;

   double t506 = 2187.0/4.0*t6;

   double t507 = 729.0/20.0*t3;

   double t508 = 729.0/80.0*xi;

   double t509 = 243.0/40.0*t1;

   double t511 = 8019.0/8.0*t24;

   double t512 = 243.0/2.0*t46;

   double t513 = -t502+10935.0/4.0*t39-t504-2187.0/8.0*t33+t506+t507+t508-t509+8019.0/80.0*t44-t511+t512;

   double t517 = 729.0*t28;

   double t524 = 2187.0*t37;

   double t525 = -81.0/10.0-6561.0/10.0*t10+2187.0/5.0*t56-t517+180.0*t5-162.0*t21+486.0/5.0*t48-486.0*t32+162.0*eta+243.0/5.0*t2+t524;

   double t527 = 2430.0*t39;

   double t529 = 243.0*t33;

   double t533 = 81.0/10.0*t1;

   double t534 = 81.0*t44;

   double t535 = 810.0*t24;

   double t537 = 486.0*t42-t527+4374.0/5.0*t22+t529-972.0*t6-243.0/5.0*t3-18.0*xi+t533-t534+t535-162.0*t46;

   double t540 = 6561.0/8.0*t37;

   double t541 = 243.0/40.0+t491-t492+15309.0/80.0*t28-t494+t495-t496+t497-t498-t499-t540;

   double t545 = 1701.0/8.0*t24;

   double t546 = -t502+3645.0/4.0*t39-t504-729.0/8.0*t33+t506+t507+t508-t509+1701.0/80.0*t44-t545+t512;

   double t548 = 2187.0/40.0*t13;

   double t549 = 729.0/20.0*t21;

   double t550 = 243.0/40.0*t12;

   double t553 = 6561.0/8.0*t35;

   double t555 = 6561.0/8.0*t42;

   double t556 = 2187.0/8.0*t26;

   double t558 = 2187.0/8.0*t46;

   double t559 = 81.0/20.0*xi;

   double t560 = 81.0/2.0*t5;

   double t561 = t548-t549+t550+729.0/4.0*t44+24057.0/8.0*t39-t553-2187.0/4.0*t33+t555+t556+6561.0/8.0*t28-t511-t558-t500-t559+t560-t509;

   double t563 = 24057.0/128.0*t5;

   double t564 = 6561.0/128.0*t12;

   double t565 = 19683.0/128.0*t21;

   double t566 = 72171.0/128.0*t26;

   double t567 = 6561.0/32.0*t48;

   double t569 = 98415.0/64.0*t42;

   double t570 = 98415.0/64.0*t35;

   double t572 = 19683.0/64.0*t30;

   double t573 = 59049.0/64.0*t22;

   double t575 = 72171.0/64.0*t6;

   double t576 = 19683.0/64.0*t3;

   double t577 = 2187.0/64.0*xi;

   double t578 = 6561.0/128.0*t1;

   double t581 = 72171.0/128.0*t46;

   double t582 = -216513.0/128.0*t28+t563-t564-t565-t566+t567+295245.0/64.0*t37-t569+t570-360855.0/64.0*t39+t572+t573+32805.0/32.0*t33-t575-t576-t577+t578-24057.0/64.0*t44+264627.0/128.0*t24+t581;

   double t585 = 405.0/4.0*t12;

   double t587 = 3645.0/8.0*t26;

   double t589 = 32805.0/8.0*t37;

   double t590 = 10935.0/8.0*t42;

   double t591 = 10935.0/8.0*t35;

   double t593 = 2187.0/4.0*t30;

   double t597 = 2187.0/4.0*t3;

   double t599 = 405.0/4.0*t1;

   double t601 = 13365.0/8.0*t24;

   double t602 = 3645.0/8.0*t46;

   double t603 = 10935.0/8.0*t28-1485.0/4.0*t5+t585+1215.0/4.0*t21+t587-729.0/2.0*t48-t589+t590-t591+40095.0/8.0*t39-t593-6561.0/4.0*t22-3645.0/4.0*t33+8019.0/4.0*t6+t597+135.0/2.0*xi-t599+1215.0/4.0*t44-t601-t602;

   double t605 = 15309.0/128.0*t26;

   double t606 = 98415.0/64.0*t37;

   double t607 = 32805.0/64.0*t42;

   double t608 = 32805.0/64.0*t35;

   double t610 = 10935.0/32.0*t33;

   double t612 = 56133.0/128.0*t24;

   double t613 = 15309.0/128.0*t46;

   double t614 = -45927.0/128.0*t28+t563-t564-t565-t605+t567+t606-t607+t608-120285.0/64.0*t39+t572+t573+t610-t575-t576-t577+t578-5103.0/64.0*t44+t612+t613;

   double t626 = -729.0/10.0*t13+243.0/5.0*t21-81.0/10.0*t12-t534+162.0*t32-54.0*eta-t527+1458.0*t35+t529-1620.0*t42-486.0*t26-t517+t535+540.0*t46+t524+27.0/5.0*xi-54.0*t5+t533;

   double t628 = 1215.0/8.0*t5;

   double t629 = 729.0/8.0*t12;

   double t630 = 2187.0/16.0*t21;

   double t632 = 729.0/8.0*t48;

   double t635 = 243.0/4.0*t2;

   double t637 = 81.0/8.0+24057.0/16.0*t28-t628+t629+t630+8019.0/8.0*t26-t632-1215.0/4.0*t32+891.0/8.0*eta-t635-t589+6075.0/2.0*t42;

   double t640 = 6561.0/8.0*t22;

   double t642 = 3645.0/4.0*t6;

   double t643 = 1215.0/2.0*t3;

   double t644 = 243.0/16.0*xi;

   double t647 = -10935.0/4.0*t35+18225.0/4.0*t39-t593-t640-3645.0/8.0*t33+t642+t643+t644-t599+2673.0/16.0*t44-t601-4455.0/4.0*t46;

   double t658 = -20.0-1215.0*t28+300.0*t5-180.0*t12-270.0*t21-810.0*t26+162.0*t48+t440-t439+108.0*t2+3645.0*t37-2700.0*t42;

   double t671 = 2430.0*t35-4050.0*t39+972.0*t30+1458.0*t22+405.0*t33-1620.0*t6-1080.0*t3-30.0*xi+200.0*t1-135.0*t44+1350.0*t24+900.0*t46;

   double t677 = 10935.0/8.0*t37;

   double t679 = 81.0/8.0+5103.0/16.0*t28-t628+t629+t630+1701.0/8.0*t26-t632-405.0/4.0*t32+189.0/8.0*eta-t635-t677+2025.0/2.0*t42;

   double t684 = 2835.0/8.0*t24;

   double t686 = -3645.0/4.0*t35+6075.0/4.0*t39-t593-t640-1215.0/8.0*t33+t642+t643+t644-t599+567.0/16.0*t44-t684-945.0/4.0*t46;

   double t692 = t548-t549+t550-243.0/4.0*t44+5103.0/8.0*t39-t553+729.0/4.0*t33+t555+t556+2187.0/8.0*t28-t545-t558-t540-t559+t560-t509;

   double t694 = 5103.0/128.0*t5;

   double t695 = 6561.0/128.0*t21;

   double t696 = 2187.0/32.0*t48;

   double t698 = 19683.0/64.0*t22;

   double t699 = 15309.0/64.0*t6;

   double t700 = 729.0/64.0*xi;

   double t702 = -72171.0/128.0*t28+t694-t564-t695-t566-t696+t606-t569+t570-76545.0/64.0*t39+t572+t698-t610-t699-t576+t700+t578+8019.0/64.0*t44+t612+t581;

   double t713 = 3645.0/8.0*t28-315.0/4.0*t5+t585+405.0/4.0*t21+t587+243.0/2.0*t48-t677+t590-t591+8505.0/8.0*t39-t593-2187.0/4.0*t22+1215.0/4.0*t33+1701.0/4.0*t6+t597-45.0/2.0*xi-t599-405.0/4.0*t44-t684-t602;

   double t715 = -t101+t694-t564-t695-t605-t696+t107-t607+t608-t451+t572+t698-t114-t699-t576+t700+t578+t120+11907.0/128.0*t24+t613;


   values[0] = t41+t73;

   values[1] = t87+t99;

   values[2] = -t101-4617.0/128.0*t5-t103+4131.0/128.0*t21-t105-t106+t107-t108+t109-40095.0/64.0*t39+t123;

   values[3] = t136+t150;

   values[4] = -t152-4383.0/128.0*t5-t103+3969.0/128.0*t21-t105+t155+t156-t108+t109-t157+t165;

   values[5] = -27.0/20.0*xi-t168+27.0/2.0*t5+t81-243.0/20.0*t21+t94;

   values[6] = t189+t203;

   values[7] = t215+t229;

   values[8] = t242+t258;

   values[9] = t272+t288;

   values[10] = t295+t301;

   values[11] = t307;

   values[12] = t315+t325;

   values[13] = -27.0/100.0*xi-t304+t227+t228-t225+t208;

   values[14] = -729.0/128.0*t28-927.0/128.0*t5-t331+4293.0/640.0*t21-t333+t155-t296+t334-t335+1701.0/64.0*t39+t341;

   values[15] = t348+t359;

   values[16] = -2187.0/128.0*t28-1161.0/128.0*t5-t331+5103.0/640.0*t21-t333-t106-t243+t334-t335+8019.0/64.0*t39+t369;

   values[17] = t376+t379;

   values[18] = t383+t390;

   values[19] = t403;

   values[20] = t409+t423;

   values[21] = t433+t447;

   values[22] = t452+t456;

   values[23] = t462+t466;

   values[24] = -243.0/100.0-t468+t469-2187.0/5.0*t28-t471-t472+t473-t474+t475+t476+t489;

   values[25] = t501+t513;

   values[26] = t525+t537;

   values[27] = t541+t546;

   values[28] = t561;

   values[29] = t582;

   values[30] = t603;

   values[31] = t614;

   values[32] = t626;

   values[33] = t637+t647;

   values[34] = t658+t671;

   values[35] = t679+t686;

   values[36] = t692;

   values[37] = t702;

   values[38] = t713;

   values[39] = t715;

   values[40] = t483+t79-t471-t472+t473-t485;

   values[41] = t79-t488-t482+t479-243.0/100.0+t475+t476-t474-t468+t469;

 }


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


 TBaseFunct2D *BF_C_Q_UL6S_2D_Obj = new TBaseFunct2D

         (42, BF_C_Q_UL6S_2D, BFUnitSquare,

          C_Q_UL6S_2D_Funct, C_Q_UL6S_2D_DeriveXi,

          C_Q_UL6S_2D_DeriveEta, C_Q_UL6S_2D_DeriveXiXi,

          C_Q_UL6S_2D_DeriveXiEta, C_Q_UL6S_2D_DeriveEtaEta, 6, 6,

          0, NULL);

TBaseFunct2D
Definition: BaseFunct2D.h:27