18 September 2013 11:15:55 AM
BURGERS_SOLUTION_PRB
C++ version
Test the BURGERS_SOLUTION library.
TEST01
Compute an analytic solution to the Burgers equation.
Viscosity NU = 0.0031831
NX = 11
NT = 11
X grid points:
0: -1
1: -0.8
2: -0.6
3: -0.4
4: -0.2
5: 0
6: 0.2
7: 0.4
8: 0.6
9: 0.8
10: 1
T grid points:
0: 0
1: 0.095493
2: 0.190986
3: 0.286479
4: 0.381972
5: 0.477465
6: 0.572958
7: 0.668451
8: 0.763944
9: 0.859437
10: 0.95493
U(X,T) at grid points:
Col: 0 1 2 3 4
Row
0: 1.22465e-16 1.51662e-16 1.31656e-16 3.73866e-17 2.2444e-16
1: 0.587785 0.467745 0.385266 0.326331 0.282259
2: 0.951057 0.843167 0.74283 0.661054 0.590027
3: 0.951057 0.997698 0.990755 0.971206 0.946876
4: 0.587785 0.748375 0.869466 0.923817 0.954709
5: -0 -2.88173e-17 -2.33948e-17 1.32611e-17 -8.4109e-17
6: -0.587785 -0.748375 -0.869466 -0.923817 -0.954709
7: -0.951057 -0.997698 -0.990755 -0.971206 -0.946876
8: -0.951057 -0.843167 -0.74283 -0.661054 -0.590027
9: -0.587785 -0.467745 -0.385266 -0.326331 -0.282259
10: -1.22465e-16 -9.66198e-17 -5.42872e-17 3.70492e-17 -1.20307e-16
Col: 5 6 7 8 9
Row
0: 1.9052e-16 2.57603e-16 1.82473e-16 1.73482e-16 1.34109e-16
1: 0.248397 0.222439 0.20304 0.188983 0.178866
2: 0.525742 0.466479 0.411552 0.360885 0.314918
3: 0.919698 0.890611 0.860186 0.828804 0.796735
4: 0.974759 0.987809 0.995661 0.999404 0.999787
5: -5.8488e-17 -2.31326e-16 -3.545e-16 -3.91835e-16 -3.22888e-16
6: -0.974759 -0.987809 -0.995661 -0.999404 -0.999787
7: -0.919698 -0.890611 -0.860186 -0.828804 -0.796735
8: -0.525742 -0.466479 -0.411552 -0.360885 -0.314918
9: -0.248397 -0.222439 -0.20304 -0.188983 -0.178866
10: -6.82723e-17 -1.62661e-16 1.53606e-16 2.01942e-16 -1.76069e-17
Col: 10
Row
0: 1.73718e-16
1: 0.171141
2: 0.274681
3: 0.764181
4: 0.997357
5: -3.97525e-16
6: -0.997357
7: -0.764181
8: -0.274681
9: -0.171141
10: -7.18165e-17
Data written to file "burgers_test01.txt".
TEST02
Compute an analytic solution to the Burgers equation.
Viscosity NU = 0.0031831
NX = 41
NT = 41
X grid points:
0: -1
1: -0.95
2: -0.9
3: -0.85
4: -0.8
5: -0.75
6: -0.7
7: -0.65
8: -0.6
9: -0.55
10: -0.5
11: -0.45
12: -0.4
13: -0.35
14: -0.3
15: -0.25
16: -0.2
17: -0.15
18: -0.1
19: -0.05
20: 0
21: 0.05
22: 0.1
23: 0.15
24: 0.2
25: 0.25
26: 0.3
27: 0.35
28: 0.4
29: 0.45
30: 0.5
31: 0.55
32: 0.6
33: 0.65
34: 0.7
35: 0.75
36: 0.8
37: 0.85
38: 0.9
39: 0.95
40: 1
T grid points:
0: 0
1: 0.0238732
2: 0.0477465
3: 0.0716197
4: 0.095493
5: 0.119366
6: 0.143239
7: 0.167113
8: 0.190986
9: 0.214859
10: 0.238732
11: 0.262606
12: 0.286479
13: 0.310352
14: 0.334225
15: 0.358099
16: 0.381972
17: 0.405845
18: 0.429718
19: 0.453592
20: 0.477465
21: 0.501338
22: 0.525211
23: 0.549085
24: 0.572958
25: 0.596831
26: 0.620704
27: 0.644578
28: 0.668451
29: 0.692324
30: 0.716197
31: 0.74007
32: 0.763944
33: 0.787817
34: 0.81169
35: 0.835563
36: 0.859437
37: 0.88331
38: 0.907183
39: 0.931056
40: 0.95493
Data written to file "burgers_test02.txt".
BURGERS_SOLUTION_PRB
Normal end of execution.
18 September 2013 11:15:56 AM