27 February 2022 11:15:14 AM TEST_INTERP_ND_TEST C++ 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: 0: 0.758051 1: 0.741949 W parameters: 0: 0.123104 1: 0.00751236 F(X) X(1) X(2) ... -0.0777207 0.218418 0.956318 -0.245706 0.829509 0.561695 0.419987 0.415307 0.0661187 0.49729 0.257578 0.109957 0.289512 0.043829 0.633966 0.405005 0.0617272 0.449539 -0.0667784 0.401306 0.754673 0.190394 0.797287 0.00183837 -0.142793 0.897504 0.350752 0.656024 0.0945448 0.0136169 Problem 2 C parameters: 0: 0.443905 1: 1.55609 W parameters: 0: 0.113664 1: 0.351629 F(X) X(1) X(2) ... 0.25253 0.218418 0.956318 0.391548 0.829509 0.561695 0.391472 0.415307 0.0661187 0.416329 0.257578 0.109957 0.399564 0.043829 0.633966 0.466075 0.0617272 0.449539 0.336953 0.401306 0.754673 0.337053 0.797287 0.00183837 0.425617 0.897504 0.350752 0.373722 0.0945448 0.0136169 Problem 3 C parameters: 0: 1.39662 1: 0.453381 W parameters: 0: 0.692066 1: 0.561662 F(X) X(1) X(2) ... 0.190276 0.218418 0.956318 0.07116 0.829509 0.561695 0.239618 0.415307 0.0661187 0.357044 0.257578 0.109957 0.407672 0.043829 0.633966 0.46581 0.0617272 0.449539 0.145191 0.401306 0.754673 0.105798 0.797287 0.00183837 0.0712198 0.897504 0.350752 0.678149 0.0945448 0.0136169 Problem 4 C parameters: 0: 4.60403 1: 2.42597 W parameters: 0: 0.911977 1: 0.597917 F(X) X(1) X(2) ... 1.75178e-05 0.218418 0.956318 0.85909 0.829509 0.561695 0.00101451 0.415307 0.0661187 2.81281e-05 0.257578 0.109957 1.14406e-07 0.043829 0.633966 1.94372e-07 0.0617272 0.449539 0.00343926 0.401306 0.754673 0.093487 0.797287 0.00183837 0.694908 0.897504 0.350752 9.46477e-08 0.0945448 0.0136169 Problem 5 C parameters: 0: 4.05565 1: 16.3444 W parameters: 0: 0.396988 1: 0.185314 F(X) X(1) X(2) ... 1.63189e-06 0.218418 0.956318 0.000368577 0.829509 0.561695 0.132329 0.415307 0.0661187 0.165784 0.257578 0.109957 0.000156069 0.043829 0.633966 0.00341942 0.0617272 0.449539 8.93167e-05 0.401306 0.754673 0.00983054 0.797287 0.00183837 0.00879216 0.897504 0.350752 0.0177231 0.0945448 0.0136169 Problem 6 C parameters: 0: 2.62353 1: 1.67647 W parameters: 0: 0.617205 1: 0.361529 F(X) X(1) X(2) ... 0 0.218418 0.956318 0 0.829509 0.561695 3.3215 0.415307 0.0661187 2.3634 0.257578 0.109957 0 0.043829 0.633966 0 0.0617272 0.449539 0 0.401306 0.754673 0 0.797287 0.00183837 0 0.897504 0.350752 1.3111 0.0945448 0.0136169 TEST02 P00_F evaluates the function. Here, we use spatial dimension M = 3 The number of points is N = 10 Problem 1 C parameters: 0: 0.49084 1: 0.398353 2: 0.610807 W parameters: 0: 0.453794 1: 0.911977 2: 0.597917 F(X) X(1) X(2) ... -0.76193 0.218418 0.956318 0.829509 -0.981776 0.561695 0.415307 0.0661187 -0.995649 0.257578 0.109957 0.043829 -0.949228 0.633966 0.0617272 0.449539 -0.768521 0.401306 0.754673 0.797287 -0.960405 0.00183837 0.897504 0.350752 -0.959309 0.0945448 0.0136169 0.859097 -0.984547 0.840847 0.123104 0.00751236 -0.963674 0.260303 0.912484 0.113664 -0.931152 0.351629 0.822887 0.267132 Problem 2 C parameters: 0: 0.420699 1: 1.69543 2: 0.883876 W parameters: 0: 0.185314 1: 0.574366 2: 0.367027 F(X) X(1) X(2) ... 0.239884 0.218418 0.956318 0.829509 0.337565 0.561695 0.415307 0.0661187 0.226626 0.257578 0.109957 0.043829 0.217469 0.633966 0.0617272 0.449539 0.314955 0.401306 0.754673 0.797287 0.303822 0.00183837 0.897504 0.350752 0.175298 0.0945448 0.0136169 0.859097 0.211614 0.840847 0.123104 0.00751236 0.284576 0.260303 0.912484 0.113664 0.333284 0.351629 0.822887 0.267132 Problem 3 C parameters: 0: 0.958181 1: 0.561256 2: 0.330564 W parameters: 0: 0.714471 1: 0.117707 2: 0.299329 F(X) X(1) X(2) ... 0.0600341 0.218418 0.956318 0.829509 0.0967226 0.561695 0.415307 0.0661187 0.3264 0.257578 0.109957 0.043829 0.0972543 0.633966 0.0617272 0.449539 0.0542926 0.401306 0.754673 0.797287 0.144677 0.00183837 0.897504 0.350752 0.273963 0.0945448 0.0136169 0.859097 0.0805198 0.840847 0.123104 0.00751236 0.0954448 0.260303 0.912484 0.113664 0.078857 0.351629 0.822887 0.267132 Problem 4 C parameters: 0: 3.38866 1: 3.38725 2: 0.254094 W parameters: 0: 0.710781 1: 0.0882833 2: 0.777994 F(X) X(1) X(2) ... 1.08761e-05 0.218418 0.956318 0.829509 0.219815 0.561695 0.415307 0.0661187 0.0908348 0.257578 0.109957 0.043829 0.920524 0.633966 0.0617272 0.449539 0.00203984 0.401306 0.754673 0.797287 1.68056e-06 0.00183837 0.897504 0.350752 0.011974 0.0945448 0.0136169 0.859097 0.78153 0.840847 0.123104 0.00751236 3.89696e-05 0.260303 0.912484 0.113664 0.000457556 0.351629 0.822887 0.267132 Problem 5 C parameters: 0: 7.78364 1: 3.22367 2: 9.39268 W parameters: 0: 0.763537 1: 0.761731 2: 0.40697 F(X) X(1) X(2) ... 0.000144957 0.218418 0.956318 0.829509 0.00276885 0.561695 0.415307 0.0661187 7.86758e-05 0.257578 0.109957 0.043829 0.0256053 0.633966 0.0617272 0.449539 0.00149094 0.401306 0.754673 0.797287 0.00101334 0.00183837 0.897504 0.350752 7.02836e-06 0.0945448 0.0136169 0.859097 0.00164098 0.840847 0.123104 0.00751236 0.000778691 0.260303 0.912484 0.113664 0.00894427 0.351629 0.822887 0.267132 Problem 6 C parameters: 0: 2.65802 1: 1.59152 2: 0.0504565 W parameters: 0: 0.501103 1: 0.0419093 2: 0.368851 F(X) X(1) X(2) ... 0 0.218418 0.956318 0.829509 0 0.561695 0.415307 0.0661187 0 0.257578 0.109957 0.043829 0 0.633966 0.0617272 0.449539 0 0.401306 0.754673 0.797287 0 0.00183837 0.897504 0.350752 1.37207 0.0945448 0.0136169 0.859097 0 0.840847 0.123104 0.00751236 0 0.260303 0.912484 0.113664 0 0.351629 0.822887 0.267132 TEST02 P00_F evaluates the function. Here, we use spatial dimension M = 4 The number of points is N = 10 Problem 1 C parameters: 0: 0.448694 1: 0.286721 2: 0.48216 3: 0.282426 W parameters: 0: 0.21293 1: 0.714471 2: 0.117707 3: 0.299329 F(X) X(1) X(2) ... -0.642591 0.218418 0.956318 0.829509 0.561695 -0.127288 0.415307 0.0661187 0.257578 0.109957 -0.124916 0.043829 0.633966 0.0617272 0.449539 -0.521378 0.401306 0.754673 0.797287 0.00183837 -0.314364 0.897504 0.350752 0.0945448 0.0136169 -0.43957 0.859097 0.840847 0.123104 0.00751236 -0.295157 0.260303 0.912484 0.113664 0.351629 -0.648197 0.822887 0.267132 0.692066 0.561662 -0.778456 0.861216 0.453794 0.911977 0.597917 -0.308818 0.188955 0.761492 0.396988 0.185314 Problem 2 C parameters: 0: 1.36234 1: 1.36178 2: 0.102154 3: 1.17373 W parameters: 0: 0.0882833 1: 0.777994 2: 0.745303 3: 0.308675 F(X) X(1) X(2) ... 0.0416258 0.218418 0.956318 0.829509 0.561695 0.0201354 0.415307 0.0661187 0.257578 0.109957 0.0459851 0.043829 0.633966 0.0617272 0.449539 0.0370209 0.401306 0.754673 0.797287 0.00183837 0.0148338 0.897504 0.350752 0.0945448 0.0136169 0.0206817 0.859097 0.840847 0.123104 0.00751236 0.0450787 0.260303 0.912484 0.113664 0.351629 0.015308 0.822887 0.267132 0.692066 0.561662 0.0176016 0.861216 0.453794 0.911977 0.597917 0.0474846 0.188955 0.761492 0.396988 0.185314 Problem 3 C parameters: 0: 0.587595 1: 0.498848 2: 0.497668 3: 0.265889 W parameters: 0: 0.938749 1: 0.562088 2: 0.01782 3: 0.501103 F(X) X(1) X(2) ... 0.0208996 0.218418 0.956318 0.829509 0.561695 0.16466 0.415307 0.0661187 0.257578 0.109957 0.135141 0.043829 0.633966 0.0617272 0.449539 0.0305148 0.401306 0.754673 0.797287 0.00183837 0.0604051 0.897504 0.350752 0.0945448 0.0136169 0.0322436 0.859097 0.840847 0.123104 0.00751236 0.0595186 0.260303 0.912484 0.113664 0.351629 0.0238798 0.822887 0.267132 0.692066 0.561662 0.0140943 0.861216 0.453794 0.911977 0.597917 0.063107 0.188955 0.761492 0.396988 0.185314 Problem 4 C parameters: 0: 0.191182 1: 1.68263 2: 1.23955 3: 3.91664 W parameters: 0: 0.0290366 1: 0.0174423 2: 0.152384 3: 0.114319 F(X) X(1) X(2) ... 0.00188895 0.218418 0.956318 0.829509 0.561695 0.970974 0.415307 0.0661187 0.257578 0.109957 0.060049 0.043829 0.633966 0.0617272 0.449539 0.0928314 0.401306 0.754673 0.797287 0.00183837 0.604831 0.897504 0.350752 0.0945448 0.0136169 0.119905 0.859097 0.840847 0.123104 0.00751236 0.0434456 0.260303 0.912484 0.113664 0.351629 0.0243101 0.822887 0.267132 0.692066 0.561662 0.00648377 0.861216 0.453794 0.911977 0.597917 0.175946 0.188955 0.761492 0.396988 0.185314 Problem 5 C parameters: 0: 8.08666 1: 2.72614 2: 4.72195 3: 4.86524 W parameters: 0: 0.612948 1: 0.809519 2: 0.58709 3: 0.215492 F(X) X(1) X(2) ... 0.00162918 0.218418 0.956318 0.829509 0.561695 0.00336523 0.415307 0.0661187 0.257578 0.109957 0.000166538 0.043829 0.633966 0.0617272 0.449539 0.020384 0.401306 0.754673 0.797287 0.00183837 0.00104921 0.897504 0.350752 0.0945448 0.0136169 0.00509896 0.859097 0.840847 0.123104 0.00751236 0.0024049 0.260303 0.912484 0.113664 0.351629 0.00471863 0.822887 0.267132 0.692066 0.561662 0.00170854 0.861216 0.453794 0.911977 0.597917 0.0100109 0.188955 0.761492 0.396988 0.185314 Problem 6 C parameters: 0: 1.18182 1: 1.11294 2: 0.689371 3: 1.31587 W parameters: 0: 0.945017 1: 0.909057 2: 0.519726 3: 0.0301946 F(X) X(1) X(2) ... 0 0.218418 0.956318 0.829509 0.561695 2.42701 0.415307 0.0661187 0.257578 0.109957 4.02069 0.043829 0.633966 0.0617272 0.449539 6.46385 0.401306 0.754673 0.797287 0.00183837 4.63737 0.897504 0.350752 0.0945448 0.0136169 7.73581 0.859097 0.840847 0.123104 0.00751236 0 0.260303 0.912484 0.113664 0.351629 12.0128 0.822887 0.267132 0.692066 0.561662 18.884 0.861216 0.453794 0.911977 0.597917 4.89556 0.188955 0.761492 0.396988 0.185314 TEST03 P00_D evaluates derivative components. Here, we use spatial dimension M = 2 The number of points is N = 2 Problem 1 C parameters: 0: 1.29399 1: 0.206009 W parameters: 0: 0.257578 1: 0.109957 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.218418 0.956318 -0.503162 -1.11826 -0.178031 0.829509 0.561695 -0.944708 -0.424314 -0.0675527 Problem 2 C parameters: 0: 0.129328 1: 1.87067 W parameters: 0: 0.0617272 1: 0.449539 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.218418 0.956318 0.0308134 -0.000161444 -0.0575597 0.829509 0.561695 0.0555153 -0.00141191 -0.0417401 Problem 3 C parameters: 0: 0.64224 1: 1.20776 W parameters: 0: 0.797287 1: 0.00183837 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.218418 0.956318 0.0826977 -0.0694188 -0.130545 0.829509 0.561695 0.0925024 -0.0806039 -0.151579 Problem 4 C parameters: 0: 5.05461 1: 1.97539 W parameters: 0: 0.0945448 1: 0.0136169 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.218418 0.956318 0.0210717 -0.133378 -0.155027 0.829509 0.561695 3.14241e-07 -1.18014e-05 -1.34413e-06 Problem 5 C parameters: 0: 10.3095 1: 10.0905 W parameters: 0: 0.123104 1: 0.00751236 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.218418 0.956318 2.60221e-05 -0.000268275 -0.000262576 0.829509 0.561695 2.56213e-06 -2.64143e-05 -2.58532e-05 Problem 6 C parameters: 0: 0.954396 1: 3.3456 W parameters: 0: 0.113664 1: 0.351629 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.218418 0.956318 0 0 0 0.829509 0.561695 0 0 0 TEST03 P00_D evaluates derivative components. Here, we use spatial dimension M = 3 The number of points is N = 2 Problem 1 C parameters: 0: 0.939234 1: 0.400947 2: 0.159818 W parameters: 0: 0.633966 1: 0.0617272 2: 0.449539 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.218418 0.956318 0.829509 -0.00791493 0.939205 0.400935 0.159813 0.561695 0.415307 0.0661187 -0.0244157 0.938954 0.400828 0.159771 Problem 2 C parameters: 0: 0.616362 1: 1.15909 2: 1.22454 W parameters: 0: 0.00183837 1: 0.897504 2: 0.350752 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.218418 0.956318 0.829509 0.55702 -0.0900574 -0.0876198 -0.5952 0.561695 0.415307 0.0661187 0.46467 -0.176629 0.45875 0.353684 Problem 3 C parameters: 0: 0.180828 1: 0.026044 2: 1.64313 W parameters: 0: 0.840847 1: 0.123104 2: 0.00751236 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.218418 0.956318 0.829509 0.0288032 -0.00858276 -0.00123614 -0.0779887 0.561695 0.415307 0.0661187 0.44988 -0.266501 -0.038383 -2.4216 Problem 4 C parameters: 0: 1.42246 1: 4.9864 2: 0.621134 W parameters: 0: 0.351629 1: 0.822887 2: 0.267132 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.218418 0.956318 0.829509 0.548484 0.295675 -3.63934 -0.238008 0.561695 0.415307 0.0661187 0.0144746 -0.0123049 0.293376 0.00224508 Problem 5 C parameters: 0: 6.67543 1: 5.4176 2: 8.30698 W parameters: 0: 0.453794 1: 0.911977 2: 0.597917 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.218418 0.956318 0.829509 0.0238664 0.159318 -0.129299 -0.198258 0.561695 0.415307 0.0661187 0.00039815 -0.00265782 0.00215701 0.00330742 Problem 6 C parameters: 0: 0.603001 1: 2.43011 2: 1.26689 W parameters: 0: 0.185314 1: 0.574366 2: 0.367027 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.218418 0.956318 0.829509 0 0 0 0 0.561695 0.415307 0.0661187 0 0 0 0 TEST03 P00_D evaluates derivative components. Here, we use spatial dimension M = 4 The number of points is N = 2 Problem 1 C parameters: 0: 0.0552903 1: 0.799748 2: 0.0778689 3: 0.567093 W parameters: 0: 0.401306 1: 0.754673 2: 0.797287 3: 0.00183837 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.218418 0.956318 0.829509 0.561695 -0.857758 0.0284222 0.411114 0.0400288 0.291517 0.415307 0.0661187 0.257578 0.109957 -0.895226 -0.024638 -0.356377 -0.0346993 -0.252703 Problem 2 C parameters: 0: 2.64669 1: 1.03435 2: 0.278807 3: 0.0401554 W parameters: 0: 0.859097 1: 0.840847 2: 0.123104 3: 0.00751236 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.218418 0.956318 0.829509 0.561695 0.000229951 0.000532602 -5.60166e-05 -2.43108e-05 -4.10763e-07 0.415307 0.0661187 0.257578 0.109957 0.000240049 0.000627197 0.000242327 -5.01147e-06 -7.93047e-08 Problem 3 C parameters: 0: 0.293979 1: 1.03053 2: 0.128369 3: 0.397119 W parameters: 0: 0.822887 1: 0.267132 2: 0.692066 3: 0.561662 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.218418 0.956318 0.829509 0.561695 0.0131154 -0.00810258 -0.0284033 -0.00353807 -0.0109453 0.415307 0.0661187 0.257578 0.109957 0.306327 -0.355393 -1.24582 -0.155186 -0.48008 Problem 4 C parameters: 0: 2.1432 1: 1.1293 2: 2.26953 3: 1.48796 W parameters: 0: 0.188955 1: 0.761492 2: 0.396988 3: 0.185314 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.218418 0.956318 0.829509 0.561695 0.264579 -0.0716141 -0.131477 -1.17886 -0.440959 0.415307 0.0661187 0.257578 0.109957 0.381099 -0.792465 0.675937 0.547312 0.127168 Problem 5 C parameters: 0: 6.10224 1: 3.8994 2: 6.55737 3: 3.84099 W parameters: 0: 0.21293 1: 0.714471 2: 0.117707 3: 0.299329 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.218418 0.956318 0.829509 0.561695 0.00129164 -0.00788187 -0.00503661 -0.00846974 -0.00496116 0.415307 0.0661187 0.257578 0.109957 0.00448192 -0.0273498 0.0174768 -0.0293896 0.017215 Problem 6 C parameters: 0: 1.46452 1: 1.46391 2: 0.109815 3: 1.26176 W parameters: 0: 0.0882833 1: 0.777994 2: 0.745303 3: 0.308675 X(1) X(2) ... F(X) dFdX(1) dFdX(2) ... 0.218418 0.956318 0.829509 0.561695 0 0 0 0 0 0.415307 0.0661187 0.257578 0.109957 0 0 0 0 0 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: 0: 0.322298 1: 0.306807 2: 0.192495 3: 0.6784 W parameters: 0: 0.480647 1: 0.239298 2: 0.889494 3: 0.729954 Exact Integral Q -0.785803 -0.78778 Problem 2 C parameters: 0: 0.665417 1: 1.55434 2: 1.75061 3: 0.0296339 W parameters: 0: 0.362054 1: 0.0453086 2: 0.50144 3: 0.703622 Exact Integral Q 0.00152845 0.00153381 Problem 3 C parameters: 0: 0.610387 1: 0.389117 2: 0.191277 3: 0.659219 W parameters: 0: 0.37722 1: 0.933423 2: 0.033472 3: 0.563336 Exact Integral Q 0.0550631 0.0557084 Problem 4 C parameters: 0: 2.69374 1: 0.994509 2: 2.21878 3: 1.12297 W parameters: 0: 0.792856 1: 0.535805 2: 0.280853 3: 0.299739 Exact Integral Q 0.263474 0.264295 Problem 5 C parameters: 0: 7.66269 1: 0.878415 2: 5.77331 3: 6.08559 W parameters: 0: 0.590769 1: 0.057818 2: 0.747028 3: 0.302164 Exact Integral Q 0.016149 0.0165698 Problem 6 C parameters: 0: 0.926777 1: 1.33768 2: 0.887673 3: 1.14787 W parameters: 0: 0.950829 1: 0.57687 2: 0.459283 3: 0.173201 Exact Integral Q 4.00448 4.01973 TEST_INTERP_ND_TEST Normal end of execution. 27 February 2022 11:15:14 AM