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

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

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

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

*W->H0[0][] = 

7.0710678118654752440084436210484930e-01
3.4045953094161120620927407034909500e-35
6.4687310878906129179762073366328060e-34
0.0000000000000000000000000000000000e+00
7.1496501497738353303947554773309960e-34
8.5114882735402801552318517587273760e-36

*W->H0[1][] = 

-6.1237243569579452454932101867647290e-01
3.5355339059327376220042218105242480e-01
-5.2664833692530483460497082757125640e-34
5.4579918554077046495424249402839300e-34
-6.3304194034455833654536897455534860e-34
6.2825422819069192895805105794106440e-34

*W->H0[2][] = 

1.5320678892372504279417333165709280e-34
-6.8465319688145764182121222850100340e-01
1.7677669529663688110021109052621270e-01
-1.0384015693719141789382859145647400e-33
3.4045953094161120620927407034909500e-34
-1.2597002644839614629743140602916520e-33

*W->H0[3][] = 

2.3385358667337133659898429576978400e-01
4.0504629365049126443537296475550040e-01
-5.2291251658379721748635751611574300e-01
8.8388347648318440550105545263106810e-02
-5.5537460984850328012887832725696130e-34
8.3838159494371759529033739823464660e-34

*W->H0[4][] = 

6.5538459706260157195285258542200790e-34
1.5309310892394863113733025466911760e-01
5.9292706128157112474979253958113590e-01
-3.5078038001005700489847644365467680e-01
4.4194173824159220275052772631554010e-02
8.5114882735402801552318517587273760e-36

*W->H0[5][] = 

-1.4657549249448217358017594104826480e-01
-2.5387620014487376126437136150831340e-01
-1.6387638252658617921741461151249080e-01
5.8170345215582140294374571666169480e-01
-2.1986323874172326037026391157239750e-01
2.2097086912079610137526386315776600e-02

*W->G0[0][] = 

1.2133977682759023389298527867241750e-31
-6.9006555934235421779612686462589100e-02
-2.6726124191242438468455348087978420e-01
-4.2163702135578391093318580592441500e-01
4.7809144373375745598752687187727570e-01
-1.3213749452868198044514027567716850e-01

*W->G0[1][] = 

-1.0762273838572729001861955871216280e-01
-1.8640805093377296209355800189985920e-01
-1.8289533680592327413571665445458950e-01
5.6948600236879883495794989227731850e-02
5.3811369192863645009309779356043410e-01
-3.5694424213603681723889654848142740e-01

*W->G0[2][] = 

3.5031583436237085062903255468570130e-31
4.1781424448965464011267160944930970e-02
1.6181876107166664196235729246938180e-01
3.1972457510904210729081841957928010e-01
2.2197354266027465795417199908591610e-01
-5.6619161382559877530671894318940420e-01

*W->G0[3][] = 

-1.0556807713333623553602933857226950e-01
-1.8284927325228855536503070724394700e-01
-2.0467675073099876253922251221279470e-01
-9.3656473141812652716006795963497140e-02
2.4527821095860462661248086402436250e-01
5.8702301176296628387694400324571330e-01

*W->G0[4][] = 

-1.2499971678521255435973497492867020e-31
3.6123507707732773099237584920748520e-02
1.3990574375864429108223108895606880e-01
3.1137042777624967381126057254524980e-01
4.6925813023049334813281412868645560e-01
4.0245091498622110281109057781601840e-01

*W->G0[5][] = 

-1.6155758191037421171369086641443490e-01
-2.7982594021673869501669952381958880e-01
-3.5265245601703738148387456022286050e-01
-3.7655532889563117919130132163481280e-01
-3.1926855472764427552943671219976500e-01
-1.6585086800094435039877444255816000e-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   -5e-33   8e-34   1e-32   8e-32   -1e-31   
-5e-33   1e+00   6e-33   2e-32   6e-32   -1e-31   
8e-34   6e-33   1e+00   2e-32   -4e-32   -7e-32   
1e-32   2e-32   2e-32   1e+00   -6e-32   -1e-32   
8e-32   6e-32   -4e-32   -6e-32   1e+00   2e-31   
-1e-31   -1e-31   -7e-32   -1e-32   2e-31   1e+00   

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

1e+00   -9e-34   7e-33   -8e-33   -3e-32   5e-31   
-9e-34   1e+00   6e-33   -1e-32   -2e-32   2e-31   
7e-33   6e-33   1e+00   -6e-33   8e-33   2e-31   
-8e-33   -1e-32   -6e-33   1e+00   7e-32   -5e-32   
-3e-32   -2e-32   8e-33   7e-32   1e+00   1e-32   
5e-31   2e-31   2e-31   -5e-32   1e-32   1e+00   

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

2e-33   -8e-33   1e-32   -3e-32   7e-32   4e-31   
-7e-34   -8e-33   1e-32   -3e-32   -2e-32   8e-31   
-6e-33   -6e-33   9e-33   -3e-32   5e-33   -4e-32   
-2e-32   -9e-33   -2e-32   5e-33   -2e-32   -2e-31   
-3e-32   -4e-32   -5e-32   8e-32   -9e-32   -2e-31   
3e-31   6e-32   9e-32   2e-32   -1e-31   -7e-32   
The size of double is: 8 bytes.
The size of long double is: 16 bytes.
