
Routine: Get_LegendreRoots():
 Read in quadrature of order: 9

Routine: Get_GaussLegendreWeights():
 Read in quadrature of order: 9

Routine: Get_GaussLegendreWeights():
 Read in quadrature of order: 17

Routine: Get_LegendreRoots():
 Read in quadrature of order: 17

Routine: Get_LegendreRoots():
 Read in quadrature of order: 9

Routine: Get_GaussLegendreWeights():
 Read in quadrature of order: 9

*W->H0[0][] = 

8.9677394624175180893515288515284830e-01
4.3287010976127198874576656897769070e-01
-6.2724794791210059653755990858418870e-02
-3.8718790548102617070477368456545000e-02
4.1160002551210260213355228122905730e-02
4.0511801368623776781748053139165770e-03
-2.0698080249425369814705905815466870e-02
-1.0948938531985342670181569207059990e-02
-1.4626487037240496162325011770816160e-03

*W->H0[1][] = 

-2.1136004387220282973367821193568580e-01
5.8030547934681617061534780455352000e-01
7.5650009359140409272053574255699720e-01
1.4934333948637162524833950098893300e-01
-1.2568865037077626622249565210839600e-01
-1.1234624093306780107423243463138240e-02
5.4692955864012516466668304107175140e-02
2.8203259630049580235036159330294550e-02
3.7223695256439681404176881794178340e-03

*W->H0[2][] = 

1.2150357152902569405991685371705730e-01
-2.2479608320948304459034178592614400e-01
1.6530208968486260342713500707368340e-01
7.6677540466705633617622509163556450e-01
5.3624754018922449931275597271196490e-01
2.9241538659563923040399801481841830e-02
-1.2045812116251189601788481993348160e-01
-5.7607630758707214462927375204565820e-02
-7.3550046128868757890882996978063750e-03

*W->H0[3][] = 

-8.1762705191507560807881587387840190e-02
1.3814411041185960759707867147302750e-01
-7.9335807276835383514073557753587120e-02
-1.3252545047955107374314845754355570e-01
4.1432906672209983259561888630018420e-01
7.0845036330578384729160824999925500e-01
4.6284791309856267638866856373007540e-01
1.5134094827031872354685066061560410e-01
1.7069790761494660232965549076828800e-02

*W->H0[4][] = 

5.7953170111686133448666439586129720e-02
-9.4453346698625495492947182902345770e-02
5.0156101161104338869822918492910680e-02
7.1488232533764191232038838443461420e-02
-1.5395180750610385523196458910559060e-01
-3.0189150262472751263564437743617680e-02
3.3141058058098439375731959092574670e-01
4.7276499121127176724493953792273020e-01
3.4948753888080918653122276530180040e-01

*W->H0[5][] = 

-4.1235062674420946550919332331761910e-02
6.6028165232462624912136902093807670e-02
-3.3842807337666314637781011809688130e-02
-4.5395084656599937229232415433427100e-02
8.8342872462696752926939587156545940e-02
1.4572736701305452539078314460107890e-02
-1.1709561232699148619466650985681970e-01
-9.0255536120355286404644731361647070e-02
-1.5472081465470267690945476702623540e-02

*W->H0[6][] = 

2.8196429121494812478449237299634550e-02
-4.4719463972866310536220616180668520e-02
2.2502795098483324023133003260612780e-02
2.9309606564265500140010872991346440e-02
-5.4604555252468865784909937731059030e-02
-8.4688928751435419570065267827047640e-03
6.2724917436840617851829387859270220e-02
4.4023421996254311677067766787700410e-02
6.9828691956696521898610848815702780e-03

*W->H0[7][] = 

-1.7192161279087256692714587121250360e-02
2.7123579143294025824576741925104850e-02
-1.3514070532490013068225390621352930e-02
-1.7335599003045865906057817675485250e-02
3.1612509716928727255424817877782920e-02
4.7675226847998779056860106735689990e-03
-3.4149939542721155259569130355864270e-02
-2.3165809847347266368791798614321310e-02
-3.5832542996633366240770847160661370e-03

*W->H0[8][] = 

