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

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

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

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

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

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

*W->H0[0][] = 

8.8822997760250100729054769937233010e-01
4.4803410146378722839208248861480170e-01
-3.5276317103128201410078360071426240e-02
-7.1540501803581243622535916473434750e-02
-1.0427280445542713862331418248616630e-02

*W->H0[1][] = 

-2.0024209800316952263033644162045020e-01
5.2978782155016172831603279357868920e-01
7.4216746870939184074967286375443670e-01
3.1289885471493376449735248003811190e-01
3.6319464148597061446266064920296890e-02

*W->H0[2][] = 

1.0337789558065772684796960868142920e-01
-1.8861186417156541187634768870500590e-01
6.7764399863710804421747584701714700e-02
4.9513287816831384267095457661356160e-01
4.5165251517833306510941866837194730e-01

*W->H0[3][] = 

-5.5662892598648435576735697838252510e-02
9.3489560822740412890141586803293060e-02
-2.7493032510730527159613030467388360e-02
-1.2516365150336655137056062228887480e-01
-3.0498456787002271676542886063885440e-02

*W->H0[4][] = 

2.2410377109909155078670265652580970e-02
-3.6621454280123461743417326756748600e-02
1.0189684546835032004751492694678920e-02
4.2499403413370174936300502899634840e-02
9.4008410359288398688684962142866200e-03

*W->G0[0][] = 

-3.2694151367399630039879221785223160e-01
4.8143355124137229544644705774080170e-01
-1.1242945484970418836324291888690560e-01
-3.7914109420005794267390216567349900e-01
-7.0303845888103186131675963606494690e-02

*W->G0[1][] = 

-2.0404257253782667421575877389826250e-01
4.2649511439457577428927382556428620e-01
-4.5421997916613510686054049126496490e-01
2.6874662409172835444432315006654530e-02
2.6349669644307164650754680228402660e-01

*W->G0[2][] = 

-8.8488193159571765757048733768284240e-02
2.3575288549718986467274674189881860e-01
-4.3599101942396084199501423821824620e-01
4.8838509116284323108824522049968810e-01
8.9343440512551942066405493321630130e-02

*W->G0[3][] = 

1.7407738929280795035553062458448850e-02
-6.0497531116082337052515861280472970e-02
1.8275051899511108129965687174572520e-01
-4.8587780307074688967393169760258070e-01
4.7598532292959485700291992129333050e-01

*W->G0[4][] = 

2.3922147966636212877208836748035440e-03
-4.3269741464359883162004604791587420e-03
-7.0161557587830761216368895697935190e-03
1.2361271678338464068523891140112300e-01
-6.9616537145343701045845503200715130e-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-33   -1e-32   -5e-33   -2e-32   
2e-33   1e+00   1e-32   1e-32   3e-32   
-1e-32   1e-32   1e+00   -1e-32   -3e-32   
-5e-33   1e-32   -1e-32   1e+00   1e-32   
-2e-32   3e-32   -3e-32   1e-32   1e+00   

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

1e+00   4e-32   -5e-32   9e-32   -2e-33   
4e-32   1e+00   -8e-32   -7e-32   -6e-32   
-5e-32   -8e-32   1e+00   6e-32   7e-32   
9e-32   -7e-32   6e-32   1e+00   -2e-32   
-2e-33   -6e-32   7e-32   -2e-32   1e+00   

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

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