Module:Temperament data/Badness testing - Revision history

CompactStar at 03:10, 22 July 2024

2024-07-22T03:10:26Z

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

	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, 70 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:55Z

2024-07-22T03:08:41Z

2024-07-22T03:08:11Z

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

	local function matinv(a)		local function matinv(a)
	local ~~dbl_identity~~ = {}		local xn = scalarmatmul(a, 0.0001)
	for i = 1, #a do
	~~dbl_identity[i]~~ = {}		for i = 1, 50 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		end

	~~local xn = scalarmatmul(a, 0.000005)~~

	~~for i = 1, 30 do~~
	~~xn = matmul(xn, matsub(dbl_identity, matmul(a, xn)))~~
	~~end~~
	return xn		return xn
	end		end

2024-07-22T02:32:03Z

2024-07-22T02:31:43Z

2024-07-22T02:31:29Z

← Older revision		Revision as of 02:31, 22 July 2024
Line 334:		Line 334:
	end		end
	end		end
	c = c/(#mapping * #mapping[1])		c = c/(#(mapping) * #(mapping[1]))
			mw.logObject(c)
	local te_generator = get_te_generator(subgroup, comma_matrix, generators)		local te_generator = get_te_generator(subgroup, comma_matrix, generators)
	local te_tuning = matmul(te_generator, mapping)		local te_tuning = matmul(te_generator, mapping)
	local e = 0.0		local e = 0.0
	for i = 1, #(te_tuning[1]) do		for i = 1, #(te_tuning[1]) do
			-- mw.logObject(te_tuning)
			-- mw.logObject((te_tuning[1][i] * 1200) - ((math.log(subgroup[i])/math.log(2)) * 1200))
	e = e + math.abs((te_tuning[1][i] * 1200) - ((math.log(subgroup[i])/math.log(2)) * 1200))		e = e + math.abs((te_tuning[1][i] * 1200) - ((math.log(subgroup[i])/math.log(2)) * 1200))
	end		end
	e = e/#(te_tuning[1])		e = e/#(te_tuning[1])
	local usb = c * e		local usb = c * e
	result = result .. "~~\n\nUSB~~: " .. u._round(usb, 7)		result = result .. "USB: " .. u._round(usb, 7)
	return result		return result
	end		end

	return p		return p

2024-07-22T02:25:42Z

← Older revision		Revision as of 02:25, 22 July 2024
Line 339:		Line 339:
	local e = 0.0		local e = 0.0
	for i = 1, #(te_tuning[1]) do		for i = 1, #(te_tuning[1]) do
	e = e + math.abs(te_tuning[1][i])		e = e + math.abs((te_tuning[1][i] * 1200) - ((math.log(subgroup[i])/math.log(2)) * 1200))
	end		end
	e = e/#(te_tuning[1])		e = e/#(te_tuning[1])

2024-07-22T02:24:38Z

← Older revision		Revision as of 02:24, 22 July 2024
Line 343:		Line 343:
	e = e/#(te_tuning[1])		e = e/#(te_tuning[1])
	local usb = c * e		local usb = c * e
	result = result .. "\n\nUSB: " + u._round(usb, 7)		result = result .. "\n\nUSB: " .. u._round(usb, 7)
	return result		return result
	end		end

	return p		return p

2024-07-22T02:24:21Z

← Older revision		Revision as of 02:24, 22 July 2024
Line 341:		Line 341:
	e = e + math.abs(te_tuning[1][i])		e = e + math.abs(te_tuning[1][i])
	end		end
	e = e/(~~#tuning_matrix~~)		e = e/#(te_tuning[1])
	local usb = c * e		local usb = c * e
	result = result .. "\n\nUSB: " + u._round(usb, 7)		result = result .. "\n\nUSB: " + u._round(usb, 7)