13 September 2021 8:49:03.114 AM hermite_integrands_test(): FORTRAN90 version Test hermite_integrands(). TEST01 P00_PROBLEM_NUM returns the number of problems. P00_TITLE returns the title of a problem. P00_PROBLEM_NUM: number of problems is 8 Problem Title 1 "cos(2*omega*x)" 2 "exp(-x*x)" 3 "exp(-px) / ( 1 + exp(-qx) )" 4 "sin(x^2)" 5 "1/( (1+x^2) sqrt(4+3x^2) )" 6 "x^m exp(-x*x)" 7 "x^2 cos ( x ) exp(-x*x)" 8 "sqrt(1+x*x/2) * exp(-x*x/2)" TEST02 P00_EXACT returns the "exact" integral. Problem EXACT 1 0.6520493321732922 2 1.772453850905516 3 1.209199576156145 4 1.253314137315500 5 1.047197551196598 6 1.329340388179137 7 0.3450971117607857 8 3.008823566113644 TEST03 P00_GAUSS_HERMITE applies a Gauss-Hermite rule to estimate an integral on (-oo,+oo). Problem Order Estimate Exact Error 1 1 1.77245 0.652049 1.12040 1 2 0.276403 0.652049 0.375646 1 4 0.641433 0.652049 0.106160E-01 1 8 0.652049 0.652049 0.542791E-06 1 16 0.652049 0.652049 0.795031E-12 1 32 0.652049 0.652049 0.205036E-11 1 64 0.652049 0.652049 0.118872E-11 2 1 1.77245 1.77245 0.551581E-11 2 2 1.77245 1.77245 0.551625E-11 2 4 1.77245 1.77245 0.551603E-11 2 8 1.77245 1.77245 0.551648E-11 2 16 1.77245 1.77245 0.693734E-11 2 32 1.77245 1.77245 0.553957E-11 2 64 1.77245 1.77245 0.676947E-11 3 1 0.886227 1.20920 0.322973 3 2 0.960529 1.20920 0.248670 3 4 1.10876 1.20920 0.100440 3 8 1.18010 1.20920 0.291040E-01 3 16 1.20387 1.20920 0.532531E-02 3 32 1.20870 1.20920 0.497008E-03 3 64 1.20918 1.20920 0.175953E-04 4 1 0.00000 1.25331 1.25331 4 2 1.40102 1.25331 0.147703 4 4 1.58047 1.25331 0.327152 4 8 2.12994 1.25331 0.876630 4 16 1.55094 1.25331 0.297630 4 32 2.29436 1.25331 1.04105 4 64 -0.498965 1.25331 1.75228 5 1 0.886227 1.04720 0.160971 5 2 0.830710 1.04720 0.216487 5 4 0.947590 1.04720 0.996079E-01 5 8 1.00288 1.04720 0.443171E-01 5 16 1.02654 1.04720 0.206547E-01 5 32 1.03729 1.04720 0.990572E-02 5 64 1.04239 1.04720 0.480589E-02 6 1 0.00000 1.32934 1.32934 6 2 0.443113 1.32934 0.886227 6 4 1.32934 1.32934 0.413714E-11 6 8 1.32934 1.32934 0.413825E-11 6 16 1.32934 1.32934 0.118137E-10 6 32 1.32934 1.32934 0.431322E-11 6 64 1.32934 1.32934 0.666533E-11 7 1 0.00000 0.345097 0.345097 7 2 0.673749 0.345097 0.328652 7 4 0.348155 0.345097 0.305786E-02 7 8 0.345097 0.345097 0.274859E-08 7 16 0.345097 0.345097 0.121647E-11 7 32 0.345097 0.345097 0.105482E-11 7 64 0.345097 0.345097 0.169864E-11 8 1 1.77245 3.00882 1.23637 8 2 2.54451 3.00882 0.464317 8 4 2.94657 3.00882 0.622519E-01 8 8 3.00785 3.00882 0.972964E-03 8 16 3.00882 3.00882 0.170381E-06 8 32 3.00882 3.00882 0.127072E-10 8 64 3.00882 3.00882 0.126792E-10 TEST04 P00_TURING applies a Turing procedure to estimate an integral on (-oo,+oo). Using a tolerance of TOL = 0.100000E-03 Problem H N Estimate Exact Error 1 1.000000 9 0.670109 0.652049 0.180601E-01 1 0.500000 15 0.652049 0.652049 0.178494E-07 1 0.250000 27 0.652047 0.652049 0.192583E-05 1 0.125000 53 0.652046 0.652049 0.376749E-05 1 0.062500 103 0.652041 0.652049 0.816511E-05 1 0.031250 203 0.652038 0.652049 0.116531E-04 1 0.015625 303 0.651546 0.652049 0.503673E-03 2 1.000000 9 1.77264 1.77245 0.183354E-03 2 0.500000 15 1.77245 1.77245 0.114154E-06 2 0.250000 27 1.77245 1.77245 0.284744E-05 2 0.125000 53 1.77245 1.77245 0.482675E-05 2 0.062500 103 1.77244 1.77245 0.935244E-05 2 0.031250 203 1.77244 1.77245 0.128587E-04 2 0.015625 405 1.77244 1.77245 0.135593E-04 3 1.000000 21 1.21422 1.20920 0.501647E-02 3 0.500000 39 1.20915 1.20920 0.507025E-04 3 0.250000 75 1.20911 1.20920 0.846013E-04 3 0.125000 149 1.20911 1.20920 0.902339E-04 3 0.062500 297 1.20911 1.20920 0.931438E-04 3 0.031250 591 1.20910 1.20920 0.976260E-04 3 0.015625 1181 1.20910 1.20920 0.983947E-04 4 1.000000 67199 34.7886 1.25331 33.5353 4 0.500000 72231 192.753 1.25331 191.500 4 0.250000 144461 71.5169 1.25331 70.2635 4 0.125000 101939 32.4151 1.25331 31.1618 P00_TURING - Warning! Number of steps exceeded N_TOO_MANY = 100000 4 0.062500 200001 7.77694 1.25331 6.52362 P00_TURING - Warning! Number of steps exceeded N_TOO_MANY = 100000 4 0.031250 200001 2.38535 1.25331 1.13204 4 0.015625 29071 3.99397 1.25331 2.74066 5 1.000000 47 1.05625 1.04720 0.904800E-02 5 0.500000 93 1.04615 1.04720 0.104531E-02 5 0.250000 183 1.04610 1.04720 0.110157E-02 5 0.125000 363 1.04608 1.04720 0.111985E-02 5 0.062500 725 1.04607 1.04720 0.112294E-02 5 0.031250 1447 1.04607 1.04720 0.112760E-02 5 0.015625 2891 1.04607 1.04720 0.112994E-02 6 1.000000 9 1.34191 1.32934 0.125689E-01 6 0.500000 17 1.32934 1.32934 0.666991E-06 6 0.250000 33 1.32934 1.32934 0.270826E-05 6 0.125000 65 1.32934 1.32934 0.481629E-05 6 0.062500 3 0.189991E-05 1.32934 1.32934 6 0.031250 3 0.595465E-07 1.32934 1.32934 6 0.015625 3 0.186219E-08 1.32934 1.32934 7 1.000000 9 0.334355 0.345097 0.107424E-01 7 0.500000 17 0.345097 0.345097 0.675187E-08 7 0.250000 31 0.345098 0.345097 0.649404E-06 7 0.125000 59 0.345101 0.345097 0.352199E-05 7 0.062500 115 0.345105 0.345097 0.748815E-05 7 0.031250 227 0.345108 0.345097 0.107085E-04 7 0.015625 451 0.345110 0.345097 0.127481E-04 8 1.000000 11 3.00865 3.00882 0.171790E-03 8 0.500000 21 3.00882 3.00882 0.115982E-05 8 0.250000 39 3.00881 3.00882 0.955875E-05 8 0.125000 77 3.00881 3.00882 0.135440E-04 8 0.062500 153 3.00881 3.00882 0.159144E-04 8 0.031250 305 3.00881 3.00882 0.171969E-04 8 0.015625 607 3.00880 3.00882 0.192433E-04 Using a tolerance of TOL = 0.100000E-06 Problem H N Estimate Exact Error 1 1.000000 9 0.670109 0.652049 0.180601E-01 1 0.500000 17 0.652049 0.652049 0.147556E-08 1 0.250000 33 0.652049 0.652049 0.512318E-08 1 0.125000 63 0.652049 0.652049 0.116741E-07 1 0.062500 125 0.652049 0.652049 0.112942E-07 1 0.031250 249 0.652049 0.652049 0.107237E-07 1 0.015625 495 0.652049 0.652049 0.936090E-08 2 1.000000 11 1.77264 1.77245 0.183354E-03 2 0.500000 19 1.77245 1.77245 0.139608E-10 2 0.250000 35 1.77245 1.77245 0.889602E-09 2 0.125000 67 1.77245 1.77245 0.537981E-08 2 0.062500 133 1.77245 1.77245 0.728341E-08 2 0.031250 265 1.77245 1.77245 0.839373E-08 2 0.015625 527 1.77245 1.77245 0.102632E-07 3 1.000000 35 1.21424 1.20920 0.504286E-02 3 0.500000 67 1.20921 1.20920 0.693830E-05 3 0.250000 131 1.20920 1.20920 0.771297E-07 3 0.125000 259 1.20920 1.20920 0.932342E-07 3 0.062500 517 1.20920 1.20920 0.962408E-07 3 0.031250 1033 1.20920 1.20920 0.977683E-07 3 0.015625 2065 1.20920 1.20920 0.985381E-07 P00_TURING - Warning! Number of steps exceeded N_TOO_MANY = 100000 4 1.000000 200001 301.392 1.25331 300.139 P00_TURING - Warning! Number of steps exceeded N_TOO_MANY = 100000 4 0.500000 200001 120.152 1.25331 118.899 P00_TURING - Warning! Number of steps exceeded N_TOO_MANY = 100000 4 0.250000 200001 27.1654 1.25331 25.9121 P00_TURING - Warning! Number of steps exceeded N_TOO_MANY = 100000 4 0.125000 200001 23.0552 1.25331 21.8018 P00_TURING - Warning! Number of steps exceeded N_TOO_MANY = 100000 4 0.062500 200001 7.77694 1.25331 6.52362 P00_TURING - Warning! Number of steps exceeded N_TOO_MANY = 100000 4 0.031250 200001 2.38535 1.25331 1.13204 P00_TURING - Warning! Number of steps exceeded N_TOO_MANY = 100000 4 0.015625 200001 7.81267 1.25331 6.55936 5 1.000000 455 1.05728 1.04720 0.100802E-01 5 0.500000 907 1.04721 1.04720 0.974568E-05 5 0.250000 1811 1.04719 1.04720 0.112661E-04 5 0.125000 3619 1.04719 1.04720 0.112848E-04 5 0.062500 7235 1.04719 1.04720 0.112942E-04 5 0.031250 14467 1.04719 1.04720 0.112989E-04 5 0.015625 28933 1.04719 1.04720 0.112997E-04 6 1.000000 11 1.34191 1.32934 0.125690E-01 6 0.500000 21 1.32934 1.32934 0.669520E-10 6 0.250000 41 1.32934 1.32934 0.442681E-09 6 0.125000 79 1.32934 1.32934 0.313102E-08 6 0.062500 155 1.32934 1.32934 0.772569E-08 6 0.031250 309 1.32934 1.32934 0.899104E-08 6 0.015625 617 1.32934 1.32934 0.967355E-08 7 1.000000 11 0.334355 0.345097 0.107424E-01 7 0.500000 19 0.345097 0.345097 0.100232E-09 7 0.250000 37 0.345097 0.345097 0.125071E-09 7 0.125000 71 0.345097 0.345097 0.183235E-08 7 0.062500 141 0.345097 0.345097 0.310895E-08 7 0.031250 279 0.345097 0.345097 0.583390E-08 7 0.015625 557 0.345097 0.345097 0.645572E-08 8 1.000000 15 3.00865 3.00882 0.171657E-03 8 0.500000 27 3.00882 3.00882 0.624212E-08 8 0.250000 51 3.00882 3.00882 0.195499E-08 8 0.125000 99 3.00882 3.00882 0.687079E-08 8 0.062500 195 3.00882 3.00882 0.124539E-07 8 0.031250 389 3.00882 3.00882 0.137587E-07 8 0.015625 777 3.00882 3.00882 0.144441E-07 TEST05 P00_GAUSS_HERMITE applies a Gauss-Hermite rule to estimate the integral x^m exp(-x*x) over (-oo,+oo). M Order Estimate Exact Error 0 1 1.77245 1.77245 0.551581E-11 0 2 1.77245 1.77245 0.551625E-11 0 3 1.77245 1.77245 0.551537E-11 1 1 0.00000 0.00000 0.00000 1 2 0.00000 0.00000 0.00000 1 3 0.00000 0.00000 0.00000 2 1 0.00000 0.886227 0.886227 2 2 0.886227 0.886227 0.275835E-11 2 3 0.886227 0.886227 0.275768E-11 2 4 0.886227 0.886227 0.275813E-11 3 1 0.00000 0.00000 0.00000 3 2 0.00000 0.00000 0.00000 3 3 0.00000 0.00000 0.00000 3 4 -0.555112E-16 0.00000 0.555112E-16 4 1 0.00000 1.32934 1.32934 4 2 0.443113 1.32934 0.886227 4 3 1.32934 1.32934 0.413669E-11 4 4 1.32934 1.32934 0.413714E-11 4 5 1.32934 1.32934 0.413669E-11 5 1 0.00000 0.00000 0.00000 5 2 0.00000 0.00000 0.00000 5 3 0.00000 0.00000 0.00000 5 4 -0.111022E-15 0.00000 0.111022E-15 5 5 0.111022E-15 0.00000 0.111022E-15 6 1 0.00000 3.32335 3.32335 6 2 0.221557 3.32335 3.10179 6 3 1.99401 3.32335 1.32934 6 4 3.32335 3.32335 0.103428E-10 6 5 3.32335 3.32335 0.103419E-10 6 6 3.32335 3.32335 0.297553E-10 TEST06 P00_MONTE_CARLO applies a Monte Carlo rule to estimate Hermite-weighted integrals on (-oo,+oo). Problem Order Estimate Exact Error 1 128 0.578806 0.652049 0.732430E-01 1 512 0.629587 0.652049 0.224621E-01 1 2048 0.621772 0.652049 0.302771E-01 1 8192 0.645132 0.652049 0.691781E-02 1 32768 0.650156 0.652049 0.189337E-02 1 131072 0.650896 0.652049 0.115379E-02 1 524288 0.652398 0.652049 0.348506E-03 2 128 1.86985 1.77245 0.973932E-01 2 512 1.72293 1.77245 0.495271E-01 2 2048 1.76568 1.77245 0.677186E-02 2 8192 1.77084 1.77245 0.161031E-02 2 32768 1.77141 1.77245 0.104023E-02 2 131072 1.77153 1.77245 0.921525E-03 2 524288 1.77264 1.77245 0.187075E-03 3 128 1.13731 1.20920 0.718858E-01 3 512 1.41220 1.20920 0.202998 3 2048 1.14319 1.20920 0.660053E-01 3 8192 1.19296 1.20920 0.162433E-01 3 32768 1.20190 1.20920 0.729819E-02 3 131072 1.20270 1.20920 0.649628E-02 3 524288 1.21234 1.20920 0.314000E-02 4 128 3.50536 1.25331 2.25204 4 512 -0.555570 1.25331 1.80888 4 2048 1.53671 1.25331 0.283396 4 8192 1.31954 1.25331 0.662222E-01 4 32768 1.12858 1.25331 0.124736 4 131072 1.57008 1.25331 0.316764 4 524288 0.471595 1.25331 0.781720 5 128 0.984539 1.04720 0.626583E-01 5 512 1.01649 1.04720 0.307070E-01 5 2048 1.02121 1.04720 0.259833E-01 5 8192 1.04027 1.04720 0.693074E-02 5 32768 1.02265 1.04720 0.245509E-01 5 131072 1.01983 1.04720 0.273660E-01 5 524288 1.04247 1.04720 0.472923E-02 6 128 1.63083 1.32934 0.301492 6 512 1.45835 1.32934 0.129013 6 2048 1.32453 1.32934 0.481049E-02 6 8192 1.35058 1.32934 0.212364E-01 6 32768 1.31678 1.32934 0.125580E-01 6 131072 1.34482 1.32934 0.154819E-01 6 524288 1.32846 1.32934 0.883695E-03 7 128 0.314339 0.345097 0.307583E-01 7 512 0.366894 0.345097 0.217973E-01 7 2048 0.359545 0.345097 0.144475E-01 7 8192 0.344629 0.345097 0.467828E-03 7 32768 0.342128 0.345097 0.296926E-02 7 131072 0.343708 0.345097 0.138924E-02 7 524288 0.345581 0.345097 0.483502E-03 8 128 2.96992 3.00882 0.389019E-01 8 512 3.03581 3.00882 0.269815E-01 8 2048 3.01708 3.00882 0.825500E-02 8 8192 3.00619 3.00882 0.263578E-02 8 32768 3.00975 3.00882 0.925143E-03 8 131072 3.00764 3.00882 0.118313E-02 8 524288 3.00990 3.00882 0.107492E-02 hermite_integrands_TEST Normal end of execution. 13 September 2021 8:49:03.655 AM