DUNE-FEM (unstable)

orthonormalbase_2d.hh
1// original implementation done by Christian Engwer
2#ifndef DUNE_ORTHONORMALBASE_2D_HH
3#define DUNE_ORTHONORMALBASE_2D_HH
4
5#define PMAX2D 8
6
7#include <stdlib.h>
8#include <stdio.h>
9#include <cassert>
10
11namespace Dune {
12
13 namespace Fem {
14
15 template <class DomainField, class RangeField>
16 class OrthonormalBase_2D
17 {
18 typedef const DomainField* DomainType;
19 typedef RangeField* JacobianRangeType;
20 typedef RangeField HessianRangeType[3];
21
22 public:
23 /* \phi_i(x,y) for triangle */
24 static RangeField
25 eval_triangle_2d ( const int i, DomainType xi )
26 {
27 const RangeField &x = xi[0];
28 const RangeField &y = xi[1];
29
30 switch (i) {
31 #if (PMAX2D>=0)
32 case ( 0 ):
33 return 1.4142135623730950488016887242097;
34 #endif
35 #if (PMAX2D>=1)
36 case ( 1 ):
37 return -2.
38 + 6.*x;
39 case ( 2 ):
40 return -3.4641016151377545870548926830117
41 + 3.4641016151377545870548926830117*x
42 + 6.928203230275509174109785366023*y;
43 #endif
44 #if (PMAX2D>=2)
45 case ( 3 ):
46 return 2.4494897427831780981972840747059
47 + x*(-19.595917942265424785578272597647
48 + 24.494897427831780981972840747059*x);
49 case ( 4 ):
50 return 4.2426406871192851464050661726291
51 - 8.485281374238570292810132345258*y
52 + x*(-25.455844122715710878430397035775
53 + 21.213203435596425732025330863145*x
54 + 42.426406871192851464050661726291*y);
55 case ( 5 ):
56 return 5.477225575051661134569697828008
57 + y*(-32.863353450309966807418186968048
58 + 32.863353450309966807418186968048*y)
59 + x*(-10.954451150103322269139395656016
60 + 5.477225575051661134569697828008*x
61 + 32.863353450309966807418186968048*y);
62 #endif
63 #if (PMAX2D>=3)
64 case ( 6 ):
65 return -2.8284271247461900976033774484194
66 + x*(42.426406871192851464050661726291
67 + x*(-127.27922061357855439215198517887
68 + 98.99494936611665341611821069468*x));
69 case ( 7 ):
70 return -4.8989794855663561963945681494118
71 + 9.797958971132712392789136298824*y
72 + x*(63.686733312362630553129385942353
73 - 117.57550765359254871346963558588*y
74 + x*(-161.66632302368975448102074893059
75 + 102.87856919689348012428593113765*x
76 + 205.75713839378696024857186227529*y));
77 case ( 8 ):
78 return -6.3245553203367586639977870888654
79 + (37.947331922020551983986722533193
80 - 37.947331922020551983986722533193*y)*y
81 + x*(56.920997883030827975980083799789
82 + y*(-303.57865537616441587189378026554
83 + 265.63132345414386388790705773235*y)
84 + x*(-94.868329805051379959966806332982
85 + 44.271887242357310647984509622058*x
86 + 265.63132345414386388790705773235*y));
87 case ( 9 ):
88 return -7.4833147735478827711674974646331
89 + y*(89.799777282574593254009969575597
90 + y*(-224.49944320643648313502492393899
91 + 149.66629547095765542334994929266*y))
92 + x*(22.449944320643648313502492393899
93 + x*(-22.449944320643648313502492393899
94 + 7.483314773547882771167497464633*x
95 + 89.7997772825745932540099695756*y)
96 + y*(-179.59955456514918650801993915119
97 + 224.49944320643648313502492393899*y));
98 #endif
99 #if (PMAX2D>=4)
100 case ( 10 ):
101 return 3.1622776601683793319988935444327
102 + x*(-75.894663844041103967973445066385
103 + x*(398.44698518121579583186058659852
104 + x*(-708.35019587771697036775215395293
105 + 398.44698518121579583186058659852*x)));
106 case ( 11 ):
107 return 5.477225575051661134569697828008
108 - 10.954451150103322269139395656016*y
109 + x*(-120.49896265113654496053335221618
110 + 230.04347415216976765192730877634*y
111 + x*(575.10868538042441912981827194084
112 - 920.1738966086790706077092351053*y
113 + x*(-920.17389660867907060770923510535
114 + 460.08694830433953530385461755267*x
115 + 920.1738966086790706077092351053*y)));
116 case ( 12 ):
117 return 7.0710678118654752440084436210485
118 + y*(-42.426406871192851464050661726291
119 + 42.426406871192851464050661726291*y)
120 + x*(-127.27922061357855439215198517887
121 + (721.24891681027847488886124934695
122 - 678.8225099390856234248105876207*y)*y
123 + x*(487.90367901871779183658260985235
124 + y*(-2206.1731573020282761306344097671
125 + 1527.3506473629426527058238221465*y)
126 + x*(-622.25396744416182147274303865227
127 + 254.55844122715710878430397035775*x
128 + 1527.3506473629426527058238221465*y)));
129 case ( 13 ):
130 return 8.3666002653407554797817202578519
131 + y*(-100.39920318408906575738064309422
132 + (250.99800796022266439345160773556
133 - 167.33200530681510959563440515704*y)*y)
134 + x*(-100.39920318408906575738064309422
135 + y*(1104.3912350249797233311870740364
136 + y*(-2509.9800796022266439345160773556
137 + 1505.9880477613359863607096464133*y))
138 + x*(250.99800796022266439345160773556
139 + x*(-234.26480742954115343388816721985
140 + 75.29940238806679931803548232067*x
141 + 903.592828656801591816425787848*y)
142 + y*(-1907.5848604976922493902322187902
143 + 2258.98207164200397954106446962*y)));
144 case ( 14 ):
145 return 9.4868329805051379959966806332982
146 + y*(-189.73665961010275991993361266596
147 + y*(853.81496824546241963970125699683
148 + y*(-1328.1566172707193194395352886617
149 + 664.0783086353596597197676443309*y)))
150 + x*(-37.947331922020551983986722533193
151 + x*(56.920997883030827975980083799789
152 + x*(-37.947331922020551983986722533193
153 + 9.486832980505137995996680633298*x
154 + 189.73665961010275991993361266596*y)
155 + y*(-569.20997883030827975980083799789
156 + 853.8149682454624196397012569968*y))
157 + y*(569.20997883030827975980083799789
158 + y*(-1707.6299364909248392794025139937
159 + 1328.1566172707193194395352886617*y)));
160 #endif
161 #if (PMAX2D>=5)
162 case ( 15 ):
163 return -3.4641016151377545870548926830117
164 + x*(121.24355652982141054692124390541
165 + x*(-969.94845223857128437536995124329
166 + x*(2909.8453567157138531261098537299
167 + x*(-3637.3066958946423164076373171623
168 + 1600.4149461936426192193604195514*x))));
169 case ( 16 ):
170 return -6.
171 + 12.*y
172 + x*(198.
173 - 384.*y
174 + x*(-1488.
175 + 2592.*y
176 + x*(4176.
177 - 5760.*y
178 + x*(-4860.
179 + 1980.*x
180 + 3960.*y))));
181 case ( 17 ):
182 return -7.7459666924148337703585307995648
183 + (46.475800154489002622151184797389
184 - 46.475800154489002622151184797389*y)*y
185 + x*(224.63303408003017934039739318738
186 + y*(-1301.3224043256920734202331743269
187 + 1254.8466041712030707980819895295*y)
188 + x*(-1471.7336715588184163681208519173
189 + (7529.079625027218424788491937177
190 - 6274.233020856015353990409947647*y)*y
191 + x*(3578.6366118956532019056412293989
192 + y*(-13942.740046346700786645355439217
193 + 7668.507025490685432654945491569*y)
194 + x*(-3601.8745119728977032167168217976
195 + 1278.0845042484475721091575819282*x
196 + 7668.507025490685432654945491569*y))));
197 case ( 18 ):
198 return -9.165151389911680013176094387456
199 + y*(109.98181667894016015811313264947
200 + y*(-274.95454169735040039528283162368
201 + 183.30302779823360026352188774912*y))
202 + x*(210.79848196796864030305017091149
203 + y*(-2419.5999669366835234784889182884
204 + (5774.0453756443584083009394640973
205 - 3666.0605559646720052704377549824*y)*y)
206 + x*(-1081.4878640095782415547791377198
207 + y*(10558.254401178255375178860734349
208 + y*(-20621.590627301280029646212371776
209 + 10081.666528902848014493703826202*y))
210 + x*(2071.3242141200396829777973315651
211 + x*(-1695.5530071336608024375774616794
212 + 504.0833264451424007246851913101*x
213 + 6048.999917341708808696222295721*y)
214 + y*(-14297.636168262220820554707244431
215 + 15122.499793354272021740555739302*y))));
216 case ( 19 ):
217 return -10.392304845413263761164678049035
218 + y*(207.8460969082652752232935609807
219 + y*(-935.30743608719373850482102441317
220 + (1454.9226783578569265630549268649
221 - 727.4613391789284632815274634325*y)*y))
222 + x*(155.88457268119895641747017073553
223 + y*(-2909.8453567157138531261098537299
224 + y*(12158.996669133518600562673317371
225 + y*(-17459.072140294283118756659122379
226 + 8002.074730968213096096802097757*y)))
227 + x*(-519.61524227066318805823390245176
228 + x*(727.46133917892846328152746343247
229 + x*(-467.65371804359686925241051220659
230 + 114.31535329954590137281145853939*x
231 + 2286.3070659909180274562291707878*y)
232 + y*(-7066.767294881019357591981073344
233 + 10288.381796959131123553031268545*y))
234 + y*(7482.4594886975499080385681953054
235 + y*(-21512.071030005455985610883561503
236 + 16004.149461936426192193604195514*y))));
237 case ( 20 ):
238 return -11.489125293076057319701222936438
239 + y*(344.67375879228171959103668809314
240 + y*(-2412.716311545972037137256816652
241 + y*(6433.9101641225920990326848444052
242 + y*(-7238.1489346379161114117704499559
243 + 2895.2595738551664445647081799823*y))))
244 + x*(57.445626465380286598506114682189
245 + x*(-114.89125293076057319701222936438
246 + x*(114.89125293076057319701222936438
247 + x*(-57.445626465380286598506114682189
248 + 11.489125293076057319701222936438*x
249 + 344.67375879228171959103668809314*y)
250 + y*(-1378.6950351691268783641467523725
251 + 2412.716311545972037137256816652*y))
252 + y*(2068.0425527536903175462201285588
253 + y*(-7238.1489346379161114117704499559
254 + 6433.910164122592099032684844405*y)))
255 + y*(-1378.6950351691268783641467523725
256 + y*(7238.1489346379161114117704499559
257 + y*(-12867.82032824518419806536968881
258 + 7238.148934637916111411770449956*y))));
259 #endif
260 #if (PMAX2D>=6)
261 case ( 21 ):
262 return 3.7416573867739413855837487323165
263 + x*(-179.59955456514918650801993915119
264 + x*(2020.4949888579283482152243154509
265 + x*(-8979.9777282574593254009969575597
266 + x*(18521.204064531009858639556224967
267 + x*(-17780.355901949769464293973975968
268 + 6420.684075704083417661712824655*x)))));
269 case ( 22 ):
270 return 6.480740698407860230965967436088
271 - 12.961481396815720461931934872176*y
272 + x*(-298.11407212676157062443450206005
273 + 583.2666628567074207869370692479*y
274 + x*(3207.9666457118908143281538808636
275 - 5832.666628567074207869370692479*y
276 + x*(-13609.555466656506485028531615785
277 + 21386.44430474593876218769253909*y
278 + x*(26733.055380932423452734615673863
279 - 32079.666457118908143281538808636*y
280 + x*(-24380.546507410370188893969494563
281 + 8340.713278850916117253200090245*x
282 + 16681.426557701832234506400180491*y)))));
283 case ( 23 ):
284 return 8.3666002653407554797817202578519
285 + y*(-50.199601592044532878690321547111
286 + 50.19960159204453287869032154711*y)
287 + x*(-351.39721114431173015083225082978
288 + (2058.1836652738258480263031834316
289 - 2007.9840636817813151476128618844*y)*y
290 + x*(3438.6727090550505021902870259771
291 + y*(-18573.852589056477165115418972431
292 + 16565.868525374695849967806110547*y)
293 + x*(-13219.228419238393658055118007406
294 + (60741.517926373884783215289072005
295 - 44175.649400999188933247482961458*y)*y
296 + x*(23468.313744280819120787725323275
297 + y*(-80068.364539311029941511062867642
298 + 35892.715138311841008263579906185*y)
299 + x*(-19326.846612937145158295773795638
300 + 5982.119189718640168043929984364*x
301 + 35892.715138311841008263579906185*y)))));
302 case ( 24 ):
303 return 9.8994949366116653416118210694679
304 + y*(-118.79393923933998409934185283361
305 + (296.98484809834996024835463208404
306 - 197.98989873223330683223642138936*y)*y)
307 + x*(-356.38181771801995229802555850084
308 + y*(4157.7878733768994434769648491765
309 + y*(-10097.484835343898648444057490857
310 + 6533.666658163699125463801905849*y))
311 + x*(2969.8484809834996024835463208404
312 + y*(-31480.393898425095786325591000908
313 + (68603.499910718840817369920011412
314 - 39201.999948982194752782811435093*y)*y)
315 + x*(-9701.5050378794320347795846480785
316 + y*(84937.666556128088631029424776034
317 + y*(-143740.66647960138076020364192867
318 + 56625.11103741872575401961651736*y))
319 + x*(14700.74998086832303229355428816
320 + x*(-10453.866653061918600742083049358
321 + 2831.2555518709362877009808258678*x
322 + 33975.066622451235452411769910414*y)
323 + y*(-91471.333214291787756493226681883
324 + 84937.66655612808863102942477603*y)))));
325 case ( 25 ):
326 return 11.22497216032182415675124619695
327 + y*(-224.49944320643648313502492393899
328 + y*(1010.2474944289641741076121577255
329 + y*(-1571.496102445055381945174467573
330 + 785.7480512225276909725872337865*y)))
331 + x*(-314.29922048901107638903489351459
332 + y*(6061.4849665737850446456729463528
333 + y*(-26266.434855153068526797916100862
334 + (39287.402561126384548629361689324
335 - 18857.953229340664583342093610875*y)*y))
336 + x*(2020.4949888579283482152243154509
337 + y*(-34348.414810584781919658813362666
338 + y*(128301.43179247845011166674403113
339 + y*(-160292.60244939564895840779569244
340 + 61288.34799535715989586180423535*y)))
341 + x*(-5163.4871937480391121055732505968
342 + x*(6342.1092705818306485644541012766
343 + x*(-3771.5906458681329166684187221751
344 + 875.5478285051022842265972033621*x
345 + 17510.956570102045684531944067241*y)
346 + y*(-57920.85634726061264883643037626
347 + 78799.30456545920558039374830259*y))
348 + y*(68921.329064376000322452651649271
349 + y*(-181844.54899721355133937018839058
350 + 122576.69599071431979172360847069*y)))));
351 case ( 26 ):
352 return 12.409673645990856596133241955449
353 + y*(-372.29020937972569788399725866348
354 + y*(2606.0314656580798851879808106444
355 + y*(-6949.4172417548796938346154950517
356 + (7818.0943969742396555639424319332
357 - 3127.2377587896958622255769727733*y)*y)))
358 + x*(-223.37412562783541873039835519809
359 + y*(6328.9335594553368640279533972792
360 + y*(-41696.50345052927816300769297031
361 + y*(104241.25862632319540751923242578
362 + y*(-109453.32155763935517789519404706
363 + 40654.090864266046208932500646052*y))))
364 + x*(930.72552344931424470999314665871
365 + x*(-1737.3543104387199234586538737629
366 + x*(1675.3059422087656404779876639857
367 + x*(-819.03846063539653534479396905967
368 + 161.32575739788113574973214542084*x
369 + 4839.7727219364340724919643626253*y)
370 + y*(-19731.381097125461987851854709165
371 + 33878.409053555038507443750538377*y))
372 + y*(30527.797169137507226487775210406
373 + y*(-104241.25862632319540751923242578
374 + 90342.42414281343601985000143567*y)))
375 + y*(-21592.832144024090477271841002482
376 + y*(109453.32155763935517789519404706
377 + y*(-187634.2655273817517335346183664
378 + 101635.22716066511552233125161513*y)))));
379 case ( 27 ):
380 return 13.490737563232041465550305611496
381 + y*(-566.61097765574574155311283568281
382 + y*(5666.1097765574574155311283568281
383 + y*(-22664.439106229829662124513427313
384 + y*(42495.823324180930616483462676211
385 + y*(-37396.324525279218942505447155066
386 + 12465.441508426406314168482385022*y)))))
387 + x*(-80.944425379392248793301833668973
388 + y*(2833.0548882787287077655641784141
389 + y*(-22664.439106229829662124513427313
390 + y*(67993.317318689488986373540281938
391 + y*(-84991.646648361861232966925352422
392 + 37396.324525279218942505447155066*y))))
393 + x*(202.36106344848062198325458417243
394 + x*(-269.81475126464082931100611222991
395 + x*(202.36106344848062198325458417243
396 + x*(-80.944425379392248793301833668973
397 + 13.490737563232041465550305611496*x
398 + 566.6109776557457415531128356828*y)
399 + y*(-2833.0548882787287077655641784141
400 + 5666.109776557457415531128356828*y))
401 + y*(5666.1097765574574155311283568281
402 + y*(-22664.439106229829662124513427313
403 + 22664.439106229829662124513427313*y)))
404 + y*(-5666.1097765574574155311283568281
405 + y*(33996.658659344744493186770140969
406 + y*(-67993.317318689488986373540281938
407 + 42495.823324180930616483462676211*y)))));
408 #endif
409 #if (PMAX2D>=7)
410 case ( 28 ):
411 return -4.
412 + x*(252.
413 + x*(-3780.
414 + x*(23100.
415 + x*(-69300.
416 + x*(108108.
417 + x*(-84084.
418 + 25740.*x))))));
419 case ( 29 ):
420 return -6.9282032302755091741097853660235
421 + 13.856406460551018348219570732047*y
422 + x*(422.62039704680605962069690732743
423 - 831.3843876330611008931742439228*y
424 + x*(-6131.4598587938256190871600489308
425 + 11431.535329954590137281145853939*y
426 + x*(36199.861878189535434723628537473
427 - 60968.18842642448073216611122101*y
428 + x*(-104789.07385791707625841050366111
429 + 148609.9592894096717846548961012*y
430 + x*(157526.55684677425209173418986728
431 - 166443.15440413883239881348363335*y
432 + x*(-117897.23436959833961582621757362
433 + 34675.657167528923416419475756948*x
434 + 69351.3143350578468328389515139*y))))));
435 case ( 30 ):
436 return -8.9442719099991587856366946749251
437 + (53.665631459994952713820168049551
438 - 53.66563145999495271382016804955*y)*y
439 + x*(509.82349886995205078129159647073
440 + y*(-3005.2753617597173519739294107748
441 + 2951.6097302997223992601092427253*y)
442 + x*(-6896.0336426093514237258915943673
443 + (38370.926493896391190381420155429
444 - 35419.316763596668791121310912703*y)*y
445 + x*(37878.991538846437457171401948308
446 + y*(-188903.02273918223355264699153442
447 + 153483.70597558556476152568062171*y)
448 + x*(-101830.53569534042277447376887402
449 + (422080.19143286030309419562170972
450 - 268596.485457274738332669941088*y)*y
451 + x*(141972.42802741664740441125457509
452 + y*(-429754.3767316395813322719057408
453 + 161157.8912743648429996019646528*y)
454 + x*(-98485.3780010007373886456450656
455 + 26859.6485457274738332669941088*x
456 + 161157.8912743648429996019646528*y))))));
457 case ( 31 ):
458 return -10.583005244258362362006463014557
459 + y*(126.99606293110034834407755617468
460 + y*(-317.49015732775087086019389043671
461 + 211.66010488516724724012926029114*y))
462 + x*(539.73326745717648046232961374241
463 + y*(-6349.8031465550174172038778087342
464 + (15557.017709059792672149500631399
465 - 10159.685034488027867526204493975*y)*y)
466 + x*(-6508.5482252188928526339747539526
467 + y*(71752.775556071696814403819238697
468 + y*(-163824.92118111944936386004746534
469 + 99056.92908625827170838049381625*y))
470 + x*(31801.930758996378897829421358744
471 + y*(-309870.39355188484995954923706623
472 + (610851.06269859267553501304520023
473 - 308177.1127128035119816282029839*y)*y)
474 + x*(-76038.892679996333571016436759592
475 + y*(602596.31860807115289264800404888
476 + y*(-895639.73382158520669660696492196
477 + 288916.04316825329248277644029741*y))
478 + x*(94516.819836471434255079721183009
479 + x*(-58746.262110878169471497876193806
480 + 14445.80215841266462413882201487*x
481 + 173349.62590095197548966586417844*y)
482 + y*(-531605.51942958605816830865014723
483 + 433374.06475237993872416466044611*y))))));
484 case ( 32 ):
485 return -12.
486 + y*(240.
487 + y*(-1080.
488 + (1680.
489 - 840.*y)*y))
490 + x*(516.
491 + y*(-10080.
492 + y*(44280.
493 + y*(-67200.
494 + 32760.*y)))
495 + x*(-5220.
496 + y*(94320.
497 + y*(-380160.
498 + (524160.
499 - 229320.*y)*y))
500 + x*(21420.
501 + y*(-334080.
502 + y*(1.1232e6
503 + y*(-1.22304e6
504 + 382200.*y)))
505 + x*(-43380.
506 + x*(46332.
507 + x*(-25116.
508 + 5460.*x
509 + 109200.*y)
510 + y*(-393120.
511 + 491400.*y))
512 + y*(533520.
513 + y*(-1.27764e6
514 + 764400.*y))))));
515 case ( 33 ):
516 return -13.266499161421599396459730946683
517 + y*(397.99497484264798189379192840048
518 + y*(-2785.9648238985358732565434988034
519 + y*(7429.2395303960956620174493301423
520 + y*(-8357.8944716956076197696304964101
521 + 3343.1577886782430479078521985641*y))))
522 + x*(437.79447232691278008317112124053
523 + y*(-12735.839194964735420601341708815
524 + y*(86364.909540854612070952848462905
525 + y*(-222877.18591188286986052347990427
526 + (242378.93967917262097331928439589
527 - 93608.41808299080534141986155979*y)*y)))
528 + x*(-3382.9572861625078460972313914041
529 + y*(88752.879389910499962315600033308
530 + y*(-534905.24618851888766525635177025
531 + y*(1.2035368039241674972468267914831e6
532 + y*(-1.1115999647355158134293608560225e6
533 + 351031.56781121552003032448084923*y))))
534 + x*(10812.196816558603508114680721546
535 + x*(-17710.776380497835194273740813821
536 + x*(15800.400501253124881183539557499
537 + x*(-7336.3740362661444662422312135156
538 + 1392.9824119492679366282717494017*x
539 + 41789.472358478038098848152482051*y)
540 + y*(-178301.74872950629588841878392342
541 + 292526.30650934626669193706737435*y))
542 + y*(295710.26630808745054708740280156
543 + y*(-955585.93459719780452699442008956
544 + 780070.1506915900445118321796649*y)))
545 + y*(-235613.02510684760528112482161309
546 + y*(1.1143859295594143493026173995214e6
547 + y*(-1.7681590082342707675601529405739e6
548 + 877578.9195280388000758112021231*y))))));
549 case ( 34 ):
550 return -14.422205101855957172476885069882
551 + y*(605.73261427795020124402917293504
552 + y*(-6057.3261427795020124402917293504
553 + y*(24229.304571118008049761166917402
554 + y*(-45429.946070846265093302187970128
555 + (39978.352542344713282105925413713
556 - 13326.117514114904427368641804571*y)*y))))
557 + x*(302.86630713897510062201458646752
558 + y*(-12114.652285559004024880583458701
559 + y*(115089.19671281053823636554285766
560 + y*(-436127.48228012414489570100451323
561 + y*(772309.08320438650658613719549218
562 + y*(-639653.64067751541251369480661941
563 + 199891.76271172356641052962706856*y)))))
564 + x*(-1514.3315356948755031100729323376
565 + y*(51487.272213625767105742479699479
566 + y*(-399783.52542344713282105925413713
567 + y*(1.1630066194136643863885360120353e6
568 + y*(-1.408328328196234217892367827074e6
569 + 599675.2881351706992315888812057*y))))
570 + x*(3533.4402499547095072568368421211
571 + x*(-4542.9946070846265093302187970128
572 + x*(3331.5293785287261068421604511427
573 + x*(-1312.4206642688921026953965413593
574 + 216.33307652783935758715327604823*x
575 + 9085.989214169253018660437594026*y)
576 + y*(-46035.678685124215294546217143063
577 + 90859.89214169253018660437594026*y))
578 + y*(93888.555213082281192824521804932
579 + y*(-369496.89470954962275885779549038
580 + 363439.56856677012074641750376103*y)))
581 + y*(-96917.218284472032199044667669607
582 + y*(569388.65742127318916938742255894
583 + y*(-1.1145480102714283702890136782005e6
584 + 681449.1910626939763995328195519*y))))));
585 case ( 35 ):
586 return -15.49193338482966754071706159913
587 + y*(867.54826955046138228015544955126
588 + y*(-11711.901638931228660782098568942
589 + y*(65066.120216284603671011658716344
590 + y*(-178931.83059478266009528206146995
591 + y*(257661.83605648703053720616851672
592 + y*(-186089.10381857396649909334392874
593 + 53168.31537673541899974095540821*y))))))
594 + x*(108.44353369380767278501943119391
595 + y*(-5205.2896173027682936809326973075
596 + y*(58559.50819465614330391049284471
597 + y*(-260264.48086513841468404663486538
598 + y*(536795.49178434798028584618440984
599 + y*(-515323.67211297406107441233703345
600 + 186089.10381857396649909334392874*y)))))
601 + x*(-325.33060108142301835505829358172
602 + x*(542.21766846903836392509715596954
603 + x*(-542.21766846903836392509715596954
604 + x*(325.33060108142301835505829358172
605 + x*(-108.44353369380767278501943119391
606 + 15.49193338482966754071706159913*x
607 + 867.5482695504613822801554495513*y)
608 + y*(-5205.2896173027682936809326973075
609 + 11711.901638931228660782098568942*y))
610 + y*(13013.224043256920734202331743269
611 + y*(-58559.50819465614330391049284471
612 + 65066.12021628460367101165871634*y)))
613 + y*(-17350.965391009227645603108991025
614 + y*(117119.01638931228660782098568942
615 + y*(-260264.48086513841468404663486538
616 + 178931.83059478266009528206146995*y))))
617 + y*(13013.224043256920734202331743269
618 + y*(-117119.01638931228660782098568942
619 + y*(390396.72129770762202606995229807
620 + y*(-536795.49178434798028584618440984
621 + 257661.83605648703053720616851672*y))))));
622 #endif
623 #if (PMAX2D>=8)
624 case ( 36 ):
625 return 4.2426406871192851464050661726291
626 + x*(-339.41125496954281171240529381033
627 + x*(6533.6666581636991254638019058488
628 + x*(-52269.33326530959300371041524679
629 + x*(212344.16639032022157757356194009
630 + x*(-475650.93271431729633376477874579
631 + x*(594563.66589289662041720597343224
632 + x*(-388286.47568515697659899165611901
633 + 103138.59510386982190910715865661*x)))))));
634 case ( 37 ):
635 return 7.3484692283495342945918522241177
636 - 14.696938456699068589183704448235*y
637 + x*(-573.18059981126367497816447348118
638 + 1131.6642611658282813671452425141*y
639 + x*(10750.810481075368672987879803884
640 - 20369.956700984909064608614365254*y
641 + x*(-83743.155326271292821168747946045
642 + 147116.35395155767657772888152684*y
643 + x*(331011.79639100477229988998343538
644 - 514907.2388304518680220510853439*y
645 + x*(-720870.1343626326152308715194815
646 + 926833.0298948133624396919536191*y
647 + x*(875342.30601176817563748684508467
648 - 823851.5821287229888352817365503*y
649 + x*(-554838.82061730323737886321032978
650 + 142913.02955294174296122234205464*x
651 + 285826.05910588348592244468410928*y)))))));
652 case ( 38 ):
653 return 9.4868329805051379959966806332982
654 + y*(-56.920997883030827975980083799789
655 + 56.92099788303082797598008379979*y)
656 + x*(-702.02564055738021170375436686406
657 + (4155.2328454612504422465461173846
658 - 4098.3118475782196142705660335848*y)*y
659 + x*(12475.185369364256464735635032787
660 + y*(-70695.879370724288346167264079338
661 + 66597.56752314606873189669804575*y)
662 + x*(-91946.385247055797457199828697926
663 + (480982.43211161049639703170810822
664 - 414384.86458846442766513501006246*y)*y
665 + x*(343470.78805918852114505982232856
666 + y*(-1.5798422962435206304733272258631e6
667 + 1.1654574316550562028081922158007e6*y)
668 + x*(-706180.87340284146214600091298145
669 + (2.6572429441735281424026782520255e6
670 - 1.4917855125184719395944860362249e6*y)*y
671 + x*(808913.78774873160150464896755943
672 + y*(-2.1962397823188614666252155533311e6
673 + 704454.2698003895270307295171062*y)
674 + x*(-483449.00868654183227599084507287
675 + 117409.0449667315878384549195177*x
676 + 704454.2698003895270307295171062*y)))))));
677 case ( 39 ):
678 return 11.22497216032182415675124619695
679 + y*(-134.6996659238618898810149543634
680 + (336.74916480965472470253738590849
681 - 224.49944320643648313502492393899*y)*y)
682 + x*(-763.29810690188404265908474139258
683 + y*(9024.8776168987466220280019423475
684 + y*(-22225.44487743721183036746746996
685 + 14592.463808418371403776620056035*y))
686 + x*(12437.26915363658116568038078622
687 + y*(-140222.3522267402273661365674923
688 + (328330.43568941335658497395126078
689 - 204294.49331785719965287268078448*y)*y)
690 + x*(-83917.891870565957395872316568396
691 + y*(866792.35022005126138433123132845
692 + y*(-1.8386504398607147968758541270604e6
693 + 1.0214724665892859982643634039224e6*y))
694 + x*(286741.91383542099808421058410108
695 + y*(-2.5741106158050007156261957778845e6
696 + (4.5966260996517869921896353176509e6
697 - 2.0429449331785719965287268078448e6*y)*y)
698 + x*(-539337.46235914300708358387727104
699 + y*(3.8979389325047153693768107493679e6
700 + y*(-5.148221231610001431252391555769e6
701 + 1.3892025545614289576395342293345e6*y))
702 + x*(565895.74649046444303845732577302
703 + x*(-310527.62984314294347236647479242
704 + 69460.12772807144788197671146672*x
705 + 833521.5327368573745837205376007*y)
706 + y*(-2.8928100253808579470846771599083e6
707 + 2.0838038318421434364593013440017e6*y)))))));
708 case ( 40 ):
709 return 12.727922061357855439215198517887
710 + y*(-254.55844122715710878430397035775
711 + y*(1145.5129855222069895293678666099
712 + y*(-1781.9090885900997614901277925042
713 + 890.9545442950498807450638962521*y)))
714 + x*(-763.67532368147132635291191107324
715 + y*(15018.948032402269418273934251107
716 + y*(-66439.753160288005392703336263372
717 + (101568.81804963568640493728417274
718 - 49893.454480522793321723578190118*y)*y))
719 + x*(10946.012972767755677725070725383
720 + y*(-203901.31142295284413622748025655
721 + y*(851116.14824299979322032032489112
722 + y*(-1.2223896347728084363822276656579e6
723 + 561301.3629058814248693902546388*y)))
724 + x*(-64912.402512925062739997512441225
725 + y*(1.0943467388355484106637227685679e6
726 + y*(-4.0734441765169680547664321336646e6
727 + (5.1140790842535863154766667644871e6
728 - 1.9957381792209117328689431276047e6*y)*y))
729 + x*(195309.96403153629171475722125698
730 + y*(-2.8118525417951774236314216565717e6
731 + y*(8.5798922615613303515749653209078e6
732 + y*(-8.2324199892862608980843904013695e6
733 + 2.12047181542221871617325207308e6*y)))
734 + x*(-325020.21775883419649579930935277
735 + x*(303815.49960461200933406678862197
736 + x*(-149680.36344156837996517073457035
737 + 30292.454506031695945332172472572*x
738 + 605849.0901206339189066434494514*y)
739 + y*(-2.3877581787107336803967712419557e6
740 + 2.7263209055428526350798955225315e6*y))
741 + y*(3.6885518133815065062845645304837e6
742 + y*(-8.018590898655448926705575066269e6
743 + 4.24094363084443743234650414616e6*y)))))));
744 case ( 41 ):
745 return 14.071247279470288663696890340633
746 + y*(-422.137418384108659910906710219
747 + y*(2954.961928688760619376346971533
748 + y*(-7879.8984765033616516702585907547
749 + (8864.885786066281858129040914599
750 - 3545.9543144265127432516163658396*y)*y)))
751 + x*(-703.56236397351443318484451703167
752 + y*(20684.733500821324335634428800731
753 + y*(-141838.17257706050973006465463358
754 + y*(370355.22839565799762850215376547
755 + y*(-407784.74615904896547393588207156
756 + 159567.94414919307344632273646278*y))))
757 + x*(8372.3921312848217548996497526769
758 + y*(-230487.03043772332831135506377958
759 + y*(1.4715710404870027884494207918234e6
760 + y*(-3.5538342129030161049032866244304e6
761 + (3.5902787433568441525422615704126e6
762 - 1.2765435531935445875705818917023e6*y)*y)))
763 + x*(-41369.467001642648671268857601462
764 + y*(1.0105969796115561318267106642643e6
765 + y*(-5.6026078167938901343375538580266e6
766 + y*(1.1386453298547357586663523663641e7
767 + y*(-9.219481217508933132454202551183e6
768 + 2.4112489338100286654110991287709e6*y))))
769 + x*(104901.14846845100198786031748942
770 + x*(-149521.07359165128734044315675957
771 + x*(121645.93273102064549765961699478
772 + x*(-52907.889770808285375500307680782
773 + 9568.448150039796291313885431631*x
774 + 287053.44450119388873941656294892*y)
775 + y*(-1.3001832486230546725255926674745e6
776 + 2.0093741115083572211759159406424e6*y))
777 + y*(2.3491947333075646924041958423687e6
778 + y*(-7.0919086288530254865032327316792e6
779 + 5.358330964022285923135775841713e6*y)))
780 + y*(-2.1364374744419739278090988604184e6
781 + y*(9.352454504299927360326138164902e6
782 + y*(-1.3553425379585782040872844776098e7
783 + 6.028122334525071663527747821927e6*y)))))));
784 case ( 42 ):
785 return 15.297058540778354490084672327068
786 + y*(-642.47645871269088858355623773687
787 + y*(6424.7645871269088858355623773687
788 + y*(-25699.058348507635543342249509475
789 + y*(48185.734403451816643766717830265
790 + y*(-42403.446275037598646514711690633
791 + 14134.482091679199548838237230211*y)))))
792 + x*(-581.2882245495774706232175484286
793 + y*(23771.628972369562877591580796264
794 + y*(-231291.52513656871989008024558527
795 + y*(899467.04219776724401697873283162
796 + y*(-1.638314969717361765888068406229e6
797 + (1.3993137270762407553349854857909e6
798 - 452303.42693373438556282359136676*y)*y))))
799 + x*(5246.8910794869755900990426081844
800 + y*(-196597.79636608341190656820874748
801 + y*(1.7346864385242653991756018418896e6
802 + y*(-6.0392787118992943526854286347266e6
803 + y*(9.6853326150938151453971102838833e6
804 + y*(-7.1237789742063165726144715640264e6
805 + 1.9222895644683711386420002633087e6*y)))))
806 + x*(-20130.929039664314508951428782422
807 + y*(648901.22329981779746939180011424
808 + y*(-4.7543257944739125755183161592528e6
809 + y*(1.2978024465996355949387836002285e7
810 + y*(-1.4648463258649352259705082220401e7
811 + 5.766868693405113415926000789926e6*y))))
812 + x*(41225.572767397665350778191921449
813 + x*(-49042.369681735404495211459480581
814 + x*(34158.331721558065576359073306344
815 + x*(-12971.905642580044607591802133354
816 + 2080.3999615458562106515154364813*x
817 + 87376.79838492596084736364833221*y)
818 + y*(-457443.23860343591267149204126865
819 + 873767.9838492596084736364833221*y))
820 + y*(977206.69370200284153558903759778
821 + y*(-3.7006644021850995182412839293644e6
822 + 3.4950719353970384338945459332886e6*y)))
823 + y*(-1.0825728329308841472632922605866e6
824 + y*(6.0714025348349288971146064466134e6
825 + y*(-1.1307585673343359639070589784169e7
826 + 6.553259878869447063552273624916e6*y)))))));
827 case ( 43 ):
828 return 16.431676725154983403709093484024
829 + y*(-920.17389660867907060770923510535
830 + y*(12422.347604217167453204074673922
831 + y*(-69013.042245650930295578192632901
832 + y*(189785.86617554005831284002974048
833 + y*(-273291.64729277768397048964282629
834 + (197377.3008225616606453536309301
835 - 56393.51452073190304152960883717*y)*y)))))
836 + x*(-394.36024140371960168901824361658
837 + y*(21163.999621999618623977312407423
838 + y*(-273291.64729277768397048964282629
839 + y*(1.4492738871586695362071420452909e6
840 + y*(-3.7957173235108011662568005948096e6
841 + y*(5.1925412985627759954393032136995e6
842 + y*(-3.5527914148061098916163653567417e6
843 + 958689.7468524423517060033502319*y))))))
844 + x*(2300.4347415216976765192730877634
845 + y*(-107660.34590321545126110198050733
846 + y*(1.1801230224006309080543870940226e6
847 + y*(-5.1069651261781688418727862548347e6
848 + y*(1.0248436773479163148893361605986e7
849 + y*(-9.5652076552472189389671374989201e6
850 + 3.3554141139835482309710117258116e6*y)))))
851 + x*(-6441.2172762607534942539646457374
852 + x*(10351.956336847639544336728894935
853 + x*(-10121.912862695469776684801586159
854 + x*(5981.1303279564139589501100281848
855 + x*(-1971.8012070185980084450912180829
856 + 279.33850432763471786305458922841*x
857 + 15642.956242347544200331056996791*y)
858 + y*(-94777.911350693944272594051215851
859 + 211179.90927169184670446926945668*y))
860 + y*(240165.3870148652374286121103625
861 + y*(-1.0683218939626764009755504219573e6
862 + 1.1732217181760658150248292747593e6*y)))
863 + y*(-326661.7332960810700657367784624
864 + y*(2.1739108307380043043107130679364e6
865 + y*(-4.7618999149499141903948952916702e6
866 + 3.2263597249841809913182805055881e6*y))))
867 + y*(253047.82156738674441712003965397
868 + y*(-2.236022568759090141576733441306e6
869 + y*(7.3153824780389986113312884190875e6
870 + y*(-9.8688650411280830322676815465049e6
871 + 4.6459580039772206274983239280469e6*y)))))));
872 case ( 44 ):
873 return 17.492855684535901412622458632637
874 + y*(-1259.4856092865849017088170215498
875 + y*(22040.998162515235779904297877122
876 + y*(-161633.98652511172905263151776556
877 + y*(606127.44946916898394736819162086
878 + y*(-1.2607450948958714866105258385714e6
879 + y*(1.4708692773785167343789468116666e6
880 + y*(-900532.21063990820472180417040814
881 + 225133.05265997705118045104260203*y)))))))
882 + x*(-139.94284547628721130097966906109
883 + y*(8816.3992650060943119617191508489
884 + y*(-132245.98897509141467942578726273
885 + y*(808169.93262555864526315758882782
886 + y*(-2.4245097978766759357894727664835e6
887 + y*(3.7822352846876144598315775157142e6
888 + y*(-2.9417385547570334687578936233333e6
889 + 900532.2106399082047218041704081*y))))))
890 + x*(489.79995916700523955342884171383
891 + x*(-979.59991833401047910685768342766
892 + x*(1224.4998979175130988835721042846
893 + x*(-979.59991833401047910685768342766
894 + x*(489.79995916700523955342884171383
895 + x*(-139.94284547628721130097966906109
896 + 17.492855684535901412622458632637*x
897 + 1259.4856092865849017088170215498*y)
898 + y*(-8816.3992650060943119617191508489
899 + 22040.998162515235779904297877122*y))
900 + y*(26449.197795018282935885157452547
901 + y*(-132245.98897509141467942578726273
902 + 161633.98652511172905263151776556*y)))
903 + y*(-44081.996325030471559808595754245
904 + y*(330614.97243772853669856446815683
905 + y*(-808169.93262555864526315758882782
906 + 606127.4494691689839473681916209*y))))
907 + y*(44081.996325030471559808595754245
908 + y*(-440819.96325030471559808595754245
909 + y*(1.6163398652511172905263151776556e6
910 + y*(-2.4245097978766759357894727664835e6
911 + 1.2607450948958714866105258385714e6*y)))))
912 + y*(-26449.197795018282935885157452547
913 + y*(330614.97243772853669856446815683
914 + y*(-1.6163398652511172905263151776556e6
915 + y*(3.6367646968150139036842091497252e6
916 + y*(-3.7822352846876144598315775157142e6
917 + 1.4708692773785167343789468116666e6*y)))))));
918 #endif
919 }
920
921 printf("Error eval_triangle_2d: There is no shape function %d.\n",i);
922 assert(0);
923 return 0.0;
924 }
925
926 /* \nabla\phi_i(x,y) for triangle */
927 static void
928 grad_triangle_2d ( const int i, DomainType xi, JacobianRangeType grad )
929 {
930 const RangeField &x = xi[0];
931 const RangeField &y = xi[1];
932
933 switch (i) {
934 #if (PMAX2D>=0)
935 case ( 0 ):
936 grad[0]= 0;
937 grad[1]= 0;
938 return;
939 #endif
940 #if (PMAX2D>=1)
941 case ( 1 ):
942 grad[0]= 6.;
943 grad[1]= 0;
944 return;
945 case ( 2 ):
946 grad[0]= 3.4641016151377545870548926830117;
947 grad[1]= 6.9282032302755091741097853660235;
948 return;
949 #endif
950 #if (PMAX2D>=2)
951 case ( 3 ):
952 grad[0]= -19.595917942265424785578272597647
953 + 48.989794855663561963945681494118*x;
954 grad[1]= 0;
955 return;
956 case ( 4 ):
957 grad[0]= -25.455844122715710878430397035775
958 + 42.426406871192851464050661726291*x
959 + 42.426406871192851464050661726291*y;
960 grad[1]= -8.4852813742385702928101323452582
961 + 42.426406871192851464050661726291*x;
962 return;
963 case ( 5 ):
964 grad[0]= -10.954451150103322269139395656016
965 + 10.954451150103322269139395656016*x
966 + 32.863353450309966807418186968048*y;
967 grad[1]= -32.863353450309966807418186968048
968 + 32.863353450309966807418186968048*x
969 + 65.7267069006199336148363739361*y;
970 return;
971 #endif
972 #if (PMAX2D>=3)
973 case ( 6 ):
974 grad[0]= 42.426406871192851464050661726291
975 + x*(-254.55844122715710878430397035775
976 + 296.98484809834996024835463208404*x);
977 grad[1]= 0;
978 return;
979 case ( 7 ):
980 grad[0]= 63.686733312362630553129385942353
981 - 117.57550765359254871346963558588*y
982 + x*(-323.33264604737950896204149786118
983 + 308.63570759068044037285779341294*x
984 + 411.51427678757392049714372455059*y);
985 grad[1]= 9.7979589711327123927891362988236
986 + x*(-117.57550765359254871346963558588
987 + 205.75713839378696024857186227529*x);
988 return;
989 case ( 8 ):
990 grad[0]= 56.920997883030827975980083799789
991 + y*(-303.57865537616441587189378026554
992 + 265.63132345414386388790705773235*y)
993 + x*(-189.73665961010275991993361266596
994 + 132.81566172707193194395352886617*x
995 + 531.2626469082877277758141154647*y);
996 grad[1]= 37.947331922020551983986722533193
997 - 75.89466384404110396797344506639*y
998 + x*(-303.57865537616441587189378026554
999 + 265.63132345414386388790705773235*x
1000 + 531.2626469082877277758141154647*y);
1001 return;
1002 case ( 9 ):
1003 grad[0]= 22.449944320643648313502492393899
1004 + x*(-44.899888641287296627004984787799
1005 + 22.449944320643648313502492393899*x
1006 + 179.59955456514918650801993915119*y)
1007 + y*(-179.59955456514918650801993915119
1008 + 224.49944320643648313502492393899*y);
1009 grad[1]= 89.799777282574593254009969575597
1010 + y*(-448.99888641287296627004984787799
1011 + 448.99888641287296627004984787799*y)
1012 + x*(-179.59955456514918650801993915119
1013 + 89.7997772825745932540099695756*x
1014 + 448.99888641287296627004984787799*y);
1015 return;
1016 #endif
1017 #if (PMAX2D>=4)
1018 case ( 10 ):
1019 grad[0]= -75.894663844041103967973445066385
1020 + x*(796.89397036243159166372117319705
1021 + x*(-2125.0505876331509111032564618588
1022 + 1593.7879407248631833274423463941*x));
1023 grad[1]= 0;
1024 return;
1025 case ( 11 ):
1026 grad[0]= -120.49896265113654496053335221618
1027 + 230.04347415216976765192730877634*y
1028 + x*(1150.2173707608488382596365438817
1029 - 1840.3477932173581412154184702107*y
1030 + x*(-2760.521689826037211823127705316
1031 + 1840.3477932173581412154184702107*x
1032 + 2760.521689826037211823127705316*y));
1033 grad[1]= -10.954451150103322269139395656016
1034 + x*(230.04347415216976765192730877634
1035 + x*(-920.17389660867907060770923510535
1036 + 920.1738966086790706077092351053*x));
1037 return;
1038 case ( 12 ):
1039 grad[0]= -127.27922061357855439215198517887
1040 + (721.24891681027847488886124934695
1041 - 678.8225099390856234248105876207*y)*y
1042 + x*(975.80735803743558367316521970469
1043 + y*(-4412.3463146040565522612688195343
1044 + 3054.7012947258853054116476442929*y)
1045 + x*(-1866.7619023324854644182291159568
1046 + 1018.233764908628435137215881431*x
1047 + 4582.0519420888279581174714664394*y));
1048 grad[1]= -42.426406871192851464050661726291
1049 + 84.85281374238570292810132345258*y
1050 + x*(721.24891681027847488886124934695
1051 - 1357.6450198781712468496211752413*y
1052 + x*(-2206.1731573020282761306344097671
1053 + 1527.3506473629426527058238221465*x
1054 + 3054.7012947258853054116476442929*y));
1055 return;
1056 case ( 13 ):
1057 grad[0]= -100.39920318408906575738064309422
1058 + y*(1104.3912350249797233311870740364
1059 + y*(-2509.9800796022266439345160773556
1060 + 1505.9880477613359863607096464133*y))
1061 + x*(501.99601592044532878690321547111
1062 + x*(-702.79442228862346030166450165956
1063 + 301.19760955226719727214192928267*x
1064 + 2710.778485970404775449277363544*y)
1065 + y*(-3815.1697209953844987804644375805
1066 + 4517.96414328400795908212893924*y));
1067 grad[1]= -100.39920318408906575738064309422
1068 + (501.99601592044532878690321547111
1069 - 501.9960159204453287869032154711*y)*y
1070 + x*(1104.3912350249797233311870740364
1071 + y*(-5019.9601592044532878690321547111
1072 + 4517.96414328400795908212893924*y)
1073 + x*(-1907.5848604976922493902322187902
1074 + 903.592828656801591816425787848*x
1075 + 4517.96414328400795908212893924*y));
1076 return;
1077 case ( 14 ):
1078 grad[0]= -37.947331922020551983986722533193
1079 + y*(569.20997883030827975980083799789
1080 + y*(-1707.6299364909248392794025139937
1081 + 1328.1566172707193194395352886617*y))
1082 + x*(113.84199576606165595196016759958
1083 + x*(-113.84199576606165595196016759958
1084 + 37.947331922020551983986722533193*x
1085 + 569.2099788303082797598008379979*y)
1086 + y*(-1138.4199576606165595196016759958
1087 + 1707.6299364909248392794025139937*y));
1088 grad[1]= -189.73665961010275991993361266596
1089 + y*(1707.6299364909248392794025139937
1090 + y*(-3984.4698518121579583186058659852
1091 + 2656.3132345414386388790705773235*y))
1092 + x*(569.20997883030827975980083799789
1093 + x*(-569.20997883030827975980083799789
1094 + 189.73665961010275991993361266596*x
1095 + 1707.6299364909248392794025139937*y)
1096 + y*(-3415.2598729818496785588050279873
1097 + 3984.4698518121579583186058659852*y));
1098 return;
1099 #endif
1100 #if (PMAX2D>=5)
1101 case ( 15 ):
1102 grad[0]= 121.24355652982141054692124390541
1103 + x*(-1939.8969044771425687507399024866
1104 + x*(8729.5360701471415593783295611896
1105 + x*(-14549.226783578569265630549268649
1106 + 8002.074730968213096096802097757*x)));
1107 grad[1]= 0;
1108 return;
1109 case ( 16 ):
1110 grad[0]= 198.
1111 - 384.*y
1112 + x*(-2976.
1113 + 5184.*y
1114 + x*(12528.
1115 - 17280.*y
1116 + x*(-19440.
1117 + 9900.*x
1118 + 15840.*y)));
1119 grad[1]= 12.
1120 + x*(-384.
1121 + x*(2592.
1122 + x*(-5760.
1123 + 3960.*x)));
1124 return;
1125 case ( 17 ):
1126 grad[0]= 224.63303408003017934039739318738
1127 + y*(-1301.3224043256920734202331743269
1128 + 1254.8466041712030707980819895295*y)
1129 + x*(-2943.4673431176368327362417038346
1130 + (15058.159250054436849576983874354
1131 - 12548.466041712030707980819895295*y)*y
1132 + x*(10735.909835686959605716923688197
1133 + y*(-41828.22013904010235993606631765
1134 + 23005.521076472056297964836474707*y)
1135 + x*(-14407.498047891590812866867287191
1136 + 6390.422521242237860545787909641*x
1137 + 30674.028101962741730619781966277*y)));
1138 grad[1]= 46.475800154489002622151184797389
1139 - 92.95160030897800524430236959478*y
1140 + x*(-1301.3224043256920734202331743269
1141 + 2509.693208342406141596163979059*y
1142 + x*(7529.079625027218424788491937177
1143 - 12548.466041712030707980819895295*y
1144 + x*(-13942.740046346700786645355439217
1145 + 7668.507025490685432654945491569*x
1146 + 15337.014050981370865309890983138*y)));
1147 return;
1148 case ( 18 ):
1149 grad[0]= 210.79848196796864030305017091149
1150 + y*(-2419.5999669366835234784889182884
1151 + (5774.0453756443584083009394640973
1152 - 3666.0605559646720052704377549824*y)*y)
1153 + x*(-2162.9757280191564831095582754396
1154 + y*(21116.508802356510750357721468699
1155 + y*(-41243.181254602560059292424743552
1156 + 20163.333057805696028987407652403*y))
1157 + x*(6213.9726423601190489333919946952
1158 + x*(-6782.2120285346432097503098467175
1159 + 2520.4166322257120036234259565504*x
1160 + 24195.999669366835234784889182884*y)
1161 + y*(-42892.908504786662461664121733294
1162 + 45367.499380062816065221667217907*y)));
1163 grad[1]= 109.98181667894016015811313264947
1164 + y*(-549.90908339470080079056566324736
1165 + 549.9090833947008007905656632474*y)
1166 + x*(-2419.5999669366835234784889182884
1167 + (11548.090751288716816601878928195
1168 - 10998.181667894016015811313264947*y)*y
1169 + x*(10558.254401178255375178860734349
1170 + y*(-41243.181254602560059292424743552
1171 + 30244.999586708544043481111478605*y)
1172 + x*(-14297.636168262220820554707244431
1173 + 6048.999917341708808696222295721*x
1174 + 30244.999586708544043481111478605*y)));
1175 return;
1176 case ( 19 ):
1177 grad[0]= 155.88457268119895641747017073553
1178 + y*(-2909.8453567157138531261098537299
1179 + y*(12158.996669133518600562673317371
1180 + y*(-17459.072140294283118756659122379
1181 + 8002.074730968213096096802097757*y)))
1182 + x*(-1039.2304845413263761164678049035
1183 + x*(2182.3840175367853898445823902974
1184 + x*(-1870.6148721743874770096420488263
1185 + 571.5767664977295068640572926969*x
1186 + 9145.228263963672109824916683151*y)
1187 + y*(-21200.301884643058072775943220032
1188 + 30865.145390877393370659093805635*y))
1189 + y*(14964.918977395099816077136390611
1190 + y*(-43024.142060010911971221767123006
1191 + 32008.298923872852384387208391029*y)));
1192 grad[1]= 207.8460969082652752232935609807
1193 + y*(-1870.6148721743874770096420488263
1194 + (4364.7680350735707796891647805948
1195 - 2909.8453567157138531261098537299*y)*y)
1196 + x*(-2909.8453567157138531261098537299
1197 + y*(24317.993338267037201125346634742
1198 + y*(-52377.216420882849356269977367138
1199 + 32008.298923872852384387208391029*y))
1200 + x*(7482.4594886975499080385681953054
1201 + x*(-7066.767294881019357591981073344
1202 + 2286.3070659909180274562291707878*x
1203 + 20576.76359391826224710606253709*y)
1204 + y*(-43024.142060010911971221767123006
1205 + 48012.448385809278576580812586543*y)));
1206 return;
1207 case ( 20 ):
1208 grad[0]= 57.445626465380286598506114682189
1209 + y*(-1378.6950351691268783641467523725
1210 + y*(7238.1489346379161114117704499559
1211 + y*(-12867.82032824518419806536968881
1212 + 7238.148934637916111411770449956*y)))
1213 + x*(-229.78250586152114639402445872876
1214 + x*(344.67375879228171959103668809314
1215 + x*(-229.78250586152114639402445872876
1216 + 57.44562646538028659850611468219*x
1217 + 1378.6950351691268783641467523725*y)
1218 + y*(-4136.0851055073806350924402571176
1219 + 7238.148934637916111411770449956*y))
1220 + y*(4136.0851055073806350924402571176
1221 + y*(-14476.297869275832222823540899912
1222 + 12867.82032824518419806536968881*y)));
1223 grad[1]= 344.67375879228171959103668809314
1224 + y*(-4825.4326230919440742745136333039
1225 + y*(19301.730492367776297098054533216
1226 + y*(-28952.595738551664445647081799823
1227 + 14476.297869275832222823540899912*y)))
1228 + x*(-1378.6950351691268783641467523725
1229 + x*(2068.0425527536903175462201285588
1230 + x*(-1378.6950351691268783641467523725
1231 + 344.67375879228171959103668809314*x
1232 + 4825.4326230919440742745136333039*y)
1233 + y*(-14476.297869275832222823540899912
1234 + 19301.730492367776297098054533216*y))
1235 + y*(14476.297869275832222823540899912
1236 + y*(-38603.460984735552594196109066431
1237 + 28952.595738551664445647081799823*y)));
1238 return;
1239 #endif
1240 #if (PMAX2D>=6)
1241 case ( 21 ):
1242 grad[0]= -179.59955456514918650801993915119
1243 + x*(4040.9899777158566964304486309019
1244 + x*(-26939.933184772377976202990872679
1245 + x*(74084.816258124039434558224899868
1246 + x*(-88901.779509748847321469869879841
1247 + 38524.104454224500505970276947931*x))));
1248 grad[1]= 0;
1249 return;
1250 case ( 22 ):
1251 grad[0]= -298.11407212676157062443450206005
1252 + 583.2666628567074207869370692479*y
1253 + x*(6415.9332914237816286563077617271
1254 - 11665.333257134148415738741384958*y
1255 + x*(-40828.666399969519455085594847354
1256 + 64159.33291423781628656307761727*y
1257 + x*(106932.22152372969381093846269545
1258 - 128318.66582847563257312615523454*y
1259 + x*(-121902.73253705185094446984747282
1260 + 50044.27967310549670351920054147*x
1261 + 83407.13278850916117253200090245*y))));
1262 grad[1]= -12.961481396815720461931934872176
1263 + x*(583.26666285670742078693706924792
1264 + x*(-5832.6666285670742078693706924792
1265 + x*(21386.44430474593876218769253909
1266 + x*(-32079.666457118908143281538808636
1267 + 16681.426557701832234506400180491*x))));
1268 return;
1269 case ( 23 ):
1270 grad[0]= -351.39721114431173015083225082978
1271 + (2058.1836652738258480263031834316
1272 - 2007.9840636817813151476128618844*y)*y
1273 + x*(6877.3454181101010043805740519542
1274 + y*(-37147.705178112954330230837944862
1275 + 33131.737050749391699935612221093*y)
1276 + x*(-39657.685257715180974165354022218
1277 + (182224.55377912165434964586721601
1278 - 132526.94820299756679974244888437*y)*y
1279 + x*(93873.254977123276483150901293098
1280 + y*(-320273.45815724411976604425147057
1281 + 143570.86055324736403305431962474*y)
1282 + x*(-96634.233064685725791478868978189
1283 + 35892.715138311841008263579906185*x
1284 + 179463.57569155920504131789953092*y))));
1285 grad[1]= -50.199601592044532878690321547111
1286 + 100.39920318408906575738064309422*y
1287 + x*(2058.1836652738258480263031834316
1288 - 4015.9681273635626302952257237689*y
1289 + x*(-18573.852589056477165115418972431
1290 + 33131.737050749391699935612221093*y
1291 + x*(60741.517926373884783215289072005
1292 - 88351.29880199837786649496592292*y
1293 + x*(-80068.364539311029941511062867642
1294 + 35892.715138311841008263579906185*x
1295 + 71785.43027662368201652715981237*y))));
1296 return;
1297 case ( 24 ):
1298 grad[0]= -356.38181771801995229802555850084
1299 + y*(4157.7878733768994434769648491765
1300 + y*(-10097.484835343898648444057490857
1301 + 6533.666658163699125463801905849*y))
1302 + x*(5939.6969619669992049670926416807
1303 + y*(-62960.787796850191572651182001816
1304 + (137206.99982143768163473984002282
1305 - 78403.99989796438950556562287019*y)*y)
1306 + x*(-29104.515113638296104338753944236
1307 + y*(254812.9996683842658930882743281
1308 + y*(-431221.99943880414228061092578602
1309 + 169875.33311225617726205884955207*y))
1310 + x*(58802.999923473292129174217152639
1311 + x*(-52269.33326530959300371041524679
1312 + 16987.533311225617726205884955207*x
1313 + 169875.33311225617726205884955207*y)
1314 + y*(-365885.33285716715102597290672753
1315 + 339750.66622451235452411769910414*y))));
1316 grad[1]= -118.79393923933998409934185283361
1317 + (593.96969619669992049670926416807
1318 - 593.9696961966999204967092641681*y)*y
1319 + x*(4157.7878733768994434769648491765
1320 + y*(-20194.969670687797296888114981714
1321 + 19600.999974491097376391405717546*y)
1322 + x*(-31480.393898425095786325591000908
1323 + (137206.99982143768163473984002282
1324 - 117605.99984694658425834843430528*y)*y
1325 + x*(84937.666556128088631029424776034
1326 + y*(-287481.33295920276152040728385735
1327 + 169875.33311225617726205884955207*y)
1328 + x*(-91471.333214291787756493226681883
1329 + 33975.066622451235452411769910414*x
1330 + 169875.33311225617726205884955207*y))));
1331 return;
1332 case ( 25 ):
1333 grad[0]= -314.29922048901107638903489351459
1334 + y*(6061.4849665737850446456729463528
1335 + y*(-26266.434855153068526797916100862
1336 + (39287.402561126384548629361689324
1337 - 18857.953229340664583342093610875*y)*y))
1338 + x*(4040.9899777158566964304486309019
1339 + y*(-68696.829621169563839317626725332
1340 + y*(256602.86358495690022333348806227
1341 + y*(-320585.20489879129791681559138488
1342 + 122576.69599071431979172360847069*y)))
1343 + x*(-15490.461581244117336316719751791
1344 + x*(25368.437082327322594257816405106
1345 + x*(-18857.953229340664583342093610875
1346 + 5253.286971030613705359583220172*x
1347 + 87554.78285051022842265972033621*y)
1348 + y*(-231683.42538904245059534572150504
1349 + 315197.21826183682232157499321035*y))
1350 + y*(206763.98719312800096735795494781
1351 + y*(-545533.64699164065401811056517175
1352 + 367730.08797214295937517082541207*y))));
1353 grad[1]= -224.49944320643648313502492393899
1354 + y*(2020.4949888579283482152243154509
1355 + y*(-4714.4883073351661458355234027189
1356 + 3142.9922048901107638903489351459*y))
1357 + x*(6061.4849665737850446456729463528
1358 + y*(-52532.869710306137053595832201724
1359 + (117862.20768337915364588808506797
1360 - 75431.8129173626583333683744435*y)*y)
1361 + x*(-34348.414810584781919658813362666
1362 + y*(256602.86358495690022333348806227
1363 + y*(-480877.80734818694687522338707732
1364 + 245153.39198142863958344721694138*y))
1365 + x*(68921.329064376000322452651649271
1366 + x*(-57920.85634726061264883643037626
1367 + 17510.956570102045684531944067241*x
1368 + 157598.60913091841116078749660517*y)
1369 + y*(-363689.09799442710267874037678117
1370 + 367730.08797214295937517082541207*y))));
1371 return;
1372 case ( 26 ):
1373 grad[0]= -223.37412562783541873039835519809
1374 + y*(6328.9335594553368640279533972792
1375 + y*(-41696.50345052927816300769297031
1376 + y*(104241.25862632319540751923242578
1377 + y*(-109453.32155763935517789519404706
1378 + 40654.090864266046208932500646052*y))))
1379 + x*(1861.4510468986284894199862933174
1380 + y*(-43185.664288048180954543682004964
1381 + y*(218906.64311527871035579038809413
1382 + y*(-375268.53105476350346706923673279
1383 + 203270.45432133023104466250323026*y)))
1384 + x*(-5212.0629313161597703759616212888
1385 + x*(6701.2237688350625619119506559427
1386 + x*(-4095.1923031769826767239698452983
1387 + 967.9545443872868144983928725251*x
1388 + 24198.863609682170362459821813126*y)
1389 + y*(-78925.524388501847951407418836659
1390 + 135513.63621422015402977500215351*y))
1391 + y*(91583.391507412521679463325631217
1392 + y*(-312723.77587896958622255769727733
1393 + 271027.27242844030805955000430702*y))));
1394 grad[1]= -372.29020937972569788399725866348
1395 + y*(5212.0629313161597703759616212888
1396 + y*(-20848.251725264639081503846485155
1397 + (31272.377587896958622255769727733
1398 - 15636.188793948479311127884863866*y)*y))
1399 + x*(6328.9335594553368640279533972792
1400 + y*(-83393.00690105855632601538594062
1401 + y*(312723.77587896958622255769727733
1402 + y*(-437813.28623055742071158077618826
1403 + 203270.45432133023104466250323026*y)))
1404 + x*(-21592.832144024090477271841002482
1405 + x*(30527.797169137507226487775210406
1406 + x*(-19731.381097125461987851854709165
1407 + 4839.7727219364340724919643626253*x
1408 + 67756.81810711007701488750107675*y)
1409 + y*(-208482.51725264639081503846485155
1410 + 271027.27242844030805955000430702*y))
1411 + y*(218906.64311527871035579038809413
1412 + y*(-562902.79658214525520060385509919
1413 + 406540.90864266046208932500646052*y))));
1414 return;
1415 case ( 27 ):
1416 grad[0]= -80.944425379392248793301833668973
1417 + y*(2833.0548882787287077655641784141
1418 + y*(-22664.439106229829662124513427313
1419 + y*(67993.317318689488986373540281938
1420 + y*(-84991.646648361861232966925352422
1421 + 37396.324525279218942505447155066*y))))
1422 + x*(404.72212689696124396650916834487
1423 + x*(-809.44425379392248793301833668973
1424 + x*(809.44425379392248793301833668973
1425 + x*(-404.72212689696124396650916834487
1426 + 80.94442537939224879330183366897*x
1427 + 2833.0548882787287077655641784141*y)
1428 + y*(-11332.219553114914831062256713656
1429 + 22664.439106229829662124513427313*y))
1430 + y*(16998.329329672372246593385070484
1431 + y*(-67993.317318689488986373540281938
1432 + 67993.31731868948898637354028194*y)))
1433 + y*(-11332.219553114914831062256713656
1434 + y*(67993.317318689488986373540281938
1435 + y*(-135986.63463737897797274708056388
1436 + 84991.64664836186123296692535242*y))));
1437 grad[1]= -566.61097765574574155311283568281
1438 + y*(11332.219553114914831062256713656
1439 + y*(-67993.317318689488986373540281938
1440 + y*(169983.29329672372246593385070484
1441 + y*(-186981.62262639609471252723577533
1442 + 74792.64905055843788501089431013*y))))
1443 + x*(2833.0548882787287077655641784141
1444 + x*(-5666.1097765574574155311283568281
1445 + x*(5666.1097765574574155311283568281
1446 + x*(-2833.0548882787287077655641784141
1447 + 566.6109776557457415531128356828*x
1448 + 11332.219553114914831062256713656*y)
1449 + y*(-45328.878212459659324249026854625
1450 + 67993.31731868948898637354028194*y))
1451 + y*(67993.317318689488986373540281938
1452 + y*(-203979.95195606846695912062084581
1453 + 169983.29329672372246593385070484*y)))
1454 + y*(-45328.878212459659324249026854625
1455 + y*(203979.95195606846695912062084581
1456 + y*(-339966.58659344744493186770140969
1457 + 186981.62262639609471252723577533*y))));
1458 return;
1459 #endif
1460 #if (PMAX2D>=7)
1461 case ( 28 ):
1462 grad[0]= 252.
1463 + x*(-7560.
1464 + x*(69300.
1465 + x*(-277200.
1466 + x*(540540.
1467 + x*(-504504.
1468 + 180180.*x)))));
1469 grad[1]= 0;
1470 return;
1471 case ( 29 ):
1472 grad[0]= 422.62039704680605962069690732743
1473 - 831.3843876330611008931742439228*y
1474 + x*(-12262.919717587651238174320097862
1475 + 22863.070659909180274562291707878*y
1476 + x*(108599.58563456860630417088561242
1477 - 182904.56527927344219649833366302*y
1478 + x*(-419156.29543166830503364201464442
1479 + 594439.8371576386871386195844048*y
1480 + x*(787632.78423387126045867094933638
1481 - 832215.7720206941619940674181667*y
1482 + x*(-707383.40621759003769495730544173
1483 + 242729.60017270246391493633029863*x
1484 + 416107.88601034708099703370908337*y)))));
1485 grad[1]= 13.856406460551018348219570732047
1486 + x*(-831.38438763306110089317424392282
1487 + x*(11431.535329954590137281145853939
1488 + x*(-60968.188426424480732166111221007
1489 + x*(148609.9592894096717846548961012
1490 + x*(-166443.15440413883239881348363335
1491 + 69351.3143350578468328389515139*x)))));
1492 return;
1493 case ( 30 ):
1494 grad[0]= 509.82349886995205078129159647073
1495 + y*(-3005.2753617597173519739294107748
1496 + 2951.6097302997223992601092427253*y)
1497 + x*(-13792.067285218702847451783188735
1498 + (76741.852987792782380762840310857
1499 - 70838.63352719333758224262182541*y)*y
1500 + x*(113636.97461653931237151420584492
1501 + y*(-566709.06821754670065794097460325
1502 + 460451.11792675669428457704186514*y)
1503 + x*(-407322.14278136169109789507549609
1504 + (1.6883207657314412123767824868389e6
1505 - 1.074385941829098953330679764352e6*y)*y
1506 + x*(709862.14013708323702205627287543
1507 + y*(-2.148771883658197906661359528704e6
1508 + 805789.456371824214998009823264*y)
1509 + x*(-590912.2680060044243318738703936
1510 + 188017.5398200923168328689587616*x
1511 + 966947.3476461890579976117879168*y)))));
1512 grad[1]= 53.665631459994952713820168049551
1513 - 107.3312629199899054276403360991*y
1514 + x*(-3005.2753617597173519739294107748
1515 + 5903.219460599444798520218485451*y
1516 + x*(38370.926493896391190381420155429
1517 - 70838.63352719333758224262182541*y
1518 + x*(-188903.02273918223355264699153442
1519 + 306967.41195117112952305136124343*y
1520 + x*(422080.19143286030309419562170972
1521 - 537192.970914549476665339882176*y
1522 + x*(-429754.3767316395813322719057408
1523 + 161157.8912743648429996019646528*x
1524 + 322315.7825487296859992039293056*y)))));
1525 return;
1526 case ( 31 ):
1527 grad[0]= 539.73326745717648046232961374241
1528 + y*(-6349.8031465550174172038778087342
1529 + (15557.017709059792672149500631399
1530 - 10159.685034488027867526204493975*y)*y)
1531 + x*(-13017.096450437785705267949507905
1532 + y*(143505.55111214339362880763847739
1533 + y*(-327649.84236223889872772009493069
1534 + 198113.85817251654341676098763251*y))
1535 + x*(95405.792276989136693488264076232
1536 + y*(-929611.18065565454987864771119869
1537 + (1.8325531880957780266050391356007e6
1538 - 924531.3381384105359448846089517*y)*y)
1539 + x*(-304155.57071998533428406574703837
1540 + y*(2.4103852744322846115705920161955e6
1541 + y*(-3.5825589352863408267864278596878e6
1542 + 1.1556641726730131699311057611896e6*y))
1543 + x*(472584.09918235717127539860591504
1544 + x*(-352477.57266526901682898725716284
1545 + 101120.61510888865236897175410409*x
1546 + 1.0400977554057118529379951850707e6*y)
1547 + y*(-2.6580275971479302908415432507361e6
1548 + 2.1668703237618996936208233022306e6*y)))));
1549 grad[1]= 126.99606293110034834407755617468
1550 + y*(-634.98031465550174172038778087342
1551 + 634.9803146555017417203877808734*y)
1552 + x*(-6349.8031465550174172038778087342
1553 + (31114.035418119585344299001262798
1554 - 30479.055103464083602578613481924*y)*y
1555 + x*(71752.775556071696814403819238697
1556 + y*(-327649.84236223889872772009493069
1557 + 297170.78725877481512514148144876*y)
1558 + x*(-309870.39355188484995954923706623
1559 + (1.2217021253971853510700260904005e6
1560 - 924531.3381384105359448846089517*y)*y
1561 + x*(602596.31860807115289264800404888
1562 + y*(-1.7912794676431704133932139298439e6
1563 + 866748.1295047598774483293208922*y)
1564 + x*(-531605.51942958605816830865014723
1565 + 173349.62590095197548966586417844*x
1566 + 866748.1295047598774483293208922*y)))));
1567 return;
1568 case ( 32 ):
1569 grad[0]= 516.
1570 + y*(-10080.
1571 + y*(44280.
1572 + y*(-67200.
1573 + 32760.*y)))
1574 + x*(-10440.
1575 + y*(188640.
1576 + y*(-760320.
1577 + (1.04832e6
1578 - 458640.*y)*y))
1579 + x*(64260.
1580 + y*(-1.00224e6
1581 + y*(3.3696e6
1582 + y*(-3.66912e6
1583 + 1.1466e6*y)))
1584 + x*(-173520.
1585 + x*(231660.
1586 + x*(-150696.
1587 + 38220.*x
1588 + 655200.*y)
1589 + y*(-1.9656e6
1590 + 2.457e6*y))
1591 + y*(2.13408e6
1592 + y*(-5.11056e6
1593 + 3.0576e6*y)))));
1594 grad[1]= 240.
1595 + y*(-2160.
1596 + (5040.
1597 - 3360.*y)*y)
1598 + x*(-10080.
1599 + y*(88560.
1600 + y*(-201600.
1601 + 131040.*y))
1602 + x*(94320.
1603 + y*(-760320.
1604 + (1.57248e6
1605 - 917280.*y)*y)
1606 + x*(-334080.
1607 + y*(2.2464e6
1608 + y*(-3.66912e6
1609 + 1.5288e6*y))
1610 + x*(533520.
1611 + x*(-393120.
1612 + 109200.*x
1613 + 982800.*y)
1614 + y*(-2.55528e6
1615 + 2.2932e6*y)))));
1616 return;
1617 case ( 33 ):
1618 grad[0]= 437.79447232691278008317112124053
1619 + y*(-12735.839194964735420601341708815
1620 + y*(86364.909540854612070952848462905
1621 + y*(-222877.18591188286986052347990427
1622 + (242378.93967917262097331928439589
1623 - 93608.41808299080534141986155979*y)*y)))
1624 + x*(-6765.9145723250156921944627828082
1625 + y*(177505.75877982099992463120006662
1626 + y*(-1.0698104923770377753305127035405e6
1627 + y*(2.4070736078483349944936535829661e6
1628 + y*(-2.2231999294710316268587217120451e6
1629 + 702063.1356224310400606489616985*y))))
1630 + x*(32436.590449675810524344042164639
1631 + y*(-706839.07532054281584337446483926
1632 + y*(3.3431577886782430479078521985641e6
1633 + y*(-5.3044770247028123026804588217216e6
1634 + 2.6327367585841164002274336063692e6*y)))
1635 + x*(-70843.105521991340777094963255286
1636 + x*(79002.002506265624405917697787496
1637 + x*(-44018.244217596866797453387281093
1638 + 9750.876883644875556397902245812*x
1639 + 250736.8341508682285930889148923*y)
1640 + y*(-891508.74364753147944209391961708
1641 + 1.4626315325467313334596853368718e6*y))
1642 + y*(1.1828410652323498021883496112062e6
1643 + y*(-3.8223437383887912181079776803582e6
1644 + 3.1202806027663601780473287186598e6*y)))));
1645 grad[1]= 397.99497484264798189379192840048
1646 + y*(-5571.9296477970717465130869976068
1647 + y*(22287.718591188286986052347990427
1648 + y*(-33431.577886782430479078521985641
1649 + 16715.78894339121523953926099282*y)))
1650 + x*(-12735.839194964735420601341708815
1651 + y*(172729.81908170922414190569692581
1652 + y*(-668631.55773564860958157043971281
1653 + (969515.75871669048389327713758358
1654 - 468042.09041495402670709930779897*y)*y))
1655 + x*(88752.879389910499962315600033308
1656 + y*(-1.0698104923770377753305127035405e6
1657 + y*(3.6106104117725024917404803744492e6
1658 + y*(-4.4463998589420632537174434240902e6
1659 + 1.7551578390560776001516224042461e6*y)))
1660 + x*(-235613.02510684760528112482161309
1661 + x*(295710.26630808745054708740280156
1662 + x*(-178301.74872950629588841878392342
1663 + 41789.472358478038098848152482051*x
1664 + 585052.6130186925333838741347487*y)
1665 + y*(-1.9111718691943956090539888401791e6
1666 + 2.3402104520747701335354965389948e6*y))
1667 + y*(2.2287718591188286986052347990427e6
1668 + y*(-5.3044770247028123026804588217216e6
1669 + 3.5103156781121552003032448084923e6*y)))));
1670 return;
1671 case ( 34 ):
1672 grad[0]= 302.86630713897510062201458646752
1673 + y*(-12114.652285559004024880583458701
1674 + y*(115089.19671281053823636554285766
1675 + y*(-436127.48228012414489570100451323
1676 + y*(772309.08320438650658613719549218
1677 + y*(-639653.64067751541251369480661941
1678 + 199891.76271172356641052962706856*y)))))
1679 + x*(-3028.6630713897510062201458646752
1680 + y*(102974.54442725153421148495939896
1681 + y*(-799567.05084689426564211850827426
1682 + y*(2.3260132388273287727770720240706e6
1683 + y*(-2.816656656392468435784735654148e6
1684 + 1.1993505762703413984631777624114e6*y))))
1685 + x*(10600.320749864128521770510526363
1686 + x*(-18171.978428338506037320875188051
1687 + x*(16657.646892643630534210802255714
1688 + x*(-7874.5239856133526161723792481556
1689 + 1514.3315356948755031100729323376*x
1690 + 54515.93528501551811196262556415*y)
1691 + y*(-230178.39342562107647273108571532
1692 + 454299.46070846265093302187970128*y))
1693 + y*(375554.22085232912477129808721973
1694 + y*(-1.4779875788381984910354311819615e6
1695 + 1.4537582742670804829856700150441e6*y)))
1696 + y*(-290751.65485341609659713400300882
1697 + y*(1.7081659722638195675081622676768e6
1698 + y*(-3.3436440308142851108670410346014e6
1699 + 2.0443475731880819291985984586558e6*y)))));
1700 grad[1]= 605.73261427795020124402917293504
1701 + y*(-12114.652285559004024880583458701
1702 + y*(72687.913713354024149283500752205
1703 + y*(-181719.78428338506037320875188051
1704 + (199891.76271172356641052962706856
1705 - 79956.70508468942656421185082743*y)*y)))
1706 + x*(-12114.652285559004024880583458701
1707 + y*(230178.39342562107647273108571532
1708 + y*(-1.3083824468403724346871030135397e6
1709 + y*(3.0892363328175460263445487819687e6
1710 + y*(-3.198268203387577062568474033097e6
1711 + 1.1993505762703413984631777624114e6*y))))
1712 + x*(51487.272213625767105742479699479
1713 + x*(-96917.218284472032199044667669607
1714 + x*(93888.555213082281192824521804932
1715 + x*(-46035.678685124215294546217143063
1716 + 9085.989214169253018660437594026*x
1717 + 181719.78428338506037320875188051*y)
1718 + y*(-738993.78941909924551771559098075
1719 + 1.0903187057003103622392525112831e6*y))
1720 + y*(1.1387773148425463783387748451179e6
1721 + y*(-3.3436440308142851108670410346014e6
1722 + 2.7257967642507759055981312782077e6*y)))
1723 + y*(-799567.05084689426564211850827426
1724 + y*(3.4890198582409931591656080361058e6
1725 + y*(-5.6333133127849368715694713082959e6
1726 + 2.9983764406758534961579444060285e6*y)))));
1727 return;
1728 case ( 35 ):
1729 grad[0]= 108.44353369380767278501943119391
1730 + y*(-5205.2896173027682936809326973075
1731 + y*(58559.50819465614330391049284471
1732 + y*(-260264.48086513841468404663486538
1733 + y*(536795.49178434798028584618440984
1734 + y*(-515323.67211297406107441233703345
1735 + 186089.10381857396649909334392874*y)))))
1736 + x*(-650.66120216284603671011658716344
1737 + y*(26026.448086513841468404663486538
1738 + y*(-234238.03277862457321564197137884
1739 + y*(780793.44259541524405213990459613
1740 + y*(-1.0735909835686959605716923688197e6
1741 + 515323.6721129740610744123370334*y))))
1742 + x*(1626.6530054071150917752914679086
1743 + x*(-2168.8706738761534557003886238781
1744 + x*(1626.6530054071150917752914679086
1745 + x*(-650.66120216284603671011658716344
1746 + 108.44353369380767278501943119391*x
1747 + 5205.289617302768293680932697308*y)
1748 + y*(-26026.448086513841468404663486538
1749 + 58559.50819465614330391049284471*y))
1750 + y*(52052.896173027682936809326973075
1751 + y*(-234238.03277862457321564197137884
1752 + 260264.48086513841468404663486538*y)))
1753 + y*(-52052.896173027682936809326973075
1754 + y*(351357.04916793685982346295706826
1755 + y*(-780793.44259541524405213990459613
1756 + 536795.4917843479802858461844098*y)))));
1757 grad[1]= 867.54826955046138228015544955126
1758 + y*(-23423.803277862457321564197137884
1759 + y*(195198.36064885381101303497614903
1760 + y*(-715727.32237913064038112824587979
1761 + y*(1.2883091802824351526860308425836e6
1762 + y*(-1.1165346229114437989945600635725e6
1763 + 372178.20763714793299818668785749*y)))))
1764 + x*(-5205.2896173027682936809326973075
1765 + y*(117119.01638931228660782098568942
1766 + y*(-780793.44259541524405213990459613
1767 + y*(2.1471819671373919211433847376394e6
1768 + y*(-2.5766183605648703053720616851672e6
1769 + 1.1165346229114437989945600635725e6*y))))
1770 + x*(13013.224043256920734202331743269
1771 + x*(-17350.965391009227645603108991025
1772 + x*(13013.224043256920734202331743269
1773 + x*(-5205.2896173027682936809326973075
1774 + 867.5482695504613822801554495513*x
1775 + 23423.803277862457321564197137884*y)
1776 + y*(-117119.01638931228660782098568942
1777 + 195198.36064885381101303497614903*y))
1778 + y*(234238.03277862457321564197137884
1779 + y*(-780793.44259541524405213990459613
1780 + 715727.3223791306403811282458798*y)))
1781 + y*(-234238.03277862457321564197137884
1782 + y*(1.1711901638931228660782098568942e6
1783 + y*(-2.1471819671373919211433847376394e6
1784 + 1.2883091802824351526860308425836e6*y)))));
1785 return;
1786 #endif
1787 #if (PMAX2D>=8)
1788 case ( 36 ):
1789 grad[0]= -339.41125496954281171240529381033
1790 + x*(13067.333316327398250927603811698
1791 + x*(-156807.99979592877901113124574037
1792 + x*(849376.66556128088631029424776034
1793 + x*(-2.378254663571586481668823893729e6
1794 + x*(3.5673819953573797225032358405934e6
1795 + x*(-2.7180053297960988361929415928331e6
1796 + 825108.7608309585752728572692529*x))))));
1797 grad[1]= 0;
1798 return;
1799 case ( 37 ):
1800 grad[0]= -573.18059981126367497816447348118
1801 + 1131.6642611658282813671452425141*y
1802 + x*(21501.620962150737345975759607768
1803 - 40739.913401969818129217228730508*y
1804 + x*(-251229.46597881387846350624383814
1805 + 441349.06185467302973318664458051*y
1806 + x*(1.3240471855640190891995599337415e6
1807 - 2.0596289553218074720882043413757e6*y
1808 + x*(-3.6043506718131630761543575974075e6
1809 + 4.6341651494740668121984597680953e6*y
1810 + x*(5.252053836070609053824921070508e6
1811 - 4.9431094927723379330116904193017e6*y
1812 + x*(-3.8838717443211226616520424723085e6
1813 + 1.1433042364235339436897787364371e6*x
1814 + 2.000782413741184401457112788765e6*y))))));
1815 grad[1]= -14.696938456699068589183704448235
1816 + x*(1131.6642611658282813671452425141
1817 + x*(-20369.956700984909064608614365254
1818 + x*(147116.35395155767657772888152684
1819 + x*(-514907.23883045186802205108534393
1820 + x*(926833.02989481336243969195361907
1821 + x*(-823851.58212872298883528173655028
1822 + 285826.05910588348592244468410928*x))))));
1823 return;
1824 case ( 38 ):
1825 grad[0]= -702.02564055738021170375436686406
1826 + (4155.2328454612504422465461173846
1827 - 4098.3118475782196142705660335848*y)*y
1828 + x*(24950.370738728512929471270065574
1829 + y*(-141391.75874144857669233452815868
1830 + 133195.13504629213746379339609151*y)
1831 + x*(-275839.15574116739237159948609378
1832 + (1.4429472963348314891910951243246e6
1833 - 1.2431545937653932829954050301874e6*y)*y
1834 + x*(1.3738831522367540845802392893142e6
1835 + y*(-6.3193691849740825218933089034526e6
1836 + 4.6618297266202248112327688632027e6*y)
1837 + x*(-3.5309043670142073107300045649072e6
1838 + (1.3286214720867640712013391260128e7
1839 - 7.458927562592359697972430181124e6*y)*y
1840 + x*(4.8534827264923896090278938053566e6
1841 + y*(-1.3177438693913168799751293319986e7
1842 + 4.2267256188023371621843771026371e6*y)
1843 + x*(-3.3841430608057928259319359155101e6
1844 + 939272.3597338527027076393561416*x
1845 + 4.9311798886027266892151066197433e6*y))))));
1846 grad[1]= -56.920997883030827975980083799789
1847 + 113.84199576606165595196016759958*y
1848 + x*(4155.2328454612504422465461173846
1849 - 8196.62369515643922854113206717*y
1850 + x*(-70695.879370724288346167264079338
1851 + 133195.13504629213746379339609151*y
1852 + x*(480982.43211161049639703170810822
1853 - 828769.7291769288553302700201249*y
1854 + x*(-1.5798422962435206304733272258631e6
1855 + 2.3309148633101124056163844316014e6*y
1856 + x*(2.6572429441735281424026782520255e6
1857 - 2.9835710250369438791889720724497e6*y
1858 + x*(-2.1962397823188614666252155533311e6
1859 + 704454.2698003895270307295171062*x
1860 + 1.4089085396007790540614590342124e6*y))))));
1861 return;
1862 case ( 39 ):
1863 grad[0]= -763.29810690188404265908474139258
1864 + y*(9024.8776168987466220280019423475
1865 + y*(-22225.44487743721183036746746996
1866 + 14592.463808418371403776620056035*y))
1867 + x*(24874.53830727316233136076157244
1868 + y*(-280444.70445348045473227313498459
1869 + (656660.87137882671316994790252155
1870 - 408588.98663571439930574536156897*y)*y)
1871 + x*(-251753.67561169787218761694970519
1872 + y*(2.6003770506601537841529936939854e6
1873 + y*(-5.5159513195821443906275623811811e6
1874 + 3.0644173997678579947930902117673e6*y))
1875 + x*(1.1469676553416839923368423364043e6
1876 + y*(-1.0296442463220002862504783111538e7
1877 + (1.8386504398607147968758541270604e7
1878 - 8.171779732714287986114907231379e6*y)*y)
1879 + x*(-2.6966873117957150354179193863552e6
1880 + y*(1.948969466252357684688405374684e7
1881 + y*(-2.5741106158050007156261957778845e7
1882 + 6.946012772807144788197671146672e6*y))
1883 + x*(3.3953744789427866582307439546381e6
1884 + x*(-2.1736934089020006043065653235469e6
1885 + 555681.0218245715830558136917338*x
1886 + 5.834650729158001622086043763205e6*y)
1887 + y*(-1.735686015228514768250806295945e7
1888 + 1.250282299105286061875580806401e7*y))))));
1889 grad[1]= -134.6996659238618898810149543634
1890 + (673.49832961930944940507477181698
1891 - 673.498329619309449405074771817*y)*y
1892 + x*(9024.8776168987466220280019423475
1893 + y*(-44450.889754874423660734934939921
1894 + 43777.391425255114211329860168104*y)
1895 + x*(-140222.3522267402273661365674923
1896 + (656660.87137882671316994790252155
1897 - 612883.4799535715989586180423535*y)*y
1898 + x*(866792.35022005126138433123132845
1899 + y*(-3.6773008797214295937517082541207e6
1900 + 3.0644173997678579947930902117673e6*y)
1901 + x*(-2.5741106158050007156261957778845e6
1902 + (9.1932521993035739843792706353018e6
1903 - 6.128834799535715989586180423535e6*y)*y
1904 + x*(3.8979389325047153693768107493679e6
1905 + y*(-1.0296442463220002862504783111538e7
1906 + 4.1676076636842868729186026880035e6*y)
1907 + x*(-2.8928100253808579470846771599083e6
1908 + 833521.5327368573745837205376007*x
1909 + 4.1676076636842868729186026880035e6*y))))));
1910 return;
1911 case ( 40 ):
1912 grad[0]= -763.67532368147132635291191107324
1913 + y*(15018.948032402269418273934251107
1914 + y*(-66439.753160288005392703336263372
1915 + (101568.81804963568640493728417274
1916 - 49893.454480522793321723578190118*y)*y))
1917 + x*(21892.025945535511355450141450766
1918 + y*(-407802.62284590568827245496051311
1919 + y*(1.7022322964859995864406406497822e6
1920 + y*(-2.4447792695456168727644553313158e6
1921 + 1.1226027258117628497387805092777e6*y)))
1922 + x*(-194737.20753877518821999253732368
1923 + y*(3.2830402165066452319911683057038e6
1924 + y*(-1.2220332529550904164299296400994e7
1925 + (1.5342237252760758946430000293461e7
1926 - 5.987214537662735198606829382814e6*y)*y))
1927 + x*(781239.85612614516685902888502792
1928 + y*(-1.1247410167180709694525686626287e7
1929 + y*(3.4319569046245321406299861283631e7
1930 + y*(-3.2929679957145043592337561605478e7
1931 + 8.48188726168887486469300829232e6*y)))
1932 + x*(-1.6251010887941709824789965467638e6
1933 + x*(1.8228929976276720560044007317318e6
1934 + x*(-1.0477625440909786597561951419925e6
1935 + 242339.63604825356756265737978057*x
1936 + 4.24094363084443743234650414616e6*y)
1937 + y*(-1.4326549072264402082380627451734e7
1938 + 1.6357925433257115810479373135189e7*y))
1939 + y*(1.8442759066907532531422822652419e7
1940 + y*(-4.0092954493277244633527875331345e7
1941 + 2.12047181542221871617325207308e7*y))))));
1942 grad[1]= -254.55844122715710878430397035775
1943 + y*(2291.0259710444139790587357332197
1944 + y*(-5345.7272657702992844703833775127
1945 + 3563.8181771801995229802555850084*y))
1946 + x*(15018.948032402269418273934251107
1947 + y*(-132879.50632057601078540667252674
1948 + (304706.45414890705921481185251822
1949 - 199573.81792209117328689431276047*y)*y)
1950 + x*(-203901.31142295284413622748025655
1951 + y*(1.7022322964859995864406406497822e6
1952 + y*(-3.6671689043184253091466829969737e6
1953 + 2.2452054516235256994775610185553e6*y))
1954 + x*(1.0943467388355484106637227685679e6
1955 + y*(-8.1468883530339361095328642673293e6
1956 + (1.5342237252760758946430000293461e7
1957 - 7.982952716883646931475772510419e6*y)*y)
1958 + x*(-2.8118525417951774236314216565717e6
1959 + y*(1.7159784523122660703149930641816e7
1960 + y*(-2.4697259967858782694253171204108e7
1961 + 8.48188726168887486469300829232e6*y))
1962 + x*(3.6885518133815065062845645304837e6
1963 + x*(-2.3877581787107336803967712419557e6
1964 + 605849.0901206339189066434494514*x
1965 + 5.452641811085705270159791045063e6*y)
1966 + y*(-1.6037181797310897853411150132538e7
1967 + 1.272283089253331229703951243848e7*y))))));
1968 return;
1969 case ( 41 ):
1970 grad[0]= -703.56236397351443318484451703167
1971 + y*(20684.733500821324335634428800731
1972 + y*(-141838.17257706050973006465463358
1973 + y*(370355.22839565799762850215376547
1974 + y*(-407784.74615904896547393588207156
1975 + 159567.94414919307344632273646278*y))))
1976 + x*(16744.784262569643509799299505354
1977 + y*(-460974.06087544665662271012755915
1978 + y*(2.9431420809740055768988415836469e6
1979 + y*(-7.1076684258060322098065732488607e6
1980 + (7.1805574867136883050845231408252e6
1981 - 2.5530871063870891751411637834045e6*y)*y)))
1982 + x*(-124108.40100492794601380657280439
1983 + y*(3.0317909388346683954801319927929e6
1984 + y*(-1.680782345038167040301266157408e7
1985 + y*(3.4159359895642072759990570990922e7
1986 + y*(-2.7658443652526799397362607653549e7
1987 + 7.233746801430085996233297386313e6*y))))
1988 + x*(419604.59387380400795144126995769
1989 + y*(-8.5457498977678957112363954416735e6
1990 + y*(3.7409818017199709441304552659608e7
1991 + y*(-5.4213701518343128163491379104392e7
1992 + 2.4112489338100286654110991287709e7*y)))
1993 + x*(-747605.36795825643670221578379785
1994 + x*(729875.59638612387298595770196865
1995 + x*(-370355.22839565799762850215376547
1996 + 76547.58520031837033051108345305*x
1997 + 2.0093741115083572211759159406424e6*y)
1998 + y*(-7.8010994917383280351535560048472e6
1999 + 1.2056244669050143327055495643855e7*y))
2000 + y*(1.1745973666537823462020979211844e7
2001 + y*(-3.5459543144265127432516163658396e7
2002 + 2.6791654820111429615678879208566e7*y))))));
2003 grad[1]= -422.137418384108659910906710219
2004 + y*(5909.923857377521238752693943066
2005 + y*(-23639.695429510084955010775772264
2006 + (35459.543144265127432516163658396
2007 - 17729.771572132563716258081829198*y)*y))
2008 + x*(20684.733500821324335634428800731
2009 + y*(-283676.34515412101946012930926717
2010 + y*(1.1110656851869739928855064612964e6
2011 + y*(-1.6311389846361958618957435282862e6
2012 + 797839.7207459653672316136823139*y)))
2013 + x*(-230487.03043772332831135506377958
2014 + y*(2.9431420809740055768988415836469e6
2015 + y*(-1.0661502638709048314709859873291e7
2016 + (1.436111497342737661016904628165e7
2017 - 6.382717765967722937852909458511e6*y)*y))
2018 + x*(1.0105969796115561318267106642643e6
2019 + y*(-1.1205215633587780268675107716053e7
2020 + y*(3.4159359895642072759990570990922e7
2021 + y*(-3.6877924870035732529816810204732e7
2022 + 1.2056244669050143327055495643855e7*y)))
2023 + x*(-2.1364374744419739278090988604184e6
2024 + x*(2.3491947333075646924041958423687e6
2025 + x*(-1.3001832486230546725255926674745e6
2026 + 287053.44450119388873941656294892*x
2027 + 4.0187482230167144423518318812849e6*y)
2028 + y*(-1.4183817257706050973006465463358e7
2029 + 1.607499289206685776940732752514e7*y))
2030 + y*(1.8704909008599854720652276329804e7
2031 + y*(-4.0660276138757346122618534328294e7
2032 + 2.4112489338100286654110991287709e7*y))))));
2033 return;
2034 case ( 42 ):
2035 grad[0]= -581.2882245495774706232175484286
2036 + y*(23771.628972369562877591580796264
2037 + y*(-231291.52513656871989008024558527
2038 + y*(899467.04219776724401697873283162
2039 + y*(-1.638314969717361765888068406229e6
2040 + (1.3993137270762407553349854857909e6
2041 - 452303.42693373438556282359136676*y)*y))))
2042 + x*(10493.782158973951180198085216369
2043 + y*(-393195.59273216682381313641749496
2044 + y*(3.4693728770485307983512036837791e6
2045 + y*(-1.2078557423798588705370857269453e7
2046 + y*(1.9370665230187630290794220567767e7
2047 + y*(-1.4247557948412633145228943128053e7
2048 + 3.8445791289367422772840005266174e6*y)))))
2049 + x*(-60392.787118992943526854286347266
2050 + y*(1.9467036698994533924081754003427e6
2051 + y*(-1.4262977383421737726554948477759e7
2052 + y*(3.8934073397989067848163508006854e7
2053 + y*(-4.3945389775948056779115246661202e7
2054 + 1.7300606080215340247778002369778e7*y))))
2055 + x*(164902.2910695906614031127676858
2056 + x*(-245211.84840867702247605729740291
2057 + x*(204949.99032934839345815443983806
2058 + x*(-90803.339498060312253142614933478
2059 + 16643.19969236684968521212349185*x
2060 + 611637.5886944817259315455383255*y)
2061 + y*(-2.7446594316206154760289522476119e6
2062 + 5.242607903095557650841818899933e6*y))
2063 + y*(4.8860334685100142076779451879889e6
2064 + y*(-1.8503322010925497591206419646822e7
2065 + 1.7475359676985192169472729666443e7*y)))
2066 + y*(-4.3302913317235365890531690423465e6
2067 + y*(2.4285610139339715588458425786454e7
2068 + y*(-4.5230342693373438556282359136676e7
2069 + 2.6213039515477788254209094499664e7*y))))));
2070 grad[1]= -642.47645871269088858355623773687
2071 + y*(12849.529174253817771671124754737
2072 + y*(-77097.175045522906630026748528424
2073 + y*(192742.93761380726657506687132106
2074 + y*(-212017.23137518799323257355845317
2075 + 84806.89255007519729302942338127*y))))
2076 + x*(23771.628972369562877591580796264
2077 + y*(-462583.05027313743978016049117055
2078 + y*(2.6984011265933017320509361984949e6
2079 + y*(-6.5532598788694470635522736249161e6
2080 + (6.9965686353812037766749274289545e6
2081 - 2.7138205616024063133769415482005e6*y)*y)))
2082 + x*(-196597.79636608341190656820874748
2083 + y*(3.4693728770485307983512036837791e6
2084 + y*(-1.811783613569788305805628590418e7
2085 + y*(3.8741330460375260581588441135533e7
2086 + y*(-3.5618894871031582863072357820132e7
2087 + 1.1533737386810226831852001579852e7*y))))
2088 + x*(648901.22329981779746939180011424
2089 + x*(-1.0825728329308841472632922605866e6
2090 + x*(977206.69370200284153558903759778
2091 + x*(-457443.23860343591267149204126865
2092 + 87376.79838492596084736364833221*x
2093 + 1.7475359676985192169472729666443e6*y)
2094 + y*(-7.4013288043701990364825678587287e6
2095 + 1.0485215806191115301683637799866e7*y))
2096 + y*(1.2142805069669857794229212893227e7
2097 + y*(-3.3922757020030078917211769352507e7
2098 + 2.6213039515477788254209094499664e7*y)))
2099 + y*(-9.5086515889478251510366323185057e6
2100 + y*(3.8934073397989067848163508006854e7
2101 + y*(-5.8593853034597409038820328881603e7
2102 + 2.8834343467025567079630003949631e7*y))))));
2103 return;
2104 case ( 43 ):
2105 grad[0]= -394.36024140371960168901824361658
2106 + y*(21163.999621999618623977312407423
2107 + y*(-273291.64729277768397048964282629
2108 + y*(1.4492738871586695362071420452909e6
2109 + y*(-3.7957173235108011662568005948096e6
2110 + y*(5.1925412985627759954393032136995e6
2111 + y*(-3.5527914148061098916163653567417e6
2112 + 958689.7468524423517060033502319*y))))))
2113 + x*(4600.8694830433953530385461755267
2114 + y*(-215320.69180643090252220396101465
2115 + y*(2.3602460448012618161087741880452e6
2116 + y*(-1.0213930252356337683745572509669e7
2117 + y*(2.0496873546958326297786723211972e7
2118 + y*(-1.913041531049443787793427499784e7
2119 + 6.710828227967096461942023451623e6*y)))))
2120 + x*(-19323.651828782260482761893937212
2121 + x*(41407.825347390558177346915579741
2122 + x*(-50609.564313477348883424007930794
2123 + x*(35886.781967738483753700660169109
2124 + x*(-13802.60844913018605911563852658
2125 + 2234.7080346210777429044367138273*x
2126 + 109500.69369643280940231739897754*y)
2127 + y*(-568667.4681041636656355643072951
2128 + 1.2670794556301510802268156167401e6*y))
2129 + y*(1.2008269350743261871430605518125e6
2130 + y*(-5.3416094698133820048777521097865e6
2131 + 5.866108590880329075124146373797e6*y)))
2132 + y*(-1.3066469331843242802629471138496e6
2133 + y*(8.6956433229520172172428522717455e6
2134 + y*(-1.9047599659799656761579581166681e7
2135 + 1.2905438899936723965273122022352e7*y))))
2136 + y*(759143.46470216023325136011896191
2137 + y*(-6.708067706277270424730200323918e6
2138 + y*(2.1946147434116995833993865257263e7
2139 + y*(-2.9606595123384249096803044639515e7
2140 + 1.3937874011931661882494971784141e7*y))))));
2141 grad[1]= -920.17389660867907060770923510535
2142 + y*(24844.695208434334906408149347844
2143 + y*(-207039.1267369527908867345778987
2144 + y*(759143.46470216023325136011896191
2145 + y*(-1.3664582364638884198524482141314e6
2146 + (1.1842638049353699638721217855806e6
2147 - 394754.60164512332129070726186019*y)*y))))
2148 + x*(21163.999621999618623977312407423
2149 + y*(-546583.29458555536794097928565258
2150 + y*(4.3478216614760086086214261358728e6
2151 + y*(-1.5182869294043204665027202379238e7
2152 + y*(2.5962706492813879977196516068497e7
2153 + y*(-2.131674848883665934969819214045e7
2154 + 6.710828227967096461942023451623e6*y)))))
2155 + x*(-107660.34590321545126110198050733
2156 + y*(2.3602460448012618161087741880452e6
2157 + y*(-1.5320895378534506525618358764504e7
2158 + y*(4.0993747093916652595573446423943e7
2159 + y*(-4.78260382762360946948356874946e7
2160 + 2.013248468390128938582607035487e7*y))))
2161 + x*(253047.82156738674441712003965397
2162 + x*(-326661.7332960810700657367784624
2163 + x*(240165.3870148652374286121103625
2164 + x*(-94777.911350693944272594051215851
2165 + 15642.956242347544200331056996791*x
2166 + 422359.81854338369340893853891335*y)
2167 + y*(-2.1366437879253528019511008439146e6
2168 + 3.519665154528197445074487824278e6*y))
2169 + y*(4.3478216614760086086214261358728e6
2170 + y*(-1.4285699744849742571184685875011e7
2171 + 1.2905438899936723965273122022352e7*y)))
2172 + y*(-4.472045137518180283153466882612e6
2173 + y*(2.1946147434116995833993865257263e7
2174 + y*(-3.9475460164512332129070726186019e7
2175 + 2.3229790019886103137491619640234e7*y))))));
2176 return;
2177 case ( 44 ):
2178 grad[0]= -139.94284547628721130097966906109
2179 + y*(8816.3992650060943119617191508489
2180 + y*(-132245.98897509141467942578726273
2181 + y*(808169.93262555864526315758882782
2182 + y*(-2.4245097978766759357894727664835e6
2183 + y*(3.7822352846876144598315775157142e6
2184 + y*(-2.9417385547570334687578936233333e6
2185 + 900532.2106399082047218041704081*y))))))
2186 + x*(979.59991833401047910685768342766
2187 + y*(-52898.395590036565871770314905094
2188 + y*(661229.94487545707339712893631367
2189 + y*(-3.2326797305022345810526303553113e6
2190 + y*(7.2735293936300278073684182994504e6
2191 + y*(-7.5644705693752289196631550314284e6
2192 + 2.9417385547570334687578936233333e6*y)))))
2193 + x*(-2938.799755002031437320573050283
2194 + x*(4897.9995916700523955342884171383
2195 + x*(-4897.9995916700523955342884171383
2196 + x*(2938.799755002031437320573050283
2197 + x*(-979.59991833401047910685768342766
2198 + 139.94284547628721130097966906109*x
2199 + 8816.399265006094311961719150849*y)
2200 + y*(-52898.395590036565871770314905094
2201 + 132245.98897509141467942578726273*y))
2202 + y*(132245.98897509141467942578726273
2203 + y*(-661229.94487545707339712893631367
2204 + 808169.9326255586452631575888278*y)))
2205 + y*(-176327.98530012188623923438301698
2206 + y*(1.3224598897509141467942578726273e6
2207 + y*(-3.2326797305022345810526303553113e6
2208 + 2.4245097978766759357894727664835e6*y))))
2209 + y*(132245.98897509141467942578726273
2210 + y*(-1.3224598897509141467942578726273e6
2211 + y*(4.8490195957533518715789455329669e6
2212 + y*(-7.2735293936300278073684182994504e6
2213 + 3.7822352846876144598315775157142e6*y))))));
2214 grad[1]= -1259.4856092865849017088170215498
2215 + y*(44081.996325030471559808595754245
2216 + y*(-484901.95957533518715789455329669
2217 + y*(2.4245097978766759357894727664835e6
2218 + y*(-6.303725474479357433052629192857e6
2219 + y*(8.8252156642711004062736808699998e6
2220 + y*(-6.303725474479357433052629192857e6
2221 + 1.8010644212798164094436083408163e6*y))))))
2222 + x*(8816.3992650060943119617191508489
2223 + y*(-264491.97795018282935885157452547
2224 + y*(2.4245097978766759357894727664835e6
2225 + y*(-9.6980391915067037431578910659338e6
2226 + y*(1.8911176423438072299157887578571e7
2227 + y*(-1.765043132854220081254736174e7
2228 + 6.303725474479357433052629192857e6*y)))))
2229 + x*(-26449.197795018282935885157452547
2230 + x*(44081.996325030471559808595754245
2231 + x*(-44081.996325030471559808595754245
2232 + x*(26449.197795018282935885157452547
2233 + x*(-8816.3992650060943119617191508489
2234 + 1259.4856092865849017088170215498*x
2235 + 44081.996325030471559808595754245*y)
2236 + y*(-264491.97795018282935885157452547
2237 + 484901.95957533518715789455329669*y))
2238 + y*(661229.94487545707339712893631367
2239 + y*(-2.4245097978766759357894727664835e6
2240 + 2.4245097978766759357894727664835e6*y)))
2241 + y*(-881639.92650060943119617191508489
2242 + y*(4.8490195957533518715789455329669e6
2243 + y*(-9.6980391915067037431578910659338e6
2244 + 6.303725474479357433052629192857e6*y))))
2245 + y*(661229.94487545707339712893631367
2246 + y*(-4.8490195957533518715789455329669e6
2247 + y*(1.4547058787260055614736836598901e7
2248 + y*(-1.8911176423438072299157887578571e7
2249 + 8.82521566427110040627368087e6*y))))));
2250 return;
2251 #endif
2252 }
2253 }
2254
2255 /* H\phi_i(x,y) for triangle */
2256 static void
2257 hess_triangle_2d ( const int i, DomainType xi, HessianRangeType &h )
2258 {
2259 const RangeField &x = xi[0];
2260 const RangeField &y = xi[1];
2261
2262 switch (i) {
2263 #if (PMAX2D>=0)
2264 case ( 0 ):
2265 h[0] = 0;
2266 h[1] = 0;
2267 h[2] = 0;
2268 return;
2269 #endif
2270 #if (PMAX2D>=1)
2271 case ( 1 ):
2272 h[0] = 0;
2273 h[1] = 0;
2274 h[2] = 0;
2275 return;
2276 case ( 2 ):
2277 h[0] = 0;
2278 h[1] = 0;
2279 h[2] = 0;
2280 return;
2281 #endif
2282 #if (PMAX2D>=2)
2283 case ( 3 ):
2284 h[0] = 48.989794855663561963945681494118;
2285 h[1] = 0;
2286 h[2] = 0;
2287 return;
2288 case ( 4 ):
2289 h[0] = 42.426406871192851464050661726290;
2290 h[1] = 42.426406871192851464050661726291;
2291 h[2] = 0;
2292 return;
2293 case ( 5 ):
2294 h[0] = 10.954451150103322269139395656016;
2295 h[1] = 32.863353450309966807418186968048;
2296 h[2] = 65.726706900619933614836373936096;
2297 return;
2298 #endif
2299 #if (PMAX2D>=3)
2300 case ( 6 ):
2301 h[0] = -254.55844122715710878430397035774+593.96969619669992049670926416808*x;
2302 h[1] = 0;
2303 h[2] = 0;
2304 return;
2305 case ( 7 ):
2306 h[0] = -323.33264604737950896204149786118+617.27141518136088074571558682590*x+411.51427678757392049714372455058*y;
2307 h[1] = -117.57550765359254871346963558588+411.51427678757392049714372455058*x;
2308 h[2] = 0;
2309 return;
2310 case ( 8 ):
2311 h[0] = -189.73665961010275991993361266596+265.63132345414386388790705773235*x+531.26264690828772777581411546470*y;
2312 h[1] = -303.57865537616441587189378026554+531.26264690828772777581411546470*y+531.26264690828772777581411546470*x;
2313 h[2] = -75.894663844041103967973445066386+531.26264690828772777581411546470*x;
2314 return;
2315 case ( 9 ):
2316 h[0] = -44.899888641287296627004984787798+44.899888641287296627004984787798*x+179.5995545651491865080199391512*y;
2317 h[1] = 179.5995545651491865080199391512*x-179.59955456514918650801993915119+448.99888641287296627004984787798*y;
2318 h[2] = -448.99888641287296627004984787798+897.99777282574593254009969575596*y+448.99888641287296627004984787798*x;
2319 return;
2320 #endif
2321 #if (PMAX2D>=4)
2322 case ( 10 ):
2323 h[0] = 796.89397036243159166372117319704+(-4250.1011752663018222065129237177+4781.3638221745895499823270391822*x)*x;
2324 h[1] = 0;
2325 h[2] = 0;
2326 return;
2327 case ( 11 ):
2328 h[0] = 1150.2173707608488382596365438817-1840.3477932173581412154184702106*y+(-5521.0433796520744236462554106321+5521.0433796520744236462554106318*y+5521.0433796520744236462554106320*x)*x;
2329 h[1] = 230.04347415216976765192730877634+(-1840.3477932173581412154184702106+2760.5216898260372118231277053159*x)*x;
2330 h[2] = 0;
2331 return;
2332 case ( 12 ):
2333 h[0] = 975.80735803743558367316521970470+(-4412.3463146040565522612688195342+3054.7012947258853054116476442930*y)*y+(-3733.5238046649709288364582319135+9164.1038841776559162349429328790*y+3054.7012947258853054116476442930*x)*x;
2334 h[1] = -1357.6450198781712468496211752414*y+721.24891681027847488886124934695+(-4412.3463146040565522612688195342+6109.4025894517706108232952885860*y+4582.0519420888279581174714664395*x)*x;
2335 h[2] = 84.852813742385702928101323452582+x*(-1357.6450198781712468496211752414+3054.7012947258853054116476442930*x);
2336 return;
2337 case ( 13 ):
2338 h[0] = 501.99601592044532878690321547112+(-3815.1697209953844987804644375804+4517.96414328400795908212893924*y)*y+(-1405.5888445772469206033290033191+5421.556971940809550898554727088*y+903.59282865680159181642578784804*x)*x;
2339 h[1] = 1104.3912350249797233311870740364+(-5019.9601592044532878690321547112+4517.9641432840079590821289392399*y)*y+(-3815.1697209953844987804644375804+9035.92828656801591816425787848*y+2710.778485970404775449277363544*x)*x;
2340 h[2] = -1003.9920318408906575738064309422*y+501.99601592044532878690321547112+x*(-5019.9601592044532878690321547112+9035.9282865680159181642578784798*y+4517.96414328400795908212893924*x);
2341 return;
2342 case ( 14 ):
2343 h[0] = 113.84199576606165595196016759958+(-1138.4199576606165595196016759958+1707.6299364909248392794025139936*y)*y+(-227.68399153212331190392033519916+1138.4199576606165595196016759958*y+113.84199576606165595196016759958*x)*x;
2344 h[1] = 569.20997883030827975980083799789+(-3415.2598729818496785588050279874+3984.4698518121579583186058659851*y)*y+(-1138.4199576606165595196016759958+3415.2598729818496785588050279872*y+569.20997883030827975980083799788*x)*x;
2345 h[2] = 1707.6299364909248392794025139937+(-7968.9397036243159166372117319702+7968.9397036243159166372117319708*y)*y+x*(1707.6299364909248392794025139936*x-3415.2598729818496785588050279874+7968.9397036243159166372117319702*y);
2346 return;
2347 #endif
2348 #if (PMAX2D>=5)
2349 case ( 15 ):
2350 h[0] = -1939.8969044771425687507399024866+(17459.072140294283118756659122380+(-43647.680350735707796891647805948+32008.298923872852384387208391028*x)*x)*x;
2351 h[1] = 0;
2352 h[2] = 0;
2353 return;
2354 case ( 16 ):
2355 h[0] = -2976.+5184.*y+(25056.-34560.*y+(-58320.+47520.*y+39600.*x)*x)*x;
2356 h[1] = -384.+(5184.+(-17280.+15840.*x)*x)*x;
2357 h[2] = 0;
2358 return;
2359 case ( 17 ):
2360 h[0] = -2943.4673431176368327362417038346+(15058.159250054436849576983874354-12548.466041712030707980819895294*y)*y+(21471.819671373919211433847376394+(-83656.440278080204719872132635302+46011.042152944112595929672949414*y)*y+(-43222.494143674772438600601861570+92022.084305888225191859345898828*y+25561.690084968951442183151638564*x)*x)*x;
2361 h[1] = -1301.3224043256920734202331743269+2509.6932083424061415961639790590*y+(-25096.932083424061415961639790588*y+15058.159250054436849576983874354+(-41828.220139040102359936066317651+46011.042152944112595929672949414*y+30674.028101962741730619781966276*x)*x)*x;
2362 h[2] = -92.951600308978005244302369594778+x*(2509.6932083424061415961639790590+x*(-12548.466041712030707980819895294+15337.014050981370865309890983138*x));
2363 return;
2364 case ( 18 ):
2365 h[0] = -2162.9757280191564831095582754396+(21116.508802356510750357721468698+(-41243.181254602560059292424743552+20163.333057805696028987407652404*y)*y)*y+(12427.945284720238097866783989391+(-85785.817009573324923328243466586+90734.998760125632130443334435812*y)*y+(-20346.636085603929629250929540152+72587.999008100505704354667548652*y+10081.666528902848014493703826202*x)*x)*x;
2366 h[1] = -2419.5999669366835234784889182884+(11548.090751288716816601878928195-10998.181667894016015811313264947*y)*y+(21116.508802356510750357721468698+(-82486.362509205120118584849487104+60489.999173417088086962222957212*y)*y+(-42892.908504786662461664121733293+90734.998760125632130443334435812*y+24195.999669366835234784889182884*x)*x)*x;
2367 h[2] = -549.90908339470080079056566324736+1099.8181667894016015811313264947*y+x*(-21996.363335788032031622626529895*y+11548.090751288716816601878928195+x*(-41243.181254602560059292424743552+60489.999173417088086962222957212*y+30244.999586708544043481111478604*x));
2368 return;
2369 case ( 19 ):
2370 h[0] = -1039.2304845413263761164678049035+(14964.918977395099816077136390611+(-43024.142060010911971221767123006+32008.298923872852384387208391028*y)*y)*y+(4364.7680350735707796891647805947+(-42400.603769286116145551886440064+61730.290781754786741318187611270*y)*y+(-5611.8446165231624310289261464792+27435.684791891016329474750049454*y+2286.3070659909180274562291707878*x)*x)*x;
2371 h[1] = -2909.8453567157138531261098537299+(24317.993338267037201125346634742+(-52377.216420882849356269977367137+32008.298923872852384387208391028*y)*y)*y+(14964.918977395099816077136390611+(-86048.284120021823942443534246012+96024.896771618557153161625173084*y)*y+(-21200.301884643058072775943220032+61730.290781754786741318187611270*y+9145.2282639636721098249166831512*x)*x)*x;
2372 h[2] = -1870.6148721743874770096420488263+(8729.5360701471415593783295611894-8729.5360701471415593783295611900*y)*y+(24317.993338267037201125346634742+(-104754.43284176569871253995473427+96024.896771618557153161625173084*y)*y+x*(20576.763593918262247106062537090*x-43024.142060010911971221767123006+96024.896771618557153161625173084*y))*x;
2373 return;
2374 case ( 20 ):
2375 h[0] = -229.78250586152114639402445872876+(4136.0851055073806350924402571176+(-14476.297869275832222823540899912+12867.820328245184198065369688810*y)*y)*y+(689.34751758456343918207337618628+(-8272.1702110147612701848805142350+14476.297869275832222823540899912*y)*y+(-689.34751758456343918207337618628+4136.0851055073806350924402571178*y+229.78250586152114639402445872876*x)*x)*x;
2376 h[1] = -1378.6950351691268783641467523725+(14476.297869275832222823540899912+(-38603.46098473555259419610906643+28952.595738551664445647081799824*y)*y)*y+(4136.0851055073806350924402571176+(-28952.595738551664445647081799824+38603.460984735552594196109066430*y)*y+(-4136.0851055073806350924402571175+14476.297869275832222823540899912*y+1378.6950351691268783641467523725*x)*x)*x;
2377 h[2] = -4825.432623091944074274513633304+(38603.460984735552594196109066430+(-86857.787215654993336941245399472+57905.191477103328891294163599647*y)*y)*y+(14476.297869275832222823540899912+(-77206.92196947110518839221813286+86857.787215654993336941245399472*y)*y+x*(4825.432623091944074274513633304*x-14476.297869275832222823540899912+38603.460984735552594196109066430*y))*x;
2378 return;
2379 #endif
2380 #if (PMAX2D>=6)
2381 case ( 21 ):
2382 h[0] = 4040.9899777158566964304486309018+(-53879.866369544755952405981745357+(222254.44877437211830367467469960+(-355607.11803899538928587947951937+192620.52227112250252985138473965*x)*x)*x)*x;
2383 h[1] = 0;
2384 h[2] = 0;
2385 return;
2386 case ( 22 ):
2387 h[0] = 6415.9332914237816286563077617272-11665.333257134148415738741384958*y+(-81657.332799939038910171189694710+128318.66582847563257312615523454*y+(320796.66457118908143281538808636-384955.99748542689771937846570362*y+(-487610.93014820740377787938989127+333628.53115403664469012800360982*y+250221.39836552748351759600270735*x)*x)*x)*x;
2388 h[1] = 583.2666628567074207869370692479+(-11665.333257134148415738741384958+(64159.33291423781628656307761727+(-128318.66582847563257312615523454+83407.132788509161172532000902455*x)*x)*x)*x;
2389 h[2] = 0;
2390 return;
2391 case ( 23 ):
2392 h[0] = 6877.3454181101010043805740519542+(-37147.705178112954330230837944862+33131.737050749391699935612221094*y)*y+(-79315.370515430361948330708044436+(364449.10755824330869929173443203-265053.89640599513359948489776875*y)*y+(281619.76493136982944945270387930+(-960820.37447173235929813275441170+430712.58165974209209916295887422*y)*y+(-386536.93225874290316591547591277+717854.30276623682016527159812371*y+179463.57569155920504131789953092*x)*x)*x)*x;
2393 h[1] = -4015.9681273635626302952257237688*y+2058.1836652738258480263031834316+(-37147.705178112954330230837944862+66263.474101498783399871224442188*y+(-265053.89640599513359948489776875*y+182224.55377912165434964586721602+(-320273.45815724411976604425147056+287141.72110649472806610863924948*y+179463.57569155920504131789953092*x)*x)*x)*x;
2394 h[2] = 100.39920318408906575738064309422+x*(-4015.9681273635626302952257237688+x*(33131.737050749391699935612221094+x*(-88351.298801998377866494965922916+71785.430276623682016527159812370*x)));
2395 return;
2396 case ( 24 ):
2397 h[0] = 5939.6969619669992049670926416808+(-62960.787796850191572651182001816+(137206.99982143768163473984002282-78403.999897964389505565622870186*y)*y)*y+(-58209.030227276592208677507888471+(509625.99933676853178617654865620+(-862443.99887760828456122185157202+339750.66622451235452411769910416*y)*y)*y+(176408.99977041987638752265145792+(-1097655.9985715014530779187201826+1019251.9986735370635723530973124*y)*y+(-209077.33306123837201484166098716+679501.33244902470904823539820829*y+84937.666556128088631029424776035*x)*x)*x)*x;
2398 h[1] = 4157.7878733768994434769648491765+(-20194.969670687797296888114981714+19600.999974491097376391405717547*y)*y+(-62960.787796850191572651182001816+(274413.99964287536326947968004564-235211.99969389316851669686861056*y)*y+(254812.99966838426589308827432810+(-862443.99887760828456122185157202+509625.99933676853178617654865624*y)*y+(-365885.33285716715102597290672753+679501.33244902470904823539820824*y+169875.33311225617726205884955206*x)*x)*x)*x;
2399 h[2] = -1187.9393923933998409934185283362*y+593.96969619669992049670926416808+x*(-20194.969670687797296888114981714+39201.999948982194752782811435094*y+x*(-235211.99969389316851669686861056*y+137206.99982143768163473984002282+x*(-287481.33295920276152040728385734+339750.66622451235452411769910416*y+169875.33311225617726205884955206*x)));
2400 return;
2401 case ( 25 ):
2402 h[0] = 4040.9899777158566964304486309018+(-68696.829621169563839317626725332+(256602.86358495690022333348806226+(-320585.20489879129791681559138488+122576.69599071431979172360847070*y)*y)*y)*y+(-30980.923162488234672633439503582+(413527.97438625600193471590989563+(-1091067.2939832813080362211303435+735460.17594428591875034165082414*y)*y)*y+(76105.311246981967782773449215318+(-695050.27616712735178603716451512+945591.65478551046696472497963108*y)*y+(-75431.812917362658333368374443500+350219.13140204091369063888134482*y+26266.434855153068526797916100863*x)*x)*x)*x;
2403 h[1] = 6061.4849665737850446456729463528+(-52532.869710306137053595832201724+(117862.20768337915364588808506797-75431.812917362658333368374443500*y)*y)*y+(-68696.829621169563839317626725332+(513205.72716991380044666697612452+(-961755.61469637389375044677415464+490306.78396285727916689443388280*y)*y)*y+(206763.98719312800096735795494781+(-1091067.2939832813080362211303435+1103190.2639164288781255124762362*y)*y+(-231683.42538904245059534572150504+630394.43652367364464314998642072*y+87554.782850510228422659720336205*x)*x)*x)*x;
2404 h[2] = 2020.4949888579283482152243154510+(-9428.976614670332291671046805438+9428.9766146703322916710468054380*y)*y+(-52532.869710306137053595832201724+(235724.41536675830729177617013595-226295.43875208797500010512333050*y)*y+(256602.86358495690022333348806226+(-961755.61469637389375044677415464+735460.17594428591875034165082420*y)*y+x*(157598.60913091841116078749660518*x-363689.09799442710267874037678116+735460.17594428591875034165082414*y))*x)*x;
2405 return;
2406 case ( 26 ):
2407 h[0] = 1861.4510468986284894199862933174+(-43185.664288048180954543682004964+(218906.64311527871035579038809412+(-375268.5310547635034670692367328+203270.45432133023104466250323026*y)*y)*y)*y+(-10424.125862632319540751923242577+(183166.78301482504335892665126244+(-625447.55175793917244511539455468+542054.54485688061611910000861402*y)*y)*y+(20103.671306505187685735851967828+(-236776.57316550554385422225650998+406540.90864266046208932500646052*y)*y+(-16380.769212707930706895879381193+96795.454438728681449839287252507*y+4839.7727219364340724919643626252*x)*x)*x)*x;
2408 h[1] = 6328.9335594553368640279533972792+(-83393.00690105855632601538594062+(312723.77587896958622255769727734+(-437813.28623055742071158077618824+203270.45432133023104466250323026*y)*y)*y)*y+(-43185.664288048180954543682004964+(437813.28623055742071158077618824+(-1125805.5931642905104012077101984+813081.81728532092417865001292104*y)*y)*y+(91583.391507412521679463325631218+(-625447.55175793917244511539455468+813081.81728532092417865001292103*y)*y+(-78925.524388501847951407418836660+271027.27242844030805955000430701*y+24198.863609682170362459821813126*x)*x)*x)*x;
2409 h[2] = 5212.0629313161597703759616212888+(-41696.503450529278163007692970309+(93817.132763690875866767309183196-62544.755175793917244511539455467*y)*y)*y+(-83393.00690105855632601538594062+(625447.55175793917244511539455468+(-1313439.8586916722621347423285647+813081.81728532092417865001292104*y)*y)*y+(218906.64311527871035579038809412+(-1125805.5931642905104012077101984+1219622.7259279813862679750193816*y)*y+x*(67756.818107110077014887501076754*x-208482.51725264639081503846485156+542054.54485688061611910000861402*y))*x)*x;
2410 return;
2411 case ( 27 ):
2412 h[0] = 404.72212689696124396650916834486+(-11332.219553114914831062256713656+(67993.317318689488986373540281938+(-135986.63463737897797274708056388+84991.646648361861232966925352422*y)*y)*y)*y+(-1618.8885075878449758660366733794+(33996.658659344744493186770140969+(-135986.63463737897797274708056388+135986.63463737897797274708056388*y)*y)*y+(2428.3327613817674637990550100692+(-33996.658659344744493186770140969+67993.317318689488986373540281936*y)*y+(-1618.8885075878449758660366733794+11332.219553114914831062256713656*y+404.72212689696124396650916834488*x)*x)*x)*x;
2413 h[1] = 2833.0548882787287077655641784141+(-45328.878212459659324249026854626+(203979.95195606846695912062084582+(-339966.58659344744493186770140968+186981.62262639609471252723577534*y)*y)*y)*y+(-11332.219553114914831062256713656+(135986.63463737897797274708056388+(-407959.90391213693391824124169164+339966.58659344744493186770140968*y)*y)*y+(16998.329329672372246593385070484+(-135986.63463737897797274708056388+203979.95195606846695912062084582*y)*y+(-11332.219553114914831062256713656+45328.878212459659324249026854624*y+2833.0548882787287077655641784140*x)*x)*x)*x;
2414 h[2] = 11332.219553114914831062256713656+(-135986.63463737897797274708056388+(509949.87989017116739780155211452+(-747926.49050558437885010894310131+373963.24525279218942505447155066*y)*y)*y)*y+(-45328.878212459659324249026854626+(407959.90391213693391824124169164+(-1019899.7597803423347956031042290+747926.49050558437885010894310131*y)*y)*y+(67993.317318689488986373540281938+(-407959.90391213693391824124169164+509949.87989017116739780155211452*y)*y+x*(11332.219553114914831062256713656*x-45328.878212459659324249026854626+135986.63463737897797274708056388*y))*x)*x;
2415 return;
2416 #endif
2417 #if (PMAX2D>=7)
2418 case ( 28 ):
2419 h[0] = -7560.+(138600.+(-831600.+(2162160.+(-2522520.+1081080.*x)*x)*x)*x)*x;
2420 h[1] = 0;
2421 h[2] = 0;
2422 return;
2423 case ( 29 ):
2424 h[0] = -12262.919717587651238174320097862+22863.070659909180274562291707878*y+(217199.17126913721260834177122484-365809.13055854688439299666732606*y+(-1257468.8862950049151009260439333+1783319.5114729160614158587532144*y+(3150531.1369354850418346837973457-3328863.0880827766479762696726670*y+(-3536917.0310879501884747865272086+2080539.4300517354049851685454170*y+1456377.6010362147834896179817918*x)*x)*x)*x)*x;
2425 h[1] = -831.3843876330611008931742439228+(22863.070659909180274562291707878+(-182904.56527927344219649833366303+(594439.8371576386871386195844048+(-832215.77202069416199406741816675+416107.8860103470809970337090834*x)*x)*x)*x)*x;
2426 h[2] = 0;
2427 return;
2428 case ( 30 ):
2429 h[0] = -13792.067285218702847451783188735+(76741.852987792782380762840310858-70838.633527193337582242621825406*y)*y+(227273.94923307862474302841168985+(-1133418.1364350934013158819492065+920902.23585351338856915408373026*y)*y+(-1221966.4283440850732936852264882+(5064962.2971943236371303474605166-3223157.825487296859992039293056*y)*y+(2839448.5605483329480882250915018+(-8595087.5346327916266454381148160+3223157.8254872968599920392930560*y)*y+(-2954561.3400300221216593693519680+4834736.7382309452899880589395840*y+1128105.2389205539009972137525696*x)*x)*x)*x)*x;
2430 h[1] = -3005.2753617597173519739294107748+5903.2194605994447985202184854506*y+(-141677.26705438667516448524365081*y+76741.852987792782380762840310858+(-566709.06821754670065794097460326+920902.23585351338856915408373026*y+(-2148771.883658197906661359528704*y+1688320.7657314412123767824868389+(-2148771.8836581979066613595287040+1611578.9127436484299960196465280*y+966947.3476461890579976117879168*x)*x)*x)*x)*x;
2431 h[2] = -107.33126291998990542764033609910+x*(5903.2194605994447985202184854506+x*(-70838.633527193337582242621825406+x*(306967.41195117112952305136124342+x*(-537192.970914549476665339882176+322315.7825487296859992039293056*x))));
2432 return;
2433 case ( 31 ):
2434 h[0] = -13017.096450437785705267949507905+(143505.55111214339362880763847739+(-327649.84236223889872772009493068+198113.85817251654341676098763250*y)*y)*y+(190811.58455397827338697652815247+(-1859222.3613113090997572954223974+(3665106.3761915560532100782712014-1849062.6762768210718897692179034*y)*y)*y+(-912466.71215995600285219724111508+(7231155.8232968538347117760485866+(-10747676.805859022480359283579064+3466992.5180190395097933172835689*y)*y)*y+(1890336.3967294286851015944236602+(-10632110.388591721163366173002945+8667481.2950475987744832932089222*y)*y+(-1762387.8633263450841449362858142+5200488.7770285592646899759253532*y+606723.69065333191421383052462454*x)*x)*x)*x)*x;
2435 h[1] = -6349.8031465550174172038778087342+(31114.035418119585344299001262798-30479.055103464083602578613481925*y)*y+(143505.55111214339362880763847739+(-655299.68472447779745544018986136+594341.57451754963025028296289750*y)*y+(-929611.18065565454987864771119869+(3665106.3761915560532100782712014-2773594.0144152316078346538268551*y)*y+(2410385.2744322846115705920161956+(-7165117.8705726816535728557193756+3466992.5180190395097933172835689*y)*y+(-2658027.5971479302908415432507361+4333740.6475237993872416466044610*y+1040097.7554057118529379951850706*x)*x)*x)*x)*x;
2436 h[2] = -634.98031465550174172038778087342+1269.9606293110034834407755617468*y+x*(-60958.110206928167205157226963850*y+31114.035418119585344299001262798+x*(-327649.84236223889872772009493068+594341.57451754963025028296289750*y+x*(-1849062.6762768210718897692179034*y+1221702.1253971853510700260904005+x*(-1791279.4676431704133932139298439+1733496.2590095197548966586417844*y+866748.12950475987744832932089222*x))));
2437 return;
2438 case ( 32 ):
2439 h[0] = -10440.+(188640.+(-760320.+(1048320.-458640.*y)*y)*y)*y+(128520.+(-2004480.+(.67392e7+(-.733824e7+2293200.*y)*y)*y)*y+(-520560.+(6402240.+(-.1533168e8+9172800.*y)*y)*y+(926640.+(-7862400.+9828000.*y)*y+(-753480.+3276000.*y+229320.*x)*x)*x)*x)*x;
2440 h[1] = -10080.+(88560.+(-201600.+131040.*y)*y)*y+(188640.+(-1520640.+(3144960.-1834560.*y)*y)*y+(-1002240.+(.67392e7+(-.1100736e8+4586400.*y)*y)*y+(2134080.+(-.1022112e8+9172800.*y)*y+(-1965600.+4914000.*y+655200.*x)*x)*x)*x)*x;
2441 h[2] = -2160.+(10080.-10080.*y)*y+(88560.+(-403200.+393120.*y)*y+(-760320.+(3144960.-2751840.*y)*y+(.22464e7+(-.733824e7+4586400.*y)*y+x*(982800.*x-.255528e7+4586400.*y))*x)*x)*x;
2442 return;
2443 case ( 33 ):
2444 h[0] = -6765.9145723250156921944627828082+(177505.75877982099992463120006662+(-1069810.4923770377753305127035405+(2407073.6078483349944936535829662+(-2223199.9294710316268587217120450+702063.13562243104006064896169846*y)*y)*y)*y)*y+(64873.180899351621048688084329276+(-1413678.1506410856316867489296785+(6686315.5773564860958157043971284+(-10608954.049405624605360917643443+5265473.5171682328004548672127386*y)*y)*y)*y+(-212529.31656597402233128488976585+(3548523.1956970494065650488336187+(-11467031.215166373654323933041075+9360841.8082990805341419861559788*y)*y)*y+(316008.01002506249762367079114998+(-3566034.9745901259177683756784684+5850526.1301869253338387413474870*y)*y+(-220091.22108798433398726693640547+1253684.1707543411429654445744615*y+58505.261301869253338387413474871*x)*x)*x)*x)*x;
2445 h[1] = -12735.839194964735420601341708815+(172729.81908170922414190569692581+(-668631.55773564860958157043971281+(969515.75871669048389327713758356-468042.09041495402670709930779895*y)*y)*y)*y+(177505.75877982099992463120006662+(-2139620.9847540755506610254070810+(7221220.8235450049834809607488986+(-8892799.7178841265074348868481800+3510315.6781121552003032448084923*y)*y)*y)*y+(-706839.07532054281584337446483927+(6686315.5773564860958157043971284+(-15913431.074108436908041376465165+10530947.034336465600909734425477*y)*y)*y+(1182841.0652323498021883496112062+(-7644687.4767775824362159553607164+9360841.8082990805341419861559788*y)*y+(-891508.74364753147944209391961710+2925263.0650934626669193706737435*y+250736.83415086822859308891489230*x)*x)*x)*x)*x;
2446 h[2] = -5571.9296477970717465130869976068+(44575.437182376573972104695980855+(-100294.73366034729143723556595692+66863.155773564860958157043971280*y)*y)*y+(172729.81908170922414190569692581+(-1337263.1154712972191631408794256+(2908547.2761500714516798314127506-1872168.3616598161068283972311958*y)*y)*y+(-1069810.4923770377753305127035405+(7221220.8235450049834809607488986+(-13339199.576826189761152330272270+7020631.3562243104006064896169847*y)*y)*y+(2228771.8591188286986052347990428+(-10608954.049405624605360917643443+10530947.034336465600909734425477*y)*y+x*(585052.61301869253338387413474870*x-1911171.8691943956090539888401791+4680420.9041495402670709930779894*y))*x)*x)*x;
2447 return;
2448 case ( 34 ):
2449 h[0] = -3028.6630713897510062201458646752+(102974.54442725153421148495939896+(-799567.05084689426564211850827426+(2326013.2388273287727770720240706+(-2816656.656392468435784735654148+1199350.5762703413984631777624114*y)*y)*y)*y)*y+(21200.641499728257043541021052726+(-581503.30970683219319426800601764+(3416331.9445276391350163245353536+(-6687288.0616285702217340820692030+4088695.1463761638583971969173114*y)*y)*y)*y+(-54515.935285015518111962625564154+(1126662.6625569873743138942616592+(-4433962.7365145954731062935458846+4361274.8228012414489570100451324*y)*y)*y+(66630.587570574522136843209022853+(-920713.57370248430589092434286126+1817197.8428338506037320875188052*y)*y+(-39372.619928066763080861896240780+272579.67642507759055981312782078*y+9085.9892141692530186604375940259*x)*x)*x)*x)*x;
2450 h[1] = -12114.652285559004024880583458701+(230178.39342562107647273108571532+(-1308382.4468403724346871030135397+(3089236.3328175460263445487819688+(-3198268.2033875770625684740330970+1199350.5762703413984631777624114*y)*y)*y)*y)*y+(102974.54442725153421148495939896+(-1599134.1016937885312842370165485+(6978039.7164819863183312160722118+(-11266626.625569873743138942616592+5996752.8813517069923158888120570*y)*y)*y)*y+(-290751.65485341609659713400300882+(3416331.9445276391350163245353536+(-10030932.092442855332601123103804+8177390.2927523277167943938346228*y)*y)*y+(375554.22085232912477129808721972+(-2955975.1576763969820708623639230+4361274.8228012414489570100451323*y)*y+(-230178.39342562107647273108571531+908598.92141692530186604375940260*y+54515.935285015518111962625564156*x)*x)*x)*x)*x;
2451 h[2] = -12114.652285559004024880583458701+(145375.82742670804829856700150441+(-545159.35285015518111962625564154+(799567.05084689426564211850827427-399783.52542344713282105925413713*y)*y)*y)*y+(230178.39342562107647273108571532+(-2616764.8936807448693742060270794+(9267708.9984526380790336463459064+(-12793072.813550308250273896132388+5996752.8813517069923158888120569*y)*y)*y)*y+(-799567.05084689426564211850827426+(6978039.7164819863183312160722118+(-16899939.938354810614708413924888+11993505.762703413984631777624114*y)*y)*y+(1138777.3148425463783387748451179+(-6687288.0616285702217340820692030+8177390.2927523277167943938346228*y)*y+x*(181719.78428338506037320875188052*x-738993.78941909924551771559098076+2180637.4114006207244785050225662*y))*x)*x)*x;
2452 return;
2453 case ( 35 ):
2454 h[0] = -650.66120216284603671011658716344+(26026.448086513841468404663486538+(-234238.03277862457321564197137884+(780793.44259541524405213990459614+(-1073590.9835686959605716923688197+515323.67211297406107441233703344*y)*y)*y)*y)*y+(3253.3060108142301835505829358173+(-104105.79234605536587361865394615+(702714.09833587371964692591413652+(-1561586.8851908304881042798091923+1073590.9835686959605716923688197*y)*y)*y)*y+(-6506.6120216284603671011658716346+(156158.68851908304881042798091923+(-702714.09833587371964692591413652+780793.44259541524405213990459608*y)*y)*y+(6506.6120216284603671011658716343+(-104105.79234605536587361865394615+234238.03277862457321564197137884*y)*y+(-3253.3060108142301835505829358173+26026.448086513841468404663486539*y+650.66120216284603671011658716346*x)*x)*x)*x)*x;
2455 h[1] = -5205.2896173027682936809326973075+(117119.01638931228660782098568942+(-780793.44259541524405213990459614+(2147181.9671373919211433847376393+(-2576618.3605648703053720616851672+1116534.6229114437989945600635724*y)*y)*y)*y)*y+(26026.448086513841468404663486538+(-468476.06555724914643128394275768+(2342380.3277862457321564197137884+(-4294363.9342747838422867694752787+2576618.3605648703053720616851672*y)*y)*y)*y+(-52052.896173027682936809326973075+(702714.09833587371964692591413652+(-2342380.3277862457321564197137884+2147181.9671373919211433847376394*y)*y)*y+(52052.896173027682936809326973076+(-468476.06555724914643128394275768+780793.44259541524405213990459608*y)*y+(-26026.448086513841468404663486538+117119.01638931228660782098568942*y+5205.2896173027682936809326973078*x)*x)*x)*x)*x;
2456 h[2] = -23423.803277862457321564197137884+(390396.72129770762202606995229807+(-2147181.9671373919211433847376394+(5153236.7211297406107441233703343+(-5582673.1145572189949728003178621+2233069.2458228875979891201271448*y)*y)*y)*y)*y+(117119.01638931228660782098568942+(-1561586.8851908304881042798091923+(6441545.9014121757634301542129182+(-10306473.442259481221488246740669+5582673.1145572189949728003178621*y)*y)*y)*y+(-234238.03277862457321564197137884+(2342380.3277862457321564197137884+(-6441545.9014121757634301542129182+5153236.7211297406107441233703343*y)*y)*y+(234238.03277862457321564197137884+(-1561586.8851908304881042798091923+2147181.9671373919211433847376394*y)*y+x*(23423.803277862457321564197137884*x-117119.01638931228660782098568942+390396.72129770762202606995229804*y))*x)*x)*x;
2457 return;
2458 #endif
2459 #if (PMAX2D>=8)
2460 case ( 36 ):
2461 h[0] = 13067.333316327398250927603811698+(-313615.99959185755802226249148074+(2548129.9966838426589308827432810+(-9513018.6542863459266752955749160+(17836909.976786898612516179202966+(-16308031.978776593017157649556998+5775761.3258167100269100008847701*x)*x)*x)*x)*x)*x;
2462 h[1] = 0;
2463 h[2] = 0;
2464 return;
2465 case ( 37 ):
2466 h[0] = 21501.620962150737345975759607768-40739.913401969818129217228730508*y+(-502458.93195762775692701248767627+882698.12370934605946637328916104*y+(3972141.5566920572675986798012246-6178886.8659654224162646130241268*y+(-14417402.687252652304617430389630+18536660.597896267248793839072382*y+(26260269.180353045269124605352540-24715547.463861689665058452096509*y+(-23303230.465926735969912254833852+12004694.482447106408742676732590*y+8003129.6549647376058284511550599*x)*x)*x)*x)*x)*x;
2467 h[1] = 1131.6642611658282813671452425141+(-40739.913401969818129217228730508+(441349.06185467302973318664458052+(-2059628.9553218074720882043413756+(4634165.1494740668121984597680955+(-4943109.4927723379330116904193018+2000782.4137411844014571127887650*x)*x)*x)*x)*x)*x;
2468 h[2] = 0;
2469 return;
2470 case ( 38 ):
2471 h[0] = 24950.370738728512929471270065574+(-141391.75874144857669233452815868+133195.13504629213746379339609150*y)*y+(-551678.31148233478474319897218755+(2885894.5926696629783821902486493-2486309.1875307865659908100603748*y)*y+(4121649.4567102622537407178679426+(-18958107.554922247565679926710357+13985489.179860674433698306589608*y)*y+(-14123617.468056829242920018259629+(53144858.883470562848053565040510-29835710.250369438791889720724498*y)*y+(24267413.632461948045139469026783+(-65887193.469565843998756466599933+21133628.094011685810921885513186*y)*y+(-20304858.364834756955591615493060+29587079.331616360135290639718460*y+6574906.5181369689189534754929912*x)*x)*x)*x)*x)*x;
2472 h[1] = -8196.6236951564392285411320671696*y+4155.2328454612504422465461173846+(-141391.75874144857669233452815868+266390.27009258427492758679218300*y+(-2486309.1875307865659908100603747*y+1442947.2963348314891910951243247+(-6319369.1849740825218933089034524+9323659.4532404496224655377264056*y+(-14917855.125184719395944860362249*y+13286214.720867640712013391260128+(-13177438.693913168799751293319987+8453451.2376046743243687542052744*y+4931179.8886027266892151066197434*x)*x)*x)*x)*x)*x;
2473 h[2] = 113.84199576606165595196016759958+x*(-8196.6236951564392285411320671696+x*(133195.13504629213746379339609150+x*(-828769.72917692885533027002012492+x*(2330914.8633101124056163844316014+x*(-2983571.0250369438791889720724498+1408908.5396007790540614590342124*x)))));
2474 return;
2475 case ( 39 ):
2476 h[0] = 24874.53830727316233136076157244+(-280444.7044534804547322731349846+(656660.87137882671316994790252156-408588.98663571439930574536156896*y)*y)*y+(-503507.35122339574437523389941037+(5200754.1013203075683059873879707+(-11031902.639164288781255124762362+6128834.7995357159895861804235344*y)*y)*y+(3440902.9660250519770105270092130+(-30889327.389660008587514349334614+(55159513.195821443906275623811811-24515339.198142863958344721694138*y)*y)*y+(-10786749.247182860141671677545421+(77958778.650094307387536214987358+(-102964424.63220002862504783111538+27784051.091228579152790684586690*y)*y)*y+(16976872.394713933291153719773190+(-86784300.761425738412540314797249+62514114.955264303093779040320051*y)*y+(-13042160.453412003625839391941282+35007904.374948009732516262579228*y+3889767.1527720010813906958421363*x)*x)*x)*x)*x)*x;
2477 h[1] = 9024.8776168987466220280019423475+(-44450.88975487442366073493493992+43777.391425255114211329860168105*y)*y+(-280444.7044534804547322731349846+(1313321.7427576534263398958050431-1225766.9599071431979172360847069*y)*y+(2600377.0506601537841529936939854+(-11031902.639164288781255124762362+9193252.1993035739843792706353016*y)*y+(-10296442.463220002862504783111538+(36773008.797214295937517082541208-24515339.198142863958344721694137*y)*y+(19489694.662523576846884053746840+(-51482212.316100014312523915557690+20838038.318421434364593013440018*y)*y+(-17356860.152285147682508062959449+25005645.982105721237511616128019*y+5834650.7291580016220860437632049*x)*x)*x)*x)*x)*x;
2478 h[2] = -1346.9966592386188988101495436339*y+673.49832961930944940507477181698+x*(-44450.88975487442366073493493992+87554.782850510228422659720336210*y+x*(-1225766.9599071431979172360847069*y+656660.87137882671316994790252156+x*(-3677300.8797214295937517082541208+6128834.7995357159895861804235344*y+x*(-12257669.599071431979172360847069*y+9193252.1993035739843792706353018+x*(-10296442.463220002862504783111538+8335215.3273685737458372053760070*y+4167607.6636842868729186026880034*x)))));
2479 return;
2480 case ( 40 ):
2481 h[0] = 21892.025945535511355450141450766+(-407802.62284590568827245496051310+(1702232.2964859995864406406497822+(-2444779.2695456168727644553313158+1122602.7258117628497387805092776*y)*y)*y)*y+(-389474.41507755037643998507464735+(6566080.4330132904639823366114074+(-24440665.059101808328598592801988+(30684474.505521517892860000586923-11974429.075325470397213658765628*y)*y)*y)*y+(2343719.5683784355005770866550838+(-33742230.501542129083577059878860+(102958707.13873596421889958385089+(-98789039.871435130777012684816434+25445661.78506662459407902487696*y)*y)*y)*y+(-6500404.3551766839299159861870553+(73771036.267630130125691290609674+(-160371817.97310897853411150132538+84818872.61688874864693008292320*y)*y)*y+(9114464.9881383602800220036586590+(-71632745.361322010411903137258671+81789627.166285579052396865675945*y)*y+(-6286575.2645458719585371708519547+25445661.785066624594079024876959*y+1696377.4523377749729386016584640*x)*x)*x)*x)*x)*x;
2482 h[1] = 15018.948032402269418273934251107+(-132879.50632057601078540667252674+(304706.45414890705921481185251822-199573.81792209117328689431276047*y)*y)*y+(-407802.62284590568827245496051310+(3404464.5929719991728812812995644+(-7334337.8086368506182933659939474+4490410.9032470513989551220371104*y)*y)*y+(3283040.2165066452319911683057037+(-24440665.059101808328598592801988+(46026711.758282276839290000880383-23948858.150650940794427317531256*y)*y)*y+(-11247410.167180709694525686626287+(68639138.092490642812599722567262+(-98789039.871435130777012684816434+33927549.04675549945877203316928*y)*y)*y+(18442759.066907532531422822652419+(-80185908.986554489267055750662690+63614154.46266656148519756219240*y)*y+(-14326549.072264402082380627451734+32715850.866514231620958746270378*y+4240943.6308444374323465041461598*x)*x)*x)*x)*x)*x;
2483 h[2] = 2291.0259710444139790587357332198+(-10691.454531540598568940766755025+10691.454531540598568940766755025*y)*y+(-132879.50632057601078540667252674+(609412.90829781411842962370503644-598721.45376627351986068293828142*y)*y+(1702232.2964859995864406406497822+(-7334337.8086368506182933659939474+6735616.3548705770984326830556656*y)*y+(-8146888.3530339361095328642673292+(30684474.505521517892860000586922-23948858.150650940794427317531256*y)*y+(17159784.523122660703149930641816+(-49394519.935717565388506342408217+25445661.78506662459407902487696*y)*y+x*(5452641.8110857052701597910450630*x-16037181.797310897853411150132538+25445661.78506662459407902487696*y))*x)*x)*x)*x;
2484 return;
2485 case ( 41 ):
2486 h[0] = 16744.784262569643509799299505354+(-460974.06087544665662271012755916+(2943142.0809740055768988415836468+(-7107668.4258060322098065732488608+(7180557.4867136883050845231408252-2553087.1063870891751411637834046*y)*y)*y)*y)*y+(-248216.80200985589202761314560877+(6063581.8776693367909602639855858+(-33615646.900763340806025323148160+(68318719.791284145519981141981846+(-55316887.305053598794725215307098+14467493.602860171992466594772625*y)*y)*y)*y)*y+(1258813.7816214120238543238098730+(-25637249.693303687133709186325021+(112229454.05159912832391365797882+(-162641104.55502938449047413731318+72337468.014300859962332973863124*y)*y)*y)*y+(-2990421.4718330257468088631351915+(46983894.666151293848083916847374+(-141838172.57706050973006465463358+107166619.28044571846271551683426*y)*y)*y+(3649377.9819306193649297885098435+(-39005497.458691640175767780024235+60281223.345250716635277478219272*y)*y+(-2222131.3703739479857710129225928+12056244.669050143327055495643854*y+535833.09640222859231357758417133*x)*x)*x)*x)*x)*x;
2487 h[1] = 20684.733500821324335634428800731+(-283676.34515412101946012930926716+(1111065.6851869739928855064612964+(-1631138.9846361958618957435282863+797839.72074596536723161368231390*y)*y)*y)*y+(-460974.06087544665662271012755916+(5886284.1619480111537976831672936+(-21323005.277418096629419719746583+(28722229.946854753220338092563301-12765435.531935445875705818917023*y)*y)*y)*y+(3031790.9388346683954801319927929+(-33615646.900763340806025323148160+(102478079.68692621827997171297277+(-110633774.61010719758945043061420+36168734.007150429981166486931564*y)*y)*y)*y+(-8545749.8977678957112363954416736+(74819636.034399418882609105319216+(-162641104.55502938449047413731317+96449957.352401146616443965150832*y)*y)*y+(11745973.666537823462020979211844+(-70919086.288530254865032327316792+80374964.460334288847036637625695*y)*y+(-7801099.4917383280351535560048470+24112489.338100286654110991287709*y+2009374.1115083572211759159406424*x)*x)*x)*x)*x)*x;
2488 h[2] = 5909.923857377521238752693943066+(-47279.390859020169910021551544527+(106378.62943279538229754849097519-70919.086288530254865032327316791*y)*y)*y+(-283676.34515412101946012930926716+(2222131.3703739479857710129225928+(-4893416.9539085875856872305848586+3191358.8829838614689264547292556*y)*y)*y+(2943142.0809740055768988415836468+(-21323005.277418096629419719746583+(43083344.920282129830507138844950-25530871.063870891751411637834046*y)*y)*y+(-11205215.633587780268675107716053+(68318719.791284145519981141981846+(-110633774.61010719758945043061420+48224978.676200573308221982575417*y)*y)*y+(18704909.008599854720652276329804+(-81320552.277514692245237068656588+72337468.014300859962332973863124*y)*y+x*(4018748.2230167144423518318812848*x-14183817.257706050973006465463358+32149985.784133715538814655050278*y))*x)*x)*x)*x;
2489 return;
2490 case ( 42 ):
2491 h[0] = 10493.782158973951180198085216369+(-393195.59273216682381313641749496+(3469372.8770485307983512036837792+(-12078557.423798588705370857269453+(19370665.230187630290794220567767+(-14247557.948412633145228943128053+3844579.1289367422772840005266174*y)*y)*y)*y)*y)*y+(-120785.57423798588705370857269453+(3893407.3397989067848163508006854+(-28525954.766843475453109896955517+(77868146.795978135696327016013710+(-87890779.551896113558230493322406+34601212.160430680495556004739556*y)*y)*y)*y)*y+(494706.87320877198420933830305740+(-12990873.995170609767159507127039+(72856830.418019146765375277359361+(-135691028.08012031566884707741003+78639118.546433364762627283498992*y)*y)*y)*y+(-980847.39363470808990422918961160+(19544133.874040056830711780751956+(-74013288.043701990364825678587288+69901438.707940768677890918665772*y)*y)*y+(1024749.9516467419672907721991903+(-13723297.158103077380144761238060+26213039.515477788254209094499663*y)*y+(-544820.03698836187351885568960087+3669825.5321668903555892732299528*y+116502.39784656794779648486444295*x)*x)*x)*x)*x)*x;
2492 h[1] = 23771.628972369562877591580796264+(-462583.05027313743978016049117054+(2698401.1265933017320509361984948+(-6553259.878869447063552273624916+(6996568.6353812037766749274289545-2713820.5616024063133769415482007*y)*y)*y)*y)*y+(-393195.59273216682381313641749496+(6938745.7540970615967024073675584+(-36235672.271395766116112571808359+(77482660.920750521163176882271067+(-71237789.742063165726144715640263+23067474.773620453663704003159704*y)*y)*y)*y)*y+(1946703.6698994533924081754003427+(-28525954.766843475453109896955516+(116802220.19396720354449052402056+(-175781559.10379222711646098664481+86503030.401076701238890011848890*y)*y)*y)*y+(-4330291.3317235365890531690423464+(48571220.278679431176916851572908+(-135691028.08012031566884707741003+104852158.06191115301683637799866*y)*y)*y+(4886033.4685100142076779451879890+(-37006644.021850995182412839293644+52426079.030955576508418188999329*y)*y+(-2744659.4316206154760289522476118+10485215.806191115301683637799865*y+611637.58869448172593154553832547*x)*x)*x)*x)*x)*x;
2493 h[2] = 12849.529174253817771671124754737+(-154194.35009104581326005349705685+(578228.81284142179972520061396318+(-848068.92550075197293029423381267+424034.46275037598646514711690633*y)*y)*y)*y+(-462583.05027313743978016049117054+(5396802.2531866034641018723969896+(-19659779.636608341190656820874748+(27986274.541524815106699709715818-13569102.808012031566884707741003*y)*y)*y)*y+(3469372.8770485307983512036837792+(-36235672.271395766116112571808359+(116223991.38112578174476532340660+(-142475579.48412633145228943128053+57668686.934051134159260007899260*y)*y)*y)*y+(-9508651.5889478251510366323185056+(77868146.795978135696327016013710+(-175781559.10379222711646098664481+115337373.86810226831852001579852*y)*y)*y+(12142805.069669857794229212893227+(-67845514.040060157834423538705014+78639118.546433364762627283498992*y)*y+x*(1747535.9676985192169472729666442*x-7401328.8043701990364825678587288+20970431.612382230603367275599731*y))*x)*x)*x)*x;
2494 return;
2495 case ( 43 ):
2496 h[0] = 4600.8694830433953530385461755268+(-215320.69180643090252220396101466+(2360246.0448012618161087741880452+(-10213930.252356337683745572509669+(20496873.546958326297786723211972+(-19130415.310494437877934274997840+6710828.2279670964619420234516232*y)*y)*y)*y)*y)*y+(-38647.303657564520965523787874425+(1518286.9294043204665027202379238+(-13416135.412554540849460400647836+(43892294.868233991667987730514525+(-59213190.246768498193606089279029+27875748.023863323764989943568281*y)*y)*y)*y)*y+(124223.47604217167453204074673922+(-3919940.7995529728407888413415488+(26086929.968856051651728556815237+(-57142798.979398970284738743500042+38716316.699810171895819366067057*y)*y)*y)*y+(-202438.25725390939553369603172318+(4803307.7402973047485722422072500+(-21366437.879253528019511008439146+23464434.363521316300496585495186*y)*y)*y+(179433.90983869241876850330084555+(-2843337.3405208183281778215364755+6335397.2781507554011340780837004*y)*y+(-82815.650694781116354693831159481+657004.16217859685641390439386522*y+15642.956242347544200331056996791*x)*x)*x)*x)*x)*x;
2497 h[1] = 21163.999621999618623977312407423+(-546583.29458555536794097928565258+(4347821.6614760086086214261358727+(-15182869.294043204665027202379239+(25962706.492813879977196516068498+(-21316748.488836659349698192140451+6710828.2279670964619420234516233*y)*y)*y)*y)*y)*y+(-215320.69180643090252220396101466+(4720492.0896025236322175483760904+(-30641790.757069013051236717529007+(81987494.187833305191146892847888+(-95652076.552472189389671374989200+40264969.367802578771652140709741*y)*y)*y)*y)*y+(759143.46470216023325136011896191+(-13416135.412554540849460400647836+(65838442.302350987501981595771787+(-118426380.49353699638721217855806+69689370.059658309412474858920705*y)*y)*y)*y+(-1306646.9331843242802629471138496+(17391286.645904034434485704543491+(-57142798.979398970284738743500043+51621755.599746895861092488089409*y)*y)*y+(1200826.9350743261871430605518125+(-10683218.939626764009755504219573+17598325.772640987225372439121390*y)*y+(-568667.46810416366563556430729510+2534158.9112603021604536312334803*y+109500.69369643280940231739897754*x)*x)*x)*x)*x)*x;
2498 h[2] = 24844.695208434334906408149347844+(-414078.25347390558177346915579740+(2277430.3941064806997540803568858+(-5465832.9458555536794097928565259+(5921319.0246768498193606089279030-2368527.6098707399277442435711611*y)*y)*y)*y)*y+(-546583.29458555536794097928565258+(8695643.3229520172172428522717454+(-45548607.882129613995081607137714+(103850825.97125551990878606427399+(-106583742.44418329674849096070225+40264969.367802578771652140709740*y)*y)*y)*y)*y+(2360246.0448012618161087741880452+(-30641790.757069013051236717529007+(122981241.28174995778672033927183+(-191304153.10494437877934274997840+100662423.41950644692913035177435*y)*y)*y)*y+(-4472045.137518180283153466882612+(43892294.868233991667987730514525+(-118426380.49353699638721217855806+92919160.079544412549966478560940*y)*y)*y+(4347821.6614760086086214261358728+(-28571399.489699485142369371750021+38716316.699810171895819366067056*y)*y+x*(422359.81854338369340893853891336*x-2136643.7879253528019511008439146+7039330.3090563948901489756485558*y))*x)*x)*x)*x;
2499 return;
2500 case ( 44 ):
2501 h[0] = 979.59991833401047910685768342766+(-52898.395590036565871770314905094+(661229.94487545707339712893631366+(-3232679.7305022345810526303553112+(7273529.3936300278073684182994504+(-7564470.5693752289196631550314284+2941738.5547570334687578936233332*y)*y)*y)*y)*y)*y+(-5877.5995100040628746411461005659+(264491.97795018282935885157452547+(-2644919.7795018282935885157452547+(9698039.1915067037431578910659336+(-14547058.787260055614736836598901+7564470.5693752289196631550314284*y)*y)*y)*y)*y+(14693.998775010157186602865251415+(-528983.95590036565871770314905094+(3967379.6692527424403827736178820+(-9698039.1915067037431578910659338+7273529.3936300278073684182994508*y)*y)*y)*y+(-19591.998366680209582137153668553+(528983.95590036565871770314905094+(-2644919.7795018282935885157452546+3232679.7305022345810526303553112*y)*y)*y+(14693.998775010157186602865251415+(-264491.97795018282935885157452547+661229.94487545707339712893631366*y)*y+(-5877.5995100040628746411461005658+52898.395590036565871770314905091*y+979.59991833401047910685768342765*x)*x)*x)*x)*x)*x;
2502 h[1] = 8816.3992650060943119617191508489+(-264491.97795018282935885157452546+(2424509.7978766759357894727664834+(-9698039.1915067037431578910659340+(18911176.423438072299157887578571+(-17650431.328542200812547361739999+6303725.4744793574330526291928567*y)*y)*y)*y)*y)*y+(-52898.395590036565871770314905094+(1322459.8897509141467942578726273+(-9698039.1915067037431578910659336+(29094117.574520111229473673197802+(-37822352.846876144598315775157142+17650431.328542200812547361739999*y)*y)*y)*y)*y+(132245.98897509141467942578726274+(-2644919.7795018282935885157452548+(14547058.787260055614736836598900+(-29094117.574520111229473673197802+18911176.423438072299157887578571*y)*y)*y)*y+(-176327.98530012188623923438301698+(2644919.7795018282935885157452546+(-9698039.1915067037431578910659337+9698039.1915067037431578910659344*y)*y)*y+(132245.98897509141467942578726274+(-1322459.8897509141467942578726273+2424509.7978766759357894727664834*y)*y+(-52898.395590036565871770314905094+264491.97795018282935885157452545*y+8816.3992650060943119617191508486*x)*x)*x)*x)*x)*x;
2503 h[2] = 44081.996325030471559808595754244+(-969803.91915067037431578910659336+(7273529.3936300278073684182994502+(-25214901.897917429732210516771428+(44126078.321355502031368404349998+(-37822352.846876144598315775157142+12607450.948958714866105258385714*y)*y)*y)*y)*y)*y+(-264491.97795018282935885157452546+(4849019.5957533518715789455329668+(-29094117.574520111229473673197802+(75644705.693752289196631550314284+(-88252156.642711004062736808700000+37822352.846876144598315775157140*y)*y)*y)*y)*y+(661229.94487545707339712893631366+(-9698039.1915067037431578910659336+(43641176.361780166844210509796702+(-75644705.693752289196631550314284+44126078.321355502031368404349998*y)*y)*y)*y+(-881639.92650060943119617191508490+(9698039.1915067037431578910659336+(-29094117.574520111229473673197802+25214901.897917429732210516771428*y)*y)*y+(661229.94487545707339712893631366+(-4849019.5957533518715789455329668+7273529.3936300278073684182994508*y)*y+x*(44081.996325030471559808595754244*x-264491.97795018282935885157452546+969803.91915067037431578910659336*y))*x)*x)*x)*x;
2504 return;
2505 #endif
2506 }
2507 printf("Error in hess_triangle_2d -- unmatched switch/case %d\n",i);
2508 }
2509
2510 /* \phi_i(x,y) for quadrilateral */
2511 static RangeField
2512 eval_quadrilateral_2d ( const int i, DomainType xi )
2513 {
2514 const RangeField &x = xi[0];
2515 const RangeField &y = xi[1];
2516
2517 switch (i) {
2518 #if (PMAX2D>=0)
2519 case ( 0 ):
2520 return 1.;
2521 #endif
2522 #if (PMAX2D>=1)
2523 case ( 1 ):
2524 return -1.7320508075688772935274463415059
2525 + 3.4641016151377545870548926830117*x;
2526 case ( 2 ):
2527 return -1.7320508075688772935274463415059
2528 + 3.4641016151377545870548926830117*y;
2529 #endif
2530 #if (PMAX2D>=2)
2531 case ( 3 ):
2532 return 2.2360679774997896964091736687313
2533 + x*(-13.416407864998738178455042012388
2534 + 13.416407864998738178455042012388*x);
2535 case ( 4 ):
2536 return 3.
2537 - 6.*y
2538 + x*(-6.
2539 + 12.*y);
2540 case ( 5 ):
2541 return 2.2360679774997896964091736687313
2542 + y*(-13.416407864998738178455042012388
2543 + 13.416407864998738178455042012388*y);
2544 #endif
2545 #if (PMAX2D>=3)
2546 case ( 6 ):
2547 return -2.6457513110645905905016157536393
2548 + x*(31.749015732775087086019389043671
2549 + x*(-79.372539331937717715048472609178
2550 + 52.91502622129181181003231507279*x));
2551 case ( 7 ):
2552 return -3.8729833462074168851792653997824
2553 + 7.745966692414833770358530799565*y
2554 + x*(23.237900077244501311075592398694
2555 - 46.475800154489002622151184797389*y
2556 + x*(-23.237900077244501311075592398694
2557 + 46.475800154489002622151184797389*y));
2558 case ( 8 ):
2559 return -3.8729833462074168851792653997824
2560 + (23.237900077244501311075592398694
2561 - 23.237900077244501311075592398694*y)*y
2562 + x*(7.7459666924148337703585307995648
2563 + y*(-46.475800154489002622151184797389
2564 + 46.475800154489002622151184797389*y));
2565 case ( 9 ):
2566 return -2.6457513110645905905016157536393
2567 + y*(31.749015732775087086019389043671
2568 + y*(-79.372539331937717715048472609178
2569 + 52.91502622129181181003231507279*y));
2570 #endif
2571 #if (PMAX2D>=4)
2572 case ( 10 ):
2573 return 3.
2574 + x*(-60.
2575 + x*(270.
2576 + x*(-420.
2577 + 210.*x)));
2578 case ( 11 ):
2579 return 4.582575694955840006588047193728
2580 - 9.165151389911680013176094387456*y
2581 + x*(-54.990908339470080079056566324736
2582 + 109.98181667894016015811313264947*y
2583 + x*(137.47727084867520019764141581184
2584 - 274.95454169735040039528283162368*y
2585 + x*(-91.65151389911680013176094387456
2586 + 183.30302779823360026352188774912*y)));
2587 case ( 12 ):
2588 return 5.
2589 + y*(-30.
2590 + 30.*y)
2591 + x*(-30.
2592 + (180.
2593 - 180.*y)*y
2594 + x*(30.
2595 + y*(-180.
2596 + 180.*y)));
2597 case ( 13 ):
2598 return 4.582575694955840006588047193728
2599 + y*(-54.990908339470080079056566324736
2600 + (137.47727084867520019764141581184
2601 - 91.65151389911680013176094387456*y)*y)
2602 + x*(-9.165151389911680013176094387456
2603 + y*(109.98181667894016015811313264947
2604 + y*(-274.95454169735040039528283162368
2605 + 183.30302779823360026352188774912*y)));
2606 case ( 14 ):
2607 return 3.
2608 + y*(-60.
2609 + y*(270.
2610 + y*(-420.
2611 + 210.*y)));
2612 #endif
2613 #if (PMAX2D>=5)
2614 case ( 15 ):
2615 return -3.3166247903553998491149327366707
2616 + x*(99.498743710661995473447982100121
2617 + x*(-696.49120597463396831413587470084
2618 + x*(1857.3098825990239155043623325356
2619 + x*(-2089.4736179239019049424076241025
2620 + 835.789447169560761976963049641*x))));
2621 case ( 16 ):
2622 return -5.1961524227066318805823390245176
2623 + 10.392304845413263761164678049035*y
2624 + x*(103.92304845413263761164678049035
2625 - 207.8460969082652752232935609807*y
2626 + x*(-467.65371804359686925241051220659
2627 + 935.3074360871937385048210244132*y
2628 + x*(727.46133917892846328152746343247
2629 - 1454.9226783578569265630549268649*y
2630 + x*(-363.73066958946423164076373171623
2631 + 727.4613391789284632815274634325*y))));
2632 case ( 17 ):
2633 return -5.9160797830996160425673282915616
2634 + (35.49647869859769625540396974937
2635 - 35.49647869859769625540396974937*y)*y
2636 + x*(70.992957397195392510807939498739
2637 + y*(-425.95774438317235506484763699244
2638 + 425.95774438317235506484763699244*y)
2639 + x*(-177.48239349298848127701984874685
2640 + (1064.8943609579308876621190924811
2641 - 1064.8943609579308876621190924811*y)*y
2642 + x*(118.32159566199232085134656583123
2643 + y*(-709.92957397195392510807939498739
2644 + 709.9295739719539251080793949874*y))));
2645 case ( 18 ):
2646 return -5.9160797830996160425673282915616
2647 + y*(70.992957397195392510807939498739
2648 + y*(-177.48239349298848127701984874685
2649 + 118.32159566199232085134656583123*y))
2650 + x*(35.49647869859769625540396974937
2651 + y*(-425.95774438317235506484763699244
2652 + (1064.8943609579308876621190924811
2653 - 709.9295739719539251080793949874*y)*y)
2654 + x*(-35.49647869859769625540396974937
2655 + y*(425.95774438317235506484763699244
2656 + y*(-1064.8943609579308876621190924811
2657 + 709.9295739719539251080793949874*y))));
2658 case ( 19 ):
2659 return -5.1961524227066318805823390245176
2660 + y*(103.92304845413263761164678049035
2661 + y*(-467.65371804359686925241051220659
2662 + (727.46133917892846328152746343247
2663 - 363.73066958946423164076373171623*y)*y))
2664 + x*(10.392304845413263761164678049035
2665 + y*(-207.8460969082652752232935609807
2666 + y*(935.30743608719373850482102441317
2667 + y*(-1454.9226783578569265630549268649
2668 + 727.4613391789284632815274634325*y))));
2669 case ( 20 ):
2670 return -3.3166247903553998491149327366707
2671 + y*(99.498743710661995473447982100121
2672 + y*(-696.49120597463396831413587470084
2673 + y*(1857.3098825990239155043623325356
2674 + y*(-2089.4736179239019049424076241025
2675 + 835.789447169560761976963049641*y))));
2676 #endif
2677 #if (PMAX2D>=6)
2678 case ( 21 ):
2679 return 3.6055512754639892931192212674705
2680 + x*(-151.43315356948755031100729323376
2681 + x*(1514.3315356948755031100729323376
2682 + x*(-6057.3261427795020124402917293504
2683 + x*(11357.486517711566273325546992532
2684 + x*(-9994.5881355861783205264813534282
2685 + 3331.5293785287261068421604511427*x)))));
2686 case ( 22 ):
2687 return 5.7445626465380286598506114682189
2688 - 11.489125293076057319701222936438*y
2689 + x*(-172.33687939614085979551834404657
2690 + 344.67375879228171959103668809314*y
2691 + x*(1206.358155772986018568628408326
2692 - 2412.716311545972037137256816652*y
2693 + x*(-3216.9550820612960495163424222026
2694 + 6433.910164122592099032684844405*y
2695 + x*(3619.0744673189580557058852249779
2696 - 7238.148934637916111411770449956*y
2697 + x*(-1447.6297869275832222823540899912
2698 + 2895.2595738551664445647081799823*y)))));
2699 case ( 23 ):
2700 return 6.7082039324993690892275210061938
2701 + y*(-40.249223594996214535365126037163
2702 + 40.249223594996214535365126037163*y)
2703 + x*(-134.16407864998738178455042012388
2704 + (804.98447189992429070730252074326
2705 - 804.9844718999242907073025207433*y)*y
2706 + x*(603.73835392494321803047689055744
2707 + y*(-3622.4301235496593081828613433447
2708 + 3622.4301235496593081828613433447*y)
2709 + x*(-939.14855054991167249185294086714
2710 + (5634.8913032994700349511176452028
2711 - 5634.891303299470034951117645203*y)*y
2712 + x*(469.57427527495583624592647043357
2713 + y*(-2817.4456516497350174755588226014
2714 + 2817.4456516497350174755588226014*y)))));
2715 case ( 24 ):
2716 return 7.
2717 + y*(-84.
2718 + (210.
2719 - 140.*y)*y)
2720 + x*(-84.
2721 + y*(1008.
2722 + y*(-2520.
2723 + 1680.*y))
2724 + x*(210.
2725 + y*(-2520.
2726 + (6300.
2727 - 4200.*y)*y)
2728 + x*(-140.
2729 + y*(1680.
2730 + y*(-4200.
2731 + 2800.*y)))));
2732 case ( 25 ):
2733 return 6.7082039324993690892275210061938
2734 + y*(-134.16407864998738178455042012388
2735 + y*(603.73835392494321803047689055744
2736 + y*(-939.14855054991167249185294086714
2737 + 469.57427527495583624592647043357*y)))
2738 + x*(-40.249223594996214535365126037163
2739 + y*(804.98447189992429070730252074326
2740 + y*(-3622.4301235496593081828613433447
2741 + (5634.8913032994700349511176452028
2742 - 2817.4456516497350174755588226014*y)*y))
2743 + x*(40.249223594996214535365126037163
2744 + y*(-804.98447189992429070730252074326
2745 + y*(3622.4301235496593081828613433447
2746 + y*(-5634.8913032994700349511176452028
2747 + 2817.4456516497350174755588226014*y)))));
2748 case ( 26 ):
2749 return 5.7445626465380286598506114682189
2750 + y*(-172.33687939614085979551834404657
2751 + y*(1206.358155772986018568628408326
2752 + y*(-3216.9550820612960495163424222026
2753 + (3619.0744673189580557058852249779
2754 - 1447.6297869275832222823540899912*y)*y)))
2755 + x*(-11.489125293076057319701222936438
2756 + y*(344.67375879228171959103668809314
2757 + y*(-2412.716311545972037137256816652
2758 + y*(6433.9101641225920990326848444052
2759 + y*(-7238.1489346379161114117704499559
2760 + 2895.2595738551664445647081799823*y)))));
2761 case ( 27 ):
2762 return 3.6055512754639892931192212674705
2763 + y*(-151.43315356948755031100729323376
2764 + y*(1514.3315356948755031100729323376
2765 + y*(-6057.3261427795020124402917293504
2766 + y*(11357.486517711566273325546992532
2767 + y*(-9994.5881355861783205264813534282
2768 + 3331.5293785287261068421604511427*y)))));
2769 #endif
2770 #if (PMAX2D>=7)
2771 case ( 28 ):
2772 return -3.8729833462074168851792653997824
2773 + x*(216.88706738761534557003886238781
2774 + x*(-2927.9754097328071651955246422355
2775 + x*(16266.530054071150917752914679086
2776 + x*(-44732.957648695665023820515367487
2777 + x*(64415.459014121757634301542129181
2778 + x*(-46522.275954643491624773335982186
2779 + 13292.078844183854749935238852053*x))))));
2780 case ( 29 ):
2781 return -6.2449979983983982058468931209398
2782 + 12.48999599679679641169378624188*y
2783 + x*(262.28991593273272464556951107947
2784 - 524.5798318654654492911390221589*y
2785 + x*(-2622.8991593273272464556951107947
2786 + 5245.798318654654492911390221589*y
2787 + x*(10491.596637309308985822780443179
2788 - 20983.193274618617971645560886358*y
2789 + x*(-19671.74369495495434841771333096
2790 + 39343.487389909908696835426661921*y
2791 + x*(17311.134451560359826607587731245
2792 - 34622.26890312071965321517546249*y
2793 + x*(-5770.3781505201199422025292437484
2794 + 11540.756301040239884405058487497*y))))));
2795 case ( 30 ):
2796 return -7.4161984870956629487113974408007
2797 + (44.497190922573977692268384644804
2798 - 44.497190922573977692268384644804*y)*y
2799 + x*(222.48595461286988846134192322402
2800 + y*(-1334.9157276772193307680515393441
2801 + 1334.9157276772193307680515393441*y)
2802 + x*(-1557.4016822900892192293934625681
2803 + (9344.4100937405353153763607754089
2804 - 9344.410093740535315376360775409*y)*y
2805 + x*(4153.0711527735712512783825668484
2806 + y*(-24918.42691664142750767029540109
2807 + 24918.42691664142750767029540109*y)
2808 + x*(-4672.2050468702676576881803877044
2809 + (28033.230281221605946129082326227
2810 - 28033.230281221605946129082326227*y)*y
2811 + x*(1868.8820187481070630752721550818
2812 + y*(-11213.292112488642378451632930491
2813 + 11213.292112488642378451632930491*y))))));
2814 case ( 31 ):
2815 return -7.9372539331937717715048472609178
2816 + y*(95.247047198325261258058167131013
2817 + y*(-238.11761799581315314514541782753
2818 + 158.74507866387543543009694521836*y))
2819 + x*(158.74507866387543543009694521836
2820 + y*(-1904.9409439665052251611633426203
2821 + (4762.3523599162630629029083565507
2822 - 3174.9015732775087086019389043671*y)*y)
2823 + x*(-714.3528539874394594354362534826
2824 + y*(8572.2342478492735132252350417912
2825 + y*(-21430.585619623183783063087604478
2826 + 14287.057079748789188708725069652*y))
2827 + x*(1111.2155506471280480106786165285
2828 + y*(-13334.586607765536576128143398342
2829 + (33336.466519413841440320358495855
2830 - 22224.31101294256096021357233057*y)*y)
2831 + x*(-555.60777532356402400533930826424
2832 + y*(6667.2933038827682880640716991709
2833 + y*(-16668.233259706920720160179247927
2834 + 11112.155506471280480106786165285*y))))));
2835 case ( 32 ):
2836 return -7.9372539331937717715048472609178
2837 + y*(158.74507866387543543009694521836
2838 + y*(-714.3528539874394594354362534826
2839 + (1111.2155506471280480106786165285
2840 - 555.6077753235640240053393082642*y)*y))
2841 + x*(95.247047198325261258058167131013
2842 + y*(-1904.9409439665052251611633426203
2843 + y*(8572.2342478492735132252350417912
2844 + y*(-13334.586607765536576128143398342
2845 + 6667.293303882768288064071699171*y)))
2846 + x*(-238.11761799581315314514541782753
2847 + y*(4762.3523599162630629029083565507
2848 + y*(-21430.585619623183783063087604478
2849 + (33336.466519413841440320358495855
2850 - 16668.233259706920720160179247927*y)*y))
2851 + x*(158.74507866387543543009694521836
2852 + y*(-3174.9015732775087086019389043671
2853 + y*(14287.057079748789188708725069652
2854 + y*(-22224.31101294256096021357233057
2855 + 11112.155506471280480106786165285*y))))));
2856 case ( 33 ):
2857 return -7.4161984870956629487113974408007
2858 + y*(222.48595461286988846134192322402
2859 + y*(-1557.4016822900892192293934625681
2860 + y*(4153.0711527735712512783825668484
2861 + y*(-4672.2050468702676576881803877044
2862 + 1868.8820187481070630752721550818*y))))
2863 + x*(44.497190922573977692268384644804
2864 + y*(-1334.9157276772193307680515393441
2865 + y*(9344.4100937405353153763607754089
2866 + y*(-24918.42691664142750767029540109
2867 + (28033.230281221605946129082326227
2868 - 11213.292112488642378451632930491*y)*y)))
2869 + x*(-44.497190922573977692268384644804
2870 + y*(1334.9157276772193307680515393441
2871 + y*(-9344.4100937405353153763607754089
2872 + y*(24918.42691664142750767029540109
2873 + y*(-28033.230281221605946129082326227
2874 + 11213.292112488642378451632930491*y))))));
2875 case ( 34 ):
2876 return -6.2449979983983982058468931209398
2877 + y*(262.28991593273272464556951107947
2878 + y*(-2622.8991593273272464556951107947
2879 + y*(10491.596637309308985822780443179
2880 + y*(-19671.74369495495434841771333096
2881 + (17311.134451560359826607587731245
2882 - 5770.378150520119942202529243748*y)*y))))
2883 + x*(12.48999599679679641169378624188
2884 + y*(-524.57983186546544929113902215894
2885 + y*(5245.7983186546544929113902215894
2886 + y*(-20983.193274618617971645560886358
2887 + y*(39343.487389909908696835426661921
2888 + y*(-34622.26890312071965321517546249
2889 + 11540.756301040239884405058487497*y))))));
2890 case ( 35 ):
2891 return -3.8729833462074168851792653997824
2892 + y*(216.88706738761534557003886238781
2893 + y*(-2927.9754097328071651955246422355
2894 + y*(16266.530054071150917752914679086
2895 + y*(-44732.957648695665023820515367487
2896 + y*(64415.459014121757634301542129181
2897 + y*(-46522.275954643491624773335982186
2898 + 13292.078844183854749935238852053*y))))));
2899 #endif
2900 #if (PMAX2D>=8)
2901 case ( 36 ):
2902 return 4.1231056256176605498214098559741
2903 + x*(-296.86360504447155958714150963013
2904 + x*(5195.1130882782522927749764185273
2905 + x*(-38097.4959807071834803498270692
2906 + x*(142865.6099276519380513118515095
2907 + x*(-297160.46864951603114672865113976
2908 + x*(346687.21342443536967118342632972
2909 + x*(-212257.47760679716510480617938555
2910 + 53064.36940169929127620154484639*x)))))));
2911 case ( 37 ):
2912 return 6.7082039324993690892275210061938
2913 - 13.416407864998738178455042012388*y
2914 + x*(-375.65942021996466899674117634685
2915 + 751.3188404399293379934823526937*y
2916 + x*(5071.4021729695230314560058806825
2917 - 10142.804345939046062912011761365*y
2918 + x*(-28174.456516497350174755588226014
2919 + 56348.91303299470034951117645203*y
2920 + x*(77479.755420367712980577867621539
2921 - 154959.51084073542596115573524308*y
2922 + x*(-111570.84780532950669203212937502
2923 + 223141.69561065901338406425875003*y
2924 + x*(80578.9456371824214998009823264
2925 - 161157.8912743648429996019646528*y
2926 + x*(-23022.555896337834714228852093257
2927 + 46045.111792675669428457704186514*y)))))));
2928 case ( 38 ):
2929 return 8.0622577482985496523666132303038
2930 + y*(-48.373546489791297914199679381823
2931 + 48.373546489791297914199679381823*y)
2932 + x*(-338.61482542853908539939775567276
2933 + (2031.6889525712345123963865340366
2934 - 2031.6889525712345123963865340366*y)*y
2935 + x*(3386.1482542853908539939775567276
2936 + y*(-20316.889525712345123963865340366
2937 + 20316.889525712345123963865340366*y)
2938 + x*(-13544.59301714156341597591022691
2939 + (81267.558102849380495855461361462
2940 - 81267.55810284938049585546136146*y)*y
2941 + x*(25396.111907140431404954831675457
2942 + y*(-152376.67144284258842972899005274
2943 + 152376.67144284258842972899005274*y)
2944 + x*(-22348.578478283579636360251874402
2945 + (134091.47086970147781816151124641
2946 - 134091.47086970147781816151124641*y)*y
2947 + x*(7449.5261594278598787867506248007
2948 + y*(-44697.156956567159272720503748804
2949 + 44697.156956567159272720503748804*y)))))));
2950 case ( 39 ):
2951 return 8.7749643873921220604063883074163
2952 + y*(-105.299572648705464724876659689
2953 + (263.24893162176366181219164922249
2954 - 175.49928774784244120812776614833*y)*y)
2955 + x*(-263.24893162176366181219164922249
2956 + y*(3158.9871794611639417462997906699
2957 + y*(-7897.4679486529098543657494766747
2958 + 5264.97863243527323624383298445*y))
2959 + x*(1842.7425213523456326853415445574
2960 + y*(-22112.910256228147592224098534689
2961 + (55282.275640570368980560246336723
2962 - 36854.850427046912653706830891149*y)*y)
2963 + x*(-4913.9800569395883538275774521531
2964 + y*(58967.760683275060245930929425838
2965 + y*(-147419.40170818765061482732356459
2966 + 98279.60113879176707655154904306*y))
2967 + x*(5528.2275640570368980560246336723
2968 + y*(-66338.730768684442776672295604067
2969 + (165846.82692171110694168073901017
2970 - 110564.55128114073796112049267345*y)*y)
2971 + x*(-2211.2910256228147592224098534689
2972 + y*(26535.492307473777110668918241627
2973 + y*(-66338.730768684442776672295604067
2974 + 44225.820512456295184448197069378*y)))))));
2975 case ( 40 ):
2976 return 9.
2977 + y*(-180.
2978 + y*(810.
2979 + y*(-1260.
2980 + 630.*y)))
2981 + x*(-180.
2982 + y*(3600.
2983 + y*(-16200.
2984 + (25200.
2985 - 12600.*y)*y))
2986 + x*(810.
2987 + y*(-16200.
2988 + y*(72900.
2989 + y*(-113400.
2990 + 56700.*y)))
2991 + x*(-1260.
2992 + y*(25200.
2993 + y*(-113400.
2994 + (176400.
2995 - 88200.*y)*y))
2996 + x*(630.
2997 + y*(-12600.
2998 + y*(56700.
2999 + y*(-88200.
3000 + 44100.*y)))))));
3001 case ( 41 ):
3002 return 8.7749643873921220604063883074163
3003 + y*(-263.24893162176366181219164922249
3004 + y*(1842.7425213523456326853415445574
3005 + y*(-4913.9800569395883538275774521531
3006 + (5528.2275640570368980560246336723
3007 - 2211.2910256228147592224098534689*y)*y)))
3008 + x*(-105.299572648705464724876659689
3009 + y*(3158.9871794611639417462997906699
3010 + y*(-22112.910256228147592224098534689
3011 + y*(58967.760683275060245930929425838
3012 + y*(-66338.730768684442776672295604067
3013 + 26535.492307473777110668918241627*y))))
3014 + x*(263.24893162176366181219164922249
3015 + y*(-7897.4679486529098543657494766747
3016 + y*(55282.275640570368980560246336723
3017 + y*(-147419.40170818765061482732356459
3018 + (165846.82692171110694168073901017
3019 - 66338.73076868444277667229560407*y)*y)))
3020 + x*(-175.49928774784244120812776614833
3021 + y*(5264.9786324352732362438329844498
3022 + y*(-36854.850427046912653706830891149
3023 + y*(98279.601138791767076551549043063
3024 + y*(-110564.55128114073796112049267345
3025 + 44225.820512456295184448197069378*y)))))));
3026 case ( 42 ):
3027 return 8.0622577482985496523666132303038
3028 + y*(-338.61482542853908539939775567276
3029 + y*(3386.1482542853908539939775567276
3030 + y*(-13544.59301714156341597591022691
3031 + y*(25396.111907140431404954831675457
3032 + y*(-22348.578478283579636360251874402
3033 + 7449.526159427859878786750624801*y)))))
3034 + x*(-48.373546489791297914199679381823
3035 + y*(2031.6889525712345123963865340366
3036 + y*(-20316.889525712345123963865340366
3037 + y*(81267.558102849380495855461361462
3038 + y*(-152376.67144284258842972899005274
3039 + (134091.47086970147781816151124641
3040 - 44697.156956567159272720503748804*y)*y))))
3041 + x*(48.373546489791297914199679381823
3042 + y*(-2031.6889525712345123963865340366
3043 + y*(20316.889525712345123963865340366
3044 + y*(-81267.558102849380495855461361462
3045 + y*(152376.67144284258842972899005274
3046 + y*(-134091.47086970147781816151124641
3047 + 44697.156956567159272720503748804*y)))))));
3048 case ( 43 ):
3049 return 6.7082039324993690892275210061938
3050 + y*(-375.65942021996466899674117634685
3051 + y*(5071.4021729695230314560058806825
3052 + y*(-28174.456516497350174755588226014
3053 + y*(77479.755420367712980577867621539
3054 + y*(-111570.84780532950669203212937502
3055 + (80578.9456371824214998009823264
3056 - 23022.555896337834714228852093257*y)*y)))))
3057 + x*(-13.416407864998738178455042012388
3058 + y*(751.31884043992933799348235269371
3059 + y*(-10142.804345939046062912011761365
3060 + y*(56348.913032994700349511176452028
3061 + y*(-154959.51084073542596115573524308
3062 + y*(223141.69561065901338406425875003
3063 + y*(-161157.8912743648429996019646528
3064 + 46045.111792675669428457704186514*y)))))));
3065 case ( 44 ):
3066 return 4.1231056256176605498214098559741
3067 + y*(-296.86360504447155958714150963013
3068 + y*(5195.1130882782522927749764185273
3069 + y*(-38097.4959807071834803498270692
3070 + y*(142865.6099276519380513118515095
3071 + y*(-297160.46864951603114672865113976
3072 + y*(346687.21342443536967118342632972
3073 + y*(-212257.47760679716510480617938555
3074 + 53064.36940169929127620154484639*y)))))));
3075 #endif
3076 }
3077
3078 printf("Error in eval_quadrilateral_2d -- unmatched switch/case %d\n",i);
3079 return -123456789.0;
3080 }
3081
3082
3083 /* \nabla\phi_i(x,y) for quadrilateral */
3084 static void
3085 grad_quadrilateral_2d ( const int i, DomainType xi, JacobianRangeType grad )
3086 {
3087 const RangeField &x = xi[0];
3088 const RangeField &y = xi[1];
3089
3090 switch (i) {
3091 #if (PMAX2D>=0)
3092 case ( 0 ):
3093 grad[0]= 0;
3094 grad[1]= 0;
3095 return;
3096 #endif
3097 #if (PMAX2D>=1)
3098 case ( 1 ):
3099 grad[0]= 3.4641016151377545870548926830117;
3100 grad[1]= 0;
3101 return;
3102 case ( 2 ):
3103 grad[0]= 0;
3104 grad[1]= 3.4641016151377545870548926830117;
3105 return;
3106 #endif
3107 #if (PMAX2D>=2)
3108 case ( 3 ):
3109 grad[0]= -13.416407864998738178455042012388
3110 + 26.832815729997476356910084024775*x;
3111 grad[1]= 0;
3112 return;
3113 case ( 4 ):
3114 grad[0]= -6.
3115 + 12.*y;
3116 grad[1]= -6.
3117 + 12.*x;
3118 return;
3119 case ( 5 ):
3120 grad[0]= 0;
3121 grad[1]= -13.416407864998738178455042012388
3122 + 26.832815729997476356910084024775*y;
3123 return;
3124 #endif
3125 #if (PMAX2D>=3)
3126 case ( 6 ):
3127 grad[0]= 31.749015732775087086019389043671
3128 + x*(-158.74507866387543543009694521836
3129 + 158.74507866387543543009694521836*x);
3130 grad[1]= 0;
3131 return;
3132 case ( 7 ):
3133 grad[0]= 23.237900077244501311075592398694
3134 - 46.475800154489002622151184797389*y
3135 + x*(-46.475800154489002622151184797389
3136 + 92.95160030897800524430236959478*y);
3137 grad[1]= 7.7459666924148337703585307995648
3138 + x*(-46.475800154489002622151184797389
3139 + 46.475800154489002622151184797389*x);
3140 return;
3141 case ( 8 ):
3142 grad[0]= 7.7459666924148337703585307995648
3143 + y*(-46.475800154489002622151184797389
3144 + 46.475800154489002622151184797389*y);
3145 grad[1]= 23.237900077244501311075592398694
3146 - 46.475800154489002622151184797389*y
3147 + x*(-46.475800154489002622151184797389
3148 + 92.95160030897800524430236959478*y);
3149 return;
3150 case ( 9 ):
3151 grad[0]= 0;
3152 grad[1]= 31.749015732775087086019389043671
3153 + y*(-158.74507866387543543009694521836
3154 + 158.74507866387543543009694521836*y);
3155 return;
3156 #endif
3157 #if (PMAX2D>=4)
3158 case ( 10 ):
3159 grad[0]= -60.
3160 + x*(540.
3161 + x*(-1260.
3162 + 840.*x));
3163 grad[1]= 0;
3164 return;
3165 case ( 11 ):
3166 grad[0]= -54.990908339470080079056566324736
3167 + 109.98181667894016015811313264947*y
3168 + x*(274.95454169735040039528283162368
3169 - 549.9090833947008007905656632474*y
3170 + x*(-274.95454169735040039528283162368
3171 + 549.9090833947008007905656632474*y));
3172 grad[1]= -9.165151389911680013176094387456
3173 + x*(109.98181667894016015811313264947
3174 + x*(-274.95454169735040039528283162368
3175 + 183.30302779823360026352188774912*x));
3176 return;
3177 case ( 12 ):
3178 grad[0]= -30.
3179 + (180.
3180 - 180.*y)*y
3181 + x*(60.
3182 + y*(-360.
3183 + 360.*y));
3184 grad[1]= -30.
3185 + 60.*y
3186 + x*(180.
3187 - 360.*y
3188 + x*(-180.
3189 + 360.*y));
3190 return;
3191 case ( 13 ):
3192 grad[0]= -9.165151389911680013176094387456
3193 + y*(109.98181667894016015811313264947
3194 + y*(-274.95454169735040039528283162368
3195 + 183.30302779823360026352188774912*y));
3196 grad[1]= -54.990908339470080079056566324736
3197 + (274.95454169735040039528283162368
3198 - 274.95454169735040039528283162368*y)*y
3199 + x*(109.98181667894016015811313264947
3200 + y*(-549.90908339470080079056566324736
3201 + 549.9090833947008007905656632474*y));
3202 return;
3203 case ( 14 ):
3204 grad[0]= 0;
3205 grad[1]= -60.
3206 + y*(540.
3207 + y*(-1260.
3208 + 840.*y));
3209 return;
3210 #endif
3211 #if (PMAX2D>=5)
3212 return;
3213 case ( 15 ):
3214 grad[0]= 99.498743710661995473447982100121
3215 + x*(-1392.9824119492679366282717494017
3216 + x*(5571.9296477970717465130869976068
3217 + x*(-8357.8944716956076197696304964101
3218 + 4178.9472358478038098848152482051*x)));
3219 grad[1]= 0;
3220 return;
3221 case ( 16 ):
3222 grad[0]= 103.92304845413263761164678049035
3223 - 207.8460969082652752232935609807*y
3224 + x*(-935.30743608719373850482102441317
3225 + 1870.6148721743874770096420488263*y
3226 + x*(2182.3840175367853898445823902974
3227 - 4364.7680350735707796891647805948*y
3228 + x*(-1454.9226783578569265630549268649
3229 + 2909.8453567157138531261098537299*y)));
3230 grad[1]= 10.392304845413263761164678049035
3231 + x*(-207.8460969082652752232935609807
3232 + x*(935.30743608719373850482102441317
3233 + x*(-1454.9226783578569265630549268649
3234 + 727.4613391789284632815274634325*x)));
3235 return;
3236 case ( 17 ):
3237 grad[0]= 70.992957397195392510807939498739
3238 + y*(-425.95774438317235506484763699244
3239 + 425.95774438317235506484763699244*y)
3240 + x*(-354.9647869859769625540396974937
3241 + (2129.7887219158617753242381849622
3242 - 2129.7887219158617753242381849622*y)*y
3243 + x*(354.9647869859769625540396974937
3244 + y*(-2129.7887219158617753242381849622
3245 + 2129.7887219158617753242381849622*y)));
3246 grad[1]= 35.49647869859769625540396974937
3247 - 70.99295739719539251080793949874*y
3248 + x*(-425.95774438317235506484763699244
3249 + 851.9154887663447101296952739849*y
3250 + x*(1064.8943609579308876621190924811
3251 - 2129.7887219158617753242381849622*y
3252 + x*(-709.92957397195392510807939498739
3253 + 1419.8591479439078502161587899748*y)));
3254 return;
3255 case ( 18 ):
3256 grad[0]= 35.49647869859769625540396974937
3257 + y*(-425.95774438317235506484763699244
3258 + (1064.8943609579308876621190924811
3259 - 709.9295739719539251080793949874*y)*y)
3260 + x*(-70.992957397195392510807939498739
3261 + y*(851.91548876634471012969527398487
3262 + y*(-2129.7887219158617753242381849622
3263 + 1419.8591479439078502161587899748*y)));
3264 grad[1]= 70.992957397195392510807939498739
3265 + y*(-354.9647869859769625540396974937
3266 + 354.9647869859769625540396974937*y)
3267 + x*(-425.95774438317235506484763699244
3268 + (2129.7887219158617753242381849622
3269 - 2129.7887219158617753242381849622*y)*y
3270 + x*(425.95774438317235506484763699244
3271 + y*(-2129.7887219158617753242381849622
3272 + 2129.7887219158617753242381849622*y)));
3273 return;
3274 case ( 19 ):
3275 grad[0]= 10.392304845413263761164678049035
3276 + y*(-207.8460969082652752232935609807
3277 + y*(935.30743608719373850482102441317
3278 + y*(-1454.9226783578569265630549268649
3279 + 727.4613391789284632815274634325*y)));
3280 grad[1]= 103.92304845413263761164678049035
3281 + y*(-935.30743608719373850482102441317
3282 + (2182.3840175367853898445823902974
3283 - 1454.9226783578569265630549268649*y)*y)
3284 + x*(-207.8460969082652752232935609807
3285 + y*(1870.6148721743874770096420488263
3286 + y*(-4364.7680350735707796891647805948
3287 + 2909.8453567157138531261098537299*y)));
3288 return;
3289 case ( 20 ):
3290 grad[0]= 0;
3291 grad[1]= 99.498743710661995473447982100121
3292 + y*(-1392.9824119492679366282717494017
3293 + y*(5571.9296477970717465130869976068
3294 + y*(-8357.8944716956076197696304964101
3295 + 4178.9472358478038098848152482051*y)));
3296 return;
3297 #endif
3298 #if (PMAX2D>=6)
3299 case ( 21 ):
3300 grad[0]= -151.43315356948755031100729323376
3301 + x*(3028.6630713897510062201458646752
3302 + x*(-18171.978428338506037320875188051
3303 + x*(45429.946070846265093302187970128
3304 + x*(-49972.940677930891602632406767141
3305 + 19989.176271172356641052962706856*x))));
3306 grad[1]= 0;
3307 return;
3308 case ( 22 ):
3309 grad[0]= -172.33687939614085979551834404657
3310 + 344.67375879228171959103668809314*y
3311 + x*(2412.716311545972037137256816652
3312 - 4825.4326230919440742745136333039*y
3313 + x*(-9650.8652461838881485490272666078
3314 + 19301.730492367776297098054533216*y
3315 + x*(14476.297869275832222823540899912
3316 - 28952.595738551664445647081799823*y
3317 + x*(-7238.1489346379161114117704499559
3318 + 14476.297869275832222823540899912*y))));
3319 grad[1]= -11.489125293076057319701222936438
3320 + x*(344.67375879228171959103668809314
3321 + x*(-2412.716311545972037137256816652
3322 + x*(6433.9101641225920990326848444052
3323 + x*(-7238.1489346379161114117704499559
3324 + 2895.2595738551664445647081799823*x))));
3325 return;
3326 case ( 23 ):
3327 grad[0]= -134.16407864998738178455042012388
3328 + (804.98447189992429070730252074326
3329 - 804.9844718999242907073025207433*y)*y
3330 + x*(1207.4767078498864360609537811149
3331 + y*(-7244.8602470993186163657226866893
3332 + 7244.860247099318616365722686689*y)
3333 + x*(-2817.4456516497350174755588226014
3334 + (16904.673909898410104853352935608
3335 - 16904.673909898410104853352935608*y)*y
3336 + x*(1878.2971010998233449837058817343
3337 + y*(-11269.782606598940069902235290406
3338 + 11269.782606598940069902235290406*y))));
3339 grad[1]= -40.249223594996214535365126037163
3340 + 80.49844718999242907073025207433*y
3341 + x*(804.98447189992429070730252074326
3342 - 1609.9689437998485814146050414865*y
3343 + x*(-3622.4301235496593081828613433447
3344 + 7244.860247099318616365722686689*y
3345 + x*(5634.8913032994700349511176452028
3346 - 11269.782606598940069902235290406*y
3347 + x*(-2817.4456516497350174755588226014
3348 + 5634.891303299470034951117645203*y))));
3349 return;
3350 case ( 24 ):
3351 grad[0]= -84.
3352 + y*(1008.
3353 + y*(-2520.
3354 + 1680.*y))
3355 + x*(420.
3356 + y*(-5040.
3357 + (12600.
3358 - 8400.*y)*y)
3359 + x*(-420.
3360 + y*(5040.
3361 + y*(-12600.
3362 + 8400.*y))));
3363 grad[1]= -84.
3364 + (420.
3365 - 420.*y)*y
3366 + x*(1008.
3367 + y*(-5040.
3368 + 5040.*y)
3369 + x*(-2520.
3370 + (12600.
3371 - 12600.*y)*y
3372 + x*(1680.
3373 + y*(-8400.
3374 + 8400.*y))));
3375 return;
3376 case ( 25 ):
3377 grad[0]= -40.249223594996214535365126037163
3378 + y*(804.98447189992429070730252074326
3379 + y*(-3622.4301235496593081828613433447
3380 + (5634.8913032994700349511176452028
3381 - 2817.4456516497350174755588226014*y)*y))
3382 + x*(80.498447189992429070730252074326
3383 + y*(-1609.9689437998485814146050414865
3384 + y*(7244.8602470993186163657226866893
3385 + y*(-11269.782606598940069902235290406
3386 + 5634.891303299470034951117645203*y))));
3387 grad[1]= -134.16407864998738178455042012388
3388 + y*(1207.4767078498864360609537811149
3389 + y*(-2817.4456516497350174755588226014
3390 + 1878.2971010998233449837058817343*y))
3391 + x*(804.98447189992429070730252074326
3392 + y*(-7244.8602470993186163657226866893
3393 + (16904.673909898410104853352935608
3394 - 11269.782606598940069902235290406*y)*y)
3395 + x*(-804.98447189992429070730252074326
3396 + y*(7244.8602470993186163657226866893
3397 + y*(-16904.673909898410104853352935608
3398 + 11269.782606598940069902235290406*y))));
3399 return;
3400 case ( 26 ):
3401 grad[0]= -11.489125293076057319701222936438
3402 + y*(344.67375879228171959103668809314
3403 + y*(-2412.716311545972037137256816652
3404 + y*(6433.9101641225920990326848444052
3405 + y*(-7238.1489346379161114117704499559
3406 + 2895.2595738551664445647081799823*y))));
3407 grad[1]= -172.33687939614085979551834404657
3408 + y*(2412.716311545972037137256816652
3409 + y*(-9650.8652461838881485490272666078
3410 + (14476.297869275832222823540899912
3411 - 7238.148934637916111411770449956*y)*y))
3412 + x*(344.67375879228171959103668809314
3413 + y*(-4825.4326230919440742745136333039
3414 + y*(19301.730492367776297098054533216
3415 + y*(-28952.595738551664445647081799823
3416 + 14476.297869275832222823540899912*y))));
3417 return;
3418 case ( 27 ):
3419 grad[0]= 0;
3420 grad[1]= -151.43315356948755031100729323376
3421 + y*(3028.6630713897510062201458646752
3422 + y*(-18171.978428338506037320875188051
3423 + y*(45429.946070846265093302187970128
3424 + y*(-49972.940677930891602632406767141
3425 + 19989.176271172356641052962706856*y))));
3426 return;
3427 #endif
3428 #if (PMAX2D>=7)
3429 case ( 28 ):
3430 grad[0]= 216.88706738761534557003886238781
3431 + x*(-5855.950819465614330391049284471
3432 + x*(48799.590162213452753258744037258
3433 + x*(-178931.83059478266009528206146995
3434 + x*(322077.2950706087881715077106459
3435 + x*(-279133.65572786094974864001589312
3436 + 93044.55190928698324954667196437*x)))));
3437 grad[1]= 0;
3438 return;
3439 case ( 29 ):
3440 grad[0]= 262.28991593273272464556951107947
3441 - 524.5798318654654492911390221589*y
3442 + x*(-5245.7983186546544929113902215894
3443 + 10491.596637309308985822780443179*y
3444 + x*(31474.789911927926957468341329537
3445 - 62949.57982385585391493668265907*y
3446 + x*(-78686.974779819817393670853323841
3447 + 157373.94955963963478734170664768*y
3448 + x*(86555.672257801799133037938656226
3449 - 173111.34451560359826607587731245*y
3450 + x*(-34622.26890312071965321517546249
3451 + 69244.53780624143930643035092498*y)))));
3452 grad[1]= 12.48999599679679641169378624188
3453 + x*(-524.57983186546544929113902215894
3454 + x*(5245.7983186546544929113902215894
3455 + x*(-20983.193274618617971645560886358
3456 + x*(39343.487389909908696835426661921
3457 + x*(-34622.26890312071965321517546249
3458 + 11540.756301040239884405058487497*x)))));
3459 return;
3460 case ( 30 ):
3461 grad[0]= 222.48595461286988846134192322402
3462 + y*(-1334.9157276772193307680515393441
3463 + 1334.9157276772193307680515393441*y)
3464 + x*(-3114.8033645801784384587869251363
3465 + (18688.820187481070630752721550818
3466 - 18688.820187481070630752721550818*y)*y
3467 + x*(12459.213458320713753835147700545
3468 + y*(-74755.280749924282523010886203271
3469 + 74755.28074992428252301088620327*y)
3470 + x*(-18688.820187481070630752721550818
3471 + (112132.92112488642378451632930491
3472 - 112132.92112488642378451632930491*y)*y
3473 + x*(9344.4100937405353153763607754089
3474 + y*(-56066.460562443211892258164652453
3475 + 56066.46056244321189225816465245*y)))));
3476 grad[1]= 44.497190922573977692268384644804
3477 - 88.99438184514795538453676928961*y
3478 + x*(-1334.9157276772193307680515393441
3479 + 2669.8314553544386615361030786883*y
3480 + x*(9344.4100937405353153763607754089
3481 - 18688.820187481070630752721550818*y
3482 + x*(-24918.42691664142750767029540109
3483 + 49836.853833282855015340590802181*y
3484 + x*(28033.230281221605946129082326227
3485 - 56066.46056244321189225816465245*y
3486 + x*(-11213.292112488642378451632930491
3487 + 22426.584224977284756903265860981*y)))));
3488 return;
3489 case ( 31 ):
3490 grad[0]= 158.74507866387543543009694521836
3491 + y*(-1904.9409439665052251611633426203
3492 + (4762.3523599162630629029083565507
3493 - 3174.9015732775087086019389043671*y)*y)
3494 + x*(-1428.7057079748789188708725069652
3495 + y*(17144.468495698547026450470083582
3496 + y*(-42861.171239246367566126175208956
3497 + 28574.114159497578377417450139304*y))
3498 + x*(3333.6466519413841440320358495855
3499 + y*(-40003.759823296609728384430195026
3500 + (100009.39955824152432096107548756
3501 - 66672.93303882768288064071699171*y)*y)
3502 + x*(-2222.431101294256096021357233057
3503 + y*(26669.173215531073152256286796684
3504 + y*(-66672.933038827682880640716991709
3505 + 44448.62202588512192042714466114*y)))));
3506 grad[1]= 95.247047198325261258058167131013
3507 + y*(-476.23523599162630629029083565507
3508 + 476.23523599162630629029083565507*y)
3509 + x*(-1904.9409439665052251611633426203
3510 + (9524.7047198325261258058167131013
3511 - 9524.704719832526125805816713101*y)*y
3512 + x*(8572.2342478492735132252350417912
3513 + y*(-42861.171239246367566126175208956
3514 + 42861.171239246367566126175208956*y)
3515 + x*(-13334.586607765536576128143398342
3516 + (66672.933038827682880640716991709
3517 - 66672.93303882768288064071699171*y)*y
3518 + x*(6667.2933038827682880640716991709
3519 + y*(-33336.466519413841440320358495855
3520 + 33336.466519413841440320358495855*y)))));
3521 return;
3522 case ( 32 ):
3523 grad[0]= 95.247047198325261258058167131013
3524 + y*(-1904.9409439665052251611633426203
3525 + y*(8572.2342478492735132252350417912
3526 + y*(-13334.586607765536576128143398342
3527 + 6667.293303882768288064071699171*y)))
3528 + x*(-476.23523599162630629029083565507
3529 + y*(9524.7047198325261258058167131013
3530 + y*(-42861.171239246367566126175208956
3531 + (66672.933038827682880640716991709
3532 - 33336.466519413841440320358495855*y)*y))
3533 + x*(476.23523599162630629029083565507
3534 + y*(-9524.7047198325261258058167131013
3535 + y*(42861.171239246367566126175208956
3536 + y*(-66672.933038827682880640716991709
3537 + 33336.466519413841440320358495855*y)))));
3538 grad[1]= 158.74507866387543543009694521836
3539 + y*(-1428.7057079748789188708725069652
3540 + (3333.6466519413841440320358495855
3541 - 2222.431101294256096021357233057*y)*y)
3542 + x*(-1904.9409439665052251611633426203
3543 + y*(17144.468495698547026450470083582
3544 + y*(-40003.759823296609728384430195026
3545 + 26669.173215531073152256286796684*y))
3546 + x*(4762.3523599162630629029083565507
3547 + y*(-42861.171239246367566126175208956
3548 + (100009.39955824152432096107548756
3549 - 66672.93303882768288064071699171*y)*y)
3550 + x*(-3174.9015732775087086019389043671
3551 + y*(28574.114159497578377417450139304
3552 + y*(-66672.933038827682880640716991709
3553 + 44448.62202588512192042714466114*y)))));
3554 return;
3555 case ( 33 ):
3556 grad[0]= 44.497190922573977692268384644804
3557 + y*(-1334.9157276772193307680515393441
3558 + y*(9344.4100937405353153763607754089
3559 + y*(-24918.42691664142750767029540109
3560 + (28033.230281221605946129082326227
3561 - 11213.292112488642378451632930491*y)*y)))
3562 + x*(-88.994381845147955384536769289609
3563 + y*(2669.8314553544386615361030786883
3564 + y*(-18688.820187481070630752721550818
3565 + y*(49836.853833282855015340590802181
3566 + y*(-56066.460562443211892258164652453
3567 + 22426.584224977284756903265860981*y)))));
3568 grad[1]= 222.48595461286988846134192322402
3569 + y*(-3114.8033645801784384587869251363
3570 + y*(12459.213458320713753835147700545
3571 + y*(-18688.820187481070630752721550818
3572 + 9344.410093740535315376360775409*y)))
3573 + x*(-1334.9157276772193307680515393441
3574 + y*(18688.820187481070630752721550818
3575 + y*(-74755.280749924282523010886203271
3576 + (112132.92112488642378451632930491
3577 - 56066.46056244321189225816465245*y)*y))
3578 + x*(1334.9157276772193307680515393441
3579 + y*(-18688.820187481070630752721550818
3580 + y*(74755.280749924282523010886203271
3581 + y*(-112132.92112488642378451632930491
3582 + 56066.46056244321189225816465245*y)))));
3583 return;
3584 case ( 34 ):
3585 grad[0]= 12.48999599679679641169378624188
3586 + y*(-524.57983186546544929113902215894
3587 + y*(5245.7983186546544929113902215894
3588 + y*(-20983.193274618617971645560886358
3589 + y*(39343.487389909908696835426661921
3590 + y*(-34622.26890312071965321517546249
3591 + 11540.756301040239884405058487497*y)))));
3592 grad[1]= 262.28991593273272464556951107947
3593 + y*(-5245.7983186546544929113902215894
3594 + y*(31474.789911927926957468341329537
3595 + y*(-78686.974779819817393670853323841
3596 + (86555.672257801799133037938656226
3597 - 34622.26890312071965321517546249*y)*y)))
3598 + x*(-524.57983186546544929113902215894
3599 + y*(10491.596637309308985822780443179
3600 + y*(-62949.579823855853914936682659073
3601 + y*(157373.94955963963478734170664768
3602 + y*(-173111.34451560359826607587731245
3603 + 69244.53780624143930643035092498*y)))));
3604 return;
3605 case ( 35 ):
3606 grad[0]= 0;
3607 grad[1]= 216.88706738761534557003886238781
3608 + y*(-5855.950819465614330391049284471
3609 + y*(48799.590162213452753258744037258
3610 + y*(-178931.83059478266009528206146995
3611 + y*(322077.2950706087881715077106459
3612 + y*(-279133.65572786094974864001589312
3613 + 93044.55190928698324954667196437*y)))));
3614 return;
3615 #endif
3616 #if (PMAX2D>=8)
3617 case ( 36 ):
3618 grad[0]= -296.86360504447155958714150963013
3619 + x*(10390.226176556504585549952837055
3620 + x*(-114292.4879421215504410494812076
3621 + x*(571462.43971060775220524740603801
3622 + x*(-1.4858023432475801557336432556988e6
3623 + x*(2.0801232805466122180271005579783e6
3624 + x*(-1.4858023432475801557336432556988e6
3625 + 424514.95521359433020961235877109*x))))));
3626 grad[1]= 0;
3627 return;
3628 case ( 37 ):
3629 grad[0]= -375.65942021996466899674117634685
3630 + 751.3188404399293379934823526937*y
3631 + x*(10142.804345939046062912011761365
3632 - 20285.60869187809212582402352273*y
3633 + x*(-84523.369549492050524266764678042
3634 + 169046.73909898410104853352935608*y
3635 + x*(309919.02168147085192231147048615
3636 - 619838.0433629417038446229409723*y
3637 + x*(-557854.23902664753346016064687508
3638 + 1.1157084780532950669203212937502e6*y
3639 + x*(483473.6738230945289988058939584
3640 - 966947.3476461890579976117879168*y
3641 + x*(-161157.8912743648429996019646528
3642 + 322315.7825487296859992039293056*y))))));
3643 grad[1]= -13.416407864998738178455042012388
3644 + x*(751.31884043992933799348235269371
3645 + x*(-10142.804345939046062912011761365
3646 + x*(56348.913032994700349511176452028
3647 + x*(-154959.51084073542596115573524308
3648 + x*(223141.69561065901338406425875003
3649 + x*(-161157.8912743648429996019646528
3650 + 46045.111792675669428457704186514*x))))));
3651 return;
3652 case ( 38 ):
3653 grad[0]= -338.61482542853908539939775567276
3654 + (2031.6889525712345123963865340366
3655 - 2031.6889525712345123963865340366*y)*y
3656 + x*(6772.2965085707817079879551134552
3657 + y*(-40633.779051424690247927730680731
3658 + 40633.779051424690247927730680731*y)
3659 + x*(-40633.779051424690247927730680731
3660 + (243802.67430854814148756638408439
3661 - 243802.67430854814148756638408439*y)*y
3662 + x*(101584.44762856172561981932670183
3663 + y*(-609506.68577137035371891596021097
3664 + 609506.685771370353718915960211*y)
3665 + x*(-111742.89239141789818180125937201
3666 + (670457.35434850738909080755623206
3667 - 670457.3543485073890908075562321*y)*y
3668 + x*(44697.156956567159272720503748804
3669 + y*(-268182.94173940295563632302249282
3670 + 268182.94173940295563632302249282*y))))));
3671 grad[1]= -48.373546489791297914199679381823
3672 + 96.74709297958259582839935876365*y
3673 + x*(2031.6889525712345123963865340366
3674 - 4063.3779051424690247927730680731*y
3675 + x*(-20316.889525712345123963865340366
3676 + 40633.779051424690247927730680731*y
3677 + x*(81267.558102849380495855461361462
3678 - 162535.11620569876099171092272292*y
3679 + x*(-152376.67144284258842972899005274
3680 + 304753.34288568517685945798010548*y
3681 + x*(134091.47086970147781816151124641
3682 - 268182.94173940295563632302249282*y
3683 + x*(-44697.156956567159272720503748804
3684 + 89394.31391313431854544100749761*y))))));
3685 return;
3686 case ( 39 ):
3687 grad[0]= -263.24893162176366181219164922249
3688 + y*(3158.9871794611639417462997906699
3689 + y*(-7897.4679486529098543657494766747
3690 + 5264.97863243527323624383298445*y))
3691 + x*(3685.4850427046912653706830891149
3692 + y*(-44225.820512456295184448197069378
3693 + (110564.55128114073796112049267345
3694 - 73709.7008540938253074136617823*y)*y)
3695 + x*(-14741.940170818765061482732356459
3696 + y*(176903.28204982518073779278827751
3697 + y*(-442258.20512456295184448197069378
3698 + 294838.80341637530122965464712919*y))
3699 + x*(22112.910256228147592224098534689
3700 + y*(-265354.92307473777110668918241627
3701 + (663387.30768684442776672295604067
3702 - 442258.20512456295184448197069378*y)*y)
3703 + x*(-11056.455128114073796112049267345
3704 + y*(132677.46153736888555334459120813
3705 + y*(-331693.65384342221388336147802034
3706 + 221129.10256228147592224098534689*y))))));
3707 grad[1]= -105.299572648705464724876659689
3708 + (526.49786324352732362438329844498
3709 - 526.497863243527323624383298445*y)*y
3710 + x*(3158.9871794611639417462997906699
3711 + y*(-15794.935897305819708731498953349
3712 + 15794.935897305819708731498953349*y)
3713 + x*(-22112.910256228147592224098534689
3714 + (110564.55128114073796112049267345
3715 - 110564.55128114073796112049267345*y)*y
3716 + x*(58967.760683275060245930929425838
3717 + y*(-294838.80341637530122965464712919
3718 + 294838.80341637530122965464712919*y)
3719 + x*(-66338.730768684442776672295604067
3720 + (331693.65384342221388336147802034
3721 - 331693.65384342221388336147802034*y)*y
3722 + x*(26535.492307473777110668918241627
3723 + y*(-132677.46153736888555334459120813
3724 + 132677.46153736888555334459120813*y))))));
3725 return;
3726 case ( 40 ):
3727 grad[0]= -180.
3728 + y*(3600.
3729 + y*(-16200.
3730 + (25200.
3731 - 12600.*y)*y))
3732 + x*(1620.
3733 + y*(-32400.
3734 + y*(145800.
3735 + y*(-226800.
3736 + 113400.*y)))
3737 + x*(-3780.
3738 + y*(75600.
3739 + y*(-340200.
3740 + (529200.
3741 - 264600.*y)*y))
3742 + x*(2520.
3743 + y*(-50400.
3744 + y*(226800.
3745 + y*(-352800.
3746 + 176400.*y))))));
3747 grad[1]= -180.
3748 + y*(1620.
3749 + y*(-3780.
3750 + 2520.*y))
3751 + x*(3600.
3752 + y*(-32400.
3753 + (75600.
3754 - 50400.*y)*y)
3755 + x*(-16200.
3756 + y*(145800.
3757 + y*(-340200.
3758 + 226800.*y))
3759 + x*(25200.
3760 + y*(-226800.
3761 + (529200.
3762 - 352800.*y)*y)
3763 + x*(-12600.
3764 + y*(113400.
3765 + y*(-264600.
3766 + 176400.*y))))));
3767 return;
3768 case ( 41 ):
3769 grad[0]= -105.299572648705464724876659689
3770 + y*(3158.9871794611639417462997906699
3771 + y*(-22112.910256228147592224098534689
3772 + y*(58967.760683275060245930929425838
3773 + y*(-66338.730768684442776672295604067
3774 + 26535.492307473777110668918241627*y))))
3775 + x*(526.49786324352732362438329844498
3776 + y*(-15794.935897305819708731498953349
3777 + y*(110564.55128114073796112049267345
3778 + y*(-294838.80341637530122965464712919
3779 + (331693.65384342221388336147802034
3780 - 132677.46153736888555334459120813*y)*y)))
3781 + x*(-526.49786324352732362438329844498
3782 + y*(15794.935897305819708731498953349
3783 + y*(-110564.55128114073796112049267345
3784 + y*(294838.80341637530122965464712919
3785 + y*(-331693.65384342221388336147802034
3786 + 132677.46153736888555334459120813*y))))));
3787 grad[1]= -263.24893162176366181219164922249
3788 + y*(3685.4850427046912653706830891149
3789 + y*(-14741.940170818765061482732356459
3790 + (22112.910256228147592224098534689
3791 - 11056.455128114073796112049267345*y)*y))
3792 + x*(3158.9871794611639417462997906699
3793 + y*(-44225.820512456295184448197069378
3794 + y*(176903.28204982518073779278827751
3795 + y*(-265354.92307473777110668918241627
3796 + 132677.46153736888555334459120813*y)))
3797 + x*(-7897.4679486529098543657494766747
3798 + y*(110564.55128114073796112049267345
3799 + y*(-442258.20512456295184448197069378
3800 + (663387.30768684442776672295604067
3801 - 331693.65384342221388336147802034*y)*y))
3802 + x*(5264.9786324352732362438329844498
3803 + y*(-73709.700854093825307413661782297
3804 + y*(294838.80341637530122965464712919
3805 + y*(-442258.20512456295184448197069378
3806 + 221129.10256228147592224098534689*y))))));
3807 return;
3808 case ( 42 ):
3809 grad[0]= -48.373546489791297914199679381823
3810 + y*(2031.6889525712345123963865340366
3811 + y*(-20316.889525712345123963865340366
3812 + y*(81267.558102849380495855461361462
3813 + y*(-152376.67144284258842972899005274
3814 + (134091.47086970147781816151124641
3815 - 44697.156956567159272720503748804*y)*y))))
3816 + x*(96.747092979582595828399358763645
3817 + y*(-4063.3779051424690247927730680731
3818 + y*(40633.779051424690247927730680731
3819 + y*(-162535.11620569876099171092272292
3820 + y*(304753.34288568517685945798010548
3821 + y*(-268182.94173940295563632302249282
3822 + 89394.31391313431854544100749761*y))))));
3823 grad[1]= -338.61482542853908539939775567276
3824 + y*(6772.2965085707817079879551134552
3825 + y*(-40633.779051424690247927730680731
3826 + y*(101584.44762856172561981932670183
3827 + y*(-111742.89239141789818180125937201
3828 + 44697.156956567159272720503748804*y))))
3829 + x*(2031.6889525712345123963865340366
3830 + y*(-40633.779051424690247927730680731
3831 + y*(243802.67430854814148756638408439
3832 + y*(-609506.68577137035371891596021097
3833 + (670457.35434850738909080755623206
3834 - 268182.94173940295563632302249282*y)*y)))
3835 + x*(-2031.6889525712345123963865340366
3836 + y*(40633.779051424690247927730680731
3837 + y*(-243802.67430854814148756638408439
3838 + y*(609506.68577137035371891596021097
3839 + y*(-670457.35434850738909080755623206
3840 + 268182.94173940295563632302249282*y))))));
3841 return;
3842 case ( 43 ):
3843 grad[0]= -13.416407864998738178455042012388
3844 + y*(751.31884043992933799348235269371
3845 + y*(-10142.804345939046062912011761365
3846 + y*(56348.913032994700349511176452028
3847 + y*(-154959.51084073542596115573524308
3848 + y*(223141.69561065901338406425875003
3849 + y*(-161157.8912743648429996019646528
3850 + 46045.111792675669428457704186514*y))))));
3851 grad[1]= -375.65942021996466899674117634685
3852 + y*(10142.804345939046062912011761365
3853 + y*(-84523.369549492050524266764678042
3854 + y*(309919.02168147085192231147048615
3855 + y*(-557854.23902664753346016064687508
3856 + (483473.6738230945289988058939584
3857 - 161157.8912743648429996019646528*y)*y))))
3858 + x*(751.31884043992933799348235269371
3859 + y*(-20285.60869187809212582402352273
3860 + y*(169046.73909898410104853352935608
3861 + y*(-619838.04336294170384462294097231
3862 + y*(1.1157084780532950669203212937502e6
3863 + y*(-966947.3476461890579976117879168
3864 + 322315.7825487296859992039293056*y))))));
3865 return;
3866 case ( 44 ):
3867 grad[0]= 0;
3868 grad[1]= -296.86360504447155958714150963013
3869 + y*(10390.226176556504585549952837055
3870 + y*(-114292.4879421215504410494812076
3871 + y*(571462.43971060775220524740603801
3872 + y*(-1.4858023432475801557336432556988e6
3873 + y*(2.0801232805466122180271005579783e6
3874 + y*(-1.4858023432475801557336432556988e6
3875 + 424514.95521359433020961235877109*y))))));
3876 return;
3877 #endif
3878 }
3879 }
3880
3881 /* H\phi_i(x,y) for quadrilateral */
3882 static void
3883 hess_quadrilateral_2d ( const int i, DomainType xi, HessianRangeType &h )
3884 {
3885 const RangeField &x = xi[0];
3886 const RangeField &y = xi[1];
3887
3888 switch( i )
3889 {
3890 #if (PMAX2D>=0)
3891 case ( 0 ):
3892 h[0] = 0;
3893 h[1] = 0;
3894 h[2] = 0;
3895 return;
3896 #endif
3897 #if (PMAX2D>=1)
3898 case ( 1 ):
3899 h[0] = 0;
3900 h[1] = 0;
3901 h[2] = 0;
3902 return;
3903 case ( 2 ):
3904 h[0] = 0;
3905 h[1] = 0;
3906 h[2] = 0;
3907 return;
3908 #endif
3909 #if (PMAX2D>=2)
3910 case ( 3 ):
3911 h[0] = 26.832815729997476356910084024776;
3912 h[1] = 0;
3913 h[2] = 0;
3914 return;
3915 case ( 4 ):
3916 h[0] = 0;
3917 h[1] = 12.;
3918 h[2] = 0;
3919 return;
3920 case ( 5 ):
3921 h[0] = 0;
3922 h[1] = 0;
3923 h[2] = 26.832815729997476356910084024776;
3924 return;
3925 #endif
3926 #if (PMAX2D>=3)
3927 case ( 6 ):
3928 h[0] = -158.74507866387543543009694521836+317.49015732775087086019389043674*x;
3929 h[1] = 0;
3930 h[2] = 0;
3931 return;
3932 case ( 7 ):
3933 h[0] = -46.475800154489002622151184797388+92.951600308978005244302369594778*y;
3934 h[1] = -46.475800154489002622151184797389+92.951600308978005244302369594778*x;
3935 h[2] = 0;
3936 return;
3937 case ( 8 ):
3938 h[0] = 0;
3939 h[1] = -46.475800154489002622151184797389+92.951600308978005244302369594778*y;
3940 h[2] = -46.475800154489002622151184797388+92.951600308978005244302369594778*x;
3941 return;
3942 case ( 9 ):
3943 h[0] = 0;
3944 h[1] = 0;
3945 h[2] = -158.74507866387543543009694521836+317.49015732775087086019389043674*y;
3946 return;
3947 #endif
3948 #if (PMAX2D>=4)
3949 case ( 10 ):
3950 h[0] = 540.+(-2520.+2520.*x)*x;
3951 h[1] = 0;
3952 h[2] = 0;
3953 return;
3954 case ( 11 ):
3955 h[0] = 274.95454169735040039528283162368-549.90908339470080079056566324736*y+(-549.90908339470080079056566324736+1099.8181667894016015811313264947*y)*x;
3956 h[1] = 109.98181667894016015811313264947+(-549.90908339470080079056566324736+549.90908339470080079056566324736*x)*x;
3957 h[2] = 0;
3958 return;
3959 case ( 12 ):
3960 h[0] = 60.+(-360.+360.*y)*y;
3961 h[1] = -360.*y+180.+(-360.+720.*y)*x;
3962 h[2] = 60.+x*(-360.+360.*x);
3963 return;
3964 case ( 13 ):
3965 h[0] = 0;
3966 h[1] = 109.98181667894016015811313264947+(-549.90908339470080079056566324736+549.90908339470080079056566324736*y)*y;
3967 h[2] = 274.95454169735040039528283162368-549.90908339470080079056566324736*y+(-549.90908339470080079056566324736+1099.8181667894016015811313264947*y)*x;
3968 return;
3969 case ( 14 ):
3970 h[0] = 0;
3971 h[1] = 0;
3972 h[2] = 540.+(-2520.+2520.*y)*y;
3973 return;
3974 #endif
3975 #if (PMAX2D>=5)
3976 case ( 15 ):
3977 h[0] = -1392.9824119492679366282717494017+(11143.859295594143493026173995214+(-25073.683415086822859308891489230+16715.788943391215239539260992820*x)*x)*x;
3978 h[1] = 0;
3979 h[2] = 0;
3980 return;
3981 case ( 16 ):
3982 h[0] = -935.30743608719373850482102441318+1870.6148721743874770096420488264*y+(4364.7680350735707796891647805947-8729.5360701471415593783295611894*y+(-4364.7680350735707796891647805948+8729.5360701471415593783295611900*y)*x)*x;
3983 h[1] = -207.8460969082652752232935609807+(1870.6148721743874770096420488264+(-4364.7680350735707796891647805947+2909.8453567157138531261098537300*x)*x)*x;
3984 h[2] = 0;
3985 return;
3986 case ( 17 ):
3987 h[0] = -354.96478698597696255403969749370+(2129.7887219158617753242381849622-2129.7887219158617753242381849622*y)*y+(709.92957397195392510807939498738+(-4259.5774438317235506484763699243+4259.5774438317235506484763699244*y)*y)*x;
3988 h[1] = -425.95774438317235506484763699244+851.91548876634471012969527398488*y+(-4259.5774438317235506484763699244*y+2129.7887219158617753242381849622+(-2129.7887219158617753242381849622+4259.5774438317235506484763699244*y)*x)*x;
3989 h[2] = -70.99295739719539251080793949874+x*(851.91548876634471012969527398488+x*(-2129.7887219158617753242381849622+1419.8591479439078502161587899748*x));
3990 return;
3991 case ( 18 ):
3992 h[0] = -70.99295739719539251080793949874+(851.91548876634471012969527398488+(-2129.7887219158617753242381849622+1419.8591479439078502161587899748*y)*y)*y;
3993 h[1] = -425.95774438317235506484763699244+(2129.7887219158617753242381849622-2129.7887219158617753242381849622*y)*y+(851.91548876634471012969527398488+(-4259.5774438317235506484763699244+4259.5774438317235506484763699244*y)*y)*x;
3994 h[2] = -354.96478698597696255403969749370+709.92957397195392510807939498738*y+(-4259.5774438317235506484763699244*y+2129.7887219158617753242381849622+(-2129.7887219158617753242381849622+4259.5774438317235506484763699244*y)*x)*x;
3995 return;
3996 case ( 19 ):
3997 h[0] = 0;
3998 h[1] = -207.8460969082652752232935609807+(1870.6148721743874770096420488263+(-4364.7680350735707796891647805947+2909.8453567157138531261098537300*y)*y)*y;
3999 h[2] = -935.30743608719373850482102441318+(4364.7680350735707796891647805947-4364.7680350735707796891647805948*y)*y+(1870.6148721743874770096420488263+(-8729.5360701471415593783295611894+8729.5360701471415593783295611900*y)*y)*x;
4000 return;
4001 case ( 20 ):
4002 h[0] = 0;
4003 h[1] = 0;
4004 h[2] = -1392.9824119492679366282717494017+(11143.859295594143493026173995214+(-25073.683415086822859308891489230+16715.788943391215239539260992820*y)*y)*y;
4005 return;
4006 #endif
4007 #if (PMAX2D>=6)
4008 case ( 21 ):
4009 h[0] = 3028.6630713897510062201458646752+(-36343.956856677012074641750376103+(136289.83821253879527990656391038+(-199891.76271172356641052962706856+99945.881355861783205264813534280*x)*x)*x)*x;
4010 h[1] = 0;
4011 h[2] = 0;
4012 return;
4013 case ( 22 ):
4014 h[0] = 2412.716311545972037137256816652-4825.432623091944074274513633304*y+(-19301.730492367776297098054533215+38603.460984735552594196109066430*y+(43428.893607827496668470622699736-86857.787215654993336941245399472*y+(-28952.595738551664445647081799824+57905.191477103328891294163599647*y)*x)*x)*x;
4015 h[1] = 344.67375879228171959103668809314+(-4825.432623091944074274513633304+(19301.730492367776297098054533215+(-28952.595738551664445647081799824+14476.297869275832222823540899911*x)*x)*x)*x;
4016 h[2] = 0;
4017 return;
4018 case ( 23 ):
4019 h[0] = 1207.4767078498864360609537811149+(-7244.8602470993186163657226866894+7244.8602470993186163657226866894*y)*y+(-5634.8913032994700349511176452029+(33809.347819796820209706705871217-33809.347819796820209706705871218*y)*y+(5634.8913032994700349511176452028+(-33809.347819796820209706705871217+33809.347819796820209706705871217*y)*y)*x)*x;
4020 h[1] = -1609.9689437998485814146050414866*y+804.98447189992429070730252074326+(-7244.8602470993186163657226866894+14489.720494198637232731445373379*y+(-33809.347819796820209706705871218*y+16904.673909898410104853352935609+(-11269.782606598940069902235290406+22539.565213197880139804470580812*y)*x)*x)*x;
4021 h[2] = 80.498447189992429070730252074326+x*(-1609.9689437998485814146050414866+x*(7244.8602470993186163657226866894+x*(-11269.782606598940069902235290406+5634.8913032994700349511176452028*x)));
4022 return;
4023 case ( 24 ):
4024 h[0] = 420.+(-5040.+(12600.-8400.*y)*y)*y+(-840.+(10080.+(-25200.+16800.*y)*y)*y)*x;
4025 h[1] = 1008.+(-5040.+5040.*y)*y+(-5040.+(25200.-25200.*y)*y+(5040.+(-25200.+25200.*y)*y)*x)*x;
4026 h[2] = -840.*y+420.+x*(-5040.+10080.*y+x*(-25200.*y+12600.+x*(-8400.+16800.*y)));
4027 return;
4028 case ( 25 ):
4029 h[0] = 80.498447189992429070730252074326+(-1609.9689437998485814146050414865+(7244.8602470993186163657226866894+(-11269.782606598940069902235290406+5634.8913032994700349511176452028*y)*y)*y)*y;
4030 h[1] = 804.98447189992429070730252074326+(-7244.8602470993186163657226866894+(16904.673909898410104853352935609-11269.782606598940069902235290406*y)*y)*y+(-1609.9689437998485814146050414865+(14489.720494198637232731445373379+(-33809.347819796820209706705871218+22539.565213197880139804470580811*y)*y)*y)*x;
4031 h[2] = 1207.4767078498864360609537811149+(-5634.8913032994700349511176452029+5634.8913032994700349511176452028*y)*y+(-7244.8602470993186163657226866894+(33809.347819796820209706705871218-33809.347819796820209706705871218*y)*y+(7244.8602470993186163657226866894+(-33809.347819796820209706705871218+33809.347819796820209706705871218*y)*y)*x)*x;
4032 return;
4033 case ( 26 ):
4034 h[0] = 0;
4035 h[1] = 344.67375879228171959103668809314+(-4825.432623091944074274513633304+(19301.730492367776297098054533215+(-28952.595738551664445647081799824+14476.297869275832222823540899911*y)*y)*y)*y;
4036 h[2] = 2412.716311545972037137256816652+(-19301.730492367776297098054533215+(43428.893607827496668470622699736-28952.595738551664445647081799824*y)*y)*y+(-4825.432623091944074274513633304+(38603.460984735552594196109066430+(-86857.787215654993336941245399472+57905.191477103328891294163599647*y)*y)*y)*x;
4037 return;
4038 case ( 27 ):
4039 h[0] = 0;
4040 h[1] = 0;
4041 h[2] = 3028.6630713897510062201458646752+(-36343.956856677012074641750376103+(136289.83821253879527990656391038+(-199891.76271172356641052962706856+99945.881355861783205264813534280*y)*y)*y)*y;
4042 return;
4043 #endif
4044 #if (PMAX2D>=7)
4045 case ( 28 ):
4046 h[0] = -5855.9508194656143303910492844710+(97599.180324426905506517488074516+(-536795.49178434798028584618440984+(1288309.1802824351526860308425836+(-1395668.2786393047487432000794656+558267.31145572189949728003178623*x)*x)*x)*x)*x;
4047 h[1] = 0;
4048 h[2] = 0;
4049 return;
4050 case ( 29 ):
4051 h[0] = -5245.7983186546544929113902215894+10491.596637309308985822780443178*y+(62949.579823855853914936682659074-125899.15964771170782987336531815*y+(-236060.92433945945218101255997152+472121.84867891890436202511994304*y+(346222.68903120719653215175462490-692445.37806241439306430350924980*y+(-173111.34451560359826607587731245+346222.68903120719653215175462490*y)*x)*x)*x)*x;
4052 h[1] = -524.5798318654654492911390221589+(10491.596637309308985822780443178+(-62949.579823855853914936682659074+(157373.94955963963478734170664768+(-173111.34451560359826607587731245+69244.537806241439306430350924982*x)*x)*x)*x)*x;
4053 h[2] = 0;
4054 return;
4055 case ( 30 ):
4056 h[0] = -3114.8033645801784384587869251362+(18688.820187481070630752721550818-18688.820187481070630752721550818*y)*y+(24918.426916641427507670295401091+(-149510.56149984856504602177240654+149510.56149984856504602177240654*y)*y+(-56066.460562443211892258164652454+(336398.76337465927135354898791472-336398.76337465927135354898791472*y)*y+(37377.640374962141261505443101637+(-224265.84224977284756903265860982+224265.84224977284756903265860982*y)*y)*x)*x)*x;
4057 h[1] = -1334.9157276772193307680515393441+2669.8314553544386615361030786882*y+(-37377.640374962141261505443101636*y+18688.820187481070630752721550818+(-74755.28074992428252301088620327+149510.56149984856504602177240654*y+(-224265.84224977284756903265860981*y+112132.92112488642378451632930491+(-56066.460562443211892258164652455+112132.92112488642378451632930491*y)*x)*x)*x)*x;
4058 h[2] = -88.994381845147955384536769289608+x*(2669.8314553544386615361030786882+x*(-18688.820187481070630752721550818+x*(49836.85383328285501534059080218+x*(-56066.460562443211892258164652454+22426.584224977284756903265860982*x))));
4059 return;
4060 case ( 31 ):
4061 h[0] = -1428.7057079748789188708725069652+(17144.468495698547026450470083582+(-42861.171239246367566126175208956+28574.114159497578377417450139304*y)*y)*y+(6667.2933038827682880640716991710+(-80007.519646593219456768860390052+(200018.79911648304864192215097513-133345.86607765536576128143398342*y)*y)*y+(-6667.2933038827682880640716991710+(80007.519646593219456768860390051+(-200018.79911648304864192215097512+133345.86607765536576128143398342*y)*y)*y)*x)*x;
4062 h[1] = -1904.9409439665052251611633426203+(9524.7047198325261258058167131014-9524.7047198325261258058167131013*y)*y+(17144.468495698547026450470083582+(-85722.342478492735132252350417912+85722.342478492735132252350417912*y)*y+(-40003.759823296609728384430195026+(200018.79911648304864192215097512-200018.79911648304864192215097513*y)*y+(26669.173215531073152256286796684+(-133345.86607765536576128143398342+133345.86607765536576128143398342*y)*y)*x)*x)*x;
4063 h[2] = -476.23523599162630629029083565506+952.47047198325261258058167131016*y+x*(-19049.409439665052251611633426202*y+9524.7047198325261258058167131014+x*(-42861.171239246367566126175208956+85722.342478492735132252350417912*y+x*(-133345.86607765536576128143398342*y+66672.933038827682880640716991710+x*(-33336.466519413841440320358495854+66672.933038827682880640716991710*y))));
4064 return;
4065 case ( 32 ):
4066 h[0] = -476.23523599162630629029083565506+(9524.7047198325261258058167131014+(-42861.171239246367566126175208956+(66672.933038827682880640716991710-33336.466519413841440320358495854*y)*y)*y)*y+(952.47047198325261258058167131016+(-19049.409439665052251611633426203+(85722.342478492735132252350417912+(-133345.86607765536576128143398342+66672.933038827682880640716991710*y)*y)*y)*y)*x;
4067 h[1] = -1904.9409439665052251611633426203+(17144.468495698547026450470083582+(-40003.759823296609728384430195026+26669.173215531073152256286796684*y)*y)*y+(9524.7047198325261258058167131014+(-85722.342478492735132252350417912+(200018.79911648304864192215097513-133345.86607765536576128143398341*y)*y)*y+(-9524.7047198325261258058167131013+(85722.342478492735132252350417912+(-200018.79911648304864192215097513+133345.86607765536576128143398342*y)*y)*y)*x)*x;
4068 h[2] = -1428.7057079748789188708725069652+(6667.2933038827682880640716991710-6667.2933038827682880640716991704*y)*y+(17144.468495698547026450470083582+(-80007.519646593219456768860390052+80007.519646593219456768860390052*y)*y+(-42861.171239246367566126175208956+(200018.79911648304864192215097513-200018.79911648304864192215097512*y)*y+(28574.114159497578377417450139304+(-133345.86607765536576128143398342+133345.86607765536576128143398342*y)*y)*x)*x)*x;
4069 return;
4070 case ( 33 ):
4071 h[0] = -88.994381845147955384536769289608+(2669.8314553544386615361030786882+(-18688.820187481070630752721550818+(49836.85383328285501534059080218+(-56066.460562443211892258164652454+22426.584224977284756903265860982*y)*y)*y)*y)*y;
4072 h[1] = -1334.9157276772193307680515393441+(18688.820187481070630752721550818+(-74755.28074992428252301088620327+(112132.92112488642378451632930491-56066.460562443211892258164652455*y)*y)*y)*y+(2669.8314553544386615361030786882+(-37377.640374962141261505443101636+(149510.56149984856504602177240654+(-224265.84224977284756903265860981+112132.92112488642378451632930491*y)*y)*y)*y)*x;
4073 h[2] = -3114.8033645801784384587869251362+(24918.426916641427507670295401091+(-56066.460562443211892258164652454+37377.640374962141261505443101637*y)*y)*y+(18688.820187481070630752721550818+(-149510.56149984856504602177240654+(336398.76337465927135354898791472-224265.84224977284756903265860982*y)*y)*y+(-18688.820187481070630752721550818+(149510.56149984856504602177240654+(-336398.76337465927135354898791472+224265.84224977284756903265860982*y)*y)*y)*x)*x;
4074 return;
4075 case ( 34 ):
4076 h[0] = 0;
4077 h[1] = -524.57983186546544929113902215894+(10491.596637309308985822780443179+(-62949.579823855853914936682659074+(157373.94955963963478734170664768+(-173111.34451560359826607587731245+69244.537806241439306430350924982*y)*y)*y)*y)*y;
4078 h[2] = -5245.7983186546544929113902215894+(62949.579823855853914936682659074+(-236060.92433945945218101255997152+(346222.68903120719653215175462490-173111.34451560359826607587731244*y)*y)*y)*y+(10491.596637309308985822780443179+(-125899.15964771170782987336531815+(472121.84867891890436202511994304+(-692445.37806241439306430350924980+346222.68903120719653215175462491*y)*y)*y)*y)*x;
4079 return;
4080 case ( 35 ):
4081 h[0] = 0;
4082 h[1] = 0;
4083 h[2] = -5855.9508194656143303910492844710+(97599.180324426905506517488074516+(-536795.49178434798028584618440984+(1288309.1802824351526860308425836+(-1395668.2786393047487432000794656+558267.31145572189949728003178623*y)*y)*y)*y)*y;
4084 return;
4085 #endif
4086 #if (PMAX2D>=8)
4087 case ( 36 ):
4088 h[0] = 10390.226176556504585549952837055+(-228584.9758842431008820989624152+(1714387.3191318232566157422181140+(-5943209.3729903206229345730227951+(10400616.402733061090135502789892+(-8914814.0594854809344018595341932+2971604.6864951603114672865113979*x)*x)*x)*x)*x)*x;
4089 h[1] = 0;
4090 h[2] = 0;
4091 return;
4092 case ( 37 ):
4093 h[0] = 10142.804345939046062912011761365-20285.608691878092125824023522730*y+(-169046.73909898410104853352935609+338093.47819796820209706705871218*y+(929757.06504441255576693441145848-1859514.1300888251115338688229170*y+(-2231416.9561065901338406425875003+4462833.9122131802676812851750007*y+(2417368.3691154726449940294697920-4834736.7382309452899880589395840*y+(-966947.34764618905799761178791677+1933894.6952923781159952235758335*y)*x)*x)*x)*x)*x;
4094 h[1] = 751.3188404399293379934823526937+(-20285.608691878092125824023522730+(169046.73909898410104853352935609+(-619838.04336294170384462294097232+(1115708.4780532950669203212937502+(-966947.3476461890579976117879168+322315.78254872968599920392930558*x)*x)*x)*x)*x)*x;
4095 h[2] = 0;
4096 return;
4097 case ( 38 ):
4098 h[0] = 6772.2965085707817079879551134552+(-40633.779051424690247927730680732+40633.779051424690247927730680732*y)*y+(-81267.55810284938049585546136146+(487605.34861709628297513276816877-487605.34861709628297513276816876*y)*y+(304753.34288568517685945798010548+(-1828520.0573141110611567478806329+1828520.0573141110611567478806329*y)*y+(-446971.56956567159272720503748803+(2681829.4173940295563632302249282-2681829.4173940295563632302249282*y)*y+(223485.78478283579636360251874402+(-1340914.7086970147781816151124641+1340914.7086970147781816151124641*y)*y)*x)*x)*x)*x;
4099 h[1] = -4063.3779051424690247927730680732*y+2031.6889525712345123963865340366+(-40633.779051424690247927730680732+81267.558102849380495855461361464*y+(-487605.34861709628297513276816876*y+243802.67430854814148756638408438+(-609506.68577137035371891596021096+1219013.3715427407074378319204219*y+(-1340914.7086970147781816151124641*y+670457.35434850738909080755623205+(-268182.94173940295563632302249281+536365.88347880591127264604498566*y)*x)*x)*x)*x)*x;
4100 h[2] = 96.747092979582595828399358763646+x*(-4063.3779051424690247927730680732+x*(40633.779051424690247927730680732+x*(-162535.11620569876099171092272292+x*(304753.34288568517685945798010548+x*(-268182.94173940295563632302249282+89394.313913134318545441007497608*x)))));
4101 return;
4102 case ( 39 ):
4103 h[0] = 3685.4850427046912653706830891148+(-44225.820512456295184448197069378+(110564.55128114073796112049267345-73709.700854093825307413661782298*y)*y)*y+(-29483.880341637530122965464712918+(353806.56409965036147558557655503+(-884516.41024912590368896394138754+589677.60683275060245930929425836*y)*y)*y+(66338.730768684442776672295604070+(-796064.76922421331332006754724880+(1990161.9230605332833001688681220-1326774.6153736888555334459120814*y)*y)*y+(-44225.820512456295184448197069377+(530709.84614947554221337836483254+(-1326774.6153736888555334459120813+884516.41024912590368896394138756*y)*y)*y)*x)*x)*x;
4104 h[1] = 3158.9871794611639417462997906699+(-15794.935897305819708731498953349+15794.93589730581970873149895335*y)*y+(-44225.820512456295184448197069378+(221129.10256228147592224098534690-221129.10256228147592224098534690*y)*y+(176903.28204982518073779278827752+(-884516.41024912590368896394138754+884516.41024912590368896394138754*y)*y+(-265354.92307473777110668918241627+(1326774.6153736888555334459120814-1326774.6153736888555334459120814*y)*y+(132677.46153736888555334459120814+(-663387.30768684442776672295604068+663387.30768684442776672295604067*y)*y)*x)*x)*x)*x;
4105 h[2] = -1052.9957264870546472487665968900*y+526.49786324352732362438329844498+x*(-15794.935897305819708731498953349+31589.87179461163941746299790670*y+x*(-221129.10256228147592224098534690*y+110564.55128114073796112049267345+x*(-294838.80341637530122965464712918+589677.60683275060245930929425836*y+x*(-663387.30768684442776672295604070*y+331693.65384342221388336147802034+x*(-132677.46153736888555334459120813+265354.92307473777110668918241627*y)))));
4106 return;
4107 case ( 40 ):
4108 h[0] = 1620.+(-32400.+(145800.+(-226800.+113400.*y)*y)*y)*y+(-7560.+(151200.+(-680400.+(1058400.-529200.*y)*y)*y)*y+(7560.+(-151200.+(680400.+(-1058400.+529200.*y)*y)*y)*y)*x)*x;
4109 h[1] = 3600.+(-32400.+(75600.-50400.*y)*y)*y+(-32400.+(291600.+(-680400.+453600.*y)*y)*y+(75600.+(-680400.+(1587600.-1058400.*y)*y)*y+(-50400.+(453600.+(-1058400.+705600.*y)*y)*y)*x)*x)*x;
4110 h[2] = 1620.+(-7560.+7560.*y)*y+(-32400.+(151200.-151200.*y)*y+(145800.+(-680400.+680400.*y)*y+(-226800.+(1058400.-1058400.*y)*y+(113400.+(-529200.+529200.*y)*y)*x)*x)*x)*x;
4111 return;
4112 case ( 41 ):
4113 h[0] = 526.49786324352732362438329844498+(-15794.935897305819708731498953349+(110564.55128114073796112049267345+(-294838.80341637530122965464712918+(331693.65384342221388336147802034-132677.46153736888555334459120814*y)*y)*y)*y)*y+(-1052.9957264870546472487665968900+(31589.871794611639417462997906699+(-221129.10256228147592224098534689+(589677.60683275060245930929425838+(-663387.30768684442776672295604070+265354.92307473777110668918241627*y)*y)*y)*y)*y)*x;
4114 h[1] = 3158.9871794611639417462997906699+(-44225.820512456295184448197069378+(176903.28204982518073779278827752+(-265354.92307473777110668918241627+132677.46153736888555334459120814*y)*y)*y)*y+(-15794.935897305819708731498953349+(221129.10256228147592224098534690+(-884516.41024912590368896394138754+(1326774.6153736888555334459120814-663387.30768684442776672295604070*y)*y)*y)*y+(15794.935897305819708731498953350+(-221129.10256228147592224098534690+(884516.41024912590368896394138758+(-1326774.6153736888555334459120814+663387.30768684442776672295604066*y)*y)*y)*y)*x)*x;
4115 h[2] = 3685.4850427046912653706830891148+(-29483.880341637530122965464712918+(66338.730768684442776672295604070-44225.820512456295184448197069377*y)*y)*y+(-44225.820512456295184448197069378+(353806.56409965036147558557655504+(-796064.76922421331332006754724878+530709.84614947554221337836483253*y)*y)*y+(110564.55128114073796112049267345+(-884516.41024912590368896394138754+(1990161.9230605332833001688681220-1326774.6153736888555334459120814*y)*y)*y+(-73709.700854093825307413661782298+(589677.60683275060245930929425839+(-1326774.6153736888555334459120814+884516.41024912590368896394138756*y)*y)*y)*x)*x)*x;
4116 return;
4117 case ( 42 ):
4118 h[0] = 96.747092979582595828399358763646+(-4063.3779051424690247927730680732+(40633.779051424690247927730680732+(-162535.11620569876099171092272292+(304753.34288568517685945798010548+(-268182.94173940295563632302249282+89394.313913134318545441007497608*y)*y)*y)*y)*y)*y;
4119 h[1] = 2031.6889525712345123963865340366+(-40633.779051424690247927730680732+(243802.67430854814148756638408438+(-609506.68577137035371891596021096+(670457.35434850738909080755623205-268182.94173940295563632302249281*y)*y)*y)*y)*y+(-4063.3779051424690247927730680732+(81267.558102849380495855461361464+(-487605.34861709628297513276816876+(1219013.3715427407074378319204219+(-1340914.7086970147781816151124641+536365.88347880591127264604498563*y)*y)*y)*y)*y)*x;
4120 h[2] = 6772.2965085707817079879551134552+(-81267.55810284938049585546136146+(304753.34288568517685945798010548+(-446971.56956567159272720503748803+223485.78478283579636360251874403*y)*y)*y)*y+(-40633.779051424690247927730680732+(487605.34861709628297513276816876+(-1828520.0573141110611567478806329+(2681829.4173940295563632302249282-1340914.7086970147781816151124641*y)*y)*y)*y+(40633.779051424690247927730680732+(-487605.34861709628297513276816876+(1828520.0573141110611567478806329+(-2681829.4173940295563632302249282+1340914.7086970147781816151124641*y)*y)*y)*y)*x)*x;
4121 return;
4122 case ( 43 ):
4123 h[0] = 0;
4124 h[1] = 751.31884043992933799348235269371+(-20285.608691878092125824023522730+(169046.73909898410104853352935609+(-619838.04336294170384462294097232+(1115708.4780532950669203212937502+(-966947.3476461890579976117879168+322315.78254872968599920392930558*y)*y)*y)*y)*y)*y;
4125 h[2] = 10142.804345939046062912011761365+(-169046.73909898410104853352935609+(929757.06504441255576693441145848+(-2231416.9561065901338406425875003+(2417368.3691154726449940294697920-966947.34764618905799761178791677*y)*y)*y)*y)*y+(-20285.608691878092125824023522730+(338093.47819796820209706705871218+(-1859514.1300888251115338688229169+(4462833.9122131802676812851750007+(-4834736.7382309452899880589395840+1933894.6952923781159952235758336*y)*y)*y)*y)*y)*x;
4126 return;
4127 case ( 44 ):
4128 h[0] = 0;
4129 h[1] = 0;
4130 h[2] = 10390.226176556504585549952837055+(-228584.9758842431008820989624152+(1714387.3191318232566157422181140+(-5943209.3729903206229345730227951+(10400616.402733061090135502789892+(-8914814.0594854809344018595341932+2971604.6864951603114672865113979*y)*y)*y)*y)*y)*y;
4131 return;
4132 #endif
4133 }
4134 printf( "Error in hess_quadrilateral_2d -- unmatched switch/case %d\n", i );
4135 }
4136
4137 }; // end class OrthonormalBase_2D
4138 } // end namespace Fem
4139
4140} // end namespace Dune
4141#endif
Dune
Dune namespace.
Definition: alignedallocator.hh:13
Creative Commons License   |  Legal Statements / Impressum  |  Hosted by TU Dresden  |  generated with Hugo v0.111.3 (Nov 13, 23:29, 2024)