Go to the source code of this file.
Definition at line 22 of file dcfode.cpp.
References C_INT.
Referenced by CInternalSolver::dstoda_().
33 double ragq, pint, xpin, fnqm1, agamq, rqfac, tsign, rq1fac;
98 for (nq = 2; nq <= 12; ++nq)
108 fnqm1 = (double) nqm1;
113 for (ib = 1; ib <= i__1; ++ib)
117 pc[i__ - 1] = pc[i__ - 2] + fnqm1 * pc[i__ - 1];
119 pc[0] = fnqm1 * pc[0];
125 for (i__ = 2; i__ <= i__1; ++i__)
128 pint += tsign * pc[i__ - 1] / i__;
130 xpin += tsign * pc[i__ - 1] / (i__ + 1);
133 elco[nq * 13 + 1] = pint * rq1fac;
134 elco[nq * 13 + 2] = 1.;
136 for (i__ = 2; i__ <= i__1; ++i__)
139 elco[i__ + 1 + nq * 13] = rq1fac * pc[i__ - 1] / i__;
141 agamq = rqfac * xpin;
143 tesco[nq * 3 + 2] = ragq;
146 tesco[nqp1 * 3 + 1] = ragq * rqfac / nqp1;
148 tesco[nqm1 * 3 + 3] = ragq;
156 for (nq = 1; nq <= 5; ++nq)
168 for (ib = 1; ib <= i__1; ++ib)
172 pc[i__ - 1] = pc[i__ - 2] + fnq * pc[i__ - 1];
177 for (i__ = 1; i__ <= i__1; ++i__)
180 elco[i__ + nq * 13] = pc[i__ - 1] / pc[1];
182 elco[nq * 13 + 2] = 1.;
183 tesco[nq * 3 + 1] = rq1fac;
184 tesco[nq * 3 + 2] = nqp1 / elco[nq * 13 + 1];
185 tesco[nq * 3 + 3] = (nq + 2) / elco[nq * 13 + 1];