| # | Line 59 | Line 59 | namespace oopse { | |
|---|---|---|
| 59 | //P_{l+1}= \frac{(2l+1)(x)P_l-l P_{l-1}{l+1} | |
| 60 | for (int i = 2; i <= maxPower; ++i) { | |
| 61 | DoublePolynomial pn; | |
| 62 | < | |
| 63 | < | pn = polyList_[i-1] * x * static_cast<double>((2*i+1)/(i+1)) - polyList_[i-2] * static_cast<double>(i/(i+1)); |
| 62 | > | RealType tmp1 = (2.0*i-1.0)/i; |
| 63 | > | RealType tmp2 = (i-1.0)/i; |
| 64 | > | pn = polyList_[i-1] * x * tmp1 - polyList_[i-2] * tmp2; |
| 65 | polyList_.push_back(pn); | |
| 66 | } | |
| 67 | } | |
| – | Removed lines |
| + | Added lines |
| < | Changed lines |
| > | Changed lines |