
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

*W->H0[0][] = 

7.0710678118654752440084436210484900e-01
1.7022976547080560310463703517454750e-35
1.0213785928248336186278222110472850e-34
0.0000000000000000000000000000000000e+00
3.1492506612099036574357851507291290e-34

*W->H0[1][] = 

-6.1237243569579452454932101867647280e-01
3.5355339059327376220042218105242450e-01
-3.8301697230931260698543332914273190e-35
1.3724774841083701750311360960947890e-34
-2.5374874415491960212784958055705990e-34

*W->H0[2][] = 

1.7022976547080560310463703517454750e-35
-6.8465319688145764182121222850100270e-01
1.7677669529663688110021109052621230e-01
-1.8725274201788616341510073869200230e-34
1.6171827719726532294940518341582010e-34

*W->H0[3][] = 

2.3385358667337133659898429576978430e-01
4.0504629365049126443537296475549980e-01
-5.2291251658379721748635751611574220e-01
8.8388347648318440550105545263106340e-02
-3.7450548403577232683020147738400450e-34

*W->H0[4][] = 

8.5114882735402801552318517587273760e-35
1.5309310892394863113733025466911820e-01
5.9292706128157112474979253958113500e-01
-3.5078038001005700489847644365467670e-01
4.4194173824159220275052772631553690e-02

*W->G0[0][] = 

-1.4664711502135328916944728550462930e-01
-2.5400025400038100063500111125199180e-01
-1.6395645894598823910561266693558530e-01
5.8198769524737792235322224715340660e-01
-2.1997067253202993375417092825695800e-01

*W->G0[1][] = 

1.2940866771090641948014507413969100e-31
7.0243935868627047439549499773730420e-02
2.7205359379125437699912985301844390e-01
4.2919753763947606516120441809722030e-01
-4.8666426339228761449343624064109210e-01

*W->G0[2][] = 

1.2201241215222784757323159651612510e-01
2.1133169700169294103900346059852280e-01
2.0890833364530163811186478450578050e-01
-5.5339628471641716520632785066290110e-02
-6.2749240535431464466233392494083350e-01

*W->G0[3][] = 

-1.9877729714025970274528466597331910e-31
5.3649692204934508145150552262029670e-02
2.0778436443886522111484327040381320e-01
4.3902459917891904649718437358542050e-01
5.1108195990158183753293870911964170e-01

*W->G0[4][] = 

-1.8417408866805574819296835120596490e-01
-3.1899887901076797087798710330489220e-01
-3.9501656636581180821642959577427120e-01
-3.8783417603566839708259476858456640e-01
-2.4055391172970546702755049953419100e-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-34   3e-33   9e-33   -4e-32   
5e-34   1e+00   2e-33   9e-33   -3e-32   
3e-33   2e-33   1e+00   2e-33   -1e-32   
9e-33   9e-33   2e-33   1e+00   1e-33   
-4e-32   -3e-32   -1e-32   1e-33   1e+00   

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

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

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

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