Skip to content
Snippets Groups Projects
Commit 0581856e authored by Ian Bell's avatar Ian Bell
Browse files

No inplace addition for MCX :(

parent ffb88608
No related branches found
No related tags found
No related merge requests found
...@@ -10,7 +10,7 @@ public: ...@@ -10,7 +10,7 @@ public:
using result = std::common_type_t<TauType, DeltaType>; using result = std::common_type_t<TauType, DeltaType>;
result r = 0.0, lntau = log(tau), lndelta = log(delta); result r = 0.0, lntau = log(tau), lndelta = log(delta);
for (auto i = 0; i < n.size(); ++i) { for (auto i = 0; i < n.size(); ++i) {
r += n[i] * exp(t[i]*lntau + d[i]*lndelta -c[i]*powi(delta, l_i[i])); r = r + n[i] * exp(t[i]*lntau + d[i]*lndelta -c[i]*powi(delta, l_i[i]));
} }
return r; return r;
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment