Module:Limits: Difference between revisions

← Older edit Newer edit →

@@ Line 2: / Line 2: @@
 local p = {}
--- multiply the ratio by a power of `equave` so that it lies within [1; equave)
+p.rat = rat
-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 11: @@
 		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 64: / Line 50: @@
 				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