
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

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

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

*W->H0[0][] = 

8.9189506634159576114849079957920480e-01
4.4171768929215586116840541671099860e-01
-4.8506755332265860778342312101326840e-02
-6.0238062513699039054147678652987980e-02
1.9912857204826944506986621459659730e-02
2.5376473704429912000140350853227650e-02

*W->H0[1][] = 

-2.0503088510365545820920981244627780e-01
5.5161904794107287241801557255778740e-01
7.5288410549159763962049399286694940e-01
2.4606259901592953072861140156895890e-01
-6.4160982587815717259090775809459500e-02
-7.5482212001576495934254922474268890e-02

*W->H0[2][] = 

1.1126158912849193978947809307307010e-01
-2.0432501830114280756511885868054480e-01
1.0726550217260871932907912327576570e-01
6.2891758437888047125609602241193960e-01
5.9007609818445820638703444743372210e-01
2.9927287357069407271555264261042920e-01

*W->H0[3][] = 

-6.7191268422254577427417371282492060e-02
1.1312640123578521976339909205421480e-01
-4.7597130051220507432689750313455540e-02
-1.4427531891402568133030236569843360e-01
1.0011251395224762002611561161632410e-01
2.2507611198043699878758912593377240e-01

*W->H0[4][] = 

3.8428121768482366135530848149973650e-02
-6.2636275602958667815974477823096010e-02
2.4635080773997013381734158316415300e-02
6.6249393278137612378542495172418000e-02
-3.7990725144119483234774933903962100e-02
-6.8147627068387099851816558047607990e-02

*W->H0[5][] = 

-1.5860844687643677860427256687848220e-02
2.5521043385045588077720850301196520e-02
-9.7890038632020852178867245446868040e-03
-2.5311543415490180490526126967150610e-02
1.3790246716460565417659053473014720e-02
2.3602903063883250333654142083562090e-02

*W->G0[0][] = 

3.0060906601136559011985955852563580e-01
-4.4914591311274779117964658353211370e-01
1.5163059343245509653307436176679700e-01
3.3178342147754405175831051689980080e-01
-1.5167667097277605166282857596569450e-01
-2.2764940152893562438456177553204410e-01

*W->G0[1][] = 

2.0331509574796142686321962478161710e-01
-4.0651876118469017484078004121141840e-01
4.1960355578775851878869305198268360e-01
-8.5363298620527162671612626165401440e-02
-3.2625764385211260983838425269820910e-01
-6.0061482605928826792196924781685640e-02

*W->G0[2][] = 

1.0518985773665471620509964971575450e-01
-2.5560736285629154136545139643504820e-01
4.1215305640558031471676974142351750e-01
-4.4411305077513545141181465380296190e-01
6.1254975192579201233435754997261410e-02
2.2965464098051981628652034849213960e-01

*W->G0[3][] = 

-3.2967631665057676887316473532905920e-02
9.5094248664088113521139294683440500e-02
-2.1415323405224231458606533613298370e-01
4.2062818069154075805776809580420220e-01
-5.1300286210618137768201649250409940e-01
-6.2518865611637408307375062447368440e-02

*W->G0[4][] = 

-1.4026930030508929369339000186968520e-03
8.1325538312749689030973463539590570e-03
-3.7291066482571648601759069413972260e-02
1.5104813412862278397165214099674810e-01
-4.7882817380761625051337323316880210e-01
4.9643661380425784640289250793323020e-01

*W->G0[5][] = 

-2.6181866882996125101989682377595060e-03
6.8741188247199188680051014383507090e-03
-1.2537054881833811005797286888658630e-02
9.7510337134616010922889537975387640e-03
9.9953349094296217255694607150426670e-02
-6.9978779582820063480379926463523130e-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-32   -5e-32   9e-32   1e-32   4e-32   
5e-32   1e+00   1e-31   -2e-31   -9e-33   -8e-32   
-5e-32   1e-31   1e+00   1e-31   5e-32   8e-32   
9e-32   -2e-31   1e-31   1e+00   -1e-31   -4e-32   
1e-32   -9e-33   5e-32   -1e-31   1e+00   2e-32   
4e-32   -8e-32   8e-32   -4e-32   2e-32   1e+00   

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

1e+00   8e-32   -9e-32   -4e-32   -5e-32   -6e-32   
8e-32   1e+00   -1e-31   3e-31   -5e-32   2e-31   
-9e-32   -1e-31   1e+00   -2e-31   -1e-31   -2e-31   
-4e-32   3e-31   -2e-31   1e+00   2e-32   2e-31   
-5e-32   -5e-32   -1e-31   2e-32   1e+00   -6e-32   
-6e-32   2e-31   -2e-31   2e-31   -6e-32   1e+00   

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

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