From 5dca98a66d4f833c19a4fd63163f466ac5f50b56 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Mustafa=20A=C4=9Fg=C3=BCl?= <33010171+maggul@users.noreply.github.com> Date: Fri, 15 Nov 2024 15:39:39 -0600 Subject: [PATCH] Feature/arkode sts (#39) Output files for LSRKStep, see LLNL/sundials#541 --- .../gcc-9.4.0/double/ark_analytic_lsrk.out | 40 ++++++++++ .../double/ark_analytic_lsrk_varjac.out | 44 +++++++++++ .../gcc-9.4.0/double/ark_analytic_ssprk.out | 36 +++++++++ .../gcc-9.4.0/double/ark_heat2D_lsrk.out | 71 +++++++++++++++++ .../double/ark_heat2D_lsrk_p_--np_2_2.out | 77 +++++++++++++++++++ .../gcc-9.4.0/double/ark_sod_lsrk.out | 48 ++++++++++++ .../gcc-9.4.0/extended/ark_analytic_lsrk.out | 40 ++++++++++ .../extended/ark_analytic_lsrk_varjac.out | 44 +++++++++++ .../gcc-9.4.0/extended/ark_analytic_ssprk.out | 36 +++++++++ .../gcc-9.4.0/extended/ark_heat2D_lsrk.out | 72 +++++++++++++++++ .../gcc-9.4.0/extended/ark_sod_lsrk.out | 48 ++++++++++++ .../gcc-9.4.0/single/ark_analytic_lsrk.out | 40 ++++++++++ .../single/ark_analytic_lsrk_varjac.out | 44 +++++++++++ .../gcc-9.4.0/single/ark_analytic_ssprk.out | 36 +++++++++ 14 files changed, 676 insertions(+) create mode 100644 linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_analytic_lsrk.out create mode 100644 linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_analytic_lsrk_varjac.out create mode 100644 linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_analytic_ssprk.out create mode 100644 linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_heat2D_lsrk.out create mode 100644 linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_heat2D_lsrk_p_--np_2_2.out create mode 100644 linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_sod_lsrk.out create mode 100644 linux-ubuntu20.04-x86_64/gcc-9.4.0/extended/ark_analytic_lsrk.out create mode 100644 linux-ubuntu20.04-x86_64/gcc-9.4.0/extended/ark_analytic_lsrk_varjac.out create mode 100644 linux-ubuntu20.04-x86_64/gcc-9.4.0/extended/ark_analytic_ssprk.out create mode 100644 linux-ubuntu20.04-x86_64/gcc-9.4.0/extended/ark_heat2D_lsrk.out create mode 100644 linux-ubuntu20.04-x86_64/gcc-9.4.0/extended/ark_sod_lsrk.out create mode 100644 linux-ubuntu20.04-x86_64/gcc-9.4.0/single/ark_analytic_lsrk.out create mode 100644 linux-ubuntu20.04-x86_64/gcc-9.4.0/single/ark_analytic_lsrk_varjac.out create mode 100644 linux-ubuntu20.04-x86_64/gcc-9.4.0/single/ark_analytic_ssprk.out diff --git a/linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_analytic_lsrk.out b/linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_analytic_lsrk.out new file mode 100644 index 0000000..3de5e33 --- /dev/null +++ b/linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_analytic_lsrk.out @@ -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 diff --git a/linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_analytic_lsrk_varjac.out b/linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_analytic_lsrk_varjac.out new file mode 100644 index 0000000..46bc651 --- /dev/null +++ b/linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_analytic_lsrk_varjac.out @@ -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 diff --git a/linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_analytic_ssprk.out b/linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_analytic_ssprk.out new file mode 100644 index 0000000..457b582 --- /dev/null +++ b/linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_analytic_ssprk.out @@ -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 diff --git a/linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_heat2D_lsrk.out b/linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_heat2D_lsrk.out new file mode 100644 index 0000000..90a997b --- /dev/null +++ b/linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_heat2D_lsrk.out @@ -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 diff --git a/linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_heat2D_lsrk_p_--np_2_2.out b/linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_heat2D_lsrk_p_--np_2_2.out new file mode 100644 index 0000000..2e427c1 --- /dev/null +++ b/linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_heat2D_lsrk_p_--np_2_2.out @@ -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 diff --git a/linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_sod_lsrk.out b/linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_sod_lsrk.out new file mode 100644 index 0000000..24abdf2 --- /dev/null +++ b/linux-ubuntu20.04-x86_64/gcc-9.4.0/double/ark_sod_lsrk.out @@ -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 diff --git a/linux-ubuntu20.04-x86_64/gcc-9.4.0/extended/ark_analytic_lsrk.out b/linux-ubuntu20.04-x86_64/gcc-9.4.0/extended/ark_analytic_lsrk.out new file mode 100644 index 0000000..0936978 --- /dev/null +++ b/linux-ubuntu20.04-x86_64/gcc-9.4.0/extended/ark_analytic_lsrk.out @@ -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 diff --git a/linux-ubuntu20.04-x86_64/gcc-9.4.0/extended/ark_analytic_lsrk_varjac.out b/linux-ubuntu20.04-x86_64/gcc-9.4.0/extended/ark_analytic_lsrk_varjac.out new file mode 100644 index 0000000..35f4788 --- /dev/null +++ b/linux-ubuntu20.04-x86_64/gcc-9.4.0/extended/ark_analytic_lsrk_varjac.out @@ -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.019761744896271672689014309299 +Steps = 1684 +Step attempts = 1685 +Stability limited steps = 0 +Accuracy limited steps = 1685 +Error test fails = 1 +NLS step fails = 0 +Inequality constraint fails = 0 +Initial step size = 2.8691124672595165773365555626549e-11 +Last step size = 0.02463588618472410449365358510454 +Current step size = 0.02463588618472410449365358510454 +RHS fn evals = 149012 +Number of dom_eig updates = 68 +Max. num. of stages used = 196 +Max. num. of stages allowed = 200 +Max. spectral radius = 1010100.5502498682042755717702676 +Min. spectral radius = 1009899 + +ACCURACY at the final time = 1.70378e-13 diff --git a/linux-ubuntu20.04-x86_64/gcc-9.4.0/extended/ark_analytic_ssprk.out b/linux-ubuntu20.04-x86_64/gcc-9.4.0/extended/ark_analytic_ssprk.out new file mode 100644 index 0000000..03578c4 --- /dev/null +++ b/linux-ubuntu20.04-x86_64/gcc-9.4.0/extended/ark_analytic_ssprk.out @@ -0,0 +1,36 @@ + +Analytical ODE test problem: + lambda = -10 + reltol = 1.0e-09 + abstol = 1.0e-09 + + 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.007995934261343870703975245107 +Steps = 2250 +Step attempts = 2251 +Stability limited steps = 0 +Accuracy limited steps = 2251 +Error test fails = 1 +NLS step fails = 0 +Inequality constraint fails = 0 +Initial step size = 5.1020836256135713498183845473047e-12 +Last step size = 0.010680653872584741312028275717766 +Current step size = 0.010680653872584741312028275717766 +RHS fn evals = 20260 +Number of stages used = 9 + +ACCURACY at the final time = 1.67376e-10 diff --git a/linux-ubuntu20.04-x86_64/gcc-9.4.0/extended/ark_heat2D_lsrk.out b/linux-ubuntu20.04-x86_64/gcc-9.4.0/extended/ark_heat2D_lsrk.out new file mode 100644 index 0000000..2c1804f --- /dev/null +++ b/linux-ubuntu20.04-x86_64/gcc-9.4.0/extended/ark_heat2D_lsrk.out @@ -0,0 +1,72 @@ + + +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.000000000000000000e+00 3.691406250000000003e-01 0.000000000000000000e+00 +5.000000000000000000e-02 3.604068953509203244e-01 8.105295805123803447e-04 +1.000000000000000000e-01 3.341701389361607421e-01 7.545034540760439733e-04 +1.500000000000000000e-01 2.933045814989621925e-01 6.655020872766008217e-04 +2.000000000000000000e-01 2.418102694051705693e-01 5.516214330740384594e-04 +2.500000000000000000e-01 1.847278700504722997e-01 4.242334054577285240e-04 +3.000000000000000000e-01 1.276450308938331190e-01 2.958580780149134301e-04 +3.500000000000000000e-01 7.614944266701188971e-02 1.791094087891994537e-04 +4.000000000000000000e-01 3.528179910109258407e-02 8.530904694562694375e-05 +4.500000000000000000e-01 9.042542522692682192e-03 2.370608965467836992e-05 +5.000000000000000001e-01 1.483315166001197479e-07 3.059603009402151958e-07 +5.500000000000000001e-01 9.039739042935413814e-03 1.741090628602193960e-05 +6.000000000000000001e-01 3.527647409317965266e-02 7.335124076643660299e-05 +6.500000000000000001e-01 7.614211328475368128e-02 1.626500505321930095e-04 +7.000000000000000001e-01 1.276364431261181661e-01 2.765661025954269808e-04 +7.500000000000000001e-01 1.847188398303575561e-01 4.039484934969580631e-04 +8.000000000000000001e-01 2.418016806916415405e-01 5.323281462855703437e-04 +8.500000000000000001e-01 2.932973330371505503e-01 6.492070722437629535e-04 +9.000000000000000001e-01 3.341649712655011408e-01 7.429941074708659417e-04 +9.500000000000000002e-01 3.604042250429179231e-01 8.045913927347891968e-04 +1.000000000000000000e+00 3.694466133676789137e-01 8.279696744365391283e-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.877680248394454148564307479148e-06 +Last step size = 0.00031970316937780904491904844283347 +Current step size = 0.00031970316937780904491904844283347 +RHS fn evals = 17225 +Number of dom_eig updates = 50 +Max. num. of stages used = 46 +Max. num. of stages allowed = 1000 +Max. spectral radius = 320695.19999999999998863131622784 +Min. spectral radius = 320695.19999999999998863131622784 + Max error = 8.279696744365391283e-04 diff --git a/linux-ubuntu20.04-x86_64/gcc-9.4.0/extended/ark_sod_lsrk.out b/linux-ubuntu20.04-x86_64/gcc-9.4.0/extended/ark_sod_lsrk.out new file mode 100644 index 0000000..156f65e --- /dev/null +++ b/linux-ubuntu20.04-x86_64/gcc-9.4.0/extended/ark_sod_lsrk.out @@ -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.28095e-02 0.00000e+00 0.00000e+00 1.76732e+00 + 2.00e-02 7.06459e-01 7.61806e-02 0.00000e+00 0.00000e+00 1.75863e+00 + 3.00e-02 7.03444e-01 9.39396e-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.10000000000000000000135525271561 +Steps = 394 +Step attempts = 401 +Stability limited steps = 0 +Accuracy limited steps = 401 +Error test fails = 7 +NLS step fails = 0 +Inequality constraint fails = 0 +Initial step size = 8.6275149528694139164424205093341e-16 +Last step size = 0.0001522612671736054756833432302361 +Current step size = 0.0001522612671736054756833432302361 +RHS fn evals = 4005 +Number of stages used = 10 \ No newline at end of file diff --git a/linux-ubuntu20.04-x86_64/gcc-9.4.0/single/ark_analytic_lsrk.out b/linux-ubuntu20.04-x86_64/gcc-9.4.0/single/ark_analytic_lsrk.out new file mode 100644 index 0000000..c3e8503 --- /dev/null +++ b/linux-ubuntu20.04-x86_64/gcc-9.4.0/single/ark_analytic_lsrk.out @@ -0,0 +1,40 @@ + +Analytical ODE test problem: + lambda = -1000 + reltol = 1.0e-04 + abstol = 1.0e-08 + + t u + --------------------- + 1.000000 0.785402 + 2.000000 1.107148 + 3.000000 1.249045 + 4.000000 1.325815 + 5.000000 1.373399 + 6.000000 1.405649 + 7.000000 1.428898 + 8.000000 1.446444 + 9.000000 1.460139 + 10.000000 1.471127 + --------------------- + +Final Statistics: +Current time = 10.01134014129639 +Steps = 149 +Step attempts = 164 +Stability limited steps = 0 +Accuracy limited steps = 164 +Error test fails = 15 +NLS step fails = 0 +Inequality constraint fails = 0 +Initial step size = 3.452669830039667e-07 +Last step size = 0.1750497967004776 +Current step size = 0.4429996311664581 +RHS fn evals = 1882 +Number of dom_eig updates = 1 +Max. num. of stages used = 24 +Max. num. of stages allowed = 200 +Max. spectral radius = 1010 +Min. spectral radius = 1010 + +ACCURACY at the final time = 5.96046e-07 diff --git a/linux-ubuntu20.04-x86_64/gcc-9.4.0/single/ark_analytic_lsrk_varjac.out b/linux-ubuntu20.04-x86_64/gcc-9.4.0/single/ark_analytic_lsrk_varjac.out new file mode 100644 index 0000000..0681384 --- /dev/null +++ b/linux-ubuntu20.04-x86_64/gcc-9.4.0/single/ark_analytic_lsrk_varjac.out @@ -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 = -1000 + alpha = 10 + reltol = 1.0e-04 + abstol = 1.0e-08 + + t u + --------------------- + 1.000000 0.785398 + 2.000000 1.107149 + 3.000000 1.249041 + 4.000000 1.325817 + 5.000000 1.373399 + 6.000000 1.405647 + 7.000000 1.428900 + 8.000000 1.446441 + 9.000000 1.460141 + 10.000000 1.471126 + --------------------- + +Final Statistics: +Current time = 10.00951671600342 +Steps = 287 +Step attempts = 304 +Stability limited steps = 0 +Accuracy limited steps = 304 +Error test fails = 17 +NLS step fails = 0 +Inequality constraint fails = 0 +Initial step size = 3.452669830039667e-07 +Last step size = 0.06377532333135605 +Current step size = 0.06377532333135605 +RHS fn evals = 2131 +Number of dom_eig updates = 21 +Max. num. of stages used = 21 +Max. num. of stages allowed = 200 +Max. spectral radius = 1020.074462890625 +Min. spectral radius = 999.8999633789062 + +ACCURACY at the final time = 1.43051e-06 diff --git a/linux-ubuntu20.04-x86_64/gcc-9.4.0/single/ark_analytic_ssprk.out b/linux-ubuntu20.04-x86_64/gcc-9.4.0/single/ark_analytic_ssprk.out new file mode 100644 index 0000000..ef2a8d8 --- /dev/null +++ b/linux-ubuntu20.04-x86_64/gcc-9.4.0/single/ark_analytic_ssprk.out @@ -0,0 +1,36 @@ + +Analytical ODE test problem: + lambda = -10 + reltol = 1.0e-04 + abstol = 1.0e-08 + + t u + --------------------- + 1.000000 0.785403 + 2.000000 1.107155 + 3.000000 1.249053 + 4.000000 1.325830 + 5.000000 1.373416 + 6.000000 1.405659 + 7.000000 1.428910 + 8.000000 1.446454 + 9.000000 1.460148 + 10.000000 1.471140 + --------------------- + +Final Statistics: +Current time = 10.56618309020996 +Steps = 73 +Step attempts = 75 +Stability limited steps = 0 +Accuracy limited steps = 75 +Error test fails = 2 +NLS step fails = 0 +Inequality constraint fails = 0 +Initial step size = 3.452669830039667e-07 +Last step size = 0.6506694555282593 +Current step size = 0.6506694555282593 +RHS fn evals = 673 +Number of stages used = 9 + +ACCURACY at the final time = 1.19209e-05