05 October 2012 09:44:26 AM
TEST_INTERP_ND_PRB
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_PRB
Normal end of execution.
05 October 2012 09:44:26 AM