Module:Temperament data: Difference between revisions

@@ Line 67: / Line 67: @@
 end
-local function matinv(a)
+function p.matinv(a)
-	xk = {}
+	dbl_identity = {}
 	for i = 1, #a do
-		xk[i] = {}
+		dbl_identity[i] = {}
-		for j = 1, #(a[1]) do
+		for j = 1, #a do
-			xk[i][j] = {math.random(), math.random()}
+			if i == j then
+				dbl_identity[i][j] = {2, 0}
+			else
+				dbl_identity[i][j] = {0, 0}
+			end
 		end
 	end
-	for i = 1, 15 do
-		xk = matsub(scalarmatmul({2,0}, xk), matmul(xk, axk))
+	xn = scalarmatmul(a, {0.0001, 0})
+	for i = 1, 30 do
+		xn = matmul(xn, matsub(dbl_identity, matmul(a, xn)))
 	end
-	return xk
+	return xn
 end
 return p