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 19: | Line 19: | ||
end | end | ||
local function matadd(a, b) | |||
local result = {} | |||
for i = 1, #a do | |||
result[i] = {} | |||
for j = 1, #(b[1]) do | |||
result[i][j] = cadd(a[i][j], b[i][j]) | |||
end | |||
end | |||
return result | |||
end | |||
local function matmul(a, b) | local function matmul(a, b) |
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 matadd(a, b)
local result = {}
for i = 1, #a do
result[i] = {}
for j = 1, #(b[1]) do
result[i][j] = cadd(a[i][j], b[i][j])
end
end
return result
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