Module:Harmonic entropy: Difference between revisions

← Older edit Newer edit →

@@ Line 12: / Line 12: @@
 		return math.sqrt(rat.benedetti_height(ratio))
 	end
-	deviation = deviation or 0.01
+	deviation = deviation or 17.22635157248405
 	ratios = ratios
 		or limits.integer_limit(200, function(ratio)
-			if math.abs(rat.as_float(ratio) * 2^(rat.as_float(c) / 1200)) > 3 * deviation then
+			if math.abs(rat.cents(ratio) - c) > 3 * deviation then
 				return 1 / 0
 			end
@@ Line 26: / Line 26: @@
 	local function weighted_gaussian(ratio)
-		return gaussian(rat.as_float(ratio) * 2^(rat.as_float(c) / 1200)) / norm(ratio)
+		return gaussian(rat.cents(ratio) - c) / norm(ratio)
 	end