I'm guessing that the rcond number in the lstsq version (default 1e-10) is the difference. Generally the lstsq version should work better than the MPL version because at*a is not as well conditioned and vandermonde matrices are notoriously ill conditioned anyway for higher degree fits. It would help if you attached the data files in ascii form unless they happen to contain thousands of data items. Just the x will suffice, zip it up if you have to.
Here are the files.
Since the two algorithms behave differently and each has it place then can both be included in numpy?
i.e. numpy.polyfit(x,y,N, mode='alg1')
numpy.polyfit (x,y,N, mode='alg2')
replacing alg1 and alg2 with meaningful names.
Thanks
Greg
2 4 0 16 1 3 14 17 12 13 8 7 5 11 23 19 24 9 18 15 44 22 20 25 27 10 41 45 26 40 21 50 29 38 28 37 53 48 57 78 54 61 39 64 65 58 43 35 55 47 60 66 46 68 52 51 42 75 93 56 87 74 59 71 84 49 72 81 73 63 94 67 36 80 30 70 106 69 62 77 76 90 98 103 6 172 79 83 82 31 85 97 105 91 121 89 95 101 116 135 92 102 96 88 120 137 136 34 389 499 489 565 479 555 412 390 443 442 456 365 404 593 457 364 375 441 367 360 453 376 482 576 550 445 409 658 363 435 381 477 549 396 460 544 601 466 384 557 391 461 511 642 547 558 542 610 608 639 438 591 551 378 599 552 397 569 578 437 643 361 385 399 538 600 395 561 602 649 454 464 447 386 535 554 577 596 686 635 459 514 628 448 583 407 408 603 379 626 410 500 570 449 582 434 401 417 501 462 566 366 394 631 485 382 563 369 784 781 824 838 783 826 825 854 867 763 801 799 879 795 759 852 872 831 802 769 873 773 762 779 772 789 836 881 839 827 796 882 851 786 797 755 726 841 896 865 890 837 857 775 853 767 785 756 750 754 887 861 894 770 749 771 840 844 782 780 877 859 774 832 871 730 722 901 794 884 846 760 848 729 876 724 720 757 830 798 850 818 764 875 787 829 899 885 874 960 815 864 777 723 820 891 747 751 725 821 776 823 866 800 833 788 758 878 1411 1350 1358 1353 1362 1378 1326 1348 1321 1355 1357 1349 1352 1346 1373 1430 1408 1351 1310 1418 1281 1336 1416 1428 1423 1436 1364 1432 1419 1377 1421 1359 1279 1401 1424 1323 1340 1322 1339 1407 1438 1427 1247 1435 1290 1311 1262 1343 1307 1437 1369 1315 1288 1356 1412 1309 1327 1276 1367 1325 1425 1420 1344 1361 1337 1368 1314 1415 1335 1278 1320 1280 1382 1282 1434 1328 1273 1332 1410 1347 1284 1338 1417 1354 1360 1329 1270 1305 1375 1414 1287 1291 1429 1379 1256 1372 1370 1289 1306 1275 1342 1274 1286 1313 1271 1248 1319 1268 1479 1440 1452 1449 1441 1542 1444 1447 1645 1569 1485 1537 1527 1522 1462 1461 1554 1651 1448 1652 1538 1529 1510 1451 1639 1570 1486 1737 1481 1623 1725 1453 1458 1459 1680 1484 1445 1671 1460 1524 1673 1666 1627 1540 1564 1626 1672 1499 1674 1561 1555 1519 1472 1632 1534 1548 1650 1470 1656 1637 1648 1717 1661 1521 1724 1755 1631 1735 1517 1667 1551 1728 1566 1575 1515 1539 1531 1541 1563 1558 1643 1454 1549 1758 1483 1761 1714 1638 1505 1654 1664 1475 1556 1473 1528 1516 1552 1443 1709 1553 1794 1571 1644 1520 1718 1742 1498 1736 2149 2052 2128 2141 2019 2108 2125 1998 2055 2087 2106 2123 2029 2036 1960 2121 2049 1945 2132 2058 2035 2144 2053 2109 1963 1949 1966 2137 2060 2153 2147 1956 2126 2115 2130 2129 1808 2021 2038 1804 2064 2104 1952 2107 1911 1958 2116 1810 2031 1940 2148 2117 2150 1814 2059 1916 2010 2026 1809 1964 2103 2032 2039 2118 2151 1973 2152 1951 2084 2011 2065 2105 2028 2112 2044 2040 2122 2033 1812 1941 1946 1987 1996 2127 2113 2048 1995 1950 2155 1959 2138 1803 1970 1961 2020 2008 1811 1968 1965 1947 1819 2037 2063 1817 2133 1967 2140 2047 2185 2161 2217 2241 2163 2164 2212 2220 2247 2162 2208 2248 2225 2255 2206 2191 2218 2165 2238 2237 2166 2231 2211 2192 2258 2233 2201 2221 2239 2251 2249 2219 2186 2197 2246 2210 2199 2184 2226 2243 2188 2240 2250 2204 2205 2193 2260 2160 2236 2187 2228 2182 2223 2215 2224 2235 2213 2198 2168 2202 2195 2183 2253 2209 2234 2288 2261 2256 2254 2200 2257 2242 2167 2280 2252 2245 2207 2263 2216 2277 2284 2244 2190 2259 2266 2267 2180 2227 2196 2229 2297 2189 2169 2296 2181 2279 2281 2179 2178 2194 2287 2222 2285 2232 2301 2264 2214 2289 2787 2714 2823 2591 2640 2717 2650 2612 2838 2630 2827 2588 2685 2684 2677 2720 2613 2872 2729 2804 2622 2841 2642 2694 2867 2693 2865 2601 2808 2870 2691 2607 2877 2620 2836 2835 2576 2596 2636 2805 2802 2603 2579 2797 2610 2577 2716 2832 2626 2676 2678 2873 2785 2632 2608 2692 2643 2726 2813 2614 2675 2791 2803 2644 2701 2801 2594 2737 2625 2826 2839 2649 2700 2799 2634 2695 2730 2592 2874 2667 2647 2683 2586 2722 2679 2669 2629 2871 2549 2733 2582 2631 2807 2606 2641 2727 2724 2652 2623 2783 2646 2821 2602 2702 2708 2589 2815 2842 3222 3230 3229 2898 3228 3214 3232 3208 3213 2931 3202 2886 2922 3205 3221 3235 3212 3204 2887 2893 3209 2907 2921 3231 2891 2914 3223 3224 2890 2911 3189 2896 3233 2909 2917 3220 3198 3207 3226 3227 2900 2889 3225 2923 2884 3203 3192 2910 2926 3210 2885 2905 2913 2895 2897 2918 2924 3200 2927 2934 3234 2912 3206 3187 3236 3196 2881 3238 3190 3219 2899 2936 3237 3199 2919 2915 2882 3215 3188 2902 2937 2906 2930 2901 2883 3197 2908 2904 3193 3201 3218 2925 2892 2916 3194 2939 3195 2888 3191 2880 2894 2935 2932 2944 2920 2928 2929 2938 3581 3576 3572 3571 3568 3575 3569 3566 3562 3574 3577 3565 3570 3559 3550 3582 3578 3585 3583 3560 3573 3563 3580 3588 3586 3584 3579 3564 3567 3561 3549 3587 3551 3590 3591 3589 3558 3515 3523 3511 3526 3552 3512 3547 3510 3508 3520 3503 3524 3509 3518 3557 3545 3548 3514 3592 3517 3519 3528 3513 3522 3488 3521 3487 3516 3544 3497 3553 3507 3525 3504 3546 3499 3490 3527 3492 3494 3501 3489 3500 3593 3491 3496 3311 3458 3306 3594 3495 3394 3543 3498 3277 3486 3502 3370 3505 3308 3460 3469 3459 3466 3483 3366 3482 3269 3323 3264 3485
3956.970544 4114.773479 4313.241696 4384.792216 4391.872310 4411.476715 4482.940096 4594.108179 4645.381743 4730.829231 4763.895598 4781.032135 4880.854414 5008.228140 5013.866252 5026.548238 5040.509458 5061.949397 5098.508781 5118.791696 5162.805554 5186.805236 5193.941710 5256.845108 5286.838670 5364.175004 5413.297293 5418.790951 5446.085519 5475.639165 5508.616723 5604.737760 5614.409364 5636.101761 5638.542797 5649.773243 5660.913105 5719.715027 5815.218454 5832.621843 5837.578731 5862.529421 5863.607227 5877.191757 5952.905334 5984.647641 6033.194005 6048.375830 6058.973718 6073.108082 6100.784624 6112.341023 6138.705305 6142.491247 6170.101238 6190.399915 6224.593338 6239.605558 6248.244072 6254.401594 6280.211422 6317.197093 6317.614065 6325.296754 6339.852487 6381.476997 6383.259524 6391.571145 6397.082360 6418.080290 6438.735101 6440.724082 6454.977105 6477.089430 6560.501879 6637.842362 6710.064249 6724.529239 6728.116095 6728.478821 6774.863744 6874.242169 6902.087390 6919.527972 6979.799498 6980.973237 6989.850615 7015.371907 7019.392350 7050.062660 7050.597597 7061.397784 7105.324281 7107.397379 7113.332893 7123.788680 7132.977337 7138.892938 7172.808848 7176.242427 7186.719367 7191.745735 7239.760194 7243.749478 7255.868073 7257.267870 7268.538062 7272.629743 7679.940589 7683.236804 7780.699438 7873.928129 7899.908299 7904.060156 7925.056466 7986.214500 7992.112331 7999.563794 8001.244799 8018.175375 8030.246998 8042.876999 8066.748960 8068.723221 8070.288644 8071.628413 8130.859172 8156.118329 8175.197378 8206.961155 8218.103232 8222.490958 8239.676758 8241.640380 8249.615006 8267.669676 8268.053374 8283.319127 8288.302563 8293.092398 8316.055561 8327.772501 8334.822552 8338.317644 8339.073386 8342.683772 8350.156751 8357.692104 8400.200111 8408.910930 8417.560026 8424.256975 8424.301520 8435.165529 8442.029865 8463.791695 8467.858939 8468.002339 8473.892687 8475.457054 8481.274131 8484.884301 8503.376563 8506.217207 8513.976462 8516.443233 8518.839920 8521.927566 8522.369277 8524.426356 8529.354446 8530.550718 8537.266324 8543.730529 8554.276036 8555.003337 8573.504421 8575.663769 8578.398385 8586.508857 8597.926278 8601.801143 8607.638604 8613.355557 8618.223021 8619.832545 8625.856397 8628.547664 8635.510731 8635.616818 8638.093330 8638.725205 8639.290086 8640.963579 8641.074684 8644.931541 8650.771065 8651.748195 8652.153167 8654.191277 8662.647644 8663.908731 8672.156947 8673.406207 8673.653387 8676.156890 8684.297953 8688.001563 8691.611803 8691.714375 8692.654871 8705.209794 8714.218292 8715.503995 8720.084210 8720.465880 8950.316961 9014.234275 9016.582013 9068.120557 9079.086925 9084.084293 9092.854202 9152.516280 9206.442055 9211.797882 9259.935823 9274.088555 9319.200628 9418.736103 9426.658011 9437.435612 9448.483510 9448.919497 9451.123274 9475.891179 9488.585768 9492.490822 9513.276633 9516.127121 9521.205075 9530.898966 9564.375640 9588.755348 9596.967298 9615.132594 9642.107437 9645.510988 9646.175388 9666.070621 9677.726960 9695.380196 9701.239340 9709.602677 9711.670031 9739.504151 9761.667030 9761.780136 9769.510284 9799.923976 9807.513067 9809.803762 9810.694795 9823.853204 9838.334226 9848.403779 9876.275178 9887.129520 9887.560756 9887.888122 9895.889715 9899.907916 9909.691312 9911.111820 9949.714270 9952.911904 9960.038486 9962.479402 9971.996514 9980.326000 9984.951467 9994.006658 10001.073320 10003.476996 10014.757757 10022.105369 10025.979253 10049.807479 10054.324996 10058.877706 10058.904792 10085.774604 10088.651241 10091.809589 10097.187658 10107.697162 10109.684042 10113.805433 10114.595009 10137.018824 10138.092392 10150.691072 10151.749667 10171.572498 10177.785950 10182.552424 10185.711350 10193.402181 10198.124474 10202.876980 10203.927967 10204.398934 10213.144657 10229.042155 10235.314449 10238.453722 10241.626409 10241.820711 10263.338102 10266.747743 10319.492303 10332.302989 10339.341659 10343.745803 8793.912077 9240.132446 9371.045485 9559.747528 9566.311321 9605.890226 9663.291742 9727.898364 9817.130726 9829.579358 9894.337029 9943.329879 9986.701787 9988.681733 10007.595675 10039.175362 10058.567031 10071.732150 10080.223263 10080.572538 10082.640448 10091.338544 10096.808553 10104.572156 10196.766192 10243.213998 10247.809490 10253.163862 10292.904202 10312.906320 10315.642277 10327.044105 10342.839692 10346.473991 10355.275639 10362.684038 10372.115530 10372.748272 10382.069791 10387.721477 10389.628136 10399.630692 10430.582283 10442.936940 10449.245767 10454.870137 10472.091864 10478.317998 10489.109327 10509.306150 10513.053502 10514.097318 10525.194284 10527.980647 10529.909410 10531.238327 10537.874168 10544.844175 10555.905185 10563.925159 10565.773122 10575.339600 10586.124716 10590.426276 10611.035077 10622.729994 10636.607153 10657.784468 10666.334794 10672.555469 10686.503366 10686.595292 10693.772247 10703.728069 10703.988774 10704.045868 10708.960756 10720.027629 10738.600374 10777.951535 10783.717838 10797.095152 10803.299357 10804.900797 10805.415021 10809.890716 10825.610894 10830.556433 10854.840855 10861.376313 10885.763384 10918.933837 10961.511982 10966.953981 10967.776642 10973.773434 10977.243562 10978.798310 10978.830169 10989.002409 11009.161106 11015.743459 11020.016124 11056.625898 11092.348893 11093.400822 11100.135715 11153.180487 9814.960894 9981.143054 10098.094627 10098.936901 10251.240937 10282.032995 10284.640897 10322.598305 10360.982393 10376.872740 10392.495710 10417.277021 10425.686178 10434.329513 10450.355048 10478.041847 10499.163586 10499.215862 10526.877232 10532.554246 10536.409254 10539.253460 10608.396814 10645.575298 10659.240945 10670.922154 10689.184454 10715.836202 10722.298297 10725.640874 10748.404608 10753.162544 10756.517743 10757.745888 10772.695671 10782.428719 10790.583373 10791.898039 10801.684370 10811.977492 10816.266996 10818.283291 10820.330119 10834.649198 10857.733858 10874.140510 10876.059901 10877.637723 10882.405065 10894.207293 10900.073600 10912.167792 10924.820170 10926.306574 10929.257567 10959.628651 10960.178273 10975.859148 10976.998078 10978.812143 10991.596620 10999.187685 11000.809537 11010.126635 11024.307375 11071.514474 11085.007968 11094.736694 11109.770745 11113.527504 11116.002506 11122.254903 11133.988453 11136.612348 11156.603051 11165.834508 11166.019724 11167.998566 11175.515681 11193.546241 11199.134090 11238.198692 11244.422619 11253.787513 11264.283982 11271.505442 11274.799495 11299.434399 11316.162315 11317.946622 11318.921151 11319.649335 11346.282844 11358.697447 11362.651974 11372.919362 11374.332828 11382.460151 11394.420472 11396.858534 11397.989178 11401.709776 11418.079306 11432.017976 11438.011867 11459.901142 11468.659319 11472.542199 10681.749907 10956.703906 11125.372074 11136.314500 11172.483691 11178.012600 11269.727750 11300.670178 11340.773673 11383.722922 11457.659862 11462.856526 11489.727744 11551.874237 11579.544500 11637.402145 11686.213732 11690.080383 11741.072475 11747.628633 11787.678534 11801.682051 11814.371802 11829.473150 11838.299686 11858.353106 11870.960722 11871.392282 11871.694919 11876.150858 11883.402109 11895.260550 11913.376985 11913.856733 11942.894620 11951.949830 11969.258031 11983.875254 11993.436462 11995.024242 12021.141064 12023.333762 12048.485972 12050.270337 12055.970075 12066.039635 12073.944508 12093.030117 12098.907173 12116.424190 12116.649970 12120.302248 12124.960350 12136.994784 12137.012910 12150.298122 12150.452929 12166.326678 12170.039350 12171.669367 12173.044335 12174.841150 12179.606648 12207.744517 12209.538641 12225.212482 12229.865951 12232.316882 12246.263509 12256.108638 12276.069698 12286.868822 12297.149059 12306.631699 12312.322071 12320.164640 12324.374077 12336.310756 12339.376427 12368.598153 12380.799010 12382.583485 12384.933702 12389.988805 12391.470566 12395.745999 12401.831477 12402.594363 12413.304896 12413.370202 12444.759649 12448.495831 12459.577196 12470.502275 12470.742632 12472.315225 12474.617021 12478.433235 12478.483165 12482.236078 12498.028231 12513.399284 12525.231137 12532.312092 12534.743170 12537.889753 12542.266646 12551.826599 11193.319247 11514.884367 11602.445732 11622.344038 11622.438404 11630.019170 11641.485017 11667.278665 11708.588012 11789.904842 11836.441818 11850.623879 11881.255290 11925.047593 11951.347144 11964.909208 11994.423115 12038.711459 12090.526264 12106.978903 12160.550605 12177.736209 12215.007376 12230.055459 12250.739418 12257.282265 12259.861311 12260.730253 12262.735060 12291.224191 12316.483640 12320.419844 12358.768148 12371.123837 12374.651193 12388.097106 12445.256592 12448.681215 12473.311313 12504.297769 12505.134285 12527.592355 12575.220793 12625.944644 12651.232120 12654.872920 12667.986158 12673.676463 12698.035009 12713.500228 12766.790350 12795.378324 12811.393790 12816.705173 12821.887367 12836.305776 12859.372179 12868.193924 12869.219546 12883.997847 12913.567482 12928.892730 12951.562632 12957.166348 13033.097686 13039.613435 13052.920297 13103.059009 13121.399000 13137.420222 13143.819581 13189.273853 13217.056695 13227.901100 13251.107604 13252.729284 13256.649463 13269.470154 13291.730794 13301.492326 13306.257808 13307.889142 13360.021884 13423.017549 13429.234543 13437.179872 13484.090893 13496.163341 13595.522328 13606.157527 13609.179758 13617.822523 13619.608609 13690.242120 13717.833712 13735.254554 13755.434901 13807.103472 13824.931396 13827.734231 13834.988156 13914.748586 13978.380023 14060.743382 14083.144765 14103.587172 14155.685069 14169.683047 12042.398439 12228.734395 12520.514566 12535.035045 12580.133137 12590.412616 12629.404482 12654.158926 12675.304854 12733.821411 12802.369761 12809.260088 12868.977918 12891.461146 12902.360057 12944.236947 12970.555105 12984.851197 12993.054529 13024.315480 13028.277354 13036.186899 13065.224983 13074.307427 13077.734413 13120.257899 13138.725590 13139.563867 13144.916995 13147.708541 13150.495330 13150.670526 13153.697961 13158.210056 13170.801478 13173.506940 13174.016751 13175.306920 13181.264005 13182.217312 13189.197422 13191.435130 13212.256657 13212.308700 13216.024874 13216.050142 13217.485968 13262.914785 13271.722595 13282.591589 13291.680677 13313.503687 13317.334527 13331.044319 13332.073092 13332.970644 13351.033630 13367.919842 13372.757479 13383.369447 13400.388651 13403.526696 13406.483139 13410.260061 13413.739986 13423.056454 13427.512376 13428.578560 13431.234935 13438.083882 13440.731763 13467.215303 13471.345252 13475.473964 13476.412131 13481.174967 13491.373425 13496.015255 13499.707775 13500.224718 13512.408367 13515.414742 13519.258531 13519.955607 13529.853371 13555.701103 13567.556673 13581.752972 13592.412676 13612.170183 13618.520577 13621.250454 13627.232255 13638.740682 13644.590236 13645.126659 13649.339579 13652.869745 13660.267656 13663.052077 13666.620509 13667.808315 13667.827996 13678.639363 13685.255280 13689.203630 13691.048025 13713.641976 11598.693262 11704.336246 11919.040078 11981.528008 12057.828308 12094.424616 12149.980041 12179.583269 12183.210365 12251.745007 12262.323134 12331.253864 12412.695494 12418.496422 12425.118997 12528.458621 12550.929713 12566.605855 12624.572874 12625.259561 12647.308193 12679.031675 12683.923804 12704.924039 12745.050992 12770.103628 12787.591668 12809.165792 12819.808330 12822.677224 12835.990692 12854.239067 12913.753810 12933.290259 12939.277223 12959.460061 12965.124068 12971.594291 12994.159714 13007.621915 13010.606602 13018.524252 13043.586740 13060.549764 13069.738691 13081.265005 13083.369723 13111.744055 13132.289855 13231.443965 13243.583985 13306.335621 13307.422095 13315.641892 13373.612750 13379.343564 13437.392770 13486.774587 13498.942593 13516.270938 13517.006573 13617.129826 13651.302768 13655.976384 13659.231345 13669.637967 13692.471342 13702.477965 13709.766221 13724.485954 13748.119265 13764.708728 13780.746631 13821.058818 13822.464110 13864.480461 13867.177653 13868.782455 13869.689934 13920.956099 13928.184482 13975.859048 13985.023508 13988.347825 13992.894227 14009.444963 14037.708338 14058.926075 14150.438533 14172.273453 14186.946281 14194.290515 14246.683826 14305.245169 14308.202607 14310.050716 14351.339722 14376.973672 14392.142221 14397.889970 14434.301413 14436.668387 14454.755269 14546.002108 14564.107635 14609.352294 14610.444061 14658.388376 6191.820693 6417.995072 6450.171771 6639.135507 6702.781151 6808.454342 6858.701101 6875.949295 6896.356613 6899.675204 6980.449533 7001.842011 7076.463311 7081.714631 7082.396339 7176.817676 7216.309297 7266.648313 7280.514286 7297.524297 7355.848621 7373.904267 7378.860143 7408.329335 7458.162090 7556.879500 7567.809539 7642.724706 7735.706825 8038.110570 8048.230497 8061.020685 8148.930937 8210.087296 8262.783633 8377.205167 8474.150854 8624.533634 8822.084093 8902.373733 9086.760508 9093.644637 9159.773432 9176.508060 9188.687554 9189.349289 9191.366199 9192.863741 9292.189993 9293.081584 9297.837340 9391.546496 9403.902790 9423.233460 9528.854294 9531.233339 9555.442442 9596.527664 9621.645865 9664.165188 9748.118478 9801.103502 9845.340170 9883.906937 9910.602414 9914.261309 9914.838154 9944.961769 9947.652699 10030.420720 10032.086427 10085.772122 10106.567350 10131.776082 10171.059604 10171.271504 10250.514131 10294.686428 10310.671898 10354.773039 10359.713965 10417.971882 10419.604148 10518.670282 10543.442886 10604.550121 10604.651083 10638.592841 10715.252457 10750.162225 10763.995006 10825.711128 10857.760357 10872.843194 10894.852725 10910.851029 10927.298566 10932.509272 10938.224926 11014.452664 11019.602796 11031.241511 11052.111308 11055.718676 11153.896765 11160.080716 11179.872948 11193.467570
------------------------------------------------------------------------- Using Tomcat but need to do more? Need to support web services, security? Get stuff done quickly with pre-integrated technology to make your job easier Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo http://sel.as-us.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642
_______________________________________________ Numpy-discussion mailing list Numpy-discussion@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/numpy-discussion
