special_functions/detail/hypergeometric_1F1_negative_b_regions.hpp

0001
0002 ///////////////////////////////////////////////////////////////////////////////
0003 //  Copyright 2019 John Maddock
0004 //  Distributed under the Boost
0005 //  Software License, Version 1.0. (See accompanying file
0006 //  LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
0007
0008 #ifndef BOOST_MATH_DETIAL_1F1_MAP_NEG_B_HPP
0009 #define BOOST_MATH_DETIAL_1F1_MAP_NEG_B_HPP
0010
0011 namespace boost {
0012    namespace math {
0013       namespace detail {
0014          //
0015          // hypergeometric_1F1_negative_b_recurrence_region maps out the domains over which
0016          // forward and backwards recurrences are stable when b < 0.
0017          // Returns -1, 0, or 1:
0018          // -1:  Backwards recurrence is stable.
0019          // +1:  Forwards recurrence is stable.
0020          //  0:  Neither recurrence is stable.
0021          //
0022          template <class T>
0023          int hypergeometric_1F1_negative_b_recurrence_region(const T& a, const T& b, const T& z)
0024          {
0025             BOOST_MATH_STD_USING
0026             static const double domain[][4] = {
0027                { 1e-300, -1000000.1, 278609.88983674475, 465687},
0028                { 1e-300, -400000.03999999998, 111530.83622048119, 111544},
0029                { 1e-300, -160000.016, 44699.113858309654, 44712},
0030                { 1e-300, -64000.006399999998, 17966.035866477272, 17980},
0031                { 1e-300, -25600.002560000001, 7273.793013139204, 7287},
0032                { 1e-300, -10240.001024000001, 2998.0791015625, 3012},
0033                { 1e-300, -4096.0004096000002, 1291.0461647727273, 1306},
0034                { 1e-300, -1638.40016384, 619.421875, 636},
0035                { 1e-300, -655.36006553599998, 374.62926136363637, 397},
0036                { 1e-300, -262.14402621440001, 335.29958677685943, 370},
0037                { 1e-300, -104.85761048576001, 412.83057851239664, 462},
0038                { 1e-300, -41.943044194304001, 513.00603693181824, 570},
0039                { 1e-300, -16.777217677721602, 584.14449896694214, 644},
0040                { 1e-300, -6.7108870710886404, 623.30417097107443, 684},
0041                { 1e-300, -2.6843548284354561, 642.31960227272725, 704},
0042                { 1e-300, -1.0737419313741825, 650.20738636363626, 711},
0043                { 1.0000000000000001e-275, -1000000.1, 278597.3203069513, 448872},
0044                { 1.0000000000000001e-275, -400000.03999999998, 111518.14741654831, 111531},
0045                { 1.0000000000000001e-275, -160000.016, 44686.621781782669, 44700},
0046                { 1.0000000000000001e-275, -64000.006399999998, 17953.944740988994, 17967},
0047                { 1.0000000000000001e-275, -25600.002560000001, 7261.1431107954559, 7274},
0048                { 1.0000000000000001e-275, -10240.001024000001, 2985.014559659091, 2999},
0049                { 1.0000000000000001e-275, -4096.0004096000002, 1277.4573863636363, 1292},
0050                { 1.0000000000000001e-275, -1638.40016384, 603.88778409090912, 620},
0051                { 1.0000000000000001e-275, -655.36006553599998, 354.10653409090912, 376},
0052                { 1.0000000000000001e-275, -262.14402621440001, 303.91367510330571, 337},
0053                { 1.0000000000000001e-275, -104.85761048576001, 367.0454545454545, 416},
0054                { 1.0000000000000001e-275, -41.943044194304001, 460.06028861053716, 516},
0055                { 1.0000000000000001e-275, -16.777217677721602, 528.36389462809916, 588},
0056                { 1.0000000000000001e-275, -6.7108870710886404, 566.33555010330565, 627},
0057                { 1.0000000000000001e-275, -2.6843548284354561, 585.32541322314046, 646},
0058                { 1.0000000000000001e-275, -1.0737419313741825, 593.0082967458676, 654},
0059                { 1.0000000000000002e-250, -1000000.1, 278584.75303113955, 431059},
0060                { 1.0000000000000002e-250, -400000.03999999998, 111505.53022349965, 111519},
0061                { 1.0000000000000002e-250, -160000.016, 44674.060190374199, 44687},
0062                { 1.0000000000000002e-250, -64000.006399999998, 17941.25452769886, 17954},
0063                { 1.0000000000000002e-250, -25600.002560000001, 7248.2002840909081, 7262},
0064                { 1.0000000000000002e-250, -10240.001024000001, 2972.017267400568, 2985},
0065                { 1.0000000000000002e-250, -4096.0004096000002, 1263.6363636363637, 1278},
0066                { 1.0000000000000002e-250, -1638.40016384, 588.35369318181824, 604},
0067                { 1.0000000000000002e-250, -655.36006553599998, 334.18039772727275, 355},
0068                { 1.0000000000000002e-250, -262.14402621440001, 273.7357954545455, 306},
0069                { 1.0000000000000002e-250, -104.85761048576001, 323.05010330578506, 370},
0070                { 1.0000000000000002e-250, -41.943044194304001, 407.41541838842977, 463},
0071                { 1.0000000000000002e-250, -16.777217677721602, 472.78376807851242, 532},
0072                { 1.0000000000000002e-250, -6.7108870710886404, 509.47443181818176, 570},
0073                { 1.0000000000000002e-250, -2.6843548284354561, 528.11466942148763, 589},
0074                { 1.0000000000000002e-250, -1.0737419313741825, 535.52815082644634, 596},
0075                { 1.0000000000000003e-225, -1000000.1, 278572.12184906006, 278585},
0076                { 1.0000000000000003e-225, -400000.03999999998, 111493.21068649805, 551165},
0077                { 1.0000000000000003e-225, -160000.016, 44661.126509232956, 44674},
0078                { 1.0000000000000003e-225, -64000.006399999998, 17928.069602272728, 17942},
0079                { 1.0000000000000003e-225, -25600.002560000001, 7235.550426136364, 7249},
0080                { 1.0000000000000003e-225, -10240.001024000001, 2959.153053977273, 2973},
0081                { 1.0000000000000003e-225, -4096.0004096000002, 1250.2649147727275, 1264},
0082                { 1.0000000000000003e-225, -1638.40016384, 573.37073863636363, 589},
0083                { 1.0000000000000003e-225, -655.36006553599998, 315.21946022727275, 335},
0084                { 1.0000000000000003e-225, -262.14402621440001, 245.11137654958674, 275},
0085                { 1.0000000000000003e-225, -104.85761048576001, 280.71022727272725, 326},
0086                { 1.0000000000000003e-225, -41.943044194304001, 355.45519111570246, 411},
0087                { 1.0000000000000003e-225, -16.777217677721602, 417.41670971074382, 476},
0088                { 1.0000000000000003e-225, -6.7108870710886404, 452.71984762396687, 513},
0089                { 1.0000000000000003e-225, -2.6843548284354561, 470.45648243801656, 532},
0090                { 1.0000000000000003e-225, -1.0737419313741825, 478.06438855888422, 539},
0091                { 1.0000000000000004e-200, -1000000.1, 278559.54284667969, 278573},
0092                { 1.0000000000000004e-200, -400000.03999999998, 111480.15247691762, 111493},
0093                { 1.0000000000000004e-200, -160000.016, 44648.944635564636, 44662},
0094                { 1.0000000000000004e-200, -64000.006399999998, 17916.076071999294, 17929},
0095                { 1.0000000000000004e-200, -25600.002560000001, 7223.06005859375, 7236},
0096                { 1.0000000000000004e-200, -10240.001024000001, 2946.222301136364, 2960},
0097                { 1.0000000000000004e-200, -4096.0004096000002, 1236.0916193181818, 1251},
0098                { 1.0000000000000004e-200, -1638.40016384, 558.40909090909076, 574},
0099                { 1.0000000000000004e-200, -655.36006553599998, 297.05220170454544, 316},
0100                { 1.0000000000000004e-200, -262.14402621440001, 218.16438533057845, 247},
0101                { 1.0000000000000004e-200, -104.85761048576001, 239.76949896694211, 283},
0102                { 1.0000000000000004e-200, -41.943044194304001, 304.40728305785115, 359},
0103                { 1.0000000000000004e-200, -16.777217677721602, 362.01349431818181, 421},
0104                { 1.0000000000000004e-200, -6.7108870710886404, 396.18123708677683, 456},
0105                { 1.0000000000000004e-200, -2.6843548284354561, 413.36518595041321, 474},
0106                { 1.0000000000000004e-200, -1.0737419313741825, 420.31249999999994, 482},
0107                { 1.0000000000000006e-175, -1000000.1, 278547.12638697342, 465687},
0108                { 1.0000000000000006e-175, -400000.03999999998, 111468.15407492899, 111481},
0109                { 1.0000000000000006e-175, -160000.016, 44636.134832208802, 44649},
0110                { 1.0000000000000006e-175, -64000.006399999998, 17903.091796875, 17917},
0111                { 1.0000000000000006e-175, -25600.002560000001, 7210.247247869318, 7223},
0112                { 1.0000000000000006e-175, -10240.001024000001, 2933.2915482954545, 2947},
0113                { 1.0000000000000006e-175, -4096.0004096000002, 1223.600497159091, 1237},
0114                { 1.0000000000000006e-175, -1638.40016384, 544.06640625, 559},
0115                { 1.0000000000000006e-175, -655.36006553599998, 279.11931818181813, 298},
0116                { 1.0000000000000006e-175, -262.14402621440001, 192.94921875, 220},
0117                { 1.0000000000000006e-175, -104.85761048576001, 200.78641528925618, 242},
0118                { 1.0000000000000006e-175, -41.943044194304001, 254.19679752066114, 308},
0119                { 1.0000000000000006e-175, -16.777217677721602, 307.10227272727263, 366},
0120                { 1.0000000000000006e-175, -6.7108870710886404, 339.36918904958679, 400},
0121                { 1.0000000000000006e-175, -2.6843548284354561, 356.1554106404958, 417},
0122                { 1.0000000000000006e-175, -1.0737419313741825, 363.13694473140492, 424},
0123                { 1.0000000000000007e-150, -1000000.1, 278534.48812689661, 448872},
0124                { 1.0000000000000007e-150, -400000.03999999998, 111455.15580610794, 111468},
0125                { 1.0000000000000007e-150, -160000.016, 44623.643454811798, 44637},
0126                { 1.0000000000000007e-150, -64000.006399999998, 17891.00373493541, 17904},
0127                { 1.0000000000000007e-150, -25600.002560000001, 7197.631392045454, 7211},
0128                { 1.0000000000000007e-150, -10240.001024000001, 2920.3607954545455, 2934},
0129                { 1.0000000000000007e-150, -4096.0004096000002, 1210.3267045454545, 1224},
0130                { 1.0000000000000007e-150, -1638.40016384, 530.0033735795455, 545},
0131                { 1.0000000000000007e-150, -655.36006553599998, 262.37215909090907, 280},
0132                { 1.0000000000000007e-150, -262.14402621440001, 169.07024793388427, 194},
0133                { 1.0000000000000007e-150, -104.85761048576001, 164.12706611570246, 203},
0134                { 1.0000000000000007e-150, -41.943044194304001, 206.02079028925615, 258},
0135                { 1.0000000000000007e-150, -16.777217677721602, 253.0191115702479, 311},
0136                { 1.0000000000000007e-150, -6.7108870710886404, 283.23056559917353, 343},
0137                { 1.0000000000000007e-150, -2.6843548284354561, 299.29041838842977, 360},
0138                { 1.0000000000000007e-150, -1.0737419313741825, 305.7925490702479, 367},
0139                { 1.0000000000000008e-125, -1000000.1, 278522.04145396838, 285024},
0140                { 1.0000000000000008e-125, -400000.03999999998, 111443.00279374557, 111456},
0141                { 1.0000000000000008e-125, -160000.016, 44611.143155184654, 44624},
0142                { 1.0000000000000008e-125, -64000.006399999998, 17878.306263316765, 17891},
0143                { 1.0000000000000008e-125, -25600.002560000001, 7185.1384943181802, 7198},
0144                { 1.0000000000000008e-125, -10240.001024000001, 2907.430042613636, 2921},
0145                { 1.0000000000000008e-125, -4096.0004096000002, 1197.0621448863635, 1211},
0146                { 1.0000000000000008e-125, -1638.40016384, 516.05983664772725, 531},
0147                { 1.0000000000000008e-125, -655.36006553599998, 246.17897727272725, 263},
0148                { 1.0000000000000008e-125, -262.14402621440001, 147.2197830578512, 171},
0149                { 1.0000000000000008e-125, -104.85761048576001, 130.599173553719, 166},
0150                { 1.0000000000000008e-125, -41.943044194304001, 159.12190082644628, 209},
0151                { 1.0000000000000008e-125, -16.777217677721602, 199.5635330578512, 256},
0152                { 1.0000000000000008e-125, -6.7108870710886404, 227.1807205578512, 287},
0153                { 1.0000000000000008e-125, -2.6843548284354561, 242.41670971074376, 303},
0154                { 1.0000000000000008e-125, -1.0737419313741825, 248.39036673553716, 309},
0155                { 1.0000000000000009e-100, -1000000.1, 278509.36086203833, 278546},
0156                { 1.0000000000000009e-100, -400000.03999999998, 111430.10303574696, 551165},
0157                { 1.0000000000000009e-100, -160000.016, 44598.652055220169, 44612},
0158                { 1.0000000000000009e-100, -64000.006399999998, 17865.125532670456, 17879},
0159                { 1.0000000000000009e-100, -25600.002560000001, 7172.046519886364, 7186},
0160                { 1.0000000000000009e-100, -10240.001024000001, 2895.015980113636, 2908},
0161                { 1.0000000000000009e-100, -4096.0004096000002, 1184.2329545454545, 1198},
0162                { 1.0000000000000009e-100, -1638.40016384, 502.04403409090901, 517},
0163                { 1.0000000000000009e-100, -655.36006553599998, 231.13849431818178, 247},
0164                { 1.0000000000000009e-100, -262.14402621440001, 128.02290482954544, 149},
0165                { 1.0000000000000009e-100, -104.85761048576001, 100.76188016528923, 133},
0166                { 1.0000000000000009e-100, -41.943044194304001, 115.10072314049584, 162},
0167                { 1.0000000000000009e-100, -16.777217677721602, 147.05255681818178, 203},
0168                { 1.0000000000000009e-100, -6.7108870710886404, 171.19963842975204, 231},
0169                { 1.0000000000000009e-100, -2.6843548284354561, 185.58367768595039, 246},
0170                { 1.0000000000000009e-100, -1.0737419313741825, 191.0575929752066, 252},
0171                { 1.0000000000000009e-75, -1000000.1, 278497.00512279174, 278510},
0172                { 1.0000000000000009e-75, -400000.03999999998, 111417.96751369131, 111431},
0173                { 1.0000000000000009e-75, -160000.016, 44586.151478160507, 44599},
0174                { 1.0000000000000009e-75, -64000.006399999998, 17853.128506747154, 17866},
0175                { 1.0000000000000009e-75, -25600.002560000001, 7160.0009918212891, 7173},
0176                { 1.0000000000000009e-75, -10240.001024000001, 2882.563210227273, 2896},
0177                { 1.0000000000000009e-75, -4096.0004096000002, 1171.4037642045455, 1185},
0178                { 1.0000000000000009e-75, -1638.40016384, 489.03941761363637, 503},
0179                { 1.0000000000000009e-75, -655.36006553599998, 216.05113636363637, 232},
0180                { 1.0000000000000009e-75, -262.14402621440001, 110.08296745867767, 129},
0181                { 1.0000000000000009e-75, -104.85761048576001, 74.400826446280988, 102},
0182                { 1.0000000000000009e-75, -41.943044194304001, 75.619834710743788, 118},
0183                { 1.0000000000000009e-75, -16.777217677721602, 96.373966942148741, 150},
0184                { 1.0000000000000009e-75, -6.7108870710886404, 116.33910123966942, 175},
0185                { 1.0000000000000009e-75, -2.6843548284354561, 129.00947507747935, 189},
0186                { 1.0000000000000009e-75, -1.0737419313741825, 134.06379132231399, 195},
0187                { 1.0000000000000011e-50, -1000000.1, 278484.14271457132, 465687},
0188                { 1.0000000000000011e-50, -400000.03999999998, 111405.16246448863, 111418},
0189                { 1.0000000000000011e-50, -160000.016, 44573.165838068184, 44587},
0190                { 1.0000000000000011e-50, -64000.006399999998, 17840.147727272728, 17854},
0191                { 1.0000000000000011e-50, -25600.002560000001, 7147.375887784091, 7160},
0192                { 1.0000000000000011e-50, -10240.001024000001, 2869.888583096591, 2883},
0193                { 1.0000000000000011e-50, -4096.0004096000002, 1158.1576704545455, 1172},
0194                { 1.0000000000000011e-50, -1638.40016384, 476.00257457386363, 490},
0195                { 1.0000000000000011e-50, -655.36006553599998, 202.20170454545453, 217},
0196                { 1.0000000000000011e-50, -262.14402621440001, 94.02682722107437, 111},
0197                { 1.0000000000000011e-50, -104.85761048576001, 53.0810950413223, 76},
0198                { 1.0000000000000011e-50, -41.943044194304001, 42.355371900826441, 78},
0199                { 1.0000000000000011e-50, -16.777217677721602, 50.051652892561975, 100},
0200                { 1.0000000000000011e-50, -6.7108870710886404, 63.174070247933884, 120},
0201                { 1.0000000000000011e-50, -2.6843548284354561, 72.956805268595033, 132},
0202                { 1.0000000000000011e-50, -1.0737419313741825, 77.104855371900811, 138},
0203                { 1.0000000000000012e-25, -1000000.1, 278471.90274041548, 448844},
0204                { 1.0000000000000012e-25, -400000.03999999998, 111392.85056374289, 111423},
0205                { 1.0000000000000012e-25, -160000.016, 44561.0361328125, 44574},
0206                { 1.0000000000000012e-25, -64000.006399999998, 17828.149314186794, 17841},
0207                { 1.0000000000000012e-25, -25600.002560000001, 7134.765625, 7148},
0208                { 1.0000000000000012e-25, -10240.001024000001, 2857.211647727273, 2870},
0209                { 1.0000000000000012e-25, -4096.0004096000002, 1146.1576704545453, 1159},
0210                { 1.0000000000000012e-25, -1638.40016384, 463.0625, 476},
0211                { 1.0000000000000012e-25, -655.36006553599998, 189.27556818181819, 203},
0212                { 1.0000000000000012e-25, -262.14402621440001, 79.163223140495859, 95},
0213                { 1.0000000000000012e-25, -104.85761048576001, 36.15702479338843, 54},
0214                { 1.0000000000000012e-25, -41.943044194304001, 19.008264462809912, 44},
0215                { 1.0000000000000012e-25, -16.777217677721602, 14.220686983471072, 53},
0216                { 1.0000000000000012e-25, -6.7108870710886404, 15.207081141998497, 66},
0217                { 1.0000000000000012e-25, -2.6843548284354561, 19.000046957175051, 76},
0218                { 1.0000000000000012e-25, -1.0737419313741825, 21.00738364892468, 81},
0219                { 1, -1000000.1, 278455, 285026},
0220                { 1, -400000.03999999998, 111376, 111392},
0221                { 1, -160000.016, 44544, 44561},
0222                { 1, -64000.006399999998, 17812, 17828},
0223                { 1, -25600.002560000001, 7119, 7135},
0224                { 1, -10240.001024000001, 2841, 2858},
0225                { 1, -4096.0004096000002, 1130, 1147},
0226                { 1, -1638.40016384, 447, 464},
0227                { 1, -655.36006553599998, 174, 190},
0228                { 1, -262.14402621440001, 64, 81},
0229                { 1, -104.85761048576001, 21, 37},
0230                { 1, -41.943044194304001, 5, 20},
0231                { 1, -16.777217677721602, 0, 16},
0232                { 1, -6.7108870710886404, 0, 17},
0233                { 1, -2.6843548284354561, 0, 20},
0234                { 1, -1.0737419313741825, 0, 21},
0235                { 2.5, -1000000.1, 278450, 278466},
0236                { 2.5, -400000.03999999998, 111372, 111388},
0237                { 2.5, -160000.016, 44540, 44557},
0238                { 2.5, -64000.006399999998, 17808, 17824},
0239                { 2.5, -25600.002560000001, 7115, 7131},
0240                { 2.5, -10240.001024000001, 2838, 2854},
0241                { 2.5, -4096.0004096000002, 1127, 1144},
0242                { 2.5, -1638.40016384, 445, 461},
0243                { 2.5, -655.36006553599998, 172, 188},
0244                { 2.5, -262.14402621440001, 63, 79},
0245                { 2.5, -104.85761048576001, 20, 35},
0246                { 2.5, -41.943044194304001, 4, 19},
0247                { 2.5, -16.777217677721602, 0, 13},
0248                { 2.5, -6.7108870710886404, 0, 13},
0249                { 2.5, -2.6843548284354561, 0, 15},
0250                { 2.5, -1.0737419313741825, 2, 16},
0251                { 6.25, -1000000.1, 278440, 278456},
0252                { 6.25, -400000.03999999998, 111362, 111378},
0253                { 6.25, -160000.016, 44531, 44548},
0254                { 6.25, -64000.006399999998, 17800, 17816},
0255                { 6.25, -25600.002560000001, 7108, 7124},
0256                { 6.25, -10240.001024000001, 2831, 2848},
0257                { 6.25, -4096.0004096000002, 1121, 1138},
0258                { 6.25, -1638.40016384, 440, 456},
0259                { 6.25, -655.36006553599998, 167, 183},
0260                { 6.25, -262.14402621440001, 59, 75},
0261                { 6.25, -104.85761048576001, 18, 32},
0262                { 6.25, -41.943044194304001, 3, 16},
0263                { 6.25, -16.777217677721602, 0, 10},
0264                { 6.25, -6.7108870710886404, 0, 9},
0265                { 6.25, -2.6843548284354561, 0, 9},
0266                { 6.25, -1.0737419313741825, 4, 9},
0267                { 15.625, -1000000.1, 278415, 278432},
0268                { 15.625, -400000.03999999998, 111339, 111425},
0269                { 15.625, -160000.016, 44510, 44527},
0270                { 15.625, -64000.006399999998, 17781, 17797},
0271                { 15.625, -25600.002560000001, 7091, 7107},
0272                { 15.625, -10240.001024000001, 2816, 2833},
0273                { 15.625, -4096.0004096000002, 1108, 1125},
0274                { 15.625, -1638.40016384, 429, 445},
0275                { 15.625, -655.36006553599998, 159, 174},
0276                { 15.625, -262.14402621440001, 53, 67},
0277                { 15.625, -104.85761048576001, 14, 27},
0278                { 15.625, -41.943044194304001, 2, 11},
0279                { 15.625, -16.777217677721602, 0, 7},
0280                { 15.625, -6.7108870710886404, 0, 5},
0281                { 15.625, -2.6843548284354561, 5, 5},
0282                { 15.625, -1.0737419313741825, 5, 5},
0283                { 39.0625, -1000000.1, 278359, 431061},
0284                { 39.0625, -400000.03999999998, 111287, 111304},
0285                { 39.0625, -160000.016, 44463, 44480},
0286                { 39.0625, -64000.006399999998, 17738, 17755},
0287                { 39.0625, -25600.002560000001, 7053, 7069},
0288                { 39.0625, -10240.001024000001, 2784, 2800},
0289                { 39.0625, -4096.0004096000002, 1081, 1097},
0290                { 39.0625, -1638.40016384, 407, 422},
0291                { 39.0625, -655.36006553599998, 143, 157},
0292                { 39.0625, -262.14402621440001, 43, 56},
0293                { 39.0625, -104.85761048576001, 10, 19},
0294                { 39.0625, -41.943044194304001, 0, 7},
0295                { 39.0625, -16.777217677721602, 0, 4},
0296                { 39.0625, -6.7108870710886404, 0, 3},
0297                { 39.0625, -2.6843548284354561, 3, 3},
0298                { 39.0625, -1.0737419313741825, 3, 3},
0299                { 97.65625, -1000000.1, 278230, 278277},
0300                { 97.65625, -400000.03999999998, 111170, 111425},
0301                { 97.65625, -160000.016, 44358, 44374},
0302                { 97.65625, -64000.006399999998, 17645, 17661},
0303                { 97.65625, -25600.002560000001, 6972, 6988},
0304                { 97.65625, -10240.001024000001, 2715, 2731},
0305                { 97.65625, -4096.0004096000002, 1026, 1041},
0306                { 97.65625, -1638.40016384, 366, 380},
0307                { 97.65625, -655.36006553599998, 117, 129},
0308                { 97.65625, -262.14402621440001, 30, 40},
0309                { 97.65625, -104.85761048576001, 5, 12},
0310                { 97.65625, -41.943044194304001, 0, 4},
0311                { 97.65625, -16.777217677721602, 0, 3},
0312                { 97.65625, -6.7108870710886404, 0, 1},
0313                { 97.65625, -2.6843548284354561, 0, 1},
0314                { 97.65625, -1.0737419313741825, 0, 1},
0315                { 244.140625, -1000000.1, 277936, 796691},
0316                { 244.140625, -400000.03999999998, 110906, 110923},
0317                { 244.140625, -160000.016, 44124, 44141},
0318                { 244.140625, -64000.006399999998, 17442, 17458},
0319                { 244.140625, -25600.002560000001, 6801, 6816},
0320                { 244.140625, -10240.001024000001, 2577, 2593},
0321                { 244.140625, -4096.0004096000002, 924, 938},
0322                { 244.140625, -1638.40016384, 300, 312},
0323                { 244.140625, -655.36006553599998, 82, 92},
0324                { 244.140625, -262.14402621440001, 17, 24},
0325                { 244.140625, -104.85761048576001, 2, 7},
0326                { 244.140625, -41.943044194304001, 0, 3},
0327                { 244.140625, -16.777217677721602, 0, 1},
0328                { 244.140625, -6.7108870710886404, 0, 1},
0329                { 244.140625, -2.6843548284354561, 0, 1},
0330                { 244.140625, -1.0737419313741825, 0, 1},
0331                { 610.3515625, -1000000.1, 277277, 277294},
0332                { 610.3515625, -400000.03999999998, 110322, 111647},
0333                { 610.3515625, -160000.016, 43617, 43633},
0334                { 610.3515625, -64000.006399999998, 17014, 17030},
0335                { 610.3515625, -25600.002560000001, 6456, 6471},
0336                { 610.3515625, -10240.001024000001, 2321, 2335},
0337                { 610.3515625, -4096.0004096000002, 757, 769},
0338                { 610.3515625, -1638.40016384, 212, 221},
0339                { 610.3515625, -655.36006553599998, 48, 55},
0340                { 610.3515625, -262.14402621440001, 8, 13},
0341                { 610.3515625, -104.85761048576001, 0, 4},
0342                { 610.3515625, -41.943044194304001, 0, 1},
0343                { 610.3515625, -16.777217677721602, 0, 1},
0344                { 610.3515625, -6.7108870710886404, 0, 1},
0345                { 610.3515625, -2.6843548284354561, 0, 1},
0346                { 610.3515625, -1.0737419313741825, 0, 1},
0347                { 1525.87890625, -1000000.1, 275817, 431061},
0348                { 1525.87890625, -400000.03999999998, 109054, 109070},
0349                { 1525.87890625, -160000.016, 42546, 42562},
0350                { 1525.87890625, -64000.006399999998, 16151, 16167},
0351                { 1525.87890625, -25600.002560000001, 5814, 5828},
0352                { 1525.87890625, -10240.001024000001, 1902, 1914},
0353                { 1525.87890625, -4096.0004096000002, 535, 545},
0354                { 1525.87890625, -1638.40016384, 125, 131},
0355                { 1525.87890625, -655.36006553599998, 23, 28},
0356                { 1525.87890625, -262.14402621440001, 3, 7},
0357                { 1525.87890625, -104.85761048576001, 0, 1},
0358                { 1525.87890625, -41.943044194304001, 0, 1},
0359                { 1525.87890625, -16.777217677721602, 0, 1},
0360                { 1525.87890625, -6.7108870710886404, 0, 1},
0361                { 1525.87890625, -2.6843548284354561, 0, 1},
0362                { 1525.87890625, -1.0737419313741825, 0, 1},
0363                { 3814.697265625, -1000000.1, 272645, 272661},
0364                { 3814.697265625, -400000.03999999998, 106377, 106393},
0365                { 3814.697265625, -160000.016, 40390, 40406},
0366                { 3814.697265625, -64000.006399999998, 14545, 14560},
0367                { 3814.697265625, -25600.002560000001, 4765, 4777},
0368                { 3814.697265625, -10240.001024000001, 1346, 1356},
0369                { 3814.697265625, -4096.0004096000002, 316, 323},
0370                { 3814.697265625, -1638.40016384, 61, 66},
0371                { 3814.697265625, -655.36006553599998, 9, 13},
0372                { 3814.697265625, -262.14402621440001, 0, 4},
0373                { 3814.697265625, -104.85761048576001, 0, 1},
0374                { 3814.697265625, -41.943044194304001, 0, 1},
0375                { 3814.697265625, -16.777217677721602, 0, 1},
0376                { 3814.697265625, -6.7108870710886404, 0, 1},
0377                { 3814.697265625, -2.6843548284354561, 0, 1},
0378                { 3814.697265625, -1.0737419313741825, 0, 1},
0379                { 9536.7431640625, -1000000.1, 265953, 266155},
0380                { 9536.7431640625, -400000.03999999998, 100987, 101002},
0381                { 9536.7431640625, -160000.016, 36374, 36388},
0382                { 9536.7431640625, -64000.006399999998, 11921, 11934},
0383                { 9536.7431640625, -25600.002560000001, 3374, 3384},
0384                { 9536.7431640625, -10240.001024000001, 795, 802},
0385                { 9536.7431640625, -4096.0004096000002, 157, 162},
0386                { 9536.7431640625, -1638.40016384, 27, 30},
0387                { 9536.7431640625, -655.36006553599998, 3, 6},
0388                { 9536.7431640625, -262.14402621440001, 0, 1},
0389                { 9536.7431640625, -104.85761048576001, 0, 1},
0390                { 9536.7431640625, -41.943044194304001, 0, 1},
0391                { 9536.7431640625, -16.777217677721602, 0, 1},
0392                { 9536.7431640625, -6.7108870710886404, 0, 1},
0393                { 9536.7431640625, -2.6843548284354561, 0, 1},
0394                { 9536.7431640625, -1.0737419313741825, 0, 1},
0395             };
0396             static const unsigned total_elements = sizeof(domain) / sizeof(domain[0]);
0397             static const unsigned stride = 16;
0398             //static const unsigned a_elements = total_elements / stride;
0399             BOOST_MATH_ASSERT(total_elements % stride == 0);
0400
0401             static const double a_max = domain[total_elements - 1][0];
0402             static const double a_min = domain[0][0];
0403             static const double b_max = domain[stride - 1][1];
0404             static const double b_min = domain[0][1];
0405
0406             if (a < a_min)
0407                return 0; // TODO: Use series?
0408             if (b < b_min)
0409             {
0410                //
0411                // This is a general heuristic that's more or less correct,
0412                // but a bit too safe in most cases.  Checked OK with
0413                // random testing.
0414                //
0415                if (z > -b)
0416                   return 1;
0417                if (a < 100)
0418                   return z < -b / (4 - 5 * (log(a)) * a / b) ? -1 : 0;
0419                else
0420                   return z < -b / (4 - 5 * sqrt(log(a)) * a / b) ? -1 : 0;
0421             }
0422             if (a > a_max)
0423             {
0424                //
0425                // Crossover point is decreasing with increasing a
0426             // upper limit is fine, lower limit is not:
0427                //
0428                if (b > b_max)
0429                   return 0;  // TODO: don't know what else to do???
0430                unsigned index = total_elements - stride;
0431                BOOST_MATH_ASSERT(domain[index][0] == a_max);
0432                while (domain[index][1] < b)
0433                   ++index;
0434                double b0 = domain[index - 1][1];
0435                double b1 = domain[index][1];
0436                double z0 = domain[index - 1][3];
0437                double z1 = domain[index][3];
0438                T upper_z_limit = static_cast<T>((z0 * (b1 - b) + z1 * (b - b0)) / (b1 - b0));
0439                if (z > upper_z_limit)
0440                   return 1;
0441                return z < -b / (4 - 5 * sqrt(log(a)) * a / b) ? -1 : 0;
0442
0443             }
0444             if (b > b_max)
0445             {
0446                return 0;  // TODO: is there a better way???
0447             }
0448             //
0449             // Bi-linear interpolation case:
0450             //
0451             unsigned index = 0;
0452             while (domain[index][0] < a)
0453                index += stride;
0454             while (domain[index][1] < b)
0455                ++index;
0456             //
0457             // We now have 4 corners:
0458             //
0459             double x1 = domain[index - stride - 1][0];
0460             double x2 = domain[index][0];
0461             double y1 = domain[index - 1][1];
0462             double y2 = domain[index][1];
0463             double f11 = domain[index - stride - 1][2];
0464             BOOST_MATH_ASSERT(domain[index - stride - 1][0] == x1);
0465             BOOST_MATH_ASSERT(domain[index - stride - 1][1] == y1);
0466             double f12 = domain[index - stride][2];
0467             BOOST_MATH_ASSERT(domain[index - stride][0] == x1);
0468             BOOST_MATH_ASSERT(domain[index - stride][1] == y2);
0469             double f21 = domain[index - 1][2];
0470             BOOST_MATH_ASSERT(domain[index - 1][0] == x2);
0471             BOOST_MATH_ASSERT(domain[index - 1][1] == y1);
0472             double f22 = domain[index][2];
0473             BOOST_MATH_ASSERT(domain[index][0] == x2);
0474             BOOST_MATH_ASSERT(domain[index][1] == y2);
0475             //
0476             // Interpolation is a crude tool and our corners have quite tight bounds,
0477             // and the "impossible region" is bounded by convex planes.  For the upper
0478             // bound this is fine - bilinear interpolation will be over-conservative
0479             // but safe.  For the lower limit though we may calculate that the "safe"
0480             // zone is in an unsafe place.  To work around this we move our coordinates
0481             // closer to the lowest corner of our bounding box by an amount
0482             // proportional to the distance from the nearest corner.  This has the effect
0483             // of making our lower bound smaller than it would otherwise be, and more so nearer
0484             // the centre of the bounding box.
0485             //
0486             T effective_x = static_cast<T>(a + (std::min)(T(a - x1), T(x2 - a)) * 0.25);
0487             T effective_y = static_cast<T>(b + (std::min)(T(b - y1), T(y2 - b)) * 0.25);
0488
0489             T lower_limit = static_cast<T>(1 / ((x2 - x1) * (y2 - y1)) * (f11 * (x2 - effective_x) * (y2 - effective_y) + f21 * (effective_x - x1) * (y2 - effective_y) + f12 * (x2 - effective_x) * (effective_y - y1) + f22 * (effective_x - x1) * (effective_y - y1)));
0490
0491             //
0492             // If one f value was zero, take the whole box as zero:
0493             //
0494             double min_f = (std::min)((std::min)(f11, f12), (std::min)(f21, f22));
0495             if (min_f == 0)
0496                lower_limit = 0;
0497             //
0498             // Now we can test!
0499             //
0500             if (z < lower_limit)
0501                return -1;
0502
0503             BOOST_MATH_ASSERT(f11 <= domain[index - stride - 1][3]);
0504             f11 = domain[index - stride - 1][3];
0505             BOOST_MATH_ASSERT(f12 <= domain[index - stride][3]);
0506             f12 = domain[index - stride][3];
0507             BOOST_MATH_ASSERT(f21 <= domain[index - 1][3]);
0508             f21 = domain[index - 1][3];
0509             BOOST_MATH_ASSERT(f22 <= domain[index][3]);
0510             f22 = domain[index][3];
0511
0512             T upper_limit = static_cast<T>(1 / ((x2 - x1) * (y2 - y1)) * (f11 * (x2 - a) * (y2 - b) + f21 * (a - x1) * (y2 - b) + f12 * (x2 - a) * (b - y1) + f22 * (a - x1) * (b - y1)));
0513             if (z > upper_limit)
0514                return 1;
0515             return 0;
0516          }
0517
0518 } } }
0519
0520
0521 #endif   // BOOST_MATH_DETIAL_1F1_MAP_NEG_B_HPP