Module:Temperament data/Badness testing: Difference between revisions

← Older edit

@@ Line 68: / Line 68: @@
 local function matinv(a)
-	local dbl_identity = {}
+	local xn = scalarmatmul(a, 1e-7)
-	for i = 1, #a do
-		dbl_identity[i] = {}
+	for i = 1, 70 do
-		for j = 1, #a do
+		xn = matsub(scalarmatmul(xn, 2), matmul(xn, matmul(a, xn)))
-			if i == j then
-				dbl_identity[i][j] = 2
-			else
-				dbl_identity[i][j] = 0
-			end
-		end
 	end
-	local xn = scalarmatmul(a, 0.000005)
-	for i = 1, 30 do
-		xn = matmul(xn, matsub(dbl_identity, matmul(a, xn)))
-	end
 	return xn
 end