29-Mar-2019 12:09:01 test_interp_nd_test MATLAB version Test test_interp_nd. The R8LIB library is also needed. TEST_INTERP_ND_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" TEST_INTERP_ND_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.758051 2: 0.741949 W parameters: 1: 0.123104 2: 0.00751236 F(X) X(1) X(2) ... -0.0777207 0.2184 0.9563 -0.245706 0.8295 0.5617 0.419987 0.4153 0.0661 0.49729 0.2576 0.1100 0.289512 0.0438 0.6340 0.405005 0.0617 0.4495 -0.0667784 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.443905 2: 1.55609 W parameters: 1: 0.113664 2: 0.351629 F(X) X(1) X(2) ... 0.25253 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.39662 2: 0.453381 W parameters: 1: 0.692066 2: 0.561662 F(X) X(1) X(2) ... 0.190276 0.2184 0.9563 0.07116 0.8295 0.5617 0.239618 0.4153 0.0661 0.357044 0.2576 0.1100 0.407672 0.0438 0.6340 0.46581 0.0617 0.4495 0.145191 0.4013 0.7547 0.105798 0.7973 0.0018 0.0712198 0.8975 0.3508 0.678149 0.0945 0.0136 Problem 4 C parameters: 1: 4.60403 2: 2.42597 W parameters: 1: 0.911977 2: 0.597917 F(X) X(1) X(2) ... 1.75178e-05 0.2184 0.9563 0.85909 0.8295 0.5617 0.00101451 0.4153 0.0661 2.81281e-05 0.2576 0.1100 1.14406e-07 0.0438 0.6340 1.94372e-07 0.0617 0.4495 0.00343926 0.4013 0.7547 0.093487 0.7973 0.0018 0.694908 0.8975 0.3508 9.46477e-08 0.0945 0.0136 Problem 5 C parameters: 1: 4.05565 2: 16.3444 W parameters: 1: 0.396988 2: 0.185314 F(X) X(1) X(2) ... 1.63189e-06 0.2184 0.9563 0.000368577 0.8295 0.5617 0.132329 0.4153 0.0661 0.165784 0.2576 0.1100 0.000156069 0.0438 0.6340 0.00341942 0.0617 0.4495 8.93167e-05 0.4013 0.7547 0.00983054 0.7973 0.0018 0.00879216 0.8975 0.3508 0.0177231 0.0945 0.0136 Problem 6 C parameters: 1: 2.62353 2: 1.67647 W parameters: 1: 0.617205 2: 0.361529 F(X) X(1) X(2) ... 0 0.2184 0.9563 0 0.8295 0.5617 3.3215 0.4153 0.0661 2.3634 0.2576 0.1100 0 0.0438 0.6340 0 0.0617 0.4495 0 0.4013 0.7547 0 0.7973 0.0018 0 0.8975 0.3508 1.3111 0.0945 0.0136 TEST_INTERP_ND_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.49084 2: 0.398353 3: 0.610807 W parameters: 1: 0.453794 2: 0.911977 3: 0.597917 F(X) X(1) X(2) ... -0.76193 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.420699 2: 1.69543 3: 0.883876 W parameters: 1: 0.185314 2: 0.574366 3: 0.367027 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.958181 2: 0.561256 3: 0.330564 W parameters: 1: 0.714471 2: 0.117707 3: 0.299329 F(X) X(1) X(2) ... 0.0600341 0.2184 0.9563 0.8295 0.0967226 0.5617 0.4153 0.0661 0.3264 0.2576 0.1100 0.0438 0.0972543 0.6340 0.0617 0.4495 0.0542926 0.4013 0.7547 0.7973 0.144677 0.0018 0.8975 0.3508 0.273963 0.0945 0.0136 0.8591 0.0805198 0.8408 0.1231 0.0075 0.0954448 0.2603 0.9125 0.1137 0.078857 0.3516 0.8229 0.2671 Problem 4 C parameters: 1: 3.38866 2: 3.38725 3: 0.254094 W parameters: 1: 0.710781 2: 0.0882833 3: 0.777994 F(X) X(1) X(2) ... 1.08761e-05 0.2184 0.9563 0.8295 0.219815 0.5617 0.4153 0.0661 0.0908348 0.2576 0.1100 0.0438 0.920524 0.6340 0.0617 0.4495 0.00203984 0.4013 0.7547 0.7973 1.68056e-06 0.0018 0.8975 0.3508 0.011974 0.0945 0.0136 0.8591 0.78153 0.8408 0.1231 0.0075 3.89696e-05 0.2603 0.9125 0.1137 0.000457556 0.3516 0.8229 0.2671 Problem 5 C parameters: 1: 7.78364 2: 3.22367 3: 9.39268 W parameters: 1: 0.763537 2: 0.761731 3: 0.40697 F(X) X(1) X(2) ... 0.000144957 0.2184 0.9563 0.8295 0.00276885 0.5617 0.4153 0.0661 7.86758e-05 0.2576 0.1100 0.0438 0.0256053 0.6340 0.0617 0.4495 0.00149094 0.4013 0.7547 0.7973 0.00101334 0.0018 0.8975 0.3508 7.02836e-06 0.0945 0.0136 0.8591 0.00164098 0.8408 0.1231 0.0075 0.000778691 0.2603 0.9125 0.1137 0.00894427 0.3516 0.8229 0.2671 Problem 6 C parameters: 1: 2.65802 2: 1.59152 3: 0.0504565 W parameters: 1: 0.501103 2: 0.0419093 3: 0.368851 F(X) X(1) X(2) ... 0 0.2184 0.9563 0.8295 0 0.5617 0.4153 0.0661 0 0.2576 0.1100 0.0438 0 0.6340 0.0617 0.4495 0 0.4013 0.7547 0.7973 0 0.0018 0.8975 0.3508 1.37207 0.0945 0.0136 0.8591 0 0.8408 0.1231 0.0075 0 0.2603 0.9125 0.1137 0 0.3516 0.8229 0.2671 TEST_INTERP_ND_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.448694 2: 0.286721 3: 0.48216 4: 0.282426 W parameters: 1: 0.21293 2: 0.714471 3: 0.117707 4: 0.299329 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.43957 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.36234 2: 1.36178 3: 0.102154 4: 1.17373 W parameters: 1: 0.0882833 2: 0.777994 3: 0.745303 4: 0.308675 F(X) X(1) X(2) ... 0.0416258 0.2184 0.9563 0.8295 0.5617 0.0201354 0.4153 0.0661 0.2576 0.1100 0.0459851 0.0438 0.6340 0.0617 0.4495 0.0370209 0.4013 0.7547 0.7973 0.0018 0.0148338 0.8975 0.3508 0.0945 0.0136 0.0206817 0.8591 0.8408 0.1231 0.0075 0.0450787 0.2603 0.9125 0.1137 0.3516 0.015308 0.8229 0.2671 0.6921 0.5617 0.0176016 0.8612 0.4538 0.9120 0.5979 0.0474846 0.1890 0.7615 0.3970 0.1853 Problem 3 C parameters: 1: 0.587595 2: 0.498848 3: 0.497668 4: 0.265889 W parameters: 1: 0.938749 2: 0.562088 3: 0.01782 4: 0.501103 F(X) X(1) X(2) ... 0.0208996 0.2184 0.9563 0.8295 0.5617 0.16466 0.4153 0.0661 0.2576 0.1100 0.135141 0.0438 0.6340 0.0617 0.4495 0.0305148 0.4013 0.7547 0.7973 0.0018 0.0604051 0.8975 0.3508 0.0945 0.0136 0.0322436 0.8591 0.8408 0.1231 0.0075 0.0595186 0.2603 0.9125 0.1137 0.3516 0.0238798 0.8229 0.2671 0.6921 0.5617 0.0140943 0.8612 0.4538 0.9120 0.5979 0.063107 0.1890 0.7615 0.3970 0.1853 Problem 4 C parameters: 1: 0.191182 2: 1.68263 3: 1.23955 4: 3.91664 W parameters: 1: 0.0290366 2: 0.0174423 3: 0.152384 4: 0.114319 F(X) X(1) X(2) ... 0.00188895 0.2184 0.9563 0.8295 0.5617 0.970974 0.4153 0.0661 0.2576 0.1100 0.060049 0.0438 0.6340 0.0617 0.4495 0.0928314 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.0434456 0.2603 0.9125 0.1137 0.3516 0.0243101 0.8229 0.2671 0.6921 0.5617 0.00648377 0.8612 0.4538 0.9120 0.5979 0.175946 0.1890 0.7615 0.3970 0.1853 Problem 5 C parameters: 1: 8.08666 2: 2.72614 3: 4.72195 4: 4.86524 W parameters: 1: 0.612948 2: 0.809519 3: 0.58709 4: 0.215492 F(X) X(1) X(2) ... 0.00162918 0.2184 0.9563 0.8295 0.5617 0.00336523 0.4153 0.0661 0.2576 0.1100 0.000166538 0.0438 0.6340 0.0617 0.4495 0.020384 0.4013 0.7547 0.7973 0.0018 0.00104921 0.8975 0.3508 0.0945 0.0136 0.00509896 0.8591 0.8408 0.1231 0.0075 0.0024049 0.2603 0.9125 0.1137 0.3516 0.00471863 0.8229 0.2671 0.6921 0.5617 0.00170854 0.8612 0.4538 0.9120 0.5979 0.0100109 0.1890 0.7615 0.3970 0.1853 Problem 6 C parameters: 1: 1.18182 2: 1.11294 3: 0.689371 4: 1.31587 W parameters: 1: 0.945017 2: 0.909057 3: 0.519726 4: 0.0301946 F(X) X(1) X(2) ... 0 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 0.2603 0.9125 0.1137 0.3516 12.0128 0.8229 0.2671 0.6921 0.5617 18.884 0.8612 0.4538 0.9120 0.5979 4.89556 0.1890 0.7615 0.3970 0.1853 TEST_INTERP_ND_TEST03 P00_F evaluates derivative components. Here, we use spatial dimension M = 2 The number of points is N = 2 Problem 1 C parameters: 1: 1.29399 2: 0.206009 W parameters: 1: 0.257578 2: 0.109957 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.129328 2: 1.87067 W parameters: 1: 0.0617272 2: 0.449539 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.0308134 -0.0002 -0.0576 0.8295 0.5617 0.0555153 -0.0014 -0.0417 Problem 3 C parameters: 1: 0.64224 2: 1.20776 W parameters: 1: 0.797287 2: 0.00183837 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.0826977 -0.0694 -0.1305 0.8295 0.5617 0.0925024 -0.0806 -0.1516 Problem 4 C parameters: 1: 5.05461 2: 1.97539 W parameters: 1: 0.0945448 2: 0.0136169 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.0210717 -0.1334 -0.1550 0.8295 0.5617 3.14241e-07 -0.0000 -0.0000 Problem 5 C parameters: 1: 10.3095 2: 10.0905 W parameters: 1: 0.123104 2: 0.00751236 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 2.60221e-05 -0.0003 -0.0003 0.8295 0.5617 2.56213e-06 -0.0000 -0.0000 Problem 6 C parameters: 1: 0.954396 2: 3.3456 W parameters: 1: 0.113664 2: 0.351629 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0 0.0000 0.0000 0.8295 0.5617 0 0.0000 0.0000 TEST_INTERP_ND_TEST03 P00_F evaluates derivative components. Here, we use spatial dimension M = 3 The number of points is N = 2 Problem 1 C parameters: 1: 0.939234 2: 0.400947 3: 0.159818 W parameters: 1: 0.633966 2: 0.0617272 3: 0.449539 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.8295 -0.00791493 0.9392 0.4009 0.1598 0.5617 0.4153 0.0661 -0.0244157 0.9390 0.4008 0.1598 Problem 2 C parameters: 1: 0.616362 2: 1.15909 3: 1.22454 W parameters: 1: 0.00183837 2: 0.897504 3: 0.350752 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.8295 0.55702 -0.0901 -0.0876 -0.5952 0.5617 0.4153 0.0661 0.46467 -0.1766 0.4587 0.3537 Problem 3 C parameters: 1: 0.180828 2: 0.026044 3: 1.64313 W parameters: 1: 0.840847 2: 0.123104 3: 0.00751236 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.8295 0.0288032 -0.0086 -0.0012 -0.0780 0.5617 0.4153 0.0661 0.44988 -0.2665 -0.0384 -2.4216 Problem 4 C parameters: 1: 1.42246 2: 4.9864 3: 0.621134 W parameters: 1: 0.351629 2: 0.822887 3: 0.267132 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.0144746 -0.0123 0.2934 0.0022 Problem 5 C parameters: 1: 6.67543 2: 5.4176 3: 8.30698 W parameters: 1: 0.453794 2: 0.911977 3: 0.597917 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.8295 0.0238664 0.1593 -0.1293 -0.1983 0.5617 0.4153 0.0661 0.00039815 -0.0027 0.0022 0.0033 Problem 6 C parameters: 1: 0.603001 2: 2.43011 3: 1.26689 W parameters: 1: 0.185314 2: 0.574366 3: 0.367027 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.8295 0 0.0000 0.0000 0.0000 0.5617 0.4153 0.0661 0 0.0000 0.0000 0.0000 TEST_INTERP_ND_TEST03 P00_F evaluates derivative components. Here, we use spatial dimension M = 4 The number of points is N = 2 Problem 1 C parameters: 1: 0.0552903 2: 0.799748 3: 0.0778689 4: 0.567093 W parameters: 1: 0.401306 2: 0.754673 3: 0.797287 4: 0.00183837 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.64669 2: 1.03435 3: 0.278807 4: 0.0401554 W parameters: 1: 0.859097 2: 0.840847 3: 0.123104 4: 0.00751236 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.8295 0.5617 0.000229951 0.0005 -0.0001 -0.0000 -0.0000 0.4153 0.0661 0.2576 0.1100 0.000240049 0.0006 0.0002 -0.0000 -0.0000 Problem 3 C parameters: 1: 0.293979 2: 1.03053 3: 0.128369 4: 0.397119 W parameters: 1: 0.822887 2: 0.267132 3: 0.692066 4: 0.561662 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.8295 0.5617 0.0131154 -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.1432 2: 1.1293 3: 2.26953 4: 1.48796 W parameters: 1: 0.188955 2: 0.761492 3: 0.396988 4: 0.185314 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.10224 2: 3.8994 3: 6.55737 4: 3.84099 W parameters: 1: 0.21293 2: 0.714471 3: 0.117707 4: 0.299329 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.8295 0.5617 0.00129164 -0.0079 -0.0050 -0.0085 -0.0050 0.4153 0.0661 0.2576 0.1100 0.00448192 -0.0273 0.0175 -0.0294 0.0172 Problem 6 C parameters: 1: 1.46452 2: 1.46391 3: 0.109815 4: 1.26176 W parameters: 1: 0.0882833 2: 0.777994 3: 0.745303 4: 0.308675 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.2184 0.9563 0.8295 0.5617 0 0.0000 0.0000 0.0000 0.0000 0.4153 0.0661 0.2576 0.1100 0 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.322298 2: 0.306807 3: 0.192495 4: 0.6784 W parameters: 1: 0.480647 2: 0.239298 3: 0.889494 4: 0.729954 Exact Integral Q -0.785803 -0.78778 Problem 2 C parameters: 1: 0.665417 2: 1.55434 3: 1.75061 4: 0.0296339 W parameters: 1: 0.362054 2: 0.0453086 3: 0.50144 4: 0.703622 Exact Integral Q 0.00152845 0.00153381 Problem 3 C parameters: 1: 0.610387 2: 0.389117 3: 0.191277 4: 0.659219 W parameters: 1: 0.37722 2: 0.933423 3: 0.033472 4: 0.563336 Exact Integral Q 0.0550631 0.0557084 Problem 4 C parameters: 1: 2.69374 2: 0.994509 3: 2.21878 4: 1.12297 W parameters: 1: 0.792856 2: 0.535805 3: 0.280853 4: 0.299739 Exact Integral Q 0.263474 0.264295 Problem 5 C parameters: 1: 7.66269 2: 0.878415 3: 5.77331 4: 6.08559 W parameters: 1: 0.590769 2: 0.057818 3: 0.747028 4: 0.302164 Exact Integral Q 0.016149 0.0165698 Problem 6 C parameters: 1: 0.926777 2: 1.33768 3: 0.887673 4: 1.14787 W parameters: 1: 0.950829 2: 0.57687 3: 0.459283 4: 0.173201 Exact Integral Q 4.00448 4.01973 test_interp_nd_test Normal end of execution. 29-Mar-2019 12:09:02