Xenharmonic Wiki

Module:Temperament data: Difference between revisions

← Older edit
CompactStar (talk | contribs)
5,020 edits
No edit summary
ArrowHead294 (talk | contribs)
14,512 edits
mNo edit summary
 
(6 intermediate revisions by one other user not shown)
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
Line 65: Line 66:
end
end


function p.matinv(a)
local function matinv(a)
-- local dbl_identity = {}
local xn = scalarmatmul(a, 1e-7)
-- 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.01)


for i = 1, 20 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)))
mw.logObject(xn)
end
end
Retrieved from "https://en.xen.wiki/w/Module:Temperament_data"