Module:Temperament data: Difference between revisions

← Older edit Newer edit →

@@ Line 65: / Line 65: @@
 end
-local function matinv(a)
+function p.matinv(a)
-	local dbl_identity = {}
+	-- local dbl_identity = {}
-	for i = 1, #a do
+	-- for i = 1, #a do
-		dbl_identity[i] = {}
+	-- 	dbl_identity[i] = {}
-		for j = 1, #a do
+	-- 	for j = 1, #a do
-			if i == j then
+	-- 		if i == j then
-				dbl_identity[i][j] = 2
+	-- 			dbl_identity[i][j] = 2
-			else
+	-- 		else
-				dbl_identity[i][j] = 0
+	-- 			dbl_identity[i][j] = 0
-			end
+	-- 		end
-		end
+	-- 	end
-	end
+	-- end
+	local xn = scalarmatmul(a, 0.01)
-	local xn = scalarmatmul(a, 2)
-	-- for i = 1, 10 do
-	-- 	xn = matmul(xn, matsub(dbl_identity, matmul(a, xn)))
-	-- end
 	for i = 1, 20 do
 		xn = matsub(scalarmatmul(xn, 2), matmul(xn, matmul(a, xn)))
+		mw.logObject(xn)
 	end