Skip to content

Commit

Permalink
reorder logging output
Browse files Browse the repository at this point in the history
  • Loading branch information
gardner48 committed Aug 6, 2024
1 parent 8afc896 commit d37e678
Show file tree
Hide file tree
Showing 3 changed files with 35 additions and 32 deletions.
48 changes: 27 additions & 21 deletions src/arkode/arkode_arkstep.c
Original file line number Diff line number Diff line change
Expand Up @@ -1985,33 +1985,39 @@ int arkStep_TakeStep_Z(ARKodeMem ark_mem, sunrealtype* dsmPtr, int* nflagPtr)
retval = step_mem->fi(ark_mem->tcur, ark_mem->ycur, step_mem->Fi[is],
ark_mem->user_data);
step_mem->nfi++;

SUNLogExtraDebugVec(ARK_LOGGER, __func__, "implicit RHS",
"Fi_%i(:) =", step_mem->Fi[is], is);
SUNLogInfoIf(retval != 0, ARK_LOGGER, __func__, "end-stage",
"status = failed implicit rhs eval, retval = %i", retval);
}
else if (step_mem->mass_type == MASS_FIXED)
{
retval = step_mem->mmult((void*)ark_mem, step_mem->zcor, ark_mem->tempv1);
if (retval != ARK_SUCCESS)
{
SUNLogInfo(ARK_LOGGER, __func__, "end-stage",
"status = failed mass mult, retval = %i", retval);
return (ARK_MASSMULT_FAIL);
}

N_VLinearSum(ONE / step_mem->gamma, ark_mem->tempv1,
-ONE / step_mem->gamma, step_mem->sdata, step_mem->Fi[is]);
if (retval < 0) { return (ARK_RHSFUNC_FAIL); }
if (retval > 0) { return (ARK_UNREC_RHSFUNC_ERR); }
}
else
{
N_VLinearSum(ONE / step_mem->gamma, step_mem->zcor,
-ONE / step_mem->gamma, step_mem->sdata, step_mem->Fi[is]);
}
if (step_mem->mass_type == MASS_FIXED)
{
retval = step_mem->mmult((void*)ark_mem, step_mem->zcor, ark_mem->tempv1);
if (retval != ARK_SUCCESS)
{
SUNLogInfo(ARK_LOGGER, __func__, "end-stage",
"status = failed mass mult, retval = %i", retval);
return (ARK_MASSMULT_FAIL);
}

SUNLogExtraDebugVec(ARK_LOGGER, __func__, "implicit RHS",
"Fi_%i(:) =", step_mem->Fi[is], is);
N_VLinearSum(ONE / step_mem->gamma, ark_mem->tempv1,
-ONE / step_mem->gamma, step_mem->sdata, step_mem->Fi[is]);
}
else
{
N_VLinearSum(ONE / step_mem->gamma, step_mem->zcor,
-ONE / step_mem->gamma, step_mem->sdata, step_mem->Fi[is]);
}

if (retval < 0) { return (ARK_RHSFUNC_FAIL); }
if (retval > 0) { return (ARK_UNREC_RHSFUNC_ERR); }
SUNLogExtraDebugVec(ARK_LOGGER, __func__, "implicit RHS",
"Fi_%i(:) =", step_mem->Fi[is], is);
}
}

/* store explicit RHS */
Expand All @@ -2020,11 +2026,11 @@ int arkStep_TakeStep_Z(ARKodeMem ark_mem, sunrealtype* dsmPtr, int* nflagPtr)
retval = step_mem->fe(ark_mem->tn + step_mem->Be->c[is] * ark_mem->h,
ark_mem->ycur, step_mem->Fe[is], ark_mem->user_data);
step_mem->nfe++;
SUNLogInfoIf(retval != 0, ARK_LOGGER, __func__, "end-stage",
"status = failed explicit rhs eval, retval = %i", retval);

SUNLogExtraDebug(ARK_LOGGER, __func__, "explicit RHS",
"Fe_%i(:) =", step_mem->Fe[is], is);
SUNLogInfoIf(retval != 0, ARK_LOGGER, __func__, "end-stage",
"status = failed explicit rhs eval, retval = %i", retval);

if (retval < 0) { return (ARK_RHSFUNC_FAIL); }
if (retval > 0) { return (ARK_UNREC_RHSFUNC_ERR); }
Expand Down
5 changes: 2 additions & 3 deletions src/arkode/arkode_erkstep.c
Original file line number Diff line number Diff line change
Expand Up @@ -698,15 +698,14 @@ int erkStep_TakeStep(ARKodeMem ark_mem, sunrealtype* dsmPtr, int* nflagPtr)
ark_mem->user_data);
step_mem->nfe++;

SUNLogExtraDebugVec(ARK_LOGGER, __func__, "stage RHS",
"F_%i(:) =", step_mem->F[is], is);
SUNLogInfoIf(retval != 0, ARK_LOGGER, __func__, "end-stage",
"status = failed rhs eval, retval = %i", retval);

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

SUNLogExtraDebugVec(ARK_LOGGER, __func__, "stage RHS",
"F_%i(:) =", step_mem->F[is], is);

SUNLogInfo(ARK_LOGGER, __func__, "end-stage", "status = success", "");

} /* loop over stages */
Expand Down
14 changes: 6 additions & 8 deletions src/arkode/arkode_mristep.c
Original file line number Diff line number Diff line change
Expand Up @@ -1583,15 +1583,14 @@ int mriStep_TakeStep(ARKodeMem ark_mem, sunrealtype* dsmPtr, int* nflagPtr)
ark_mem->user_data);
step_mem->nfse++;

SUNLogExtraDebugVec(ARK_LOGGER, __func__, "slow explicit RHS",
"Fse_%i(:) =", step_mem->Fse[step_mem->stage_map[is]],
is);
SUNLogInfoIf(retval != 0, ARK_LOGGER, __func__, "end-stage",
"status = failed explicit rhs eval, retval = %i", retval);

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

SUNLogExtraDebugVec(ARK_LOGGER, __func__, "slow explicit RHS",
"Fse_%i(:) =", step_mem->Fse[step_mem->stage_map[is]],
is);
}

/* store implicit slow rhs */
Expand All @@ -1612,15 +1611,14 @@ int mriStep_TakeStep(ARKodeMem ark_mem, sunrealtype* dsmPtr, int* nflagPtr)
step_mem->Fsi[step_mem->stage_map[is]]);
}

SUNLogExtraDebugVec(ARK_LOGGER, __func__, "slow implicit RHS",
"Fsi_%i(:) =", step_mem->Fsi[step_mem->stage_map[is]],
is);
SUNLogInfoIf(retval != 0, ARK_LOGGER, __func__, "end-stage",
"status = failed implicit rhs eval, retval = %i", retval);

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

SUNLogExtraDebugVec(ARK_LOGGER, __func__, "slow implicit RHS",
"Fsi_%i(:) =", step_mem->Fsi[step_mem->stage_map[is]],
is);
}

} /* compute slow RHS */
Expand Down

0 comments on commit d37e678

Please sign in to comment.