Module:Temperament data/Badness testing: Difference between revisions

CompactStar (talk | contribs)
No edit summary
CompactStar (talk | contribs)
No edit summary
 
(5 intermediate revisions by the same user not shown)
Line 68: Line 68:


local function matinv(a)
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
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
Line 346: Line 335:
e = e/#(te_tuning[1])
e = e/#(te_tuning[1])
local usb = c * e
local usb = c * e
result = result .. "USB: " .. u._round(usb, 7)
result = u._round(usb, 7)
return result
return result
end
end


return p
return p