5 October 2012 9:44:19.273 AM TEST_INTERP_ND_PRB FORTRAN77 version Test the TEST_INTERP_ND library. The R8LIB library is also needed. TEST01 P00_TITLE returns the problem title. There are a total of 6 problems. 1 "Oscillatory" 2 "Product Peak" 3 "Corner Peak" 4 "Gaussian" 5 "Continuous" 6 "Discontinuous" TEST02 P00_F evaluates the function. Here, we use spatial dimension M = 2 The number of points is N = 10 Problem 1 C parameters: 1: 0.75805147 2: 0.74194853 W parameters: 1: 0.12310392 2: 0.75123641E-02 F(X) X(1) X(2) ... -0.777207E-01 0.2184 0.9563 -0.245706 0.8295 0.5617 0.419987 0.4153 0.0661 0.497290 0.2576 0.1100 0.289512 0.0438 0.6340 0.405005 0.0617 0.4495 -0.667784E-01 0.4013 0.7547 0.190394 0.7973 0.0018 -0.142793 0.8975 0.3508 0.656024 0.0945 0.0136 Problem 2 C parameters: 1: 0.44390510 2: 1.5560949 W parameters: 1: 0.11366405 2: 0.35162866 F(X) X(1) X(2) ... 0.252530 0.2184 0.9563 0.391548 0.8295 0.5617 0.391472 0.4153 0.0661 0.416329 0.2576 0.1100 0.399564 0.0438 0.6340 0.466075 0.0617 0.4495 0.336953 0.4013 0.7547 0.337053 0.7973 0.0018 0.425617 0.8975 0.3508 0.373722 0.0945 0.0136 Problem 3 C parameters: 1: 1.3966185 2: 0.45338149 W parameters: 1: 0.69206650 2: 0.56166247 F(X) X(1) X(2) ... 0.190276 0.2184 0.9563 0.711600E-01 0.8295 0.5617 0.239618 0.4153 0.0661 0.357044 0.2576 0.1100 0.407672 0.0438 0.6340 0.465810 0.0617 0.4495 0.145191 0.4013 0.7547 0.105798 0.7973 0.0018 0.712198E-01 0.8975 0.3508 0.678149 0.0945 0.0136 Problem 4 C parameters: 1: 4.6040327 2: 2.4259673 W parameters: 1: 0.91197703 2: 0.59791688 F(X) X(1) X(2) ... 0.175178E-04 0.2184 0.9563 0.859090 0.8295 0.5617 0.101451E-02 0.4153 0.0661 0.281281E-04 0.2576 0.1100 0.114406E-06 0.0438 0.6340 0.194372E-06 0.0617 0.4495 0.343926E-02 0.4013 0.7547 0.934870E-01 0.7973 0.0018 0.694908 0.8975 0.3508 0.946477E-07 0.0945 0.0136 Problem 5 C parameters: 1: 4.0556462 2: 16.344354 W parameters: 1: 0.39698848 2: 0.18531412 F(X) X(1) X(2) ... 0.163189E-05 0.2184 0.9563 0.368577E-03 0.8295 0.5617 0.132329 0.4153 0.0661 0.165784 0.2576 0.1100 0.156069E-03 0.0438 0.6340 0.341942E-02 0.0617 0.4495 0.893167E-04 0.4013 0.7547 0.983054E-02 0.7973 0.0018 0.879216E-02 0.8975 0.3508 0.177231E-01 0.0945 0.0136 Problem 6 C parameters: 1: 2.6235318 2: 1.6764682 W parameters: 1: 0.61720483 2: 0.36152870 F(X) X(1) X(2) ... 0.00000 0.2184 0.9563 0.00000 0.8295 0.5617 3.32150 0.4153 0.0661 2.36340 0.2576 0.1100 0.00000 0.0438 0.6340 0.00000 0.0617 0.4495 0.00000 0.4013 0.7547 0.00000 0.7973 0.0018 0.00000 0.8975 0.3508 1.31110 0.0945 0.0136 TEST02 P00_F evaluates the function. Here, we use spatial dimension M = 3 The number of points is N = 10 Problem 1 C parameters: 1: 0.49084012 2: 0.39835258 3: 0.61080729 W parameters: 1: 0.45379378 2: 0.91197703 3: 0.59791688 F(X) X(1) X(2) ... -0.761930 0.2184 0.9563 0.8295 -0.981776 0.5617 0.4153 0.0661 -0.995649 0.2576 0.1100 0.0438 -0.949228 0.6340 0.0617 0.4495 -0.768521 0.4013 0.7547 0.7973 -0.960405 0.0018 0.8975 0.3508 -0.959309 0.0945 0.0136 0.8591 -0.984547 0.8408 0.1231 0.0075 -0.963674 0.2603 0.9125 0.1137 -0.931152 0.3516 0.8229 0.2671 Problem 2 C parameters: 1: 0.42069857 2: 1.6954256 3: 0.88387583 W parameters: 1: 0.18531412 2: 0.57436586 3: 0.36702667 F(X) X(1) X(2) ... 0.239884 0.2184 0.9563 0.8295 0.337565 0.5617 0.4153 0.0661 0.226626 0.2576 0.1100 0.0438 0.217469 0.6340 0.0617 0.4495 0.314955 0.4013 0.7547 0.7973 0.303822 0.0018 0.8975 0.3508 0.175298 0.0945 0.0136 0.8591 0.211614 0.8408 0.1231 0.0075 0.284576 0.2603 0.9125 0.1137 0.333284 0.3516 0.8229 0.2671 Problem 3 C parameters: 1: 0.95818065 2: 0.56125583 3: 0.33056352 W parameters: 1: 0.71447121 2: 0.11770687 3: 0.29932915 F(X) X(1) X(2) ... 0.600341E-01 0.2184 0.9563 0.8295 0.967226E-01 0.5617 0.4153 0.0661 0.326400 0.2576 0.1100 0.0438 0.972543E-01 0.6340 0.0617 0.4495 0.542926E-01 0.4013 0.7547 0.7973 0.144677 0.0018 0.8975 0.3508 0.273963 0.0945 0.0136 0.8591 0.805198E-01 0.8408 0.1231 0.0075 0.954448E-01 0.2603 0.9125 0.1137 0.788570E-01 0.3516 0.8229 0.2671 Problem 4 C parameters: 1: 3.3886572 2: 3.3872488 3: 0.25409404 W parameters: 1: 0.71078052 2: 0.88283334E-01 3: 0.77799401 F(X) X(1) X(2) ... 0.108761E-04 0.2184 0.9563 0.8295 0.219815 0.5617 0.4153 0.0661 0.908348E-01 0.2576 0.1100 0.0438 0.920524 0.6340 0.0617 0.4495 0.203984E-02 0.4013 0.7547 0.7973 0.168056E-05 0.0018 0.8975 0.3508 0.119740E-01 0.0945 0.0136 0.8591 0.781530 0.8408 0.1231 0.0075 0.389696E-04 0.2603 0.9125 0.1137 0.457556E-03 0.3516 0.8229 0.2671 Problem 5 C parameters: 1: 7.7836405 2: 3.2236747 3: 9.3926848 W parameters: 1: 0.76353672 2: 0.76173065 3: 0.40696964 F(X) X(1) X(2) ... 0.144957E-03 0.2184 0.9563 0.8295 0.276885E-02 0.5617 0.4153 0.0661 0.786758E-04 0.2576 0.1100 0.0438 0.256053E-01 0.6340 0.0617 0.4495 0.149094E-02 0.4013 0.7547 0.7973 0.101334E-02 0.0018 0.8975 0.3508 0.702836E-05 0.0945 0.0136 0.8591 0.164098E-02 0.8408 0.1231 0.0075 0.778691E-03 0.2603 0.9125 0.1137 0.894427E-02 0.3516 0.8229 0.2671 Problem 6 C parameters: 1: 2.6580200 2: 1.5915236 3: 0.50456459E-01 W parameters: 1: 0.50110323 2: 0.41909255E-01 3: 0.36885058 F(X) X(1) X(2) ... 0.00000 0.2184 0.9563 0.8295 0.00000 0.5617 0.4153 0.0661 0.00000 0.2576 0.1100 0.0438 0.00000 0.6340 0.0617 0.4495 0.00000 0.4013 0.7547 0.7973 0.00000 0.0018 0.8975 0.3508 1.37207 0.0945 0.0136 0.8591 0.00000 0.8408 0.1231 0.0075 0.00000 0.2603 0.9125 0.1137 0.00000 0.3516 0.8229 0.2671 TEST02 P00_F evaluates the function. Here, we use spatial dimension M = 4 The number of points is N = 10 Problem 1 C parameters: 1: 0.44869387 2: 0.28672076 3: 0.48215961 4: 0.28242576 W parameters: 1: 0.21293000 2: 0.71447121 3: 0.11770687 4: 0.29932915 F(X) X(1) X(2) ... -0.642591 0.2184 0.9563 0.8295 0.5617 -0.127288 0.4153 0.0661 0.2576 0.1100 -0.124916 0.0438 0.6340 0.0617 0.4495 -0.521378 0.4013 0.7547 0.7973 0.0018 -0.314364 0.8975 0.3508 0.0945 0.0136 -0.439570 0.8591 0.8408 0.1231 0.0075 -0.295157 0.2603 0.9125 0.1137 0.3516 -0.648197 0.8229 0.2671 0.6921 0.5617 -0.778456 0.8612 0.4538 0.9120 0.5979 -0.308818 0.1890 0.7615 0.3970 0.1853 Problem 2 C parameters: 1: 1.3623435 2: 1.3617773 3: 0.10215355 4: 1.1737257 W parameters: 1: 0.88283334E-01 2: 0.77799401 3: 0.74530307 4: 0.30867492 F(X) X(1) X(2) ... 0.416258E-01 0.2184 0.9563 0.8295 0.5617 0.201354E-01 0.4153 0.0661 0.2576 0.1100 0.459851E-01 0.0438 0.6340 0.0617 0.4495 0.370209E-01 0.4013 0.7547 0.7973 0.0018 0.148338E-01 0.8975 0.3508 0.0945 0.0136 0.206817E-01 0.8591 0.8408 0.1231 0.0075 0.450787E-01 0.2603 0.9125 0.1137 0.3516 0.153080E-01 0.8229 0.2671 0.6921 0.5617 0.176016E-01 0.8612 0.4538 0.9120 0.5979 0.474846E-01 0.1890 0.7615 0.3970 0.1853 Problem 3 C parameters: 1: 0.58759510 2: 0.49884797 3: 0.49766799 4: 0.26588895 W parameters: 1: 0.93874945 2: 0.56208829 3: 0.17820022E-01 4: 0.50110323 F(X) X(1) X(2) ... 0.208996E-01 0.2184 0.9563 0.8295 0.5617 0.164660 0.4153 0.0661 0.2576 0.1100 0.135141 0.0438 0.6340 0.0617 0.4495 0.305148E-01 0.4013 0.7547 0.7973 0.0018 0.604051E-01 0.8975 0.3508 0.0945 0.0136 0.322436E-01 0.8591 0.8408 0.1231 0.0075 0.595186E-01 0.2603 0.9125 0.1137 0.3516 0.238798E-01 0.8229 0.2671 0.6921 0.5617 0.140943E-01 0.8612 0.4538 0.9120 0.5979 0.631070E-01 0.1890 0.7615 0.3970 0.1853 Problem 4 C parameters: 1: 0.19118193 2: 1.6826252 3: 1.2395506 4: 3.9166423 W parameters: 1: 0.29036559E-01 2: 0.17442279E-01 3: 0.15238379 4: 0.11431867 F(X) X(1) X(2) ... 0.188895E-02 0.2184 0.9563 0.8295 0.5617 0.970974 0.4153 0.0661 0.2576 0.1100 0.600490E-01 0.0438 0.6340 0.0617 0.4495 0.928314E-01 0.4013 0.7547 0.7973 0.0018 0.604831 0.8975 0.3508 0.0945 0.0136 0.119905 0.8591 0.8408 0.1231 0.0075 0.434456E-01 0.2603 0.9125 0.1137 0.3516 0.243101E-01 0.8229 0.2671 0.6921 0.5617 0.648377E-02 0.8612 0.4538 0.9120 0.5979 0.175946 0.1890 0.7615 0.3970 0.1853 Problem 5 C parameters: 1: 8.0866648 2: 2.7261447 3: 4.7219478 4: 4.8652426 W parameters: 1: 0.61294755 2: 0.80951907 3: 0.58708963 4: 0.21549164 F(X) X(1) X(2) ... 0.162918E-02 0.2184 0.9563 0.8295 0.5617 0.336523E-02 0.4153 0.0661 0.2576 0.1100 0.166538E-03 0.0438 0.6340 0.0617 0.4495 0.203840E-01 0.4013 0.7547 0.7973 0.0018 0.104921E-02 0.8975 0.3508 0.0945 0.0136 0.509896E-02 0.8591 0.8408 0.1231 0.0075 0.240490E-02 0.2603 0.9125 0.1137 0.3516 0.471863E-02 0.8229 0.2671 0.6921 0.5617 0.170854E-02 0.8612 0.4538 0.9120 0.5979 0.100109E-01 0.1890 0.7615 0.3970 0.1853 Problem 6 C parameters: 1: 1.1818162 2: 1.1129441 3: 0.68937142 4: 1.3158683 W parameters: 1: 0.94501750 2: 0.90905692 3: 0.51972572 4: 0.30194625E-01 F(X) X(1) X(2) ... 0.00000 0.2184 0.9563 0.8295 0.5617 2.42701 0.4153 0.0661 0.2576 0.1100 4.02069 0.0438 0.6340 0.0617 0.4495 6.46385 0.4013 0.7547 0.7973 0.0018 4.63737 0.8975 0.3508 0.0945 0.0136 7.73581 0.8591 0.8408 0.1231 0.0075 0.00000 0.2603 0.9125 0.1137 0.3516 12.0128 0.8229 0.2671 0.6921 0.5617 18.8840 0.8612 0.4538 0.9120 0.5979 4.89556 0.1890 0.7615 0.3970 0.1853 TEST03 P00_D evaluates derivative components. Here, we use spatial dimension M = 2 The number of points is N = 2 Problem 1 C parameters: 1: 1.2939909 2: 0.20600911 W parameters: 1: 0.25757779 2: 0.10995679 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 -0.503162 -1.1183 -0.1780 0.8295 0.5617 -0.944708 -0.4243 -0.0676 Problem 2 C parameters: 1: 0.12932824 2: 1.8706718 W parameters: 1: 0.61727229E-01 2: 0.44953896 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.308134E-01 -0.0002 -0.0576 0.8295 0.5617 0.555153E-01 -0.0014 -0.0417 Problem 3 C parameters: 1: 0.64224015 2: 1.2077598 W parameters: 1: 0.79728695 2: 0.18383712E-02 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.826977E-01 -0.0694 -0.1305 0.8295 0.5617 0.925024E-01 -0.0806 -0.1516 Problem 4 C parameters: 1: 5.0546134 2: 1.9753866 W parameters: 1: 0.94544750E-01 2: 0.13616892E-01 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.210717E-01 -0.1334 -0.1550 0.8295 0.5617 0.314241E-06 -0.0000 -0.0000 Problem 5 C parameters: 1: 10.309500 2: 10.090500 W parameters: 1: 0.12310392 2: 0.75123641E-02 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.260221E-04 -10.3092 -10.0902 0.8295 0.5617 0.256213E-05 -10.3092 -10.0902 Problem 6 C parameters: 1: 0.95439596 2: 3.3456040 W parameters: 1: 0.11366405 2: 0.35162866 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.00000 0.0000 0.0000 0.8295 0.5617 0.00000 0.0000 0.0000 TEST03 P00_D evaluates derivative components. Here, we use spatial dimension M = 3 The number of points is N = 2 Problem 1 C parameters: 1: 0.93923406 2: 0.40094748 3: 0.15981847 W parameters: 1: 0.63396571 2: 0.61727229E-01 3: 0.44953896 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.8295 -0.791493E-02 0.9392 0.4009 0.1598 0.5617 0.4153 0.0661 -0.244157E-01 0.9390 0.4008 0.1598 Problem 2 C parameters: 1: 0.61636172 2: 1.1590944 3: 1.2245439 W parameters: 1: 0.18383712E-02 2: 0.89750406 3: 0.35075234 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.8295 0.557020 -0.0901 -0.0876 -0.5952 0.5617 0.4153 0.0661 0.464670 -0.1766 0.4587 0.3537 Problem 3 C parameters: 1: 0.18082838 2: 0.26043968E-01 3: 1.6431277 W parameters: 1: 0.84084745 2: 0.12310392 3: 0.75123641E-02 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.8295 0.288032E-01 -0.0086 -0.0012 -0.0780 0.5617 0.4153 0.0661 0.449880 -0.2665 -0.0384 -2.4216 Problem 4 C parameters: 1: 1.4224641 2: 4.9864019 3: 0.62113396 W parameters: 1: 0.35162866 2: 0.82288732 3: 0.26713227 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.8295 0.548484 0.2957 -3.6393 -0.2380 0.5617 0.4153 0.0661 0.144746E-01 -0.0123 0.2934 0.0022 Problem 5 C parameters: 1: 6.6754257 2: 5.4175951 3: 8.3069792 W parameters: 1: 0.45379378 2: 0.91197703 3: 0.59791688 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.8295 0.238664E-01 6.5180 -5.2898 -8.1111 0.5617 0.4153 0.0661 0.398150E-03 -6.5180 5.2898 8.1111 Problem 6 C parameters: 1: 0.60300129 2: 2.4301100 3: 1.2668887 W parameters: 1: 0.18531412 2: 0.57436586 3: 0.36702667 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.8295 0.00000 0.0000 0.0000 0.0000 0.5617 0.4153 0.0661 0.00000 0.0000 0.0000 0.0000 TEST03 P00_D evaluates derivative components. Here, we use spatial dimension M = 4 The number of points is N = 2 Problem 1 C parameters: 1: 0.55290269E-01 2: 0.79974757 3: 0.77868883E-01 4: 0.56709328 W parameters: 1: 0.40130628 2: 0.75467349 3: 0.79728695 4: 0.18383712E-02 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.8295 0.5617 -0.857758 0.0284 0.4111 0.0400 0.2915 0.4153 0.0661 0.2576 0.1100 -0.895226 -0.0246 -0.3564 -0.0347 -0.2527 Problem 2 C parameters: 1: 2.6466887 2: 1.0343488 3: 0.27880712 4: 0.40155442E-01 W parameters: 1: 0.85909686 2: 0.84084745 3: 0.12310392 4: 0.75123641E-02 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.8295 0.5617 0.229951E-03 0.0005 -0.0001 -0.0000 -0.0000 0.4153 0.0661 0.2576 0.1100 0.240049E-03 0.0006 0.0002 -0.0000 -0.0000 Problem 3 C parameters: 1: 0.29397876 2: 1.0305330 3: 0.12836892 4: 0.39711934 W parameters: 1: 0.82288732 2: 0.26713227 3: 0.69206650 4: 0.56166247 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.8295 0.5617 0.131154E-01 -0.0081 -0.0284 -0.0035 -0.0109 0.4153 0.0661 0.2576 0.1100 0.306327 -0.3554 -1.2458 -0.1552 -0.4801 Problem 4 C parameters: 1: 2.1432049 2: 1.1293024 3: 2.2695283 4: 1.4879643 W parameters: 1: 0.18895469 2: 0.76149206 3: 0.39698848 4: 0.18531412 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.8295 0.5617 0.264579 -0.0716 -0.1315 -1.1789 -0.4410 0.4153 0.0661 0.2576 0.1100 0.381099 -0.7925 0.6759 0.5473 0.1272 Problem 5 C parameters: 1: 6.1022366 2: 3.8994023 3: 6.5573708 4: 3.8409903 W parameters: 1: 0.21293000 2: 0.71447121 3: 0.11770687 4: 0.29932915 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.8295 0.5617 0.129164E-02 -6.0944 -3.8944 -6.5489 -3.8360 0.4153 0.0661 0.2576 0.1100 0.448192E-02 -6.0944 3.8944 -6.5489 3.8360 Problem 6 C parameters: 1: 1.4645192 2: 1.4639105 3: 0.10981506 4: 1.2617552 W parameters: 1: 0.88283334E-01 2: 0.77799401 3: 0.74530307 4: 0.30867492 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.8295 0.5617 0.00000 0.0000 0.0000 0.0000 0.0000 0.4153 0.0661 0.2576 0.1100 0.00000 0.0000 0.0000 0.0000 0.0000 TEST04 P00_Q returns the integral of F over [0,1]^m. Here, we use spatial dimension M = 4 The number of sample points is N = 10000 Problem 1 C parameters: 1: 0.32229788 2: 0.30680667 3: 0.19249538 4: 0.67840007 W parameters: 1: 0.48064731 2: 0.23929848 3: 0.88949426 4: 0.72995386 Exact Integral Q -0.785803 -0.787780 Problem 2 C parameters: 1: 0.66541693 2: 1.5543428 3: 1.7506064 4: 0.29633872E-01 W parameters: 1: 0.36205422 2: 0.45308588E-01 3: 0.50144009 4: 0.70362217 Exact Integral Q 0.152845E-02 0.153381E-02 Problem 3 C parameters: 1: 0.61038748 2: 0.38911689 3: 0.19127706 4: 0.65921856 W parameters: 1: 0.37721981 2: 0.93342259 3: 0.33471966E-01 4: 0.56333555 Exact Integral Q 0.550631E-01 0.557084E-01 Problem 4 C parameters: 1: 2.6937360 2: 0.99450918 3: 2.2187831 4: 1.1229717 W parameters: 1: 0.79285630 2: 0.53580537 3: 0.28085320 4: 0.29973860 Exact Integral Q 0.263474 0.264295 Problem 5 C parameters: 1: 7.6626898 2: 0.87841454 3: 5.7733081 4: 6.0855876 W parameters: 1: 0.59076919 2: 0.57817994E-01 3: 0.74702815 4: 0.30216421 Exact Integral Q 0.161490E-01 0.165698E-01 Problem 6 C parameters: 1: 0.92677675 2: 1.3376846 3: 0.88767315 4: 1.1478655 W parameters: 1: 0.95082864 2: 0.57687031 3: 0.45928323 4: 0.17320148 Exact Integral Q 4.00448 4.01973 TEST_INTERP_ND_PRB Normal end of execution. 5 October 2012 9:44:19.280 AM