| 2344 | ******************************************************************/ |
| 2345 | |
| 2346 | static void CVChooseEta(CVodeMem cv_mem, real etaqm1, real etaq, real etaqp1) |
| 2347 | { |
| 2348 | real etam; |
| 2349 | |
| 2350 | etam = MAX(etaqm1, MAX(etaq, etaqp1)); |
| 2351 | |
| 2352 | if (etam < THRESH) { |
| 2353 | eta = ONE; |
| 2354 | qprime = q; |
| 2355 | return; |
| 2356 | } |
| 2357 | |
| 2358 | if (etam == etaq) { |
| 2359 | eta = etaq; |
| 2360 | qprime = q; |
| 2361 | } else if (etam == etaqm1) { |
| 2362 | eta = etaqm1; |
| 2363 | qprime = q - 1; |
| 2364 | } else { |
| 2365 | eta = etaqp1; |
| 2366 | qprime = q + 1; |
| 2367 | N_VScale(ONE, acor, zn[qmax]); |
| 2368 | } |
| 2369 | } |
| 2370 | |
| 2371 | /****************** CVHandleFailure ****************************** |
| 2372 |
no test coverage detected