|
Line
Link Here
|
| 0 |
-- a/src/calf/orfanidis_eq.h |
0 |
++ b/src/calf/orfanidis_eq.h |
|
Lines 748-754
private:
Link Here
|
| 748 |
std::vector<eq_double_t> v = landen(k, tol); |
748 |
std::vector<eq_double_t> v = landen(k, tol); |
| 749 |
|
749 |
|
| 750 |
std::transform(v.begin(), v.end(), v.begin(), |
750 |
std::transform(v.begin(), v.end(), v.begin(), |
| 751 |
bind2nd(std::plus<eq_double_t>(), 1.0)); |
751 |
bind(std::plus<eq_double_t>(), 1.0, std::placeholders::_1)); |
| 752 |
|
752 |
|
| 753 |
K = std::accumulate(v.begin(), v.end(), |
753 |
K = std::accumulate(v.begin(), v.end(), |
| 754 |
1, std::multiplies<eq_double_t>()) * M_PI/2.0; |
754 |
1, std::multiplies<eq_double_t>()) * M_PI/2.0; |
|
Lines 764-770
private:
Link Here
|
| 764 |
std::vector<eq_double_t> vp = landen(kp, tol); |
764 |
std::vector<eq_double_t> vp = landen(kp, tol); |
| 765 |
|
765 |
|
| 766 |
std::transform(vp.begin(), vp.end(), vp.begin(), |
766 |
std::transform(vp.begin(), vp.end(), vp.begin(), |
| 767 |
bind2nd(std::plus<eq_double_t>(), 1.0)); |
767 |
bind(std::plus<eq_double_t>(), 1.0, std::placeholders::_1)); |
| 768 |
|
768 |
|
| 769 |
Kprime = std::accumulate(vp.begin(), vp.end(), |
769 |
Kprime = std::accumulate(vp.begin(), vp.end(), |
| 770 |
1.0, std::multiplies<eq_double_t>()) * M_PI/2.0; |
770 |
1.0, std::multiplies<eq_double_t>()) * M_PI/2.0; |
| 771 |
- |
|
|
