MCPcopy Create free account
hub / github.com/boutproject/BOUT-dev / CVComputeEtaqp1

Function CVComputeEtaqp1

externalpackages/PVODE/source/cvode.cpp:2319–2331  ·  view source on GitHub ↗

Source from the content-addressed store, hash-verified

2317******************************************************************/
2318
2319static real CVComputeEtaqp1(CVodeMem cv_mem)
2320{
2321 real etaqp1, dup, cquot;
2322
2323 etaqp1 = ZERO;
2324 if (q != qmax) {
2325 cquot = (tq[5] / saved_tq5) * RPowerI(h/tau[2], L);
2326 N_VLinearSum(-cquot, zn[qmax], ONE, acor, tempv);
2327 dup = N_VWrmsNorm(tempv, ewt) /tq[3];
2328 etaqp1 = ONE / (RPowerR(BIAS3*dup, ONE/(L+1)) + ADDON);
2329 }
2330 return(etaqp1);
2331}
2332
2333/******************* CVChooseEta **********************************
2334

Callers 1

CVPrepareNextStepFunction · 0.85

Calls 4

RPowerIFunction · 0.85
N_VLinearSumFunction · 0.85
N_VWrmsNormFunction · 0.85
RPowerRFunction · 0.85

Tested by

no test coverage detected