Module:Temperament data: Difference between revisions
Jump to navigation
Jump to search
CompactStar (talk | contribs) No edit summary |
CompactStar (talk | contribs) No edit summary |
||
Line 15: | Line 15: | ||
end | end | ||
function | local function cdiv(a, b) | ||
return {(a[1] * b[1] + a[2] * b[2])/(b[1]*b[1] + b[2]*b[2]), (a[2] * b[1] - a[1] * b[2])/(b[1]*b[1] + b[2]*b[2])} | return {(a[1] * b[1] + a[2] * b[2])/(b[1]*b[1] + b[2]*b[2]), (a[2] * b[1] - a[1] * b[2])/(b[1]*b[1] + b[2]*b[2])} | ||
end | end | ||
Line 25: | Line 25: | ||
result[i] = {} | result[i] = {} | ||
for j = 1, #(b[1]) do | for j = 1, #(b[1]) do | ||
result[i][j] = 0 | result[i][j] = {0, 0} | ||
for k = 1, #(a[1]) do | for k = 1, #(a[1]) do | ||
result[i][j] = result[i][j] | result[i][j] = cadd(result[i][j], cmul(a[i][k], b[k][j])) | ||
end | end | ||
end | end |
Revision as of 03:09, 14 October 2023
Note: Do not invoke this module directly; use the corresponding template instead: Template:Temperament data.
local rat = require('Module:Rational')
local p = {}
-- Complex number functions (1st element = real and 2nd element = imaginary)
local function cadd(a, b)
return {(a[1] + b[1]), (a[2] + b[2])}
end
local function csub(a, b)
return {(a[1] - b[1]), (a[2] - b[2])}
end
local function cmul(a, b)
return {(a[1] * b[1] - a[2] * b[2]), (a[1] * b[2] + a[2] * b[1])}
end
local function cdiv(a, b)
return {(a[1] * b[1] + a[2] * b[2])/(b[1]*b[1] + b[2]*b[2]), (a[2] * b[1] - a[1] * b[2])/(b[1]*b[1] + b[2]*b[2])}
end
local function matmul(a, b)
local result = {}
for i = 1, #a do
result[i] = {}
for j = 1, #(b[1]) do
result[i][j] = {0, 0}
for k = 1, #(a[1]) do
result[i][j] = cadd(result[i][j], cmul(a[i][k], b[k][j]))
end
end
end
return result
end
local function matinv(a)
end
return p