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

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

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

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

*W->H0[0][] = 

7.0710678118654752440084436210484890e-01
0.0000000000000000000000000000000000e+00

*W->H0[1][] = 

-6.1237243569579452454932101867647260e-01
3.5355339059327376220042218105242440e-01

*W->G0[0][] = 

6.8091906188322241241854814069819010e-35
-7.0710678118654752440084436210484930e-01

*W->G0[1][] = 

-3.5355339059327376220042218105242460e-01
-6.1237243569579452454932101867647290e-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-34   
2e-34   1e+00   

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

1e+00   -1e-34   
-1e-34   1e+00   

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

5e-35   3e-34   
2e-35   -2e-34   
The size of double is: 8 bytes.
The size of long double is: 16 bytes.
