| 1680 | ******************************************************************/ |
| 1681 | |
| 1682 | static void CVAdamsFinish(CVodeMem cv_mem, real m[], real M[], real hsum) |
| 1683 | { |
| 1684 | int i; |
| 1685 | real M0_inv, xi, xi_inv; |
| 1686 | |
| 1687 | M0_inv = ONE / M[0]; |
| 1688 | |
| 1689 | l[0] = ONE; |
| 1690 | for (i=1; i <= q; i++) l[i] = M0_inv * (m[i-1] / i); |
| 1691 | xi = hsum / h; |
| 1692 | xi_inv = ONE / xi; |
| 1693 | |
| 1694 | tq[2] = xi * M[0] / M[1]; |
| 1695 | tq[5] = xi / l[q]; |
| 1696 | |
| 1697 | if (qwait == 1) { |
| 1698 | for (i=q; i >= 1; i--) m[i] += m[i-1] * xi_inv; |
| 1699 | M[2] = CVAltSum(q, m, 2); |
| 1700 | tq[3] = L * M[0] / M[2]; |
| 1701 | } |
| 1702 | |
| 1703 | tq[4] = CORTES * tq[2]; |
| 1704 | } |
| 1705 | |
| 1706 | /****************** CVAltSum ************************************** |
| 1707 |
no test coverage detected