diff --git a/src/arkode/arkode_lsrkstep.c b/src/arkode/arkode_lsrkstep.c index 3680f8b7b2..d57eaa8cd2 100644 --- a/src/arkode/arkode_lsrkstep.c +++ b/src/arkode/arkode_lsrkstep.c @@ -541,7 +541,8 @@ int lsrkStep_TakeStepRKC(ARKodeMem ark_mem, sunrealtype* dsmPtr, int* nflagPtr) #endif /* Compute forcing function, if necessary. */ - if ((!ark_mem->fn_is_current && ark_mem->initsetup) || (ark_mem->tn != step_mem->tnext)) + if ((!ark_mem->fn_is_current && ark_mem->initsetup) || + (ark_mem->tn != step_mem->tnext)) { retval = step_mem->fe(ark_mem->tn, ark_mem->yn, ark_mem->fn, ark_mem->user_data); @@ -824,7 +825,8 @@ int lsrkStep_TakeStepRKL(ARKodeMem ark_mem, sunrealtype* dsmPtr, int* nflagPtr) #endif /* Compute forcing function, if necessary. */ - if ((!ark_mem->fn_is_current && ark_mem->initsetup) || (ark_mem->tn != step_mem->tnext)) + if ((!ark_mem->fn_is_current && ark_mem->initsetup) || + (ark_mem->tn != step_mem->tnext)) { retval = step_mem->fe(ark_mem->tn, ark_mem->yn, ark_mem->fn, ark_mem->user_data); @@ -1056,7 +1058,7 @@ int lsrkStep_TakeStepSSPs2(ARKodeMem ark_mem, sunrealtype* dsmPtr, int* nflagPtr /* The method is not FSAL, and fn ​is used as storage in each step. Therefore, fn ​is computed at the beginning regardless. */ retval = step_mem->fe(ark_mem->tn, ark_mem->yn, ark_mem->fn, - ark_mem->user_data); + ark_mem->user_data); step_mem->nfe++; if (retval != ARK_SUCCESS) { return (ARK_RHSFUNC_FAIL); } @@ -1251,7 +1253,6 @@ int lsrkStep_TakeStepSSPs3(ARKodeMem ark_mem, sunrealtype* dsmPtr, int* nflagPtr ark_mem->nst, ark_mem->h, ark_mem->tcur); #endif - /* The method is not FSAL, and fn ​is used as storage in each step. Therefore, fn ​is computed at the beginning regardless. */ retval = step_mem->fe(ark_mem->tn, ark_mem->yn, ark_mem->fn, ark_mem->user_data); diff --git a/src/arkode/arkode_lsrkstep_io.c b/src/arkode/arkode_lsrkstep_io.c index 4f3e83752d..9606501ad5 100644 --- a/src/arkode/arkode_lsrkstep_io.c +++ b/src/arkode/arkode_lsrkstep_io.c @@ -54,7 +54,7 @@ int LSRKStepSetMethod(void* arkode_mem, ARKODE_LSRKMethodType method) step_mem->nfusedopvecs = 5; step_mem->q = ark_mem->hadapt_mem->q = 2; step_mem->p = ark_mem->hadapt_mem->p = 2; - step_mem->tnext = ZERO; + step_mem->tnext = ZERO; break; case ARKODE_LSRK_RKL_2: ark_mem->step = lsrkStep_TakeStepRKL; @@ -62,7 +62,7 @@ int LSRKStepSetMethod(void* arkode_mem, ARKODE_LSRKMethodType method) step_mem->nfusedopvecs = 5; step_mem->q = ark_mem->hadapt_mem->q = 2; step_mem->p = ark_mem->hadapt_mem->p = 2; - step_mem->tnext = ZERO; + step_mem->tnext = ZERO; break; case ARKODE_LSRK_SSP_S_2: ark_mem->step = lsrkStep_TakeStepSSPs2; @@ -446,7 +446,7 @@ int lsrkStep_SetDefaults(ARKodeMem ark_mem) (overwrite some adaptivity params for LSRKStep use) */ step_mem->req_stages = 0; /* no stages */ -/* Spectral info */ + /* Spectral info */ step_mem->lambdaR = ZERO; step_mem->lambdaI = ZERO; step_mem->spectral_radius = ZERO;