Module:Limits: Difference between revisions

← Older edit Newer edit →

@@ Line 1: / Line 1: @@
 local rat = require('Module:Rational')
 local p = {}
--- multiply the ratio by a power of `equave` so that it lies within [1; equave)
-function p.canonical(a, equave)
-	equave = equave or 2
-	if type(a) == 'number' then
-		a = rat.new(a)
-	end
-	local b = rat.copy(a)
-	while rat.lt(b, 1) do
-		b = rat.mul(b, equave)
-	end
-	while rat.geq(b, equave) do
-		b = rat.div(b, equave)
-	end
-	return b
-end
 -- compute all positive ratios n/m with n and m <= q modulo powers of equave
@@ Line 25: / Line 9: @@
 		for m = 1, q, 2 do
 			local a = rat.new(n, m)
-			a = p.canonical(a, equave)
+			a = rat.modulo_mul(a, equave)
 			local key = rat.as_ratio(a)
 			ratios[key] = a
@@ Line 44: / Line 28: @@
 		local approx_set = {}
 		for a_key, a in pairs(ratios) do
-			local a_approx = approximate(a)
+			local a_approx = approximate(a) % size
 			if approx_set[a_approx] then
 				return false
@@ Line 64: / Line 48: @@
 				local c_approx = approximate(c)
-				c = p.canonical(c, equave)
+				c = rat.modulo_mul(c, equave)
 				local c_key = rat.as_ratio(c)
 				if ratios[c_key] then