From f95edbabde134eb13bcee50cc2388fc98bf3de31 Mon Sep 17 00:00:00 2001 From: "David J. Gardner" Date: Thu, 21 Nov 2024 17:40:50 -0800 Subject: [PATCH] add logging to SSP104 --- src/arkode/arkode_lsrkstep.c | 146 ++++++------ test/unit_tests/logging/CMakeLists.txt | 1 + .../logging/test_logging_arkode_lsrkstep.cpp | 2 +- .../test_logging_arkode_lsrkstep_5.out | 22 ++ .../test_logging_arkode_lsrkstep_lvl3_5.out | 94 ++++++++ .../test_logging_arkode_lsrkstep_lvl4_5.out | 103 ++++++++ .../test_logging_arkode_lsrkstep_lvl5_5.out | 220 ++++++++++++++++++ 7 files changed, 510 insertions(+), 78 deletions(-) create mode 100644 test/unit_tests/logging/test_logging_arkode_lsrkstep_5.out create mode 100644 test/unit_tests/logging/test_logging_arkode_lsrkstep_lvl3_5.out create mode 100644 test/unit_tests/logging/test_logging_arkode_lsrkstep_lvl4_5.out create mode 100644 test/unit_tests/logging/test_logging_arkode_lsrkstep_lvl5_5.out diff --git a/src/arkode/arkode_lsrkstep.c b/src/arkode/arkode_lsrkstep.c index 2072a2faea..dec509507f 100644 --- a/src/arkode/arkode_lsrkstep.c +++ b/src/arkode/arkode_lsrkstep.c @@ -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. */ @@ -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); @@ -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); @@ -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); @@ -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); @@ -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; @@ -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) @@ -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; } diff --git a/test/unit_tests/logging/CMakeLists.txt b/test/unit_tests/logging/CMakeLists.txt index 1fd73ca943..9c2d2e6aed 100644 --- a/test/unit_tests/logging/CMakeLists.txt +++ b/test/unit_tests/logging/CMakeLists.txt @@ -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 diff --git a/test/unit_tests/logging/test_logging_arkode_lsrkstep.cpp b/test/unit_tests/logging/test_logging_arkode_lsrkstep.cpp index 12df07d61e..d2f7062374 100644 --- a/test/unit_tests/logging/test_logging_arkode_lsrkstep.cpp +++ b/test/unit_tests/logging/test_logging_arkode_lsrkstep.cpp @@ -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 { diff --git a/test/unit_tests/logging/test_logging_arkode_lsrkstep_5.out b/test/unit_tests/logging/test_logging_arkode_lsrkstep_5.out new file mode 100644 index 0000000000..74c20b7c01 --- /dev/null +++ b/test/unit_tests/logging/test_logging_arkode_lsrkstep_5.out @@ -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 diff --git a/test/unit_tests/logging/test_logging_arkode_lsrkstep_lvl3_5.out b/test/unit_tests/logging/test_logging_arkode_lsrkstep_lvl3_5.out new file mode 100644 index 0000000000..87dad7e353 --- /dev/null +++ b/test/unit_tests/logging/test_logging_arkode_lsrkstep_lvl3_5.out @@ -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 diff --git a/test/unit_tests/logging/test_logging_arkode_lsrkstep_lvl4_5.out b/test/unit_tests/logging/test_logging_arkode_lsrkstep_lvl4_5.out new file mode 100644 index 0000000000..eb0354adac --- /dev/null +++ b/test/unit_tests/logging/test_logging_arkode_lsrkstep_lvl4_5.out @@ -0,0 +1,103 @@ +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 +[DEBUG][rank 0][arkAdapt][new-step-before-bounds] h_acc = 5.234960318547938e-10, h_cfl = 6.103515625000001e+18 +[DEBUG][rank 0][arkAdapt][new-step-after-max-min-bounds] h_acc = 5.025561905806021e-10, h_cfl = 3.0517578125e+18 +[DEBUG][rank 0][arkAdapt][new-step-eta] eta = 82.33880626472583 +[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 +[DEBUG][rank 0][arkAdapt][new-step-before-bounds] h_acc = 8.600386333723065e-09, h_cfl = 5.02556190580602e+20 +[DEBUG][rank 0][arkAdapt][new-step-after-max-min-bounds] h_acc = 8.256370880374143e-09, h_cfl = 2.51278095290301e+20 +[DEBUG][rank 0][arkAdapt][new-step-eta] eta = 16.4287517199531 +[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 +[DEBUG][rank 0][arkAdapt][new-step-before-bounds] h_acc = 3.04407858682386e-07, h_cfl = 8.256370880374143e+21 +[DEBUG][rank 0][arkAdapt][new-step-after-max-min-bounds] h_acc = 1.651274176074829e-07, h_cfl = 4.128185440187071e+21 +[DEBUG][rank 0][arkAdapt][new-step-eta] eta = 20 +[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 diff --git a/test/unit_tests/logging/test_logging_arkode_lsrkstep_lvl5_5.out b/test/unit_tests/logging/test_logging_arkode_lsrkstep_lvl5_5.out new file mode 100644 index 0000000000..6930759601 --- /dev/null +++ b/test/unit_tests/logging/test_logging_arkode_lsrkstep_lvl5_5.out @@ -0,0 +1,220 @@ +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 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage] z_0(:) = +0.0000000000000000e+00 + +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_0(:) = +1.0000000000000000e+00 + +[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success +[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 1, tcur = 1.017252604166667e-12 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_1(:) = +1.0000000000000000e+00 + +[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success +[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 2, tcur = 2.034505208333333e-12 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_2(:) = +1.0000000000000000e+00 + +[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success +[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 3, tcur = 3.0517578125e-12 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_3(:) = +1.0000000000000000e+00 + +[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success +[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 4, tcur = 4.069010416666667e-12 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_4(:) = +1.0000000000000000e+00 + +[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 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_6(:) = +1.0000000000000000e+00 + +[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success +[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 7, tcur = 3.0517578125e-12 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_7(:) = +1.0000000000000000e+00 + +[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success +[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 8, tcur = 4.069010416666667e-12 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_8(:) = +1.0000000000000000e+00 + +[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success +[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 9, tcur = 5.086263020833334e-12 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_9(:) = +1.0000000000000000e+00 + +[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success +[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 10, tcur = 6.103515625000001e-12 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_10(:) = +1.0000000000000000e+00 + +[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][updated solution] ycur(:) = +6.1035156250000013e-12 + +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][embedded solution] y_embedded(:) = +6.1035156250000013e-12 + +[DEBUG][rank 0][arkAdapt][new-step-before-bounds] h_acc = 5.234960318547938e-10, h_cfl = 6.103515625000001e+18 +[DEBUG][rank 0][arkAdapt][new-step-after-max-min-bounds] h_acc = 5.025561905806021e-10, h_cfl = 3.0517578125e+18 +[DEBUG][rank 0][arkAdapt][new-step-eta] eta = 82.33880626472583 +[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 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage] z_0(:) = +6.1035156250000013e-12 + +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_0(:) = +1.0000000000000000e+00 + +[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success +[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 1, tcur = 8.986288072176701e-11 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_1(:) = +1.0000000000000000e+00 + +[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success +[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 2, tcur = 1.73622245818534e-10 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_2(:) = +1.0000000000000000e+00 + +[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success +[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 3, tcur = 2.57381610915301e-10 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_3(:) = +1.0000000000000000e+00 + +[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success +[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 4, tcur = 3.41140976012068e-10 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_4(:) = +1.0000000000000000e+00 + +[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 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_6(:) = +1.0000000000000000e+00 + +[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success +[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 7, tcur = 2.57381610915301e-10 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_7(:) = +1.0000000000000000e+00 + +[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success +[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 8, tcur = 3.41140976012068e-10 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_8(:) = +1.0000000000000000e+00 + +[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success +[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 9, tcur = 4.24900341108835e-10 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_9(:) = +1.0000000000000000e+00 + +[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success +[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 10, tcur = 5.086597062056021e-10 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_10(:) = +1.0000000000000000e+00 + +[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][updated solution] ycur(:) = +5.0865970620560179e-10 + +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][embedded solution] y_embedded(:) = +5.0865970620560200e-10 + +[DEBUG][rank 0][arkAdapt][new-step-before-bounds] h_acc = 8.600386333723065e-09, h_cfl = 5.02556190580602e+20 +[DEBUG][rank 0][arkAdapt][new-step-after-max-min-bounds] h_acc = 8.256370880374143e-09, h_cfl = 2.51278095290301e+20 +[DEBUG][rank 0][arkAdapt][new-step-eta] eta = 16.4287517199531 +[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 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage] z_0(:) = +5.0865970620560179e-10 + +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_0(:) = +1.0000000000000000e+00 + +[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success +[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 1, tcur = 1.884721519601293e-09 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_1(:) = +1.0000000000000000e+00 + +[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success +[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 2, tcur = 3.260783332996983e-09 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_2(:) = +1.0000000000000000e+00 + +[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success +[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 3, tcur = 4.636845146392673e-09 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_3(:) = +1.0000000000000000e+00 + +[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success +[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 4, tcur = 6.012906959788364e-09 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_4(:) = +1.0000000000000000e+00 + +[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 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_6(:) = +1.0000000000000000e+00 + +[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success +[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 7, tcur = 4.636845146392673e-09 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_7(:) = +1.0000000000000000e+00 + +[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success +[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 8, tcur = 6.012906959788364e-09 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_8(:) = +1.0000000000000000e+00 + +[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success +[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 9, tcur = 7.388968773184054e-09 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_9(:) = +1.0000000000000000e+00 + +[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success +[INFO][rank 0][lsrkStep_TakeStepSSP104][begin-stage] stage = 10, tcur = 8.765030586579745e-09 +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][stage RHS] F_10(:) = +1.0000000000000000e+00 + +[INFO][rank 0][lsrkStep_TakeStepSSP104][end-stage] status = success +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][updated solution] ycur(:) = +8.7650305865797416e-09 + +[DEBUG][rank 0][lsrkStep_TakeStepSSP104][embedded solution] y_embedded(:) = +8.7650305865797449e-09 + +[DEBUG][rank 0][arkAdapt][new-step-before-bounds] h_acc = 3.04407858682386e-07, h_cfl = 8.256370880374143e+21 +[DEBUG][rank 0][arkAdapt][new-step-after-max-min-bounds] h_acc = 1.651274176074829e-07, h_cfl = 4.128185440187071e+21 +[DEBUG][rank 0][arkAdapt][new-step-eta] eta = 20 +[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