
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

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

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

*W->H0[0][] = 

8.3651630373780790557529378091687300e-01
4.8296291314453414337487159986444860e-01

*W->H0[1][] = 

-1.2940952255126038117444941881202420e-01
2.2414386804201338102597276224040030e-01

*W->G0[0][] = 

5.0000000000000000000000000000000020e-01
-5.0000000000000000000000000000000010e-01

*W->G0[1][] = 

1.8301270189221932338186158537646810e-01
-6.8301270189221932338186158537646830e-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   -3e-34   
-3e-34   1e+00   

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

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

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

-2e-34   1e-34   
-1e-34   1e-34   
The size of double is: 8 bytes.
The size of long double is: 16 bytes.
