We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent 38b2ac3 commit e52617cCopy full SHA for e52617c
include/gemmi.hpp
@@ -104,9 +104,7 @@ struct MatrixSplit {
104
// to shift correctly.
105
frexp(this->powersVector[i], this->scalingExponents.data() + i);
106
const auto largest_log = log2(this->powersVector[i]);
107
- auto temp = ceil(largest_log) +
108
- (ceil(largest_log) == floor(largest_log) ? 1 : 0);
109
- this->powersVector[i] = std::ldexp(1.0, temp);
+ this->powersVector[i] = std::ldexp(1.0, floor(largest_log) + 1);
110
}
111
112
0 commit comments