29 May 2014 3:09:38.988 PM QWV_2D_PRB: FORTRAN77 version Test the QWV_2D library. TEST01: Compute the weights associated with an interpolatory quadrature rule defined by N=(T+1)*(T+2)/2 points, exact for polynomials of total degree T or less. Degree T = 1 Number of points N = 3 X Interval = [ -1.0000, 1.0000] Y Interval = [ -1.0000, 1.0000] Abscissas: 1: 1.00000 -1.00000 2: 1.00000 1.00000 3: -1.00000 1.00000 Weights: 1: 2.0000000 2: 0.0000000 3: 2.0000000 TEST02: Compute the weights associated with an interpolatory quadrature rule defined by N=(T+1)*(T+2)/2 points, exact for polynomials of total degree T or less. X Interval = [ -1.0000, 1.0000] Y Interval = [ -1.0000, 1.0000] Degree T = 0 Number of points N = 1 Weights: 1: 4.000000000000000 Degree T = 1 Number of points N = 3 Weights: 1: 1.000000000000000 2: 1.000000000000000 3: 2.000000000000000 Degree T = 2 Number of points N = 6 Weights: 1: 0.6661338147750939E-15 2: 0.6666666666666663 3: 2.222222222222222 4: 0.4444444444444444 5: 0.5551115123125783E-16 6: 0.6666666666666664 Degree T = 3 Number of points N = 10 Weights: 1: -0.5555555555555480E-01 2: 0.3333333333333331 3: -0.5555555555555580E-01 4: 0.8888888888888886 5: 0.8888888888888893 6: 0.2222222222222224 7: 1.333333333333333 8: 0.2222222222222220 9: 0.1111111111111109 10: 0.1111111111111112 Degree T = 4 Number of points N = 15 Weights: 1: -0.8888888888888932E-02 2: 0.8104919101110961E-01 3: 0.6117303121111219E-01 4: 0.3874097078666789 5: 0.6259236254666545 6: 0.5333333333333362E-01 7: 0.7111111111111067E-01 8: 0.9830822022444241 9: 0.5458066866444642 10: 0.3874097078666780 11: 0.6259236254666568 12: 0.5333333333333383E-01 13: -0.8888888888888703E-02 14: 0.8104919101110968E-01 15: 0.6117303121111135E-01 Degree T = 5 Number of points N = 21 Weights: 1: -0.1037037037037093E-01 2: 0.5037037037036911E-01 3: 0.5037037037037081E-01 4: -0.1037037037036947E-01 5: 0.1876963678740801 6: 0.3460933466518654 7: 0.1876963678740763 8: 0.4514390511851724E-01 9: 0.5541130536814713 10: 0.5541130536814728 11: 0.4514390511851834E-01 12: 0.2804517802740705 13: 0.6376103570518378 14: 0.2804517802740683 15: 0.3189313191851883E-01 16: 0.3288499092814910 17: 0.3288499092814925 18: 0.3189313191851956E-01 19: 0.2074074074074123E-01 20: 0.3851851851851849E-01 21: 0.2074074074074051E-01 Degree T = 6 Number of points N = 28 Weights: 1: -0.3023431594858565E-02 2: 0.1957267632451884E-01 3: 0.2633929313290840E-01 4: 0.1425431928029237E-01 5: 0.1006383046329639 6: 0.2208900184526934 7: 0.1743144584714012 8: 0.1209372637943976E-01 9: 0.1934996220710680E-01 10: 0.3245064820875231 11: 0.4027058473592984 12: 0.1677234226317961 13: 0.1953319357827178 14: 0.4489633053035124 15: 0.3721824611057551 16: 0.2479213907785274E-01 17: 0.1934996220710561E-01 18: 0.3245064820875153 19: 0.4027058473592959 20: 0.1677234226317944 21: 0.1006383046329745 22: 0.2208900184526933 23: 0.1743144584714027 24: 0.1209372637944051E-01 25: -0.3023431594861990E-02 26: 0.1957267632451757E-01 27: 0.2633929313290797E-01 28: 0.1425431928029198E-01 Degree T = 7 Number of points N = 36 Weights: 1: -0.3287981859413765E-02 2: 0.1337868480725671E-01 3: 0.2063492063491996E-01 4: 0.1337868480725546E-01 5: -0.3287981859408898E-02 6: 0.5949324721885513E-01 7: 0.1306477599993571 8: 0.1306477599993581 9: 0.5949324721885061E-01 10: 0.1263869091685831E-01 11: 0.1979944935601103 12: 0.2832184784823740 13: 0.1979944935601143 14: 0.1263869091685747E-01 15: 0.1221817987389771 16: 0.3150266070593529 17: 0.3150266070593440 18: 0.1221817987389802 19: 0.1771365352315134E-01 20: 0.2490926964598258 21: 0.3408041116306980 22: 0.2490926964598291 23: 0.1771365352314976E-01 24: 0.9646986307476696E-01 25: 0.2557725606433917 26: 0.2557725606433927 27: 0.9646986307476431E-01 28: 0.8649923133686802E-02 29: 0.1062007918394705 30: 0.1505805844901012 31: 0.1062007918394705 32: 0.8649923133690016E-02 33: 0.6355881462931014E-02 34: 0.1405228180237514E-01 35: 0.1405228180237651E-01 36: 0.6355881462928496E-02 Degree T = 8 Number of points N = 45 Weights: 1: -0.1269841269835311E-02 2: 0.6706089639041270E-02 3: 0.1111455441352989E-01 4: 0.1026455026455282E-01 5: 0.4930678698742625E-02 6: 0.3633146869162523E-01 7: 0.8838322767333079E-01 8: 0.9965911758463214E-01 9: 0.6400185533755555E-01 10: 0.4061629144893127E-02 11: 0.6772486772485166E-02 12: 0.1258344472781388 13: 0.1927501398511116 14: 0.1699470899470907 15: 0.6342599488133535E-01 16: 0.8376332474107638E-01 17: 0.2170841444607031 18: 0.2477307250801775 19: 0.1648098048612226 20: 0.1004771829779292E-01 21: 0.1015873015872910E-01 22: 0.1784328991205164 23: 0.2729409493576765 24: 0.2364021164021134 25: 0.8936689226256009E-01 26: 0.8376332474107701E-01 27: 0.2170841444607054 28: 0.2477307250801761 29: 0.1648098048612200 30: 0.1004771829779330E-01 31: 0.6772486772485237E-02 32: 0.1258344472781358 33: 0.1927501398511135 34: 0.1699470899470926 35: 0.6342599488133838E-01 36: 0.3633146869162453E-01 37: 0.8838322767332588E-01 38: 0.9965911758463601E-01 39: 0.6400185533755502E-01 40: 0.4061629144888279E-02 41: -0.1269841269836355E-02 42: 0.6706089639046927E-02 43: 0.1111455441352761E-01 44: 0.1026455026454956E-01 45: 0.4930678698747173E-02 Degree T = 9 Number of points N = 55 Weights: 1: -0.1368606701945113E-02 2: 0.4837977417140975E-02 3: 0.8876308297144902E-02 4: 0.8876308297143068E-02 5: 0.4837977417150492E-02 6: -0.1368606701935084E-02 7: 0.2425285860992349E-01 8: 0.5727330842923516E-01 9: 0.7008257906578071E-01 10: 0.5727330842922034E-01 11: 0.2425285860989794E-01 12: 0.4659404339099723E-02 13: 0.8354521980498550E-01 14: 0.1370796991940044 15: 0.1370796991940248 16: 0.8354521980500107E-01 17: 0.4659404339109654E-02 18: 0.5564545640233619E-01 19: 0.1524391996823315 20: 0.1877107583774149 21: 0.1524391996823176 22: 0.5564545640232402E-01 23: 0.8186176158691754E-02 24: 0.1295355639606716 25: 0.2061407656847711 26: 0.2061407656847630 27: 0.1295355639606894 28: 0.8186176158692687E-02 29: 0.6234969028097752E-01 30: 0.1730419031522391 31: 0.2169418247419051 32: 0.1730419031522361 33: 0.6234969028097048E-01 34: 0.7506172839505762E-02 35: 0.1142161960569350 36: 0.1802176663769002 37: 0.1802176663769038 38: 0.1142161960569279 39: 0.7506172839512260E-02 40: 0.4031900987631698E-01 41: 0.1142976211857364 42: 0.1413353845521477 43: 0.1142976211857414 44: 0.4031900987631700E-01 45: 0.3239075586856897E-02 46: 0.4317587564913915E-01 47: 0.7015250533601934E-01 48: 0.7015250533601930E-01 49: 0.4317587564913908E-01 50: 0.3239075586852207E-02 51: 0.2550690557469151E-02 52: 0.6084230077461027E-02 53: 0.7421516754852508E-02 54: 0.6084230077458821E-02 55: 0.2550690557473353E-02 Degree T = 10 Number of points N = 66 Weights: 1: -0.6240762604463766E-03 2: 0.2843227149025789E-02 3: 0.5250031948150784E-02 4: 0.5891746241568810E-02 5: 0.4705736485964679E-02 6: 0.2135354637732944E-02 7: 0.1610939653924566E-01 8: 0.4099595211758227E-01 9: 0.5326500934654063E-01 10: 0.4863338516658277E-01 11: 0.2843474741781434E-01 12: 0.1719619179693151E-02 13: 0.2883769745121509E-02 14: 0.5724711668876453E-01 15: 0.9659872841640438E-01 16: 0.1053210323353631 17: 0.8066212502628711E-01 18: 0.2855765663647366E-01 19: 0.3981286043310814E-01 20: 0.1090390674981577 21: 0.1430169021081585 22: 0.1313686303763064 23: 0.7932850918298831E-01 24: 0.4610696968783255E-02 25: 0.5086495679684716E-02 26: 0.9311356395361167E-01 27: 0.1562320334111262 28: 0.1696057154254139 29: 0.1283581371975154 30: 0.4603059518094556E-01 31: 0.4894888812994630E-01 32: 0.1347281473526573 33: 0.1764193542601264 34: 0.1635037456303485 35: 0.9822749154565460E-01 36: 0.5704840613923174E-02 37: 0.5086495679679268E-02 38: 0.9311356395362781E-01 39: 0.1562320334111511 40: 0.1696057154253968 41: 0.1283581371975113 42: 0.4603059518094044E-01 43: 0.3981286043311782E-01 44: 0.1090390674981293 45: 0.1430169021081508 46: 0.1313686303763217 47: 0.7932850918299997E-01 48: 0.4610696968790496E-02 49: 0.2883769745110260E-02 50: 0.5724711668875122E-01 51: 0.9659872841642343E-01 52: 0.1053210323353932 53: 0.8066212502626474E-01 54: 0.2855765663644533E-01 55: 0.1610939653928420E-01 56: 0.4099595211758404E-01 57: 0.5326500934649123E-01 58: 0.4863338516656233E-01 59: 0.2843474741784810E-01 60: 0.1719619179720036E-02 61: -0.6240762604606350E-03 62: 0.2843227149011163E-02 63: 0.5250031948172295E-02 64: 0.5891746241587802E-02 65: 0.4705736485965663E-02 66: 0.2135354637703863E-02 QWV_2D_PRB: Normal end of execution. 29 May 2014 3:09:38.992 PM