Skip to content

Commit

Permalink
Feature/arkode sts (#39)
Browse files Browse the repository at this point in the history 
Output files for LSRKStep, see LLNL/sundials#541
  • Loading branch information
@maggul
maggul authored Nov 15, 2024
1 parent 71ee007 commit 5dca98a
Show file tree
Hide file tree
Showing 14 changed files with 676 additions and 0 deletions.
40 changes: 40 additions & 0 deletions linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_analytic_lsrk.out
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,40 @@

Analytical ODE test problem:
lambda = -1e+06
reltol = 1.0e-08
abstol = 1.0e-08

t u
---------------------
1.000000 0.785398
2.000000 1.107149
3.000000 1.249046
4.000000 1.325818
5.000000 1.373401
6.000000 1.405648
7.000000 1.428899
8.000000 1.446441
9.000000 1.460139
10.000000 1.471128
---------------------

Final Statistics:
Current time = 10.00984671802051
Steps = 1481
Step attempts = 1482
Stability limited steps = 34
Accuracy limited steps = 1482
Error test fails = 1
NLS step fails = 0
Inequality constraint fails = 0
Initial step size = 1.930101110942615e-10
Last step size = 0.019104
Current step size = 0.03021704965932211
RHS fn evals = 161751
Number of dom_eig updates = 1
Max. num. of stages used = 196
Max. num. of stages allowed = 200
Max. spectral radius = 1010000
Min. spectral radius = 1010000

ACCURACY at the final time = 1.14353e-13
44 changes: 44 additions & 0 deletions linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_analytic_lsrk_varjac.out
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,44 @@

Analytical ODE test problem with a variable Jacobian:
The stiffness of the problem is directly proportional to
"lambda - alpha*cos((10 - t)/10*pi)"

lambda = -1e+06
alpha = 100
reltol = 1.0e-08
abstol = 1.0e-08

t u
---------------------
1.000000 0.785398
2.000000 1.107149
3.000000 1.249046
4.000000 1.325818
5.000000 1.373401
6.000000 1.405648
7.000000 1.428899
8.000000 1.446441
9.000000 1.460139
10.000000 1.471128
---------------------

Final Statistics:
Current time = 10.00171383959044
Steps = 1607
Step attempts = 1610
Stability limited steps = 0
Accuracy limited steps = 1610
Error test fails = 3
NLS step fails = 0
Inequality constraint fails = 0
Initial step size = 1.930101110942615e-10
Last step size = 0.02250442118061891
Current step size = 0.02250442118061891
RHS fn evals = 146958
Number of dom_eig updates = 67
Max. num. of stages used = 188
Max. num. of stages allowed = 200
Max. spectral radius = 1010099.676417927
Min. spectral radius = 1009899

ACCURACY at the final time = 1.27232e-13
36 changes: 36 additions & 0 deletions linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_analytic_ssprk.out
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,36 @@

Analytical ODE test problem:
lambda = -10
reltol = 1.0e-08
abstol = 1.0e-08

t u
---------------------
1.000000 0.785398
2.000000 1.107149
3.000000 1.249046
4.000000 1.325818
5.000000 1.373401
6.000000 1.405648
7.000000 1.428899
8.000000 1.446441
9.000000 1.460139
10.000000 1.471128
---------------------

Final Statistics:
Current time = 10.01679801682566
Steps = 1047
Step attempts = 1048
Stability limited steps = 0
Accuracy limited steps = 1048
Error test fails = 1
NLS step fails = 0
Inequality constraint fails = 0
Initial step size = 1.930101110942615e-10
Last step size = 0.02305133855414821
Current step size = 0.02305133855414821
RHS fn evals = 9433
Number of stages used = 9

ACCURACY at the final time = 1.69688e-09
71 changes: 71 additions & 0 deletions linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_heat2D_lsrk.out
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,71 @@

2D Heat PDE test problem:
---------------------------------
kx = 10
ky = 10
forcing = 1
tf = 1
xu = 1
yu = 1
nx = 64
ny = 64
dx = 0.015873
dy = 0.015873
---------------------------------
rtol = 1e-05
atol = 1e-10
fixed h = 0
controller = 0
method = 0
eigfrequency = 25
stage_max_limit = 1000
eigsafety = 1.01
---------------------------------
output = 1
max steps = 0
---------------------------------

t ||u||_rms max error
-----------------------------------------------------------------------
0.000000000000000e+00 3.691406250000004e-01 0.000000000000000e+00
5.000000000000000e-02 3.604068953520301e-01 8.105295841017712e-04
1.000000000000000e-01 3.341701389350745e-01 7.545034505144210e-04
1.500000000000000e-01 2.933045814980530e-01 6.655020848065352e-04
2.000000000000000e-01 2.418102694052765e-01 5.516214332939029e-04
2.500000000000000e-01 1.847278700505845e-01 4.242334056734909e-04
3.000000000000000e-01 1.276450308939409e-01 2.958580782215225e-04
3.500000000000000e-01 7.614944266700012e-02 1.791094087538192e-04
4.000000000000000e-01 3.528179910107998e-02 8.530904691973762e-05
4.500000000000000e-01 9.042542522688570e-03 2.370608964642751e-05
4.999999999999999e-01 1.483315165963637e-07 3.059603009336159e-07
5.499999999999999e-01 9.039739042933443e-03 1.741090628159039e-05
6.000000000000000e-01 3.527647409317574e-02 7.335124075792676e-05
6.500000000000000e-01 7.614211328474822e-02 1.626500505197415e-04
7.000000000000001e-01 1.276364431261110e-01 2.765661025802557e-04
7.500000000000001e-01 1.847188398303503e-01 4.039484934804660e-04
8.000000000000002e-01 2.418016806916353e-01 5.323281462729090e-04
8.500000000000002e-01 2.932973330370502e-01 6.492070720421905e-04
9.000000000000002e-01 3.341649712654264e-01 7.429941073312518e-04
9.500000000000003e-01 3.604042250428761e-01 8.045913926612025e-04
1.000000000000000e+00 3.694466133676554e-01 8.279696748304533e-04
-----------------------------------------------------------------------

Final integrator statistics:
Current time = 1
Steps = 617
Step attempts = 617
Stability limited steps = 0
Accuracy limited steps = 617
Error test fails = 0
NLS step fails = 0
Inequality constraint fails = 0
Initial step size = 2.877680248425025e-06
Last step size = 0.0003199001977636226
Current step size = 0.0003199001977636226
RHS fn evals = 17225
Number of dom_eig updates = 25
Max. num. of stages used = 46
Max. num. of stages allowed = 1000
Max. spectral radius = 320695.2
Min. spectral radius = 320695.2
Max error = 8.279696748304533e-04
77 changes: 77 additions & 0 deletions linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_heat2D_lsrk_p_--np_2_2.out
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,77 @@

2D Heat PDE test problem:
---------------------------------
nprocs = 4
npx = 2
npy = 2
---------------------------------
kx = 10
ky = 10
forcing = 1
tf = 1
xu = 1
yu = 1
nx = 64
ny = 64
nxl (proc 0) = 32
nyl (proc 0) = 32
dx = 0.015873
dy = 0.015873
---------------------------------
rtol = 1e-05
atol = 1e-10
fixed h = 0
controller = 0
method = 0
eigfrequency = 25
stage_max_limit = 1000
eigsafety = 1.01
---------------------------------
output = 1
max steps = 0
---------------------------------

t ||u||_rms max error
-----------------------------------------------------------------------
0.000000000000000e+00 3.691406249999997e-01 0.000000000000000e+00
5.000000000000000e-02 3.604068953517967e-01 8.105295833209514e-04
1.000000000000000e-01 3.341701389364218e-01 7.545034548108731e-04
1.500000000000000e-01 2.933045814986309e-01 6.655020864086980e-04
2.000000000000000e-01 2.418102694055522e-01 5.516214337898395e-04
2.500000000000000e-01 1.847278700508768e-01 4.242334062354858e-04
3.000000000000000e-01 1.276450308942233e-01 2.958580787605358e-04
3.500000000000000e-01 7.614944266702617e-02 1.791094087855438e-04
4.000000000000000e-01 3.528179910109314e-02 8.530904694552255e-05
4.500000000000000e-01 9.042542522694921e-03 2.370608965900772e-05
4.999999999999999e-01 1.483315166225389e-07 3.059603009683608e-07
5.499999999999999e-01 9.039739042784675e-03 1.741090593756697e-05
6.000000000000000e-01 3.527647409287470e-02 7.335124006428717e-05
6.500000000000000e-01 7.614211328432417e-02 1.626500495444660e-04
7.000000000000001e-01 1.276364431256062e-01 2.765661014183518e-04
7.500000000000001e-01 1.847188398298133e-01 4.039484922465642e-04
8.000000000000002e-01 2.418016806911184e-01 5.323281450848594e-04
8.500000000000002e-01 2.932973330369171e-01 6.492070716110909e-04
9.000000000000002e-01 3.341649712653251e-01 7.429941071179780e-04
9.500000000000003e-01 3.604042250428204e-01 8.045913925563974e-04
1.000000000000000e+00 3.694466133673749e-01 8.279696807845793e-04
-----------------------------------------------------------------------

Final integrator statistics:
Current time = 1
Steps = 617
Step attempts = 617
Stability limited steps = 0
Accuracy limited steps = 617
Error test fails = 0
NLS step fails = 0
Inequality constraint fails = 0
Initial step size = 2.877680248425019e-06
Last step size = 0.0003228665778641779
Current step size = 0.0003228665778641779
RHS fn evals = 17225
Number of dom_eig updates = 25
Max. num. of stages used = 46
Max. num. of stages allowed = 1000
Max. spectral radius = 320695.2
Min. spectral radius = 320695.2
Max error = 8.279696807845793e-04
48 changes: 48 additions & 0 deletions linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_sod_lsrk.out
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,48 @@

Problem parameters and options:
---------------------------------
gamma = 1.4
---------------------------------
tf = 0.1
xl = 0
xr = 1
nx = 512
dx = 0.00195312
---------------------------------
integrator = ARKODE_LSRK_SSP_10_4
rtol = 0.0001
atol = 1e-11
fixed h = 0
---------------------------------
output = 1
---------------------------------

t ||rho|| ||mx|| ||my|| ||mz|| ||et||
--------------------------------------------------------------------------
0.00e+00 7.12610e-01 0.00000e+00 0.00000e+00 0.00000e+00 1.77658e+00
1.00e-02 7.09460e-01 5.28096e-02 0.00000e+00 0.00000e+00 1.76732e+00
2.00e-02 7.06459e-01 7.61804e-02 0.00000e+00 0.00000e+00 1.75863e+00
3.00e-02 7.03444e-01 9.39395e-02 0.00000e+00 0.00000e+00 1.74990e+00
4.00e-02 7.00418e-01 1.08850e-01 0.00000e+00 0.00000e+00 1.74113e+00
5.00e-02 6.97379e-01 1.21955e-01 0.00000e+00 0.00000e+00 1.73232e+00
6.00e-02 6.94326e-01 1.33784e-01 0.00000e+00 0.00000e+00 1.72346e+00
7.00e-02 6.91261e-01 1.44650e-01 0.00000e+00 0.00000e+00 1.71456e+00
8.00e-02 6.88182e-01 1.54756e-01 0.00000e+00 0.00000e+00 1.70561e+00
9.00e-02 6.85089e-01 1.64242e-01 0.00000e+00 0.00000e+00 1.69661e+00
1.00e-01 6.81982e-01 1.73210e-01 0.00000e+00 0.00000e+00 1.68756e+00
--------------------------------------------------------------------------

Final integrator statistics:
Current time = 0.1
Steps = 382
Step attempts = 397
Stability limited steps = 0
Accuracy limited steps = 397
Error test fails = 15
NLS step fails = 0
Inequality constraint fails = 0
Initial step size = 1.122284259749542e-14
Last step size = 0.0002334416078592081
Current step size = 0.0002334416078592081
RHS fn evals = 3957
Number of stages used = 10
40 changes: 40 additions & 0 deletions linux-ubuntu20.04-x86_64/gcc-9.4.0/extended/ark_analytic_lsrk.out
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,40 @@

Analytical ODE test problem:
lambda = -1e+06
reltol = 1.0e-08
abstol = 1.0e-08

t u
---------------------
1.000000 0.785398
2.000000 1.107149
3.000000 1.249046
4.000000 1.325818
5.000000 1.373401
6.000000 1.405648
7.000000 1.428899
8.000000 1.446441
9.000000 1.460139
10.000000 1.471128
---------------------

Final Statistics:
Current time = 10.003415913614472770865049255917
Steps = 1510
Step attempts = 1511
Stability limited steps = 0
Accuracy limited steps = 1511
Error test fails = 1
NLS step fails = 0
Inequality constraint fails = 0
Initial step size = 2.8691124672595165773365555626549e-11
Last step size = 0.019103999999999999999481534521117
Current step size = 0.030441780662772963579203221520797
RHS fn evals = 162376
Number of dom_eig updates = 1
Max. num. of stages used = 196
Max. num. of stages allowed = 200
Max. spectral radius = 1010000
Min. spectral radius = 1010000

ACCURACY at the final time = 1.06895e-13
Loading

0 comments on commit 5dca98a

Please sign in to comment.