7.3128975970306274105866786340362280e-03
-1.1508013062175357628568033423240560e-02
5.7066678353819191123758478101190060e-03
7.2681978194179660047123731361450710e-03
-1.3124913843853625584697190212373390e-02
-1.9550845540340733142206906008489390e-03
1.3808158040086009482052983935221000e-02
9.2393785530239950730171495291928860e-03
1.4153249305106395098593720842618380e-03

*W->G0[0][] = 

-2.4859093314530531664023399924848980e-01
3.7795928504103152994749905873188580e-01
-1.7636230472866835569415338213074500e-01
-2.0632568576521048401167936041477890e-01
3.3553082552213826990843985263743460e-01
4.4461322478774175743256317962993900e-02
-2.7925949147149921935891866372985960e-01
-1.6894852556863133657062325127772340e-01
-2.4252087181572630979246783849170540e-02

*W->G0[1][] = 

1.9353237994978105228611101306250880e-01
-3.5958659587388866218310982931477150e-01
3.3783221624828429140601780440293170e-01
-8.8646165000192755228268250007729860e-02
-2.4757551635068360171236451980590000e-01
2.8762730632148416302555941811236020e-01
1.3297228298399626906520796779759530e-01
-1.5632433253639860312606091424290260e-01
-1.5846741132279009007172399506014040e-01

*W->G0[2][] = 

-1.3319883883870905123124696256213260e-01
2.8438453533799581611918386004201680e-01
-3.7210465761318370788189393349136100e-01
3.2797566669560081128158016267710750e-01
-9.1825267821185521847820551281341050e-02
-2.3661205988716950831658582433554300e-01
2.2639542595814095770559343511618970e-01
1.9713507732349130836393835015882140e-01
2.8623712725193252563392491296115730e-02

*W->G0[3][] = 

-7.4751394037559663504856142540221370e-02
1.7645486228831040950948180244409280e-01
-2.8359828442253959317601282535124080e-01
3.7206632282252146367498478940865270e-01
-3.6638395991602604901612630606051190e-01
1.4665317863340390524573036223565130e-01
2.1574380217493415586216536778766670e-01
-1.1967970754082728692059416616973670e-01
-1.6673809911711883582019838859205050e-01

*W->G0[4][] = 

-2.9304542004345794897786591889674880e-02
7.5386425949587609565807747497632620e-02
-1.4258190216245648782506493504650000e-01
2.4316131401823728189286103435029530e-01
-3.7007104120303128622071815988724210e-01
4.3186043558099376908046370381632720e-01
-2.0113721468156848330838133755547590e-01
-2.1780388661610328492870978118053670e-01
-5.1481876703182364041823387169638070e-02

*W->G0[5][] = 

4.6110752220046976052019375422661460e-03
-1.3953252989674132503145910988574390e-02
3.4045371884240490069614697316745620e-02
-8.0147150656717956773612886138757020e-02
1.8153475778567219228359882737689120e-01
-3.6408058856125798595231844102876130e-01
5.1034792448394780632660054283796660e-01
-1.8284128123920607280223153778546430e-01
-1.8112401886413036495746729785662660e-01

*W->G0[6][] = 

2.4700024627234099672120898775372290e-03
-5.6187136090233206552736806473720450e-03
7.9221003774293677679553368555533730e-03
-5.2089251387237807823846461213973810e-03
-1.8734934147559608008332880433863420e-02
1.1498037336802982577432827828076110e-01
-3.7211458583443269352110960444582460e-01
5.8977006771100599292009575199279570e-01
1.7189513611831977020888624726087840e-03

*W->G0[7][] = 

-1.9438804980078217646746678833007910e-03
4.9518377665329218138971794830807360e-03
-9.3273521670253767885269124441006600e-03
1.6225075423252757573996014452824970e-02
-2.5796615870562012324548282068065620e-02
2.6024154304232522274718129888902020e-02
6.1736337845568780680404169689711230e-02
-4.4981070598571721518266358104902270e-01
5.4049762204179036368296800344526920e-01

*W->G0[8][] = 

5.8483857527948985064600437871869290e-04
-1.5222355116067128552574446445161050e-03
3.0229044774672154019881579595191830e-03
-5.9184529621039552727678789972657020e-03
1.2319078501569952217083110139081170e-02
-2.6903739636159471620282745404979210e-02
4.8442481133737514728673423660982330e-02
3.8142501282230080250379238863321920e-02
-7.0375853774000354969282709929634560e-01

Checking the orthogonality conditions on the filters:
(see: Alpert, Beylkin, Gines, Vozovoi).
OBS: These filters should really be computed using extended precision.

The matrix identity: Id = (H0^T)H0+(G0^T)G0, has righthand side equal:

1e+00   -2e-29   2e-29   -1e-29   -8e-29   -3e-28   -3e-28   -3e-28   -2e-28   
-2e-29   1e+00   -2e-29   4e-30   2e-28   7e-28   7e-28   6e-28   3e-28   
2e-29   -2e-29   1e+00   -4e-30   -3e-28   -9e-28   -1e-27   -9e-28   -5e-28   
-1e-29   4e-30   -4e-30   1e+00   4e-28   1e-27   1e-27   1e-27   5e-28   
-8e-29   2e-28   -3e-28   4e-28   1e+00   -6e-28   -2e-27   -7e-28   -7e-28   
-3e-28   7e-28   -9e-28   1e-27   -6e-28   1e+00   3e-28   1e-27   2e-28   
-3e-28   7e-28   -1e-27   1e-27   -2e-27   3e-28   1e+00   -5e-29   -6e-28   
-3e-28   6e-28   -9e-28   1e-27   -7e-28   1e-27   -5e-29   1e+00   4e-29   
-2e-28   3e-28   -5e-28   5e-28   -7e-28   2e-28   -6e-28   4e-29   1e+00   

The matrix identity: Id = (H1^T)H1+(G1^T)G1, has righthand side equal:

1e+00   -2e-27   5e-27   -6e-27   4e-27   -3e-27   3e-27   -3e-27   1e-27   
-2e-27   1e+00   4e-27   -4e-27   6e-27   -6e-27   3e-27   -1e-27   2e-28   
5e-27   4e-27   1e+00   -4e-27   7e-28   -2e-28   2e-27   -2e-27   1e-27   
-6e-27   -4e-27   -4e-27   1e+00   -4e-28   -6e-29   -2e-27   2e-27   -1e-27   
4e-27   6e-27   7e-28   -4e-28   1e+00   3e-27   -7e-28   -7e-28   7e-28   
-3e-27   -6e-27   -2e-28   -6e-29   3e-27   1e+00   8e-28   5e-28   -6e-28   
3e-27   3e-27   2e-27   -2e-27   -7e-28   8e-28   1e+00   -1e-27   7e-28   
-3e-27   -1e-27   -2e-27   2e-27   -7e-28   5e-28   -1e-27   1e+00   -7e-28   
1e-27   2e-28   1e-27   -1e-27   7e-28   -6e-28   7e-28   -7e-28   1e+00   

The matrix identity: 0 = (H0^T)H1+(G0^T)G1, has righthand side equal:

-3e-28   2e-28   -6e-28   6e-28   -5e-28   4e-28   -3e-28   3e-28   -1e-28   
5e-28   -4e-29   6e-28   -7e-28   4e-28   -3e-28   3e-28   -3e-28   2e-28   
-3e-28   -8e-28   4e-28   -3e-28   6e-28   -6e-28   3e-28   -7e-29   -1e-29   
2e-28   1e-27   -8e-28   8e-28   -1e-27   1e-27   -7e-28   3e-28   -7e-29   
-2e-29   -1e-27   1e-27   -9e-28   1e-27   -1e-27   6e-28   -2e-28   4e-29   
3e-27   -5e-27   7e-27   -7e-27   7e-27   -6e-27   4e-27   -3e-27   1e-27   
3e-27   -8e-27   1e-26   -1e-26   1e-26   -9e-27   7e-27   -4e-27   2e-27   
3e-27   -6e-27   9e-27   -9e-27   9e-27   -7e-27   6e-27   -4e-27   2e-27   
2e-27   -4e-27   5e-27   -5e-27   5e-27   -4e-27   3e-27   -2e-27   1e-27   
The size of double is: 8 bytes.
The size of long double is: 16 bytes.
