floor ratio

jacobdenobel · jacobdenobel · commit 1e1b572a2b0f · 2025-06-12T14:16:05.000+02:00
diff --git a/src/weights.cpp b/src/weights.cpp
@@ -156,7 +156,7 @@ namespace parameters
 	void Weights::weights_default(const size_t mu, const size_t lambda)
 	{
 		const Float ratio = static_cast<Float>(lambda) / static_cast<Float>(mu);
-		const Float base = std::log((static_cast<Float>(lambda) + 1.) / std::round(ratio));
+		const Float base = std::log((static_cast<Float>(lambda) + 1.) / std::floor(ratio));
 
 		for (auto i = 0; i < positive.size(); ++i)
 			positive(i) = base - std::log(static_cast<Float>(i + 1));

Original file line number	Diff line number	Diff line change
`@@ -156,7 +156,7 @@ namespace parameters`
`156`	`156`	`void Weights::weights_default(const size_t mu, const size_t lambda)`
`157`	`157`	`{`
`158`	`158`	`const Float ratio = static_cast<Float>(lambda) / static_cast<Float>(mu);`
`159`		`- const Float base = std::log((static_cast<Float>(lambda) + 1.) / std::round(ratio));`
	`159`	`+ const Float base = std::log((static_cast<Float>(lambda) + 1.) / std::floor(ratio));`
`160`	`160`
`161`	`161`	`for (auto i = 0; i < positive.size(); ++i)`
`162`	`162`	`positive(i) = base - std::log(static_cast<Float>(i + 1));`