Skip to content

Commit

Permalink
add logging to SSP104
Browse files Browse the repository at this point in the history
  • Loading branch information
@gardner48
gardner48 committed Nov 22, 2024
1 parent cdba63a commit f95edba
Show file tree
Hide file tree
Showing 7 changed files with 510 additions and 78 deletions.
146 changes: 69 additions & 77 deletions src/arkode/arkode_lsrkstep.c
View file
Original file line number Diff line number Diff line change
Expand Up @@ -1808,12 +1808,9 @@ int lsrkStep_TakeStepSSP104(ARKodeMem ark_mem, sunrealtype* dsmPtr, int* nflagPt
sunrealtype* cvals = step_mem->cvals;
N_Vector* Xvecs = step_mem->Xvecs;

#if SUNDIALS_LOGGING_LEVEL >= SUNDIALS_LOGGING_DEBUG
SUNLogger_QueueMsg(ARK_LOGGER, SUN_LOGLEVEL_DEBUG,
"ARKODE::lsrkStep_TakeStepSSP104", "start-stage",
"step = %li, stage = 0, h = %" RSYM ", tcur = %" RSYM,
ark_mem->nst, ark_mem->h, ark_mem->tcur);
#endif
SUNLogInfo(ARK_LOGGER, "begin-stage", "stage = %i, tcur = %" RSYM, 0,
ark_mem->tcur);
SUNLogExtraDebugVec(ARK_LOGGER, "stage", ark_mem->yn, "z_0(:) =");

/* The method is not FSAL. Therefore, fn ​is computed at the beginning
of the step unless ARKODE updated fn. */
Expand All @@ -1822,19 +1819,20 @@ int lsrkStep_TakeStepSSP104(ARKodeMem ark_mem, sunrealtype* dsmPtr, int* nflagPt
retval = step_mem->fe(ark_mem->tn, ark_mem->yn, ark_mem->fn,
ark_mem->user_data);
step_mem->nfe++;
if (retval != ARK_SUCCESS)
{
SUNLogExtraDebugVec(ARK_LOGGER, "stage RHS", ark_mem->fn, "F_0(:) =");
SUNLogInfo(ARK_LOGGER, "end-stage",
"status = failed rhs eval, retval = %i", retval);
return (ARK_RHSFUNC_FAIL);
}
ark_mem->fn_is_current = SUNTRUE;
if (retval != ARK_SUCCESS) { return (ARK_RHSFUNC_FAIL); }
}

#ifdef SUNDIALS_LOGGING_EXTRA_DEBUG
SUNLogger_QueueMsg(ARK_LOGGER, SUN_LOGLEVEL_DEBUG,
"ARKODE::lsrkStep_TakeStepSSP104", "stage", "z_0(:) =", "");
N_VPrintFile(ark_mem->ycur, ARK_LOGGER->debug_fp);
SUNLogger_QueueMsg(ARK_LOGGER, SUN_LOGLEVEL_DEBUG,
"ARKODE::lsrkStep_TakeStepSSP104", "stage RHS",
"F(:) =", "");
N_VPrintFile(ark_mem->fn, ARK_LOGGER->debug_fp);
#endif
SUNLogExtraDebugVec(ARK_LOGGER, "stage RHS", ark_mem->fn, "F_0(:) =");
SUNLogInfo(ARK_LOGGER, "end-stage", "status = success");
SUNLogInfo(ARK_LOGGER, "begin-stage", "stage = %i, tcur = %" RSYM, 1,
ark_mem->tn + onesixth * ark_mem->h);

N_VScale(ONE, ark_mem->yn, ark_mem->tempv2);

Expand All @@ -1849,37 +1847,35 @@ int lsrkStep_TakeStepSSP104(ARKodeMem ark_mem, sunrealtype* dsmPtr, int* nflagPt
/* Evaluate stages j = 2,...,step_mem->req_stages */
for (int j = 2; j <= 5; j++)
{
#if SUNDIALS_LOGGING_LEVEL >= SUNDIALS_LOGGING_DEBUG
SUNLogger_QueueMsg(ARK_LOGGER, SUN_LOGLEVEL_DEBUG,
"ARKODE::lsrkStep_TakeStepSSP104", "start-stage",
"step = %li, stage = %i, h = %" RSYM ", tcur = %" RSYM,
ark_mem->nst, j, ark_mem->h,
ark_mem->tcur +
((sunrealtype)j - ONE) * onesixth * ark_mem->h);
#endif

/* apply user-supplied stage postprocessing function (if supplied) */
if (ark_mem->ProcessStage != NULL)
{
retval = ark_mem->ProcessStage(ark_mem->tcur + ((sunrealtype)j - ONE) *
onesixth * ark_mem->h,
ark_mem->ycur, ark_mem->user_data);
if (retval != 0) { return ARK_POSTPROCESS_STAGE_FAIL; }
if (retval != 0)
{
SUNLogInfo(ARK_LOGGER, "end-stage",
"status = failed postprocess stage, retval = %i", retval);
return ARK_POSTPROCESS_STAGE_FAIL;
}
}

retval = step_mem->fe(ark_mem->tcur +
((sunrealtype)j - ONE) * onesixth * ark_mem->h,
ark_mem->ycur, ark_mem->tempv3, ark_mem->user_data);
step_mem->nfe++;

SUNLogExtraDebugVec(ARK_LOGGER, "stage RHS", ark_mem->fn, "F_%i(:) =", j - 1);
SUNLogInfoIf(retval != 0, ARK_LOGGER, "end-stage",
"status = failed rhs eval, retval = %i", retval);

if (retval < 0) { return ARK_RHSFUNC_FAIL; }
if (retval > 0) { return RHSFUNC_RECVR; }

#ifdef SUNDIALS_LOGGING_EXTRA_DEBUG
SUNLogger_QueueMsg(ARK_LOGGER, SUN_LOGLEVEL_DEBUG,
"ARKODE::lsrkStep_TakeStepSSP104", "stage RHS",
"F_%i(:) =", j);
N_VPrintFile(ark_mem->tempv3, ARK_LOGGER->debug_fp);
#endif
SUNLogInfo(ARK_LOGGER, "end-stage", "status = success");
SUNLogInfo(ARK_LOGGER, "begin-stage", "stage = %i, tcur = %" RSYM, j,
ark_mem->tn + j * onesixth * ark_mem->h);

N_VLinearSum(ONE, ark_mem->ycur, onesixth * ark_mem->h, ark_mem->tempv3,
ark_mem->ycur);
Expand All @@ -1889,6 +1885,11 @@ int lsrkStep_TakeStepSSP104(ARKodeMem ark_mem, sunrealtype* dsmPtr, int* nflagPt
ark_mem->tempv3, ark_mem->tempv1);
}
}

SUNLogInfo(ARK_LOGGER, "end-stage", "status = success");
SUNLogInfo(ARK_LOGGER, "begin-stage", "stage = %i, tcur = %" RSYM, 6,
ark_mem->tn + TWO * onesixth * ark_mem->h);

N_VLinearSum(SUN_RCONST(1.0) / SUN_RCONST(25.0), ark_mem->tempv2,
SUN_RCONST(9.0) / SUN_RCONST(25.0), ark_mem->ycur,
ark_mem->tempv2);
Expand All @@ -1900,33 +1901,31 @@ int lsrkStep_TakeStepSSP104(ARKodeMem ark_mem, sunrealtype* dsmPtr, int* nflagPt
{
retval = ark_mem->ProcessStage(ark_mem->tcur + TWO * onesixth * ark_mem->h,
ark_mem->ycur, ark_mem->user_data);
if (retval != 0) { return ARK_POSTPROCESS_STAGE_FAIL; }
if (retval != 0)
{
SUNLogInfo(ARK_LOGGER, "end-stage",
"status = failed postprocess stage, retval = %i", retval);
return ARK_POSTPROCESS_STAGE_FAIL;
}
}

for (int j = 6; j <= 9; j++)
{
#if SUNDIALS_LOGGING_LEVEL >= SUNDIALS_LOGGING_DEBUG
SUNLogger_QueueMsg(ARK_LOGGER, SUN_LOGLEVEL_DEBUG,
"ARKODE::lsrkStep_TakeStepSSP104", "start-stage",
"step = %li, stage = %i, h = %" RSYM ", tcur = %" RSYM,
ark_mem->nst, j, ark_mem->h,
ark_mem->tcur +
((sunrealtype)j - FOUR) * onesixth * ark_mem->h);
#endif

retval = step_mem->fe(ark_mem->tcur +
((sunrealtype)j - FOUR) * onesixth * ark_mem->h,
ark_mem->ycur, ark_mem->tempv3, ark_mem->user_data);
step_mem->nfe++;

SUNLogExtraDebugVec(ARK_LOGGER, "stage RHS", ark_mem->fn, "F_%i(:) =", j);
SUNLogInfoIf(retval != 0, ARK_LOGGER, "end-stage",
"status = failed rhs eval, retval = %i", retval);

if (retval < 0) { return ARK_RHSFUNC_FAIL; }
if (retval > 0) { return RHSFUNC_RECVR; }

#ifdef SUNDIALS_LOGGING_EXTRA_DEBUG
SUNLogger_QueueMsg(ARK_LOGGER, SUN_LOGLEVEL_DEBUG,
"ARKODE::lsrkStep_TakeStepSSP104", "stage RHS",
"F_%i(:) =", j);
N_VPrintFile(ark_mem->tempv3, ARK_LOGGER->debug_fp);
#endif
SUNLogInfo(ARK_LOGGER, "end-stage", "status = success");
SUNLogInfo(ARK_LOGGER, "begin-stage", "stage = %i, tcur = %" RSYM, j + 1,
ark_mem->tn + (j - 3) * onesixth * ark_mem->h);

N_VLinearSum(ONE, ark_mem->ycur, onesixth * ark_mem->h, ark_mem->tempv3,
ark_mem->ycur);
Expand All @@ -1948,30 +1947,26 @@ int lsrkStep_TakeStepSSP104(ARKodeMem ark_mem, sunrealtype* dsmPtr, int* nflagPt
retval = ark_mem->ProcessStage(ark_mem->tcur + ((sunrealtype)j - THREE) *
onesixth * ark_mem->h,
ark_mem->ycur, ark_mem->user_data);
if (retval != 0) { return ARK_POSTPROCESS_STAGE_FAIL; }
if (retval != 0)
{
SUNLogInfo(ARK_LOGGER, "end-stage",
"status = failed postprocess stage, retval = %i", retval);
return ARK_POSTPROCESS_STAGE_FAIL;
}
}
}

#if SUNDIALS_LOGGING_LEVEL >= SUNDIALS_LOGGING_DEBUG
SUNLogger_QueueMsg(ARK_LOGGER, SUN_LOGLEVEL_DEBUG,
"ARKODE::lsrkStep_TakeStepSSP104", "start-stage",
"step = %li, stage = %i, h = %" RSYM ", tcur = %" RSYM,
ark_mem->nst, 10, ark_mem->h, ark_mem->tcur + ark_mem->h);
#endif

retval = step_mem->fe(ark_mem->tcur + ark_mem->h, ark_mem->ycur,
ark_mem->tempv3, ark_mem->user_data);
step_mem->nfe++;

SUNLogExtraDebugVec(ARK_LOGGER, "stage RHS", ark_mem->fn, "F_10(:) =", 9);
SUNLogInfoIf(retval != 0, ARK_LOGGER, "end-stage",
"status = failed rhs eval, retval = %i", retval);

if (retval < 0) { return ARK_RHSFUNC_FAIL; }
if (retval > 0) { return RHSFUNC_RECVR; }

#ifdef SUNDIALS_LOGGING_EXTRA_DEBUG
SUNLogger_QueueMsg(ARK_LOGGER, SUN_LOGLEVEL_DEBUG,
"ARKODE::lsrkStep_TakeStepSSP104", "stage RHS",
"F_%i(:) =", 10);
N_VPrintFile(ark_mem->tempv3, ARK_LOGGER->debug_fp);
#endif

cvals[0] = SUN_RCONST(0.6);
Xvecs[0] = ark_mem->ycur;
cvals[1] = ONE;
Expand All @@ -1980,6 +1975,17 @@ int lsrkStep_TakeStepSSP104(ARKodeMem ark_mem, sunrealtype* dsmPtr, int* nflagPt
Xvecs[2] = ark_mem->tempv3;

retval = N_VLinearCombination(3, cvals, Xvecs, ark_mem->ycur);
if (retval != 0)
{
SUNLogInfo(ARK_LOGGER, "end-stage",
"status = failed vector op, retval = %i", retval);
return ARK_VECTOROP_ERR;
}

SUNLogInfo(ARK_LOGGER, "end-stage", "status = success");
SUNLogExtraDebugVec(ARK_LOGGER, "updated solution", ark_mem->ycur, "ycur(:) =");
SUNLogExtraDebugVec(ARK_LOGGER, "embedded solution", ark_mem->tempv1,
"y_embedded(:) =");

/* Compute yerr (if step adaptivity enabled) */
if (!ark_mem->fixedstep)
Expand All @@ -1989,20 +1995,6 @@ int lsrkStep_TakeStepSSP104(ARKodeMem ark_mem, sunrealtype* dsmPtr, int* nflagPt
*dsmPtr = N_VWrmsNorm(ark_mem->tempv1, ark_mem->ewt);
}

#ifdef SUNDIALS_LOGGING_EXTRA_DEBUG
SUNLogger_QueueMsg(ARK_LOGGER, SUN_LOGLEVEL_DEBUG,
"ARKODE::lsrkStep_TakeStepSSP104", "updated solution",
"ycur(:) =", "");
N_VPrintFile(ark_mem->ycur, ARK_LOGGER->debug_fp);
#endif

#if SUNDIALS_LOGGING_LEVEL >= SUNDIALS_LOGGING_DEBUG
SUNLogger_QueueMsg(ARK_LOGGER, SUN_LOGLEVEL_DEBUG,
"ARKODE::lsrkStep_TakeStepSSP104", "error-test",
"step = %li, h = %" RSYM ", dsm = %" RSYM, ark_mem->nst,
ark_mem->h, *dsmPtr);
#endif

return ARK_SUCCESS;
}

Expand Down
1 change: 1 addition & 0 deletions test/unit_tests/logging/CMakeLists.txt
View file
Original file line number Diff line number Diff line change
Expand Up @@ -38,6 +38,7 @@ if(BUILD_ARKODE)
list(APPEND unit_tests "test_logging_arkode_lsrkstep.cpp\;2") # SSPs2
list(APPEND unit_tests "test_logging_arkode_lsrkstep.cpp\;3") # SSPs3
list(APPEND unit_tests "test_logging_arkode_lsrkstep.cpp\;4") # SSP43
list(APPEND unit_tests "test_logging_arkode_lsrkstep.cpp\;5") # SSP104
# MRIStep
list(APPEND unit_tests "test_logging_arkode_mristep.cpp\;0") # Ex-MRI
list(APPEND unit_tests "test_logging_arkode_mristep.cpp\;1 1 1") # Im-MRI
Expand Down
2 changes: 1 addition & 1 deletion test/unit_tests/logging/test_logging_arkode_lsrkstep.cpp
View file
Original file line number Diff line number Diff line change
Expand Up @@ -88,7 +88,7 @@ int main(int argc, char* argv[])
}
else if (method == 5)
{
flag = LSRKStepSetSSPMethodByName(arkode_mem, "ARKODE_LSRK_SSP_S_10_4");
flag = LSRKStepSetSSPMethodByName(arkode_mem, "ARKODE_LSRK_SSP_10_4");
}
else
{
Expand Down
22 changes: 22 additions & 0 deletions test/unit_tests/logging/test_logging_arkode_lsrkstep_5.out
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,22 @@
Start LSRKStep Logging test
t y y err
---------------------------------------------------------------------
0.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
6.103515625000001e-12 6.103515625000001e-12 8.077935669463161e-28
5.086597062056021e-10 5.086597062056018e-10 3.101927297073854e-25
8.765030586579745e-09 8.765030586579742e-09 3.308722450212111e-24
---------------------------------------------------------------------
Current time = 8.765030586579745e-09
Steps = 3
Step attempts = 3
Stability limited steps = 0
Accuracy limited steps = 3
Error test fails = 0
NLS step fails = 0
Inequality constraint fails = 0
Initial step size = 6.103515625000001e-12
Last step size = 8.256370880374143e-09
Current step size = 1.651274176074829e-07
RHS fn evals = 32
Number of stages used = 10
End LSRKStep Logging test
94 changes: 94 additions & 0 deletions test/unit_tests/logging/test_logging_arkode_lsrkstep_lvl3_5.out
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,94 @@
Start LSRKStep Logging test
t y y err
---------------------------------------------------------------------
0.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
[INFO][rank 0][ARKodeEvolve][begin-step-attempt] step = 1, tn = 0, h = 6.103515625000001e-12
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 0, tcur = 0
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 1, tcur = 1.017252604166667e-12
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 2, tcur = 2.034505208333333e-12
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 3, tcur = 3.0517578125e-12
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 4, tcur = 4.069010416666667e-12
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 5, tcur = 5.086263020833334e-12
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 6, tcur = 2.034505208333333e-12
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 7, tcur = 3.0517578125e-12
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 8, tcur = 4.069010416666667e-12
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 9, tcur = 5.086263020833334e-12
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 10, tcur = 6.103515625000001e-12
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][ARKodeEvolve][end-step-attempt] status = success, dsm = 0
6.103515625000001e-12 6.103515625000001e-12 8.077935669463161e-28
[INFO][rank 0][ARKodeEvolve][begin-step-attempt] step = 2, tn = 6.103515625000001e-12, h = 5.025561905806021e-10
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 0, tcur = 6.103515625000001e-12
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 1, tcur = 8.986288072176701e-11
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 2, tcur = 1.73622245818534e-10
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 3, tcur = 2.57381610915301e-10
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 4, tcur = 3.41140976012068e-10
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 5, tcur = 4.24900341108835e-10
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 6, tcur = 1.73622245818534e-10
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 7, tcur = 2.57381610915301e-10
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 8, tcur = 3.41140976012068e-10
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 9, tcur = 4.24900341108835e-10
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 10, tcur = 5.086597062056021e-10
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][ARKodeEvolve][end-step-attempt] status = success, dsm = 2.067951405164832e-15
5.086597062056021e-10 5.086597062056018e-10 3.101927297073854e-25
[INFO][rank 0][ARKodeEvolve][begin-step-attempt] step = 3, tn = 5.086597062056021e-10, h = 8.256370880374143e-09
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 0, tcur = 5.086597062056021e-10
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 1, tcur = 1.884721519601293e-09
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 2, tcur = 3.260783332996983e-09
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 3, tcur = 4.636845146392673e-09
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 4, tcur = 6.012906959788364e-09
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 5, tcur = 7.388968773184054e-09
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 6, tcur = 3.260783332996983e-09
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 7, tcur = 4.636845146392673e-09
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 8, tcur = 6.012906959788364e-09
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 9, tcur = 7.388968773184054e-09
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 10, tcur = 8.765030586579745e-09
[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success
[INFO][rank 0][ARKodeEvolve][end-step-attempt] status = success, dsm = 3.308705620159824e-14
8.765030586579745e-09 8.765030586579742e-09 3.308722450212111e-24
---------------------------------------------------------------------
Current time = 8.765030586579745e-09
Steps = 3
Step attempts = 3
Stability limited steps = 0
Accuracy limited steps = 3
Error test fails = 0
NLS step fails = 0
Inequality constraint fails = 0
Initial step size = 6.103515625000001e-12
Last step size = 8.256370880374143e-09
Current step size = 1.651274176074829e-07
RHS fn evals = 32
Number of stages used = 10
End LSRKStep Logging test
Loading

0 comments on commit f95edba

Please sign in to comment.