Skip to content

Commit

Permalink
Add full RHS mode back
Browse files Browse the repository at this point in the history
  • Loading branch information
Steven-Roberts committed Oct 30, 2024
1 parent b8a0055 commit 09bd9da
Show file tree
Hide file tree
Showing 2 changed files with 4 additions and 3 deletions.
4 changes: 2 additions & 2 deletions src/arkode/arkode_forcingstep.c
Original file line number Diff line number Diff line change
Expand Up @@ -134,15 +134,15 @@ static int forcingStep_FullRHS(const ARKodeMem ark_mem, const sunrealtype t,
if (retval != ARK_SUCCESS) { return retval; }

SUNErrCode err = SUNStepper_FullRhs(step_mem->stepper[0], t, y,
ark_mem->tempv1);
ark_mem->tempv1, SUN_FULLRHS_OTHER);
if (err != SUN_SUCCESS)
{
arkProcessError(ark_mem, ARK_RHSFUNC_FAIL, __LINE__, __func__, __FILE__,
MSG_ARK_RHSFUNC_FAILED, t);
return ARK_RHSFUNC_FAIL;
}

err = SUNStepper_FullRhs(step_mem->stepper[1], t, y, f);
err = SUNStepper_FullRhs(step_mem->stepper[1], t, y, f, SUN_FULLRHS_OTHER);
if (err != SUN_SUCCESS)
{
arkProcessError(ark_mem, ARK_RHSFUNC_FAIL, __LINE__, __func__, __FILE__,
Expand Down
3 changes: 2 additions & 1 deletion src/arkode/arkode_splittingstep.c
Original file line number Diff line number Diff line change
Expand Up @@ -188,7 +188,8 @@ static int splittingStep_FullRHS(const ARKodeMem ark_mem, const sunrealtype t,
for (int i = 0; i < step_mem->partitions; i++)
{
const SUNErrCode err = SUNStepper_FullRhs(step_mem->steppers[i], t, y,
i == 0 ? f : ark_mem->tempv1);
i == 0 ? f : ark_mem->tempv1,
SUN_FULLRHS_OTHER);
if (err != SUN_SUCCESS)
{
arkProcessError(ark_mem, ARK_RHSFUNC_FAIL, __LINE__, __func__, __FILE__,
Expand Down

0 comments on commit 09bd9da

Please sign in to comment.