Module:Temperament data - Revision history

ArrowHead294 at 18:22, 19 May 2025

2025-05-19T18:22:14Z

← Older revision		Revision as of 18:22, 19 May 2025
Line 1:		Line 1:
			local p = {}

	local rat = require("Module:Rational")		local rat = require("Module:Rational")
	~~local p = {}~~
	local u = require("Module:Utils")		local u = require("Module:Utils")

	local function gcd(a,b)		local function gcd(a, b)
	if type(a) == "number" and type(b) == "number" and a == math.floor(a) and b == math.floor(b) then		if type(a) == "number" and type(b) == "number" and a == math.floor(a) and b == math.floor(b) then
	if b == 0 then		if b == 0 then

2024-07-22T03:10:16Z

← Older revision		Revision as of 03:10, 22 July 2024
Line 66:		Line 66:

	local function matinv(a)		local function matinv(a)
	local xn = scalarmatmul(a, ~~0.0002~~)		local xn = scalarmatmul(a, 1e-7)

	for i = 1, 50 do		for i = 1, 75 do
	xn = matsub(scalarmatmul(xn, 2), matmul(xn, matmul(a, xn)))		xn = matsub(scalarmatmul(xn, 2), matmul(xn, matmul(a, xn)))
	end		end

2024-07-22T03:09:41Z

2024-07-22T03:08:37Z

2024-07-22T03:08:05Z

← Older revision		Revision as of 03:08, 22 July 2024
Line 66:		Line 66:

	local function matinv(a)		local function matinv(a)
	~~-- local dbl_identity = {}~~
	~~-- for i = 1, #a do~~
	~~-- dbl_identity[i] = {}~~
	~~-- for j = 1, #a do~~
	~~-- if i == j then~~
	~~-- dbl_identity[i][j] = 2~~
	~~-- else~~
	~~-- dbl_identity[i][j] = 0~~
	~~-- end~~
	~~-- end~~
	~~-- end~~

	local xn = scalarmatmul(a, 0.0001)		local xn = scalarmatmul(a, 0.0001)

	for i = 1, 50 do		for i = 1, 50 do
	xn = matsub(scalarmatmul(xn, 2), matmul(xn, matmul(a, xn)))		xn = matsub(scalarmatmul(xn, 2), matmul(xn, matmul(a, xn)))
	~~mw.logObject(xn)~~
	end		end

2024-07-22T03:07:36Z

← Older revision		Revision as of 03:07, 22 July 2024
Line 78:		Line 78:
	-- end		-- end

	local xn = scalarmatmul(a, 0.01)		local xn = scalarmatmul(a, 0.0001)

	for i = 1, 20 do		for i = 1, 50 do
	xn = matsub(scalarmatmul(xn, 2), matmul(xn, matmul(a, xn)))		xn = matsub(scalarmatmul(xn, 2), matmul(xn, matmul(a, xn)))
	mw.logObject(xn)		mw.logObject(xn)

2024-07-22T03:07:14Z

2024-07-22T03:06:42Z

← Older revision		Revision as of 03:06, 22 July 2024
Line 65:		Line 65:
	end		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		for i = 1, 20 do
	xn = matsub(scalarmatmul(xn, 2), matmul(xn, matmul(a, xn)))		xn = matsub(scalarmatmul(xn, 2), matmul(xn, matmul(a, xn)))
			mw.logObject(xn)
	end		end

2024-07-22T03:04:53Z

← Older revision		Revision as of 03:04, 22 July 2024
Line 78:		Line 78:
	end		end

	local xn = scalarmatmul(a, ~~0.0005~~)		local xn = scalarmatmul(a, 2)

	for i = 1, 10 do		-- for i = 1, 10 do
	xn = matmul(xn, matsub(dbl_identity, matmul(a, xn)))		-- xn = matmul(xn, matsub(dbl_identity, matmul(a, xn)))
	end		-- end

	for i = 1, 20 do		for i = 1, 20 do

2024-07-22T03:04:31Z

← Older revision		Revision as of 03:04, 22 July 2024
Line 78:		Line 78:
	end		end

	local xn = scalarmatmul(a, 0.~~00005~~)		local xn = scalarmatmul(a, 0.0005)

	for i = 1, ~~100~~ do		for i = 1, 10 do
	xn = matmul(xn, matsub(dbl_identity, matmul(a, xn)))		xn = matmul(xn, matsub(dbl_identity, matmul(a, xn)))
	end		end

			for i = 1, 20 do
			xn = matsub(scalarmatmul(xn, 2), matmul(xn, matmul(a, xn)))
			end

	return xn		return xn
	end		end