WdgG4Basis.h

Go to the documentation of this file.

#define size_tri ((order + 1) * (order + 2) / 2)
#define size_lin (order + 1)
 
WORKSPACE(RealT, term_tri, size_tri);
WORKSPACE(RealT, term_lin, size_lin);
 
term_tri[0] = 94.87247638777873248727701138705 * ss * rr - 16.180417124823243568698671879247 * rr -
  16.18041712482329330669017508626 * ss - 164.15486234178916902237688191235 * ss * rr * rr -
  164.15486234178956692630890756845 * ss * ss * rr +
  87.321823800852058639065944589674 * ss * rr * rr * rr +
  87.321823800852484964707400649786 * ss * ss * ss * rr +
  44.201433912052038976980838924646 * ss * ss - 49.062462841081497799677890725434 * ss * ss * ss +
  19.135649775259295779505919199437 * ss * ss * ss * ss +
  44.201433912051953711852547712624 * rr * rr - 49.0624628410813699019854539074 * rr * rr * rr +
  19.135649775259253146941773593426 * rr * rr * rr * rr +
  136.34281406991490825930668506771 * ss * ss * rr * rr + 1.9440916731691306740259506113944;
term_tri[1] = 0.28135944494482556521575133956503 * ss - 1.5776612770734916857406915369211 * rr -
  4.8723350260304272296707495115697 * ss * rr + 15.370200280900945699613657779992 * ss * rr * rr +
  4.6500435436287119728149264119565 * ss * ss * rr -
  10.779224699815372900957299862057 * ss * rr * rr * rr -
  0.059067962543096985683632738073356 * ss * ss * ss * rr -
  0.31089342621638227370794993476011 * ss * ss + 0.05906796254311164062755779013969 * ss * ss * ss -
  0.029533981271555598269173970038537 * ss * ss * ss * ss +
  11.827944040363041366958896105643 * rr * rr - 27.480136259955177280289717600681 * rr * rr * rr +
  19.13564977525900090427057875786 * rr * rr * rr * rr -
  10.808758681086919395397671905812 * ss * ss * rr * rr + 0.038295394575720613505875888904484;
term_tri[2] = 0.28135944494481890387760358862579 * rr - 1.5776612770735174429148628405528 * ss -
  4.8723350260302566994141670875251 * ss * rr + 4.6500435436285414425583439879119 * ss * rr * rr +
  15.370200280900235156877897679806 * ss * ss * rr -
  0.059067962543061014457634883001447 * ss * rr * rr * rr -
  10.779224699814754728777188574895 * ss * ss * ss * rr +
  11.82794404036310353944827511441 * ss * ss - 27.480136259955202149285469204187 * ss * ss * ss +
  19.13564977525900090427057875786 * ss * ss * ss * ss -
  0.31089342621634585839274222962558 * rr * rr + 0.05906796254308444016345447380445 * rr * rr * rr -
  0.029533981271544051949717868410517 * rr * rr * rr * rr -
  10.808758681086317210429115220904 * ss * ss * rr * rr + 0.038295394575721113106236970224927;
term_tri[3] = 2.2634203558048850979389499116223 * ss + 26.432035442978936146118940087035 * rr -
  113.79487188356213778206438291818 * ss * rr + 294.91435990471728700867970474064 * ss * rr * rr +
  150.06440227676208110096922609955 * ss * ss * rr -
  190.14225760864377434700145386159 * ss * rr * rr * rr -
  62.05472947454963872360167442821 * ss * ss * ss * rr +
  1.0468615827391740680241127847694 * ss * ss - 6.6134644939322981826990144327283 * ss * ss * ss +
  4.1103664064101437247700232546777 * ss * ss * ss * ss -
  100.6694346217531261800104402937 * rr * rr + 130.5911837615829540482081938535 * rr * rr * rr -
  55.723424011090003205026732757688 * rr * rr * rr * rr -
  200.5946495724491285272961249575 * ss * ss * rr * rr - 0.7801482275236110153215918217029;
term_tri[4] = 0.55087578907791667770510457557975 * ss - 13.270388538115096110914237215184 * rr +
  15.832433567271408492160844616592 * ss * rr - 143.9873799440699144724931102246 * ss * rr * rr +
  20.890350719475399898783507524058 * ss * ss * rr +
  141.4253349149136056439601816237 * ss * rr * rr * rr -
  23.452395748631651883897575316951 * ss * ss * ss * rr -
  3.907930324714739889202519407263 * ss * ss + 4.3356570981709197099007724318653 * ss * ss * ss -
  1.3044026988552248624841922719497 * ss * ss * ss * ss +
  83.983055995572300389540032483637 * rr * rr - 141.42533491491462882549967616796 * rr * rr * rr +
  70.712667457457442310442274902016 * rr * rr * rr * rr +
  47.260271708825015934962721075863 * ss * ss * rr * rr + 0.34447455870034116021116687988979;
term_tri[5] = 6.0269785601381933304310223320499 * rr - 0.73237067154559598591845315240789 * ss +
  7.9200700929723435450569013482891 * ss * rr + 1.3951239271137367836672638077289 * ss * rr * rr -
  22.992392125964467908261212869547 * ss * ss * rr -
  32.751438435714490537975507322699 * ss * rr * rr * rr +
  8.3985071112246174607207649387419 * ss * ss * ss * rr +
  1.6731117835679034655527175345924 * ss * ss - 0.75319773878421580892705833321088 * ss * ss * ss -
  0.01072009393497452833798888605088 * ss * ss * ss * ss -
  43.236427403543473246827488765121 * rr * rr + 92.302512282776262964034685865045 * rr * rr * rr -
  55.723424011089711882505071116611 * rr * rr * rr * rr +
  35.491579186944818502524867653847 * ss * ss * rr * rr - 0.1497876558048262529698746448048;
term_tri[6] = 11.567588832493493455899624677841 * ss * rr - 0.9582156651269206903975828026887 * rr -
  0.31137930349765952087182085961103 * ss - 65.583789542012084439193131402135 * ss * rr * rr +
  4.5913635512818338213492097565904 * ss * ss * rr +
  78.496195100190291782382701057941 * ss * rr * rr * rr -
  8.441387486964590181059975293465 * ss * ss * ss * rr -
  0.65080199639461766558667932258686 * ss * ss + 0.79607811452415833119999888367602 * ss * ss * ss -
  0.010720093934997843021506014338229 * ss * ss * ss * ss +
  5.8686665394031791720408364199102 * rr * rr - 9.8280011317083335597999393939972 * rr * rr * rr +
  4.1103664064101721464794536586851 * rr * rr * rr * rr +
  10.231737289660882694874999288004 * ss * ss * rr * rr + 0.027035623498290362043405821168562;
term_tri[7] = 15.289301793006224272630788618699 * ss * rr -
  0.52437563874030790334046514544752 * rr - 0.52437563874031511979012520896504 * ss -
  46.820975434669804826626204885542 * ss * rr * rr -
  46.820975434669563242096046451479 * ss * ss * rr +
  18.234784953211025992914073867723 * ss * rr * rr * rr +
  18.234784953210777302956557832658 * ss * ss * ss * rr +
  1.2726247766667313499056035652757 * ss * ss + 0.88195369724987138226879324065521 * ss * ss * ss -
  1.3044026988551640222624428133713 * ss * ss * ss * ss +
  1.2726247766667293515041592399939 * rr * rr + 0.8819536972498970284206620817713 * rr * rr * rr -
  1.3044026988551844503660959162517 * rr * rr * rr * rr +
  109.79104276158858510825666598976 * ss * ss * rr * rr + 0.018674422379214489220444761485851;
term_tri[8] = 11.567588832493424177982888068072 * ss * rr -
  0.31137930349764952886459923320217 * rr - 0.95821566512691980221916310256347 * ss +
  4.5913635512819652717553253751248 * ss * rr * rr -
  65.583789542011899698081833776087 * ss * ss * rr -
  8.4413874869647003151840181089938 * ss * rr * rr * rr +
  78.496195100190163884690264239907 * ss * ss * ss * rr + 5.8686665394031400921903696144 * ss * ss -
  9.8280011317082731636673997854814 * ss * ss * ss +
  4.1103664064101455011268626549281 * ss * ss * ss * ss -
  0.65080199639465485805800426533096 * rr * rr + 0.79607811452421461950734737911262 * rr * rr * rr -
  0.010720093935022267928047767782118 * rr * rr * rr * rr +
  10.23173728966071394097525626421 * ss * ss * rr * rr + 0.027035623498290139998800896137254;
term_tri[9] = 6.0269785601382670492398574424442 * ss - 0.73237067154558022075150347518502 * rr +
  7.9200700929718284015734752756543 * ss * rr - 22.99239212596405934618815081194 * ss * rr * rr +
  1.3951239271160726929110751370899 * ss * ss * rr +
  8.398507111224613908007086138241 * ss * rr * rr * rr -
  32.751438435716295316524337977171 * ss * ss * ss * rr -
  43.236427403543650882511428790167 * ss * ss + 92.302512282776305596598831471056 * ss * ss * ss -
  55.723424011089662144513567909598 * ss * ss * ss * ss +
  1.6731117835678310790115119743859 * rr * rr - 0.75319773878418505574927621637471 * rr * rr * rr -
  0.010720093934989183281913938117214 * rr * rr * rr * rr +
  35.491579186943113199959043413401 * ss * ss * rr * rr - 0.14978765580482777952653350439505;
term_tri[10] = 0.55087578907791312499142577507882 * rr - 13.270388538115206245038280030712 * ss +
  15.832433567272017782556758902501 * ss * rr + 20.890350719474827911881220643409 * ss * rr * rr -
  143.98737994407287033027387224138 * ss * ss * rr -
  23.452395748631580829623999306932 * ss * rr * rr * rr +
  141.42533491491587938071461394429 * ss * ss * ss * rr +
  83.983055995572527763215475715697 * ss * ss - 141.42533491491474251233739778399 * ss * ss * ss +
  70.712667457457371256168698891997 * ss * ss * ss * ss -
  3.9079303247147025857088920020033 * rr * rr + 4.3356570981709809942117317405064 * rr * rr * rr -
  1.3044026988552579471303261016146 * rr * rr * rr * rr +
  47.260271708827190195734146982431 * ss * ss * rr * rr + 0.34447455870034321412376243642939;
term_tri[11] = 26.432035442979064043811376905069 * ss + 2.2634203558048815452252711111214 * rr -
  113.79487188356277727052656700835 * ss * rr + 150.06440227676344534302188549191 * ss * rr * rr +
  294.91435990471950390201527625322 * ss * ss * rr -
  62.054729474550427426038368139416 * ss * rr * rr * rr -
  190.14225760864539438443898689002 * ss * ss * ss * rr -
  100.66943462175336776454059872776 * ss * ss + 130.59118376158318142188363708556 * ss * ss * ss -
  55.723424011090031626736163161695 * ss * ss * ss * ss +
  1.0468615827391294370585228534765 * rr * rr - 6.6134644939324047641093784477562 * rr * rr * rr +
  4.1103664064102245490062159660738 * rr * rr * rr * rr -
  200.59464957245134542063169647008 * ss * ss * rr * rr - 0.78014822752361312474533860950032;
term_tri[12] = 180.08166136297191428639052901417 * ss * rr -
  7.1063986771366964845242364390288 * rr - 7.1063986771367328998394441441633 * ss -
  374.39092011020147765520960092545 * ss * rr * rr -
  374.39092011019999972631921991706 * ss * ss * rr +
  200.0489288196101540506788296625 * ss * rr * rr * rr +
  200.04892881960878980862617027014 * ss * ss * ss * rr +
  13.607190891937248622411971155088 * ss * ss - 4.1638228815068343280358931224328 * ss * ss * ss -
  2.7132868540380958854996151785599 * ss * ss * ss * ss +
  13.607190891937246846055131754838 * rr * rr - 4.1638228815066895549534820020199 * rr * rr * rr -
  2.713286854038212680961805745028 * rr * rr * rr * rr +
  407.8529328236072615254670381546 * ss * ss * rr * rr + 0.27670985301620926177790238398302;
term_tri[13] = 1.8699043491790923177120475884294 * ss + 3.2366329539351892385923292749794 * rr -
  61.774605355023041397544147912413 * ss * rr + 270.80677724160619845861219801009 * ss * rr * rr +
  55.24769805321965776556680793874 * ss * ss * rr -
  210.90207623576242212948272936046 * ss * rr * rr * rr +
  4.6570029526241807360520397196524 * ss * ss * ss * rr +
  0.45859712713296718611388769204495 * ss * ss - 4.6570029526241043527079455088824 * ss * ss * ss +
  2.3285014763120508440863432042534 * ss * ss * ss * ss -
  15.163998876812193827845476334915 * rr * rr + 15.016970297659707256343608605675 * rr * rr * rr -
  2.7132868540383254796211076609325 * rr * rr * rr * rr -
  208.57357475944971270109817851335 * ss * ss * rr * rr - 0.099607667728191781608693133875931;
term_tri[14] = 3.2366329539351887945031194249168 * ss + 1.8699043491790510174155315326061 * rr -
  61.774605355022742969595128670335 * ss * rr + 55.247698053219608027575304731727 * ss * rr * rr +
  270.80677724160500474681612104177 * ss * ss * rr +
  4.657002952624026193007011897862 * ss * rr * rr * rr -
  210.90207623576145579136209562421 * ss * ss * ss * rr -
  15.163998876812087246435112319887 * ss * ss + 15.016970297659621991215317393653 * ss * ss * ss -
  2.7132868540383103805879727588035 * ss * ss * ss * ss +
  0.45859712713307665410411573247984 * rr * rr - 4.657002952624301528317118936684 * rr * rr * rr +
  2.3285014763121538727830284187803 * rr * rr * rr * rr -
  208.57357475944908742349070962518 * ss * ss * rr * rr - 0.099607667728191073841514935338637;
 
term_lin[0] = 0.42148179237593254831600381883471 * tt * (tt - 0.53846931010568310771446931539685) *
  (tt - 0.90617984593866396370032134655048) *
  (2.7195304527641102139675158885874 * tt + 1.4643836867112864494512896092315);
term_lin[1] = -1.193673917625192347641621838069 * tt * (tt - 0.53846931010568310771446931539685) *
  (tt - 0.90617984593866396370032134655048) *
  (2.7195304527641102139675158885874 * tt + 2.4643836867112864494512896092315);
term_lin[2] = 3.8059553529423886561929567469294 * (tt - 0.53846931010568310771446931539685) *
  (tt + 0.53846931010568310771446931539685) * (tt - 0.90617984593866396370032134655048) *
  (1.1035337019266332074210473420733 * tt + 1.0);
term_lin[3] = -4.6896671420713595859789534540298 * tt * (tt + 0.53846931010568310771446931539685) *
  (tt - 0.90617984593866396370032134655048) *
  (0.69220959000048211835749389066225 * tt + 0.62726637962390263339969233772114);
term_lin[4] = 2.0773857067541632774276185451694 * tt * (tt - 0.53846931010568310771446931539685) *
  (tt + 0.53846931010568310771446931539685) * (0.55176685096331660371052367103666 * tt + 0.5);
 
for (int ii = 0; ii < size_lin; ++ii)
{
  for (int jj = 0; jj < size_tri; ++jj)
  {
    int i_sp = ii * size_tri + jj;
    basis[i_sp] = term_tri[jj] * term_lin[ii];
  }
}