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

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

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

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

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

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

*W->H0[0][] = 

8.8180459710597623300949110496160680e-01
4.5830053603093937793641819230668070e-01
-7.1466178063564077259940208048995680e-03
-6.7606192722111045266445529162741230e-02

*W->H0[1][] = 

-1.9177464443924776540187693744706570e-01
4.9105798874070495197491544753781000e-01
7.0417329623313415224114737876400010e-01
3.6597162457362714747303230668678440e-01

*W->H0[2][] = 

8.9142265022079756717406767160571930e-02
-1.6045381492838664376761139895770190e-01
1.0480312529707648601740083611219080e-02
2.4184224038446340515694304653892220e-01

*W->H0[3][] = 

-3.4171550341044260456562724328026400e-02
5.7214213751572047481060129896756160e-02
-3.1866464352324737459791184593243780e-03
-5.7517683651612225334563457164044190e-02

*W->G0[0][] = 

-3.6200849931389677878707463713657130e-01
5.1796785638459995001816650331024080e-01
-2.2443289316917131627372226678683820e-02
-3.1647345061075816820532620106682330e-01

*W->G0[1][] = 

-2.0253177774962747228976183350225350e-01
4.5462311379313475241267902293462000e-01
-4.8695150854106675403040020955714240e-01
-1.2319469037336529815897047631690060e-01

*W->G0[2][] = 

-6.6650227217054207960060470711570700e-02
2.1019517420630039668474569218614740e-01
-4.9404252886853518942344359519493740e-01
4.5529958886592534332382144253475960e-01

*W->G0[3][] = 

4.1919738764389198603920980638792430e-03
-2.5473324972679092599695210395168880e-02
1.4921598380625366883970594046165240e-01
-6.9070118520716857115439587151396680e-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-33   -2e-33   
-2e-33   1e+00   3e-33   2e-33   
-1e-33   3e-33   1e+00   -4e-34   
-2e-33   2e-33   -4e-34   1e+00   

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

1e+00   -1e-33   -7e-33   -6e-33   
-1e-33   1e+00   9e-33   1e-32   
-7e-33   9e-33   1e+00   -6e-33   
-6e-33   1e-32   -6e-33   1e+00   

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

-2e-34   2e-33   -3e-33   3e-33   
-6e-33   6e-33   -3e-33   8e-34   
-6e-33   4e-33   9e-34   -2e-33   
-3e-33   -2e-33   5e-33   -2e-33   
The size of double is: 8 bytes.
The size of long double is: 16 bytes.
