diff options
Diffstat (limited to 'Dragon/src/NXTQRN.f')
| -rw-r--r-- | Dragon/src/NXTQRN.f | 976 |
1 files changed, 976 insertions, 0 deletions
diff --git a/Dragon/src/NXTQRN.f b/Dragon/src/NXTQRN.f new file mode 100644 index 0000000..f25b00d --- /dev/null +++ b/Dragon/src/NXTQRN.f @@ -0,0 +1,976 @@ +*DECK NXTQRN + SUBROUTINE NXTQRN(NDIM ,ORDRE ,NQUAD ,NBANGL,DQUAD , + > DANGLT,DDENWT) +* +*----------------------------------------------------------------------- +* +*Purpose: +* To define quadruple range (QR) quadrature angles. +* +*Copyright: +* Copyright (C) 2008 Ecole Polytechnique de Montreal +* This library is free software; you can redistribute it and/or +* modify it under the terms of the GNU Lesser General Public +* License as published by the Free Software Foundation; either +* version 2.1 of the License, or (at your option) any later version. +* +*Author(s): +* N. Martin +* +*Parameters: input +* NDIM number of dimensions for geometry. +* ORDRE quadrature order. +* NQUAD number of quadrant (in 3-D) and quarter (in 2-D). +* NBANGL number of angles. +* DQUAD relative density of each quadrant. +* +*Parameters: output +* DANGLT director cosines of angles. +* DDENWT angular density for each angle. +* +*References: +* I. K. Abu-Shumays, Angular quadratures for improved transport +* computations, Transport Theory and Statistical Physics, 30:2, +* 169 - 204 (2001). +* \\\\ +* E. M. Baker, Quadruple range quadrature verification and extention, +* Los Alamos documentation LA-UR-07-8050 (2006). +* +*---------- +* + IMPLICIT NONE +*---- +* Subroutine arguments +*---- + INTEGER NDIM,ORDRE,NQUAD,NBANGL + DOUBLE PRECISION DQUAD(NQUAD),DANGLT(NDIM,NQUAD,NBANGL), + > DDENWT(NQUAD,NBANGL) +*---- +* Local parameters +*---- + INTEGER IOUT + CHARACTER NAMSBR*6 + PARAMETER (IOUT=6,NAMSBR='NXTQRN') + DOUBLE PRECISION DZERO,DONE,DTWO + PARAMETER (DZERO=0.0D0,DONE=1.0D0,DTWO=2.0D0) +*---- +* Functions +*---- + DOUBLE PRECISION XDRCST,PI +*---- +* Local variables +*---- + INTEGER NBPT,ILEVEL,JLEVEL,NLEVEL,KLEVEL,IPOINT, + > IPT,IROT,IDIR1,IDIR2,IDIR3 + DOUBLE PRECISION ALPHA,XI,ROTXI,COSA,SINA,WGTINV +*---- +* Data +*---- + INTEGER IS,IA + DOUBLE PRECISION SYST(2,703) + SAVE SYST +*---- +* Order = 2 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA= 1, 1)/ + > 0.7071067811865475D0,1.570796326794897D0/ +*---- +* Order = 4 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA= 2, 3)/ + > 3.288613193063390D-1,7.853981633974483D-1,9.443782254288237D-1, + > 0.7853981633974483D0/ +*---- +* Order = 6 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA= 4, 6)/ + > 1.797057505503707D-1, 4.414924088048230D-1,7.071067811865475D-1, + > 6.878115091852507D-1, 9.837204090691257D-1, 4.414924088048230D-1/ +*---- +* Order = 8 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA= 7, 10)/ + > 1.119194178021025D-1, 2.790495092691071D-1, + > 4.990088153604911D-1, 5.063486541283412D-1, + > 8.665969087139183D-1, 5.063486541283412D-1, + > 9.937172857099943D-1, 2.790495092691071D-1/ +*---- +* Order = 10 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA= 11, 15)/ + > 7.608393410010900D-2, 1.913207247027131D-1, + > 3.599152552241242D-1, 3.748555766138863D-1, + > 7.071067811865475D-1, 4.384437241616979D-1, + > 9.329849993740272D-1, 3.748555766138863D-1, + > 9.971014165930416D-1, 1.913207247027131D-1/ +*---- +* Order = 12 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA= 16, 21)/ + > 5.498238339818787D-2, 1.389912271015858D-1, + > 2.686450605979331D-1, 2.849847975948065D-1, + > 5.658713393447228D-1, 3.614221387010560D-1, + > 8.244935580756284D-1, 3.614221387010560D-1, + > 9.632392389309796D-1, 2.849847975948065D-1, + > 9.984873246645922D-1, 1.389912271015858D-1/ +*---- +* Order = 14 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA= 22, 28)/ + > 4.155046739519842D-2, 1.054032627117009D-1, + > 2.070312681580549D-1, 2.225734192425536D-1, + > 4.548121603136229D-1, 2.966594661924317D-1, + > 7.071067811865475D-1, 3.215240305015243D-1, + > 8.905873897775868D-1, 2.966594661924317D-1, + > 9.783343262938634D-1, 2.225734192425536D-1, + > 9.991364064326955D-1, 1.054032627117009D-1/ +*---- +* Order = 16 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA= 29, 36)/ + > 3.248697884632087D-2, 8.261058789185801D-2, + > 1.639512569400686D-1, 1.780108196819361D-1, + > 3.700379385471330D-1, 2.453650131486301D-1, + > 6.006186663088213D-1, 2.794117426750241D-1, + > 7.995356262740346D-1, 2.794117426750241D-1, + > 9.290166435730784D-1, 2.453650131486301D-1, + > 9.864684411311756D-1, 1.780108196819361D-1, + > 9.994721587945503D-1, 8.261058789185801D-2/ +*---- +* Order = 18 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA= 37, 45)/ + > 2.608861776207381D-2, 6.645559945847256D-2, + > 1.328246569072209D-1, 1.452965934855972D-1, + > 3.053252347988094D-1, 2.051444864132550D-1, + > 5.104236819705835D-1, 2.416159520924252D-1, + > 7.071067811865475D-1, 2.537710638953967D-1, + > 8.599230575368895D-1, 2.416159520924252D-1, + > 9.522481299509345D-1, 2.051444864132550D-1, + > 9.911395514847942D-1, 1.452965934855972D-1, + > 9.996596340871549D-1, 6.645559945847256D-2/ +*---- +* Order = 20 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA= 46, 55)/ + > 2.140636649836879D-2, 5.459883427719590D-2, + > 1.096782080789725D-1, 1.206689353273136D-1, + > 2.554062770706209D-1, 1.734506495359569D-1, + > 4.359965790249476D-1, 2.093517367248734D-1, + > 6.217587837057514D-1, 2.273280075321085D-1, + > 7.832087939271013D-1, 2.273280075321085D-1, + > 8.999483224488741D-1, 2.093517367248734D-1, + > 9.668338190366146D-1, 1.734506495359569D-1, + > 9.939671476827519D-1, 1.206689353273136D-1, + > 9.997708574835223D-1, 5.459883427719590D-2/ +*---- +* Order = 22 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA= 56, 66)/ + > 1.787818127983761D-2, 4.564471367918457D-2, + > 9.203377352741595D-2, 1.017137272660963D-1, + > 2.163635206471488D-1, 1.482265633964675D-1, + > 3.750373673684837D-1, 1.822612531818470D-1, + > 5.464247324681954D-1, 2.027611569667827D-1, + > 7.071067811865475D-1, 2.095814978141405D-1, + > 8.375082159280953D-1, 2.027611569667827D-1, + > 9.270096941657714D-1, 1.822612531818470D-1, + > 9.763128734853243D-1, 1.482265633964675D-1, + > 9.957558860133865D-1, 1.017137272660963D-1, + > 9.998401725446549D-1, 4.564471367918457D-2/ +*---- +* Order = 24 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA= 67, 78)/ + > 1.515427633810212D-2, 3.871995573990722D-2, + > 7.829331784823644D-2, 8.683898824818842D-2, + > 1.853856824160965D-1, 1.279220161079748D-1, + > 3.250513900479526D-1, 1.595962888003237D-1, + > 4.813341017203090D-1, 1.808538704864877D-1, + > 6.359378348511013D-1, 1.914670440145664D-1, + > 7.717402867577236D-1, 1.914670440145664D-1, + > 8.765372111445716D-1, 1.808538704864877D-1, + > 9.456963539254520D-1, 1.595962888003237D-1, + > 9.826658377877590D-1, 1.279220161079748D-1, + > 9.969303668663700D-1, 8.683898824818842D-2, + > 9.998851673610667D-1, 3.871995573990722D-2/ +*---- +* Order = 26 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA= 79, 91)/ + > 1.300784260340052D-2, 3.325591619718422D-2, + > 6.739381484979049D-2, 7.496658920041428D-2, + > 1.604653826023372D-1, 1.113898480184530D-1, + > 2.838506265695340D-1, 1.405952903244943D-1, + > 4.256053496742027D-1, 1.616563247494470D-1, + > 5.714795978613439D-1, 1.742889475411106D-1, + > 7.071067811865475D-1, 1.784904947326901D-1, + > 8.206162740454497D-1, 1.742889475411106D-1, + > 9.049088828874981D-1, 1.616563247494470D-1, + > 9.588685112131292D-1, 1.405952903244943D-1, + > 9.870414687267630D-1, 1.113898480184530D-1, + > 9.977264523505389D-1, 7.496658920041428D-2, + > 9.999153944363519D-1, 3.325591619718422D-2/ +*---- +* Order = 28 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA= 92,105)/ + > 1.128664599103795D-2, 2.886965398857025D-2, + > 5.860802825933166D-2, 6.534816595720066D-2, + > 1.401592977980301D-1, 9.778151915986106D-2, + > 2.496535423349043D-1, 1.245928602207950D-1, + > 3.780145294131805D-1, 1.449498625863828D-1, + > 5.141044083621162D-1, 1.585386998315389D-1, + > 6.460952816893911D-1, 1.653174016530996D-1, + > 7.632567634673841D-1, 1.653174016530996D-1, + > 8.577276125336286D-1, 1.585386998315389D-1, + > 9.257996627524400D-1, 1.449498625863828D-1, + > 9.683352254253865D-1, 1.245928602207950D-1, + > 9.901289669738802D-1, 9.778151915986106D-2, + > 9.982810721553091D-1, 6.534816595720066D-2, + > 9.999363037825324D-1, 2.886965398857025D-2/ +*---- +* Order = 30 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA=106,120)/ + > 9.885454799950707D-3, 2.529575708608393D-2, + > 5.142584271456661D-2, 5.745247946520937D-2, + > 1.234174526540014D-1, 8.646511772502708D-2, + > 2.210520064292830D-1, 1.110388326617575D-1, + > 3.373313406305213D-1, 1.304366407688838D-1, + > 4.634907863833664D-1, 1.443355921311845D-1, + > 5.898656904335448D-1, 1.526595716146173D-1, + > 7.071067811865475D-1, 1.554283438893696D-1, + > 8.075013729086518D-1, 1.526595716146173D-1, + > 8.861017384802541D-1, 1.443355921311845D-1, + > 9.413859817462841D-1, 1.304366407688838D-1, + > 9.752620214350542D-1, 1.110388326617575D-1, + > 9.923548419695433D-1, 8.646511772502708D-2, + > 9.986768159425233D-1, 5.745247946520937D-2, + > 9.999511376979368D-1, 2.529575708608393D-2/ +*---- +* Order = 32 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA=121,136)/ + > 8.729657082861998D-3, 2.234563433970334D-2, + > 4.548138640446708D-2, 5.089471201185089D-2, + > 1.094657179495637D-1, 7.696546183192337D-2, + > 1.969448911359501D-1, 9.948832189094168D-2, + > 3.024482742232713D-1, 1.178163725528086D-1, + > 4.190211144820449D-1, 1.316318495967909D-1, + > 5.387770134090762D-1, 1.408332013024988D-1, + > 6.537245138927594D-1, 1.454226098709307D-1, + > 7.567326211652801D-1, 1.454226098709307D-1, + > 8.424484137453142D-1, 1.408332013024988D-1, + > 9.079764895735049D-1, 1.316318495967909D-1, + > 9.531657995434818D-1, 1.178163725528086D-1, + > 9.804145602016775D-1, 9.948832189094168D-2, + > 9.939905716825420D-1, 7.696546183192337D-2, + > 9.989651863256935D-1, 5.089471201185089D-2, + > 9.999618958176435D-1, 2.234563433970334D-2/ +*---- +* Order = 34 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA=137,153)/ + > 7.765167550594897D-3, 1.988233444736077D-2, + > 4.050698458787639D-2, 4.539097865607717D-2, + > 9.772547366646427D-2, 6.892134946532859D-2, + > 1.764725073896891D-1, 8.958447192452450D-2, + > 2.724164548901113D-1, 1.068152712507604D-1, + > 3.799892988394926D-1, 1.203101631042851D-1, + > 4.927420326833888D-1, 1.299509173542970D-1, + > 6.039089138656826D-1, 1.357218046104536D-1, + > 7.071067811865475D-1, 1.376417451687221D-1, + > 7.970533380857090D-1, 1.357218046104536D-1, + > 8.701754358903969D-1, 1.299509173542970D-1, + > 9.249908825320770D-1, 1.203101631042851D-1, + > 9.621794401800030D-1, 1.068152712507604D-1, + > 9.843055694933338D-1, 8.958447192452450D-2, + > 9.952134101772671D-1, 6.892134946532859D-2, + > 9.991792552888583D-1, 4.539097865607717D-2, + > 9.999698506319634D-1, 1.988233444736077D-2/ +*---- +* Order = 36 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA=154,171)/ + > 6.952010366745025D-3, 1.780448256851573D-2, + > 3.630329822321955D-2, 4.072837739441321D-2, + > 8.775842996404620D-2, 6.205518747222860D-2, + > 1.589616059026264D-1, 8.104138647658412D-2, + > 2.464422434563339D-1, 9.719456818914889D-2, + > 3.457022371739597D-1, 1.102303994594047D-1, + > 4.514464196113868D-1, 1.200213906455903D-1, + > 5.579698061333865D-1, 1.265362696526889D-1, + > 6.596626820591664D-1, 1.297861015388739D-1, + > 7.515618044435914D-1, 1.297861015388739D-1, + > 8.298612507181373D-1, 1.265362696526889D-1, + > 8.922982294166338D-1, 1.200213906455903D-1, + > 9.383442668940432D-1, 1.102303994594047D-1, + > 9.691574797937686D-1, 9.719456818914889D-2, + > 9.872847653280476D-1, 8.104138647658412D-2, + > 9.961417860777880D-1, 6.205518747222860D-2, + > 9.993408180086091D-1, 4.072837739441321D-2, + > 9.999758344839443D-1, 1.780448256851573D-2/ +*---- +* Order = 38 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA=172,190)/ + > 6.260124726824094D-3, 1.603575186530219D-2, + > 3.271949450064389D-2, 3.674484447429247D-2, + > 7.922836042016190D-2, 5.615133571572357D-2, + > 1.438820502560538D-1, 7.362948392130560D-2, + > 2.238694254490736D-1, 8.875049346501127D-2, + > 3.155245603429847D-1, 1.012516526819094D-1, + > 4.144771888783680D-1, 1.110019513021974D-1, + > 5.159200804329141D-1, 1.179578982613787D-1, + > 6.149740744496668D-1, 1.221213481438653D-1, + > 7.071067811865475D-1, 1.235068071329247D-1, + > 7.885473275300447D-1, 1.221213481438653D-1, + > 8.566367203232041D-1, 1.179578982613787D-1, + > 9.100597012831003D-1, 1.110019513021974D-1, + > 9.489174104316804D-1, 1.012516526819094D-1, + > 9.746191463074700D-1, 8.875049346501127D-2, + > 9.895948441731669D-1, 7.362948392130560D-2, + > 9.968564926332842D-1, 5.615133571572357D-2, + > 9.994645739993101D-1, 3.674484447429247D-2, + > 9.999804052272247D-1, 1.603575186530219D-2/ +*---- +* Order = 40 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA=191,210)/ + > 5.666555435801935D-3, 1.451779509317350D-2, + > 2.963979274695100D-2, 3.331533618587398D-2, + > 7.187424763201411D-2, 5.104054858395153D-2, + > 1.308135974586322D-1, 6.716362589103700D-2, + > 2.041579376098809D-1, 8.131097620559026D-2, + > 2.888941786344192D-1, 9.324211470049753D-2, + > 3.813947484262091D-1, 1.028266819400935D-1, + > 4.776033421571356D-1, 1.100123318316945D-1, + > 5.732994237606388D-1, 1.147920152181985D-1, + > 6.644147983089062D-1, 1.171767377473380D-1, + > 7.473640182589308D-1, 1.171767377473380D-1, + > 8.193459408062747D-1, 1.147920152181985D-1, + > 8.785755787411428D-1, 1.100123318316945D-1, + > 9.244122705118689D-1, 1.028266819400935D-1, + > 9.573610361567596D-1, 9.324211470049753D-2, + > 9.789379635660678D-1, 8.131097620559026D-2, + > 9.914069813754243D-1, 6.716362589103700D-2, + > 9.974137017944620D-1, 5.104054858395153D-2, + > 9.995606448264746D-1, 3.331533618587398D-2, + > 9.999839449458641D-1, 1.451779509317350D-2/ +*---- +* Order = 42 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA=211,231)/ + > 5.153529715794118D-3, 1.320539632418366D-2, + > 2.697415577591218D-2, 3.034217800698840D-2, + > 6.549119792630967D-2, 4.658869227256107D-2, + > 1.194207057461457D-1, 6.149373560761037D-2, + > 1.868636355236370D-1, 7.473138541437220D-2, + > 2.653234848709176D-1, 8.608289338331400D-2, + > 3.517712327776513D-1, 9.542142388267645D-2, + > 4.427751437470710D-1, 1.026887640907876D-1, + > 5.346931294496370D-1, 1.078700892228989D-1, + > 6.239117782235315D-1, 1.109716362868715D-1, + > 7.071067811865475D-1, 1.120039378103684D-1, + > 7.814947811687234D-1, 1.109716362868715D-1, + > 8.450463048374058D-1, 1.078700892228989D-1, + > 8.966326851502574D-1, 1.026887640907876D-1, + > 9.360860002104996D-1, 9.542142388267645D-2, + > 9.641594517381199D-1, 8.608289338331400D-2, + > 9.823858619294607D-1, 7.473138541437220D-2, + > 9.928437415016990D-1, 6.149373560761037D-2, + > 9.978531470082043D-1, 4.658869227256107D-2, + > 9.996361312598584D-1, 3.034217800698840D-2, + > 9.999867204775614D-1, 1.320539632418366D-2/ +*---- +* Order = 44 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA=232,253)/ + > 4.707114717103505D-3, 1.206308084642689D-2, + > 2.465173376735076D-2, 2.774821193462954D-2, + > 5.991670409568439D-2, 4.268848001641819D-2, + > 1.094334822143285D-1, 5.649746567994538D-2, + > 1.716206775350476D-1, 6.889030401909524D-2, + > 2.443940481368066D-1, 7.966885558279603D-2, + > 3.252092377487526D-1, 8.871167329713092D-2, + > 4.111529714803756D-1, 9.595714699102763D-2, + > 4.990835442262655D-1, 1.013841131358227D-1, + > 5.858114200649882D-1, 1.049937517313658D-1, + > 6.683033272342810D-1, 1.067950801627900D-1, + > 7.438888779969691D-1, 1.067950801627900D-1, + > 8.104473950364958D-1, 1.049937517313658D-1, + > 8.665538736181088D-1, 1.013841131358227D-1, + > 9.115663629395599D-1, 9.595714699102763D-2, + > 9.456420843442170D-1, 8.871167329713092D-2, + > 9.696760021962513D-1, 7.966885558279603D-2, + > 9.851631047914915D-1, 6.889030401909524D-2, + > 9.939941211951136D-1, 5.649746567994538D-2, + > 9.982033803640970D-1, 4.268848001641819D-2, + > 9.996960998334762D-1, 2.774821193462954D-2, + > 9.999889214741531D-1, 1.206308084642689D-2/ +*---- +* Order = 46 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA=254,276)/ + > 4.316263740076508D-3, 1.106270311640170D-2, + > 2.261617115926903D-2, 2.547185620354161D-2, + > 5.502063421109040D-2, 3.925335576844175D-2, + > 1.006332490304859D-1, 5.207448910676365D-2, + > 1.581268223158908D-1, 6.368572226347834D-2, + > 2.257488308236949D-1, 7.390806352099788D-2, + > 3.013494595288169D-1, 8.262656871717886D-2, + > 3.824459703834832D-1, 8.977817435489181D-2, + > 4.663257595606548D-1, 9.533683737852415D-2, + > 5.501830693329847D-1, 9.929896064532817D-2, + > 6.312785125450332D-1, 1.016710684725573D-1, + > 7.071067811865475D-1, 1.024607276986862D-1, + > 7.755562130490158D-1, 1.016710684725573D-1, + > 8.350440648369020D-1, 9.929896064532817D-2, + > 8.846130713312902D-1, 9.533683737852415D-2, + > 9.239778567354501D-1, 8.977817435489181D-2, + > 9.535137666765435D-1, 8.262656871717886D-2, + > 9.741855395055578D-1, 7.390806352099788D-2, + > 9.874188108722047D-1, 6.368572226347834D-2, + > 9.949235896236294D-1, 5.207448910676365D-2, + > 9.984852176226794D-1, 3.925335576844175D-2, + > 9.997442216897753D-1, 2.547185620354161D-2, + > 9.999906848902774D-1, 1.106270311640170D-2/ +*---- +* Order = 48 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA=277,300)/ + > 3.972127715469235D-3, 1.018170331312922D-2, + > 2.082219684210920D-2, 2.346349353303599D-2, + > 5.069784852121022D-2, 3.621292193062666D-2, + > 9.284156235089204D-2, 4.814205277021821D-2, + > 1.461314233398479D-1, 5.903176230890425D-2, + > 2.090839313632731D-1, 6.872060753743357D-2, + > 2.798722573569018D-1, 7.710072080607362D-2, + > 3.563718475936031D-1, 8.410905852056142D-2, + > 4.362363542057792D-1, 8.971603507554962D-2, + > 5.170013989086162D-1, 9.391384633073917D-2, + > 5.962091754811982D-1, 9.670612366917621D-2, + > 6.715437180389010D-1, 9.809983760200038D-2, + > 7.409649349075090D-1, 9.809983760200038D-2, + > 8.028291344190430D-1, 9.670612366917621D-2, + > 8.559845521541460D-1, 9.391384633073917D-2, + > 8.998321194918805D-1, 8.971603507554962D-2, + > 9.343442118634447D-1, 8.410905852056142D-2, + > 9.600372490491982D-1, 7.710072080607362D-2, + > 9.778976989673706D-1, 6.872060753743357D-2, + > 9.892651854344568D-1, 5.903176230890425D-2, + > 9.956808948153246D-1, 4.814205277021821D-2, + > 9.987140372275341D-1, 3.621292193062666D-2, + > 9.997831945570342D-1, 2.346349353303599D-2, + > 9.999921110695884D-1, 1.018170331312922D-2/ +*---- +* Order = 50 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA=301,325)/ + > 3.667550894412352D-3, 9.401828916276707D-3, + > 1.923311271050877D-2, 2.168279634366668D-2, + > 4.686269256303738D-2, 3.350948773339766D-2, + > 8.591181096521539D-2, 4.463150509293552D-2, + > 1.354257349069718D-1, 5.485592665859271D-2, + > 1.941407846601355D-1, 6.403724294009948D-2, + > 2.604960566070634D-1, 7.207485785563070D-2, + > 3.326658203805228D-1, 7.890682566302415D-2, + > 4.086156970191086D-1, 8.450117205318058D-2, + > 4.861816572020138D-1, 8.884657138646223D-2, + > 5.631674054600955D-1, 9.194373283617630D-2, + > 6.374531035792810D-1, 9.379831926574879D-2, + > 7.071067811865475D-1, 9.441579330451362D-2, + > 7.704891567940086D-1, 9.379831926574879D-2, + > 8.263428304446916D-1, 9.194373283617630D-2, + > 8.738577665732013D-1, 8.884657138646223D-2, + > 9.127065312298296D-1, 8.450117205318058D-2, + > 9.430447772775977D-1, 7.890682566302415D-2, + > 9.654749113737599D-1, 7.207485785563070D-2, + > 9.809736773897386D-1, 6.403724294009948D-2, + > 9.907875000851124D-1, 5.485592665859271D-2, + > 9.963027455230047D-1, 4.463150509293552D-2, + > 9.989013404965188D-1, 3.350948773339766D-2, + > 9.998150265801494D-1, 2.168279634366668D-2, + > 9.999932745126024D-1, 9.401828916276707D-3/ +*---- +* Order = 52 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA=326,351)/ + > 3.396696521564986D-3, 8.708184858650910D-3, + > 1.781892050929334D-2, 2.009672209559701D-2, + > 4.344484549656114D-2, 3.109543899569998D-2, + > 7.972274758288012D-2, 4.148557359431433D-2, + > 1.258351250930473D-1, 5.109681624694140D-2, + > 1.806991980219270D-1, 5.979803209497304D-2, + > 2.429742598458005D-1, 6.749583467970308D-2, + > 3.110847255339362D-1, 7.413013546424966D-2, + > 3.832618562348258D-1, 7.966752788841578D-2, + > 4.576040998675997D-1, 8.409386224353523D-2, + > 5.321533124724104D-1, 8.740711539780842D-2, + > 6.049820056372407D-1, 8.961127788931235D-2, + > 6.742853437025943D-1, 9.071164194824712D-2, + > 7.384709034673433D-1, 9.071164194824712D-2, + > 7.962391430061333D-1, 8.961127788931235D-2, + > 8.466480095202735D-1, 8.740711539780842D-2, + > 8.891560536735741D-1, 8.409386224353523D-2, + > 9.236397293076104D-1, 7.966752788841578D-2, + > 9.503821828819580D-1, 7.413013546424966D-2, + > 9.700327360725438D-1, 6.749583467970308D-2, + > 9.835384079100482D-1, 5.979803209497304D-2, + > 9.920511686867856D-1, 5.109681624694140D-2, + > 9.968170762571418D-1, 4.148557359431433D-2, + > 9.990558269685333D-1, 3.109543899569998D-2, + > 9.998412304320539D-1, 2.009672209559701D-2, + > 9.999942312097307D-1, 8.708184858650910D-3/ +*---- +* Order = 54 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA=352,378)/ + > 3.154765578452475D-3, 8.088519222518445D-3, + > 1.655490845630564D-2, 1.867799415496943D-2, + > 4.038615839476834D-2, 2.893121574263062D-2, + > 7.417347458492733D-2, 3.865622019474338D-2, + > 1.172128302266212D-1, 4.770225064114276D-2, + > 1.685713218507331D-1, 5.595108446790947D-2, + > 2.270915986436383D-1, 6.331631488632194D-2, + > 2.914082582079124D-1, 6.974022040259096D-2, + > 3.599789431577641D-1, 7.518873114594296D-2, + > 4.311305838436490D-1, 7.964552607244304D-2, + > 5.031191739976250D-1, 8.310615532185532D-2, + > 5.741997608797084D-1, 8.557281096997747D-2, + > 6.427021081394999D-1, 8.705010633284512D-2, + > 7.071067811865475D-1, 8.754202768311478D-2, + > 7.661161793050990D-1, 8.705010633284512D-2, + > 8.187152341355850D-1, 8.557281096997747D-2, + > 8.642170426206298D-1, 8.310615532185532D-2, + > 9.022895431482260D-1, 7.964552607244304D-2, + > 9.329604281442055D-1, 7.518873114594296D-2, + > 9.565987806014759D-1, 6.974022040259096D-2, + > 9.738734033874612D-1, 6.331631488632194D-2, + > 9.856894589319176D-1, 5.595108446790947D-2, + > 9.931068182377288D-1, 4.770225064114276D-2, + > 9.972453537961454D-1, 3.865622019474338D-2, + > 9.991841462964235D-1, 2.893121574263062D-2, + > 9.998629581127622D-1, 1.867799415496943D-2, + > 9.999950237146908D-1, 8.088519222518445D-3/ +*---- +* Order = 56 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA=379,406)/ + > 2.937783068161103D-3, 7.532680050608221D-3, + > 1.542057328035316D-2, 1.740393870836662D-2, + > 3.763822344160284D-2, 2.698374387900945D-2, + > 6.917953054276372D-2, 3.610294286409177D-2, + > 1.094349422308658D-1, 4.462772409676748D-2, + > 1.575965218033132D-1, 5.245143079823476D-2, + > 2.126604305088603D-1, 5.949430183433835D-2, + > 2.734385687346029D-1, 6.570130151997985D-2, + > 3.385818199466888D-1, 7.103829899589742D-2, + > 4.066156050266695D-1, 7.548738262548303D-2, + > 4.759868689664975D-1, 7.904202969199018D-2, + > 5.451202571507631D-1, 8.170266224145589D-2, + > 6.124802047142074D-1, 8.347291855738130D-2, + > 6.766350014909858D-1, 8.435680753384398D-2, + > 7.363185959605350D-1, 8.435680753384398D-2, + > 7.904859257654386D-1, 8.347291855738130D-2, + > 8.383578622783268D-1, 8.170266224145589D-2, + > 8.794523867563669D-1, 7.904202969199018D-2, + > 9.135993376468678D-1, 7.548738262548303D-2, + > 9.409369538930799D-1, 7.103829899589742D-2, + > 9.618894682490145D-1, 6.570130151997985D-2, + > 9.771261644720124D-1, 5.949430183433835D-2, + > 9.875035880013286D-1, 5.245143079823476D-2, + > 9.939939604539492D-1, 4.462772409676748D-2, + > 9.976042264113975D-1, 3.610294286409177D-2, + > 9.992914310330896D-1, 2.698374387900945D-2, + > 9.998810958907590D-1, 1.740393870836662D-2, + > 9.999956847060113D-1, 7.532680050608221D-3/ +*---- +* Order = 58 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA=407,435)/ + > 2.742433926349367D-3, 7.032198002138893D-3, + > 1.439878984638582D-2, 1.625558572095192D-2, + > 3.516048827839076D-2, 2.522520908746962D-2, + > 6.466981500411214D-2, 3.379142344706885D-2, + > 1.023963769013716D-1, 4.183513829563056D-2, + > 1.476370627578143D-1, 4.926003874337995D-2, + > 1.995172380704506D-1, 5.599261639569382D-2, + > 2.569989752795515D-1, 6.198069665978488D-2, + > 3.188984690691561D-1, 6.719054793752873D-2, + > 3.839134581171505D-1, 7.160316710825960D-2, + > 4.506602512451940D-1, 7.521031896392870D-2, + > 5.177188670869931D-1, 7.801077736827083D-2, + > 5.836839333174925D-1, 8.000706366604473D-2, + > 6.472183977592858D-1, 8.120284618700015D-2, + > 7.071067811865475D-1, 8.160107162859419D-2, + > 7.623046278240104D-1, 8.120284618700015D-2, + > 8.119809517390300D-1, 8.000706366604473D-2, + > 8.555508019177823D-1, 7.801077736827083D-2, + > 8.926955460556631D-1, 7.521031896392870D-2, + > 9.233690793374722D-1, 7.160316710825960D-2, + > 9.477888828348582D-1, 6.719054793752873D-2, + > 9.664116755840962D-1, 6.198069665978488D-2, + > 9.798943166039586D-1, 5.599261639569382D-2, + > 9.890416056467216D-1, 4.926003874337995D-2, + > 9.947436765305283D-1, 4.183513829563056D-2, + > 9.979067165959621D-1, 3.379142344706885D-2, + > 9.993816788715036D-1, 2.522520908746962D-2, + > 9.998963320519581D-1, 1.625558572095192D-2, + > 9.999962395210092D-1, 7.032198002138893D-3/ +*---- +* Order = 60 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA=436,465)/ + > 2.565935816040464D-3, 6.579962619041378D-3, + > 1.347516554782726D-2, 1.521696776845086D-2, + > 3.291878087349994D-2, 2.363209088631991D-2, + > 6.058415580193666D-2, 3.169244521291083D-2, + > 9.600762130781930D-2, 3.929176131897425D-2, + > 1.385744973107413D-1, 4.634296332807871D-2, + > 1.875194470293309D-1, 5.277837509680457D-2, + > 2.419322530851564D-1, 5.854875269926412D-2, + > 3.007708877829240D-1, 6.362110627094798D-2, + > 3.628827241348177D-1, 6.797575914341677D-2, + > 4.270337881301431D-1, 7.160310605662001D-2, + > 4.919450854898168D-1, 7.450044532753710D-2, + > 5.563343213452209D-1, 7.666914603657626D-2, + > 6.189608105599492D-1, 7.811230591489950D-2, + > 6.786710589083211D-1, 7.883297571760606D-2, + > 7.344423692845193D-1, 7.883297571760606D-2, + > 7.854218707108751D-1, 7.811230591489950D-2, + > 8.309585566641411D-1, 7.666914603657626D-2, + > 8.706262302862330D-1, 7.450044532753710D-2, + > 9.042356682830091D-1, 7.160310605662001D-2, + > 9.318348182615274D-1, 6.797575914341677D-2, + > 9.536964260509063D-1, 6.362110627094798D-2, + > 9.702931438061076D-1, 5.854875269926412D-2, + > 9.822608904897995D-1, 5.277837509680457D-2, + > 9.903520125162948D-1, 4.634296332807871D-2, + > 9.953805988921094D-1, 3.929176131897425D-2, + > 9.981630929190713D-1, 3.169244521291083D-2, + > 9.994580300671972D-1, 2.363209088631991D-2, + > 9.999092058349391D-1, 1.521696776845086D-2, + > 9.999967079812753D-1, 6.579962619041378D-3/ +*---- +* Order = 62 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA=466,496)/ + > 2.405939628907848D-3, 6.169968963546275D-3, + > 1.263753399803356D-2, 1.427456861312436D-2, + > 3.088414635685555D-2, 2.218439608809088D-2, + > 5.687137210380912D-2, 2.978102161065714D-2, + > 9.019210028251855D-2, 3.696937092634360D-2, + > 1.303066527232758D-1, 4.367061943331126D-2, + > 1.765426012600548D-1, 4.982249645060048D-2, + > 2.280987746820670D-1, 5.537868201403226D-2, + > 2.840550769938310D-1, 6.030717484633422D-2, + > 3.433889379342239D-1, 6.458796732100688D-2, + > 4.049984831103632D-1, 6.821039566414624D-2, + > 4.677318207524097D-1, 7.117047743682441D-2, + > 5.304212533007296D-1, 7.346846419828167D-2, + > 5.919207725381885D-1, 7.510675391195046D-2, + > 6.511448995883448D-1, 7.608824036260866D-2, + > 7.071067811865475D-1, 7.641513111317905D-2, + > 7.589534371356944D-1, 7.608824036260866D-2, + > 8.059961532400723D-1, 7.510675391195046D-2, + > 8.477342119124858D-1, 7.346846419828167D-2, + > 8.838704338621333D-1, 7.117047743682441D-2, + > 9.143173566537523D-1, 6.821039566414624D-2, + > 9.391932907045310D-1, 6.458796732100688D-2, + > 9.588079647322652D-1, 6.030717484633422D-2, + > 9.736379968902917D-1, 5.537868201403226D-2, + > 9.842930000463954D-1, 4.982249645060048D-2, + > 9.914737395695640D-1, 4.367061943331126D-2, + > 9.959243872135214D-1, 3.696937092634360D-2, + > 9.983815137686745D-1, 2.978102161065714D-2, + > 9.995229709735581D-1, 2.218439608809088D-2, + > 9.999201431786683D-1, 1.427456861312436D-2, + > 9.999971057230626D-1, 6.169968963546275D-3/ +*---- +* Order = 64 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA=497,528)/ + > 2.260451006401776D-3, 5.797117680017662D-3, + > 1.187555474628660D-2, 1.341688614356055D-2, + > 2.903192296767898D-2, 2.086504628721960D-2, + > 5.348772246647195D-2, 2.803569126620853D-2, + > 8.488403572313128D-2, 3.484354779039823D-2, + > 1.227451197406277D-1, 4.121716160665880D-2, + > 1.664778916763791D-1, 4.709924845926102D-2, + > 2.153746593423551D-1, 5.244635195862705D-2, + > 2.686204966167571D-1, 5.722762992459373D-2, + > 3.253060318502849D-1, 6.142301931103583D-2, + > 3.844458467127314D-1, 6.502107279693194D-2, + > 4.450021333238041D-1, 6.801672562809563D-2, + > 5.059127742515483D-1, 7.040918982162950D-2, + > 5.661226261485058D-1, 7.220010760651431D-2, + > 6.246165176514237D-1, 7.339204014218695D-2, + > 6.804523156165483D-1, 7.398732697450898D-2, + > 7.327923622500970D-1, 7.398732697450898D-2, + > 7.809316268899658D-1, 7.339204014218695D-2, + > 8.243210370739783D-1, 7.220010760651431D-2, + > 8.625846421361221D-1, 7.040918982162950D-2, + > 8.955295089148449D-1, 6.801672562809563D-2, + > 9.231475455989314D-1, 6.502107279693194D-2, + > 9.456087910133986D-1, 6.142301931103583D-2, + > 9.632460894274976D-1, 5.722762992459373D-2, + > 9.765314926376745D-1, 5.244635195862705D-2, + > 9.860451873940605D-1, 4.709924845926102D-2, + > 9.924382275889311D-1, 4.121716160665880D-2, + > 9.963908372116607D-1, 3.484354779039823D-2, + > 9.985685071868378D-1, 2.803569126620853D-2, + > 9.995784848869041D-1, 2.086504628721960D-2, + > 9.999294831134183D-1, 1.341688614356055D-2, + > 9.999974451773603D-1, 5.797117680017662D-3/ +*---- +* Order = 66 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA=529,561)/ + > 2.127767955855947D-3, 5.457055989066269D-3, + > 1.118039441421181D-2, 1.263408341121632D-2, + > 2.734100277100129D-2, 1.965938521651337D-2, + > 5.039565382935335D-2, 2.643794451890077D-2, + > 8.002669926187768D-2, 3.289309077657504D-2, + > 1.158131591183225D-1, 3.895995661650402D-2, + > 1.572300146684183D-1, 4.458584184375271D-2, + > 2.036500183342073D-1, 4.973005671396442D-2, + > 2.543492281175076D-1, 5.436302857337294D-2, + > 3.085169689930697D-1, 5.846485650012917D-2, + > 3.652704783006284D-1, 6.202354684819099D-2, + > 4.236740543125242D-1, 6.503314324601964D-2, + > 4.827621314188725D-1, 6.749192029021972D-2, + > 5.415653801336422D-1, 6.940075957307684D-2, + > 5.991386896375901D-1, 7.076178092149795D-2, + > 6.545897375326029D-1, 7.157726645471736D-2, + > 7.071067811865475D-1, 7.184889180746158D-2, + > 7.559843090415290D-1, 7.157726645471736D-2, + > 8.006452588876996D-1, 7.076178092149795D-2, + > 8.406586340725370D-1, 6.940075957307684D-2, + > 8.757515198204952D-1, 6.749192029021972D-2, + > 9.058147137811288D-1, 6.503314324601964D-2, + > 9.309014328499179D-1, 6.202354684819099D-2, + > 9.512188390918934D-1, 5.846485650012917D-2, + > 9.671124392520386D-1, 5.436302857337294D-2, + > 9.790437528693378D-1, 4.973005671396442D-2, + > 9.875620094390878D-1, 4.458584184375271D-2, + > 9.932710164778967D-1, 3.895995661650402D-2, + > 9.967927203813484D-1, 3.289309077657504D-2, + > 9.987293317386408D-1, 2.643794451890077D-2, + > 9.996261649073998D-1, 1.965938521651337D-2, + > 9.999374974370861D-1, 1.263408341121632D-2, + > 9.999977362992008D-1, 5.457055989066269D-3/ +*---- +* Order = 68 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA=562,595)/ + > 2.006430897540039D-3, 5.146050313292597D-3, + > 1.054447088264912D-2, 1.191770807990189D-2, + > 2.579323625075468D-2, 1.855478004998946D-2, + > 4.756278722072234D-2, 2.497175458640582D-2, + > 7.557097996402496D-2, 3.109953168579275D-2, + > 1.094439543367078D-1, 3.687913529406050D-2, + > 1.487153270030891D-1, 4.226206890226890D-2, + > 1.928273387808641D-1, 4.721028879193994D-2, + > 2.411349977536837D-1, 5.169555561856539D-2, + > 2.929138510155265D-1, 5.569829899210566D-2, + > 3.473716688709153D-1, 5.920617987876196D-2, + > 4.036639760910661D-1, 6.221252672611230D-2, + > 4.609130411977627D-1, 6.471478978336624D-2, + > 5.182296606197990D-1, 6.671311918267942D-2, + > 5.747368634369850D-1, 6.820913521902009D-2, + > 6.295945185532131D-1, 6.920492910010646D-2, + > 6.820237474599884D-1, 6.970231119307893D-2, + > 7.313300266641826D-1, 6.970231119307893D-2, + > 7.769238973076757D-1, 6.920492910010646D-2, + > 8.183382783461961D-1, 6.820913521902009D-2, + > 8.552414973876606D-1, 6.671311918267942D-2, + > 8.874453044857635D-1, 6.471478978336624D-2, + > 9.149073146534305D-1, 6.221252672611230D-2, + > 9.377275316774245D-1, 5.920617987876196D-2, + > 9.561388371378154D-1, 5.569829899210566D-2, + > 9.704915830950472D-1, 5.169555561856539D-2, + > 9.812327029908297D-1, 4.721028879193994D-2, + > 9.888800491032086D-1, 4.226206890226890D-2, + > 9.939929682141341D-1, 3.687913529406050D-2, + > 9.971404249087873D-1, 3.109953168579275D-2, + > 9.988682502071013D-1, 2.497175458640582D-2, + > 9.996672991369242D-1, 1.855478004998946D-2, + > 9.999444055215295D-1, 1.191770807990189D-2, + > 9.999979871155008D-1, 5.146050313292597D-3/ +*---- +* Order = 70 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA=596,630)/ + > 1.895182390808144D-3, 4.860883544830628D-3, + > 9.961246688281479D-3, 1.126046540783217D-2, + > 2.437294975502797D-2, 1.754029683264231D-2, + > 4.496109079691220D-2, 2.362319238448706D-2, + > 7.147420525157438D-2, 2.944673122343525D-2, + > 1.035791506531498D-1, 3.495721171268495D-2, + > 1.408602617289225D-1, 4.010998546152185D-2, + > 1.828200227142416D-1, 4.486951999487671D-2, + > 2.288821565254979D-1, 4.920893733923423D-2, + > 2.783976951621312D-1, 5.310912506931291D-2, + > 3.306543210696288D-1, 5.655756629702146D-2, + > 3.848889902593347D-1, 5.954703320530382D-2, + > 4.403035834865983D-1, 6.207426689302348D-2, + > 4.960831016113237D-1, 6.413873664427702D-2, + > 5.514157377272589D-1, 6.574154187441527D-2, + > 6.055140267651983D-1, 6.688449453259325D-2, + > 6.576361927484648D-1, 6.756940076441969D-2, + > 7.071067811865475D-1, 6.779754843107255D-2, + > 7.533356741767272D-1, 6.756940076441969D-2, + > 7.958346331937275D-1, 6.688449453259325D-2, + > 8.342305941325233D-1, 6.574154187441527D-2, + > 8.682750464545719D-1, 6.413873664427702D-2, + > 8.978489596635172D-1, 6.207426689302348D-2, + > 9.229628731304146D-1, 5.954703320530382D-2, + > 9.437519377240943D-1, 5.655756629702146D-2, + > 9.604658886854926D-1, 5.310912506931291D-2, + > 9.734541378124793D-1, 4.920893733923423D-2, + > 9.831463976920040D-1, 4.486951999487671D-2, + > 9.900294877758235D-1, 4.010998546152185D-2, + > 9.946212141061401D-1, 3.495721171268495D-2, + > 9.974424484568801D-1, 2.944673122343525D-2, + > 9.989887388326016D-1, 2.362319238448706D-2, + > 9.997029355364717D-1, 1.754029683264231D-2, + > 9.999503855514108D-1, 1.126046540783217D-2, + > 9.999982041402402D-1, 4.860883544830628D-3/ +*---- +* Order = 72 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA=631,666)/ + > 1.792934454202931D-3, 4.598771655133457D-3, + > 9.425061112313890D-3, 1.065603341978431D-2, + > 2.306655207082563D-2, 1.660643477465490D-2, + > 4.256620201285560D-2, 2.238010857286012D-2, + > 6.769916617188922D-2, 2.792054146952203D-2, + > 9.816763086953637D-2, 3.317875933222189D-2, + > 1.335999707374397D-1, 3.811363223988677D-2, + > 1.735510826588126D-1, 4.269199698797202D-2, + > 2.175046751209386D-1, 4.688833809218218D-2, + > 2.648780127840390D-1, 5.068409479675100D-2, + > 3.150293990097726D-1, 5.406671002171079D-2, + > 3.672684559193769D-1, 5.702853994383097D-2, + > 4.208690454062535D-1, 5.956572815765240D-2, + > 4.750844817568741D-1, 6.167712597050399D-2, + > 5.291645281141019D-1, 6.336331655665649D-2, + > 5.823735501979260D-1, 6.462577941116331D-2, + > 6.340091223112186D-1, 6.546621488177129D-2, + > 6.834203406787824D-1, 6.588603711319039D-2, + > 7.300250940525928D-1, 6.588603711319039D-2, + > 7.733255671618250D-1, 6.546621488177129D-2, + > 8.129213049427747D-1, 6.462577941116331D-2, + > 8.485192409048718D-1, 6.336331655665649D-2, + > 8.799401884183961D-1, 6.167712597050399D-2, + > 9.071214067691430D-1, 5.956572815765240D-2, + > 9.301149828309383D-1, 5.702853994383097D-2, + > 9.490819130926168D-1, 5.406671002171079D-2, + > 9.642819288691345D-1, 5.068409479675100D-2, + > 9.760592790914571D-1, 4.688833809218218D-2, + > 9.848248685466640D-1, 4.269199698797202D-2, + > 9.910353413571865D-1, 3.811363223988677D-2, + > 9.951698931586820D-1, 3.317875933222189D-2, + > 9.977057797264838D-1, 2.792054146952203D-2, + > 9.990936484865674D-1, 2.238010857286012D-2, + > 9.997339316916096D-1, 1.660643477465490D-2, + > 9.999555831250852D-1, 1.065603341978431D-2, + > 9.999983926917297D-1, 4.598771655133457D-3/ +*---- +* Order = 74 +*---- + DATA ((SYST(IS,IA),IS=1,2),IA=667,703)/ + > 1.698741886354622D-3, 4.357295595609470D-3, + > 8.930992904023735D-3, 1.009891155595077D-2, + > 2.186221185368514D-2, 1.574490757594847D-2, + > 4.035686918939289D-2, 2.123186989304243D-2, + > 6.421330823399199D-2, 2.650852294936087D-2, + > 9.316448689887876D-2, 3.153013523379523D-2, + > 1.268771624175196D-1, 3.625879153931750D-2, + > 1.649519866208174D-1, 4.066355379567878D-2, + > 2.069251875729949D-1, 4.472025016442803D-2, + > 2.522722781263537D-1, 4.841093781051414D-2, + > 3.004140595738183D-1, 5.172313096383708D-2, + > 3.507250096820162D-1, 5.464889158975430D-2, + > 4.025439676980595D-1, 5.718387043152167D-2, + > 4.551868677602700D-1, 5.932636945044961D-2, + > 5.079611361117360D-1, 6.107647781254220D-2, + > 5.601812617687321D-1, 6.243531594168673D-2, + > 6.111849763541508D-1, 6.340440772608052D-2, + > 6.603494357038005D-1, 6.398519049059308D-2, + > 7.071067811865475D-1, 6.417866575467489D-2, + > 7.509584694013859D-1, 6.398519049059308D-2, + > 7.914877918698279D-1, 6.340440772608052D-2, + > 8.283700585868560D-1, 6.243531594168673D-2, + > 8.613799882746722D-1, 6.107647781254220D-2, + > 8.903959318295398D-1, 5.932636945044961D-2, + > 9.154006522118627D-1, 5.718387043152167D-2, + > 9.364784928569100D-1, 5.464889158975430D-2, + > 9.538088869424411D-1, 5.172313096383708D-2, + > 9.676562910914906D-1, 4.841093781051414D-2, + > 9.783567686421354D-1, 4.472025016442803D-2, + > 9.863015979454995D-1, 4.066355379567878D-2, + > 9.919184369981628D-1, 3.625879153931750D-2, + > 9.956507311205414D-1, 3.153013523379523D-2, + > 9.979361958790984D-1, 2.650852294936087D-2, + > 9.991853297107750D-1, 2.123186989304243D-2, + > 9.997609932843272D-1, 1.574490757594847D-2, + > 9.999601178875827D-1, 1.009891155595077D-2, + > 9.999985571369608D-1, 4.357295595609470D-3/ +*---- +* Start processing +*---- + PI=XDRCST('Pi',' ') + NBPT=(ORDRE*(ORDRE+2))/8 + IF(3*NBPT .NE. NBANGL) CALL XABORT(NAMSBR// + >': Number of quadrature points is invalid') + ILEVEL=1 + JLEVEL=0 + NLEVEL=ORDRE/2 + KLEVEL=(NLEVEL*(NLEVEL-1))/2 + IPOINT=0 + DO IPT=1,NBPT + JLEVEL = JLEVEL + 1 + ALPHA =PI*(DBLE(NLEVEL-ILEVEL+1-JLEVEL)+DONE/DTWO) + > /(DTWO*DBLE(NLEVEL-ILEVEL+1)) + XI=SYST(1,KLEVEL+ILEVEL) + ROTXI=SQRT(DONE-XI**2) + COSA=COS(ALPHA) + SINA=SIN(ALPHA) + WGTINV=3.0D0*PI*DBLE(ORDRE+2.0D0)/SYST(2,KLEVEL+ILEVEL) + DO IROT=1,3 +*---- +* \xi direction (3) +* \mu direction (2) +* \eta direction (1) +*---- + IDIR1=MOD(IROT-1,3)+1 + IDIR2=MOD(IROT,3)+1 + IDIR3=MOD(IROT+1,3)+1 + IPOINT=IPOINT+1 + DDENWT(1,IPOINT)=DQUAD(1)*WGTINV + DANGLT(IDIR1,1,IPOINT)=COSA*ROTXI + DANGLT(IDIR2,1,IPOINT)=SINA*ROTXI + DANGLT(IDIR3,1,IPOINT)=XI + DDENWT(2,IPOINT)=DQUAD(2)*WGTINV + DANGLT(1,2,IPOINT)=-DANGLT(1,1,IPOINT) + DANGLT(2,2,IPOINT)=DANGLT(2,1,IPOINT) + DANGLT(3,2,IPOINT)=DANGLT(3,1,IPOINT) + DDENWT(3,IPOINT)=DQUAD(3)*WGTINV + DANGLT(1,3,IPOINT)=DANGLT(1,1,IPOINT) + DANGLT(2,3,IPOINT)=-DANGLT(2,1,IPOINT) + DANGLT(3,3,IPOINT)=DANGLT(3,1,IPOINT) + DDENWT(4,IPOINT)=DQUAD(4)*WGTINV + DANGLT(1,4,IPOINT)=-DANGLT(1,1,IPOINT) + DANGLT(2,4,IPOINT)=-DANGLT(2,1,IPOINT) + DANGLT(3,4,IPOINT)=DANGLT(3,1,IPOINT) + ENDDO + IF(JLEVEL .EQ. NLEVEL-ILEVEL+1) THEN + ILEVEL = ILEVEL + 1 + JLEVEL = 0 + ENDIF + ENDDO + RETURN + END |