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

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

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

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

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

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

*W->H0[0][] = 

8.9418319374422234776698144015892960e-01
4.3762643576669227288667185863146750e-01
-5.5655620266007608796782812697294740e-02
-5.0584278194349521972919862945223460e-02
3.4089746972711339293483649364505210e-02
2.4557618730022902350929511210736430e-02
3.3950286022255865860012492704061690e-03

*W->H0[1][] = 

-2.0800560607032045104794459768427860e-01
5.6513011021657367144204783740396130e-01
7.5598585817786550009377960964231240e-01
2.0036399312160585044548312172890950e-01
-1.0653485590654771151699495542940940e-01
-7.0048325477367611912717711056721510e-02
-9.3286418056916648551238439217739550e-03

*W->H0[2][] = 

1.1610002486514829183258318329282360e-01
-2.1399164831238919067595378118514140e-01
1.3383726406565543209328043529961670e-01
7.0016285236830713771818835664556180e-01
5.9301696632030100022267427803616950e-01
2.1122219448237929016211490394410530e-01
2.3921722827926783377697558815408570e-02

*W->H0[3][] = 

-7.4130650234003775944523368732868820e-02
1.2501313287251263592702058705986890e-01
-6.1890022670319169150205356381158030e-02
-1.4426241026886259761213655066156950e-01
2.3756308453060147054827422534526850e-01
4.9466578701846226565754746410019530e-01
3.9128491589157979204731530611055990e-01

*W->H0[4][] = 

4.7831895874561059336923793312444630e-02
-7.7932206076028619473993277559851720e-02
3.5859656368881885423872709343779520e-02
7.2775934469226233781278774156954880e-02
-9.1761325883003117325863585905081170e-02
-1.0919185281118470196051523310963570e-01
-2.1098974697189851358115625124958600e-02

*W->H0[5][] = 

-2.8244544190270326480610637521914160e-02
4.5298269998032899633931999442357280e-02
-2.0206644116017817548399269684666100e-02
-3.8979409757416872383492773887151160e-02
4.5589663867721537109811598127050940e-02
4.9255924185223636558538894742822610e-02
8.7096481272843898481053222457592240e-03

*W->H0[6][] = 

1.1828183066500455615191331414773910e-02
-1.8839040171993420640045943894769260e-02
8.2937922218398584083942891765426850e-03
1.5677093047198944068601576859285670e-02
-1.7836475604155029081144857106139530e-02
-1.8671381602335310636632841996488260e-02
-3.2177371427541291899229399868697450e-03

*W->G0[0][] = 

-2.7980958149128311507684655858780660e-01
4.2161473687309451608329802668589680e-01
-1.6849499965596451658558503439302040e-01
-2.7937148209865063821019756951533930e-01
2.7334270568812110200554837501661150e-01
2.4759901493906275310885640329692390e-01
3.8563128586087252527752030710383210e-02

*W->G0[1][] = 

-2.0095672361162168282199722577716270e-01
3.8947971242443413700203942382178130e-01
-3.8874621195558146393361134058048600e-01
9.9763751545583827808233781868147120e-02
2.9875220923973958954716044236365130e-01
-1.3417759282356529219878822908931420e-01
-1.9896988835310709865812737762608090e-01

*W->G0[2][] = 

1.1774062076860076314513030025609000e-01
-2.6999546158174786668424559480429540e-01
3.9773568407258437011888719494356860e-01
-4.0083827574611893191600583126729660e-01
1.0948789078845846848824706451652470e-01
2.8197030210889325506215595078682590e-01
5.3663299915204512276821848471888880e-02

*W->G0[3][] = 

-4.8399101850027568495550558833538720e-02
1.2662897628871086342993340894531240e-01
-2.4251365878328103423564054309275710e-01
3.9584946420689195570063700345920790e-01
-4.6078093169952135749624975598846970e-01
1.0260167043049393841254780520676810e-01
2.0800669821387964526049285636581760e-01

*W->G0[4][] = 

-8.8374035761620038434875379791642880e-03
2.7775727275402984647677146174465260e-02
-7.2358493934097014894643782555711150e-02
1.8242418448800197675545765180123460e-01
-4.0659568772468516649839945732792790e-01
5.4196169994477976298800648702603890e-01
3.9919131484614610025656174887935290e-02

*W->G0[5][] = 

-2.5213782470983222389275918310928340e-03
5.4551959220607197645494077148654680e-03
-5.2453202513487294658462585033390220e-03
-1.3777546511693239049861376362130600e-02
1.1989153497342792807830065014216150e-01
-4.7046407809956366705622990679554090e-01
5.1384445153802028591044592783954310e-01

*W->G0[6][] = 

1.7815130167586809085822150344606990e-03
-4.7907228802318516707174671857814970e-03
9.9014781450664637577606314682014980e-03
-1.8883307841974680336792600166498440e-02
2.4420441257644664279249617224360100e-02
7.7982602179102474764755314172461630e-02
-7.0202678895274480784532118284367010e-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   -6e-31   3e-31   -8e-31   -6e-32   -9e-31   -2e-31   
-6e-31   1e+00   -1e-31   1e-30   7e-31   2e-30   7e-31   
3e-31   -1e-31   1e+00   -1e-31   -3e-30   -2e-30   -2e-30   
-8e-31   1e-30   -1e-31   1e+00   1e-30   3e-30   1e-30   
-6e-32   7e-31   -3e-30   1e-30   1e+00   -9e-31   -3e-30   
-9e-31   2e-30   -2e-30   3e-30   -9e-31   1e+00   4e-31   
-2e-31   7e-31   -2e-30   1e-30   -3e-30   4e-31   1e+00   

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

1e+00   -4e-31   7e-31   -7e-32   5e-31   3e-30   3e-30   
-4e-31   1e+00   3e-31   -2e-31   -1e-30   -4e-30   -6e-30   
7e-31   3e-31   1e+00   1e-31   2e-30   4e-30   7e-30   
-7e-32   -2e-31   1e-31   1e+00   -9e-31   -6e-30   -7e-30   
5e-31   -1e-30   2e-30   -9e-31   1e+00   3e-30   7e-30   
3e-30   -4e-30   4e-30   -6e-30   3e-30   1e+00   -3e-30   
3e-30   -6e-30   7e-30   -7e-30   7e-30   -3e-30   1e+00   

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

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