Commit 253133c0 authored by Gauthier Quesnel's avatar Gauthier Quesnel
Browse files

core: remove useless code in QSS3 transition

parent c34c0c63
...@@ -3515,7 +3515,6 @@ struct qss3_integrator ...@@ -3515,7 +3515,6 @@ struct qss3_integrator
auto w = c - b * a / 3. + 2. * a * a * a / 27.; auto w = c - b * a / 3. + 2. * a * a * a / 27.;
auto i1 = -w / 2.; auto i1 = -w / 2.;
auto i2 = i1 * i1 + v * v * v / 27.; auto i2 = i1 * i1 + v * v * v / 27.;
auto s = 0.;
if (i2 > 0) { if (i2 > 0) {
i2 = std::sqrt(i2); i2 = std::sqrt(i2);
...@@ -3531,7 +3530,7 @@ struct qss3_integrator ...@@ -3531,7 +3530,7 @@ struct qss3_integrator
else else
B = -std::pow(std::abs(B), 1. / 3.); B = -std::pow(std::abs(B), 1. / 3.);
auto s = A + B - a / 3.; s = A + B - a / 3.;
if (s < 0.) if (s < 0.)
s = time_domain<time>::infinity; s = time_domain<time>::infinity;
} else if (i2 == 0.) { } else if (i2 == 0.) {
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment