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

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

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

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

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

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

*W->H0[0][] = 

8.6891442027979115044796736356183010e-01
4.7501413271979535695533886264525160e-01
5.8925565098878960366737030175404070e-02

*W->H0[1][] = 

-1.7450492011853362811745351541945030e-01
4.1247895569215272256715921122782860e-01
5.4718291636849857751898246020799630e-01

*W->H0[2][] = 

5.8925565098878960366737030175404220e-02
-1.0233613646983040755380991785670560e-01
-4.3956508895485705313648941106173150e-02

*W->G0[0][] = 

-4.1290113334480119872175070518049320e-01
5.4997194092287029675621228163710490e-01
1.6444913584482456578739807532146310e-01

*W->G0[1][] = 

-1.9774861218395140709827211664853360e-01
5.1031036307982877045776751556372760e-01
-4.4774861218395140709827211664853340e-01

*W->G0[2][] = 

3.8362498927159704012838994462577170e-02
-1.7568209223157662955549408580181730e-01
6.8385972329506251820938322775964440e-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   -3e-35   
2e-34   1e+00   -1e-34   
-3e-35   -1e-34   1e+00   

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

1e+00   -2e-34   -3e-34   
-2e-34   1e+00   -5e-35   
-3e-34   -5e-35   1e+00   

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

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