Xenharmonic Wiki

Module:Interval table: Difference between revisions

← Older edit
CompactStar (talk | contribs)
autopatrolled
5,084 edits
This is the final adjustment probably
ArrowHead294 (talk | contribs)
15,213 edits
mNo edit summary
 
(88 intermediate revisions by 6 users not shown)
Line 1: Line 1:
local p = {}
local p = {}
local u = require('Module:Utils')
local iv = require('Module:Interval')
local rat = require('Module:Rational')
local ud = require('Module:Ups and downs notation')
local ET = require('Module:ET')


-- Note: some of these helper functions may be moved to utils module
local ET = require("Module:ET")
local rat = require("Module:Rational")
local ud = require("Module:Ups and downs notation")
local utils = require("Module:Utils")
local yesno = require("Module:Yesno")


-- Get prime limit of n/d
-- Gets mapping of a monzo in a val
local function get_limit(n,d)
local result = 2
local monzo = u.get_monzo(n, d)
for k,_ in pairs(monzo) do
if k > result then
result = k
end
end
return result
end
 
--  Gets mapping of a monzo in a val
local function mapping(monzo, val)
local function mapping(monzo, val)
local result = 0
local result = 0
for k,v in pairs(val) do
for k, _ in pairs(val) do
result = result + val[k] * (monzo[k] or 0)
result = result + val[k] * (monzo[k] or 0)
end
end
Line 35: Line 22:
for i = 1, max_nd do
for i = 1, max_nd do
for j = 1, max_nd do
for j = 1, max_nd do
t = i/u._gcd(i, j) .. "/" .. j/u._gcd(i, j)
local t = i / utils._gcd(i, j) .. "/" .. j / utils._gcd(i, j)
if (i/j) >= 1 and (i/j) <= max_size and get_limit(i, j) <= prime_limit and not u.table_contains(ratio_strings, t) then
if
ratios[#ratios + 1] = {i/u._gcd(i, j),j/u._gcd(i, j)}
(i / j) >= 1
and (i / j) <= max_size
and rat.max_prime(rat.new (i, j)) <= prime_limit
and not utils.table_contains(ratio_strings, t)
then
ratios[#ratios + 1] = { i / utils._gcd(i, j), j / utils._gcd(i, j) }
ratio_strings[#ratio_strings + 1] = t
ratio_strings[#ratio_strings + 1] = t
end
end
end
end
end
end
 
return ratios
return ratios
end
end


-- Utility fuunction to get specific note name with ud.get_note_names_table
-- Utility function to get specific note name with ud.get_note_names_table
-- (this is essentially what "Template:Ups and downs note name" does)
-- (this is essentially what "Template:Ups and downs note name" does)
local function ud_note(et, fifth, step)
local function ud_note(et, fifth, step)
return table.concat(ud.get_note_names_table(et, fifth)[step], ", "):sub(0, -1)
return table.concat(ud.get_note_names_table(et, fifth)[step], ", "):sub(0, -1)
end
end


local function mysplit (inputstr, sep)
-- ???
        if sep == nil then
local function mysplit(inputstr, sep)
                sep = "%s"
if sep == nil then
        end
sep = "%s"
        local t={}
end
        for str in string.gmatch(inputstr, "([^"..sep.."]+)") do
local t = {}
                table.insert(t, str)
for str in string.gmatch(inputstr, "([^" .. sep .. "]+)") do
        end
table.insert(t, str)
        return t
end
return t
end
end


function p.interval_table(frame)
function p.interval_table(frame)
local additional = frame.args['additional'] or ''
local additional = frame.args["additional"] or ""
local additional_split = mysplit(additional, '\n')
local additional_split = mysplit(additional, "\n")
for i=1,#additional_split do
for i = 1, #additional_split do
additional_split[i] = mysplit(additional_split[i], ';')
additional_split[i] = mysplit(additional_split[i], ";")
end
end
local tuning = frame.args['tuning']
local tuning = frame.args["tuning"]
local et = ET.parse(tuning) or ET.parse('12edo')
local et = ET.parse(tuning) or ET.parse("12edo")
local debugg = yesno(frame.args["debug"])
local patent_val = {
local patent_val = {
[2] = ET.approximate(et, 2),  
[2] = ET.approximate(et, 2),
[3] = ET.approximate(et, 3),  
[3] = ET.approximate(et, 3),
[5] = ET.approximate(et, 5),  
[5] = ET.approximate(et, 5),
[7] = ET.approximate(et, 7),  
[7] = ET.approximate(et, 7),
[11] = ET.approximate(et, 11),  
[11] = ET.approximate(et, 11),
[13] = ET.approximate(et, 13),  
[13] = ET.approximate(et, 13),
[17] = ET.approximate(et, 17),
[17] = ET.approximate(et, 17),
[19] = ET.approximate(et, 19),  
[19] = ET.approximate(et, 19),
[23] = ET.approximate(et, 23),
[23] = ET.approximate(et, 23),
[29] = ET.approximate(et, 29),
[29] = ET.approximate(et, 29),
Line 86: Line 80:
[41] = ET.approximate(et, 41),
[41] = ET.approximate(et, 41),
[43] = ET.approximate(et, 43),
[43] = ET.approximate(et, 43),
[47] = ET.approximate(et, 47)
[47] = ET.approximate(et, 47),
}
} -- NOTE: indices are prime numbers
 
local t_head = "{| class=\"wikitable center-1 right-2"
if et.size > 20 then
t_head = t_head .. " mw-collapsible mw-collapsed"
end
t_head = t_head .. "\"\n"
local fifth = patent_val[3] - patent_val[2]
local fifth_error = ET.cents(et, fifth) - rat.cents(rat.new(3, 2))
local dual_fifth = (math.abs(fifth_error) > (400 / et.size))
local dual_flat_fifth = ET.approximate(et, 3 / 2, -1)
local dual_sharp_fifth = ET.approximate(et, 3 / 2, 1)
 
-- List of all 47-limit ratios with numerator and denominator <= certain cents value and in the range 1/1 - 12/1
local ratios_list = get_ratios_list((20 + (et.size / 4)), 12, 47)
 
result = t_head ..
"|-\n" ..
"! Steps\n" ..
"! Cents\n"
local wikitext = '{|class="wikitable"\n'
if rat.eq(et.equave, 3) then
local octave = patent_val[2]
result = result .. "! [[Hekt]]s\n"
local fifth = patent_val[3] - patent_val[2]
end
local fifth_error = ET.cents(et, fifth) - iv._to_cents(3/2)
result = result ..
local dual_fifth = math.abs(fifth_error) > (400 / et.size)
"! Approximate ratios\n"
local dual_flat_fifth = ET.approximate(et, 3/2, -1)
local dual_sharp_fifth = ET.approximate(et, 3/2, 1)
-- List of all 13-limit ratios with numerator and denominator <= 81 and in the range 1/1 - 5/1
local ratios_list = get_ratios_list(81, 5, 13)
wikitext = wikitext .. '!Steps\n'
wikitext = wikitext .. '!Cents\n'
if rat.eq(et.equave, 2) then
if rat.eq(et.equave, 2) then
if dual_fifth then
if dual_fifth then
wikitext = wikitext .. '![[Ups and downs notation]]<br>(dual flat fifth ' .. dual_flat_fifth  .. '\\' .. et.size .. ')\n'
result = result
wikitext = wikitext .. '![[Ups and downs notation]]<br>(dual sharp fifth ' .. dual_sharp_fifth .. '\\' .. et.size .. ')\n'
.. string.format("! [[Ups and downs notation]]<br />(Dual flat fifth %s\\%s)\n", dual_flat_fifth, et.size)
.. string.format("! [[Ups and downs notation]]<br />(Dual sharp fifth %s\\%s)\n", dual_sharp_fifth, et.size)
else
else
wikitext = wikitext .. '![[Ups and downs notation]]\n'
result = result .. "! [[Ups and downs notation]]\n"
end
end
end
end
wikitext = wikitext .. '!Approximate ratios\n'
if #additional_split > 0 then
if #additional_split > 0 then
for i=1,#additional_split[1] do
for i = 1, #additional_split[1] do
wikitext = wikitext .. '!' .. additional_split[1][i] .. '\n'
result = result .. string.format("! %s\n", additional_split[1][i])
end
end
end
end


for i=0,et.size do
    local decimalPlaces = 3 - math.floor(math.log(1200 / et.size, 10))
wikitext = wikitext .. '|-\n'
    if decimalPlaces > 3 then decimalPlaces = 3 end
wikitext = wikitext .. '|' .. i .. '\n'
    if decimalPlaces < 1 then decimalPlaces = 1 end
wikitext = wikitext .. '|' .. u._round(ET.cents(et, i), 6) .. '\n'
for i = 0, et.size do
result = result .. "|-\n"
.. string.format("| %s\n", i)
.. string.format("| %s\n", utils._round_dec(ET.cents(et, i), decimalPlaces))
if rat.eq(et.equave, 3) then
result = result .. string.format("| %s\n", utils._round_dec(ET.hekts(et, i), decimalPlaces))
end
result = result .. "| "
for j = 1, #ratios_list do
-- Go through all ratios.
-- If they are mapped to the current step in the ET's patent val and < specified edostep error,
-- add to approximate ratios column.
local t = mapping(rat.new(ratios_list[j][1], ratios_list[j][2]), patent_val)
if
t == i
and math.abs(ET.cents(et, i) - utils._log(ratios_list[j][1] / ratios_list[j][2], 2) * 1200)
< (1200 / (math.log(et.size) * (et.size + 5)))
then
result = result .. string.format("[[%s/%s]], ", ratios_list[j][1], ratios_list[j][2])
end
end
-- ???
if result:sub(-2, -1) == "| " then
result = result .. "\n"
else
result = result:sub(0, -3) .. "\n"
end
if rat.eq(et.equave, 2) then
if rat.eq(et.equave, 2) then
if dual_fifth then
if dual_fifth then
wikitext = wikitext .. '|' .. ud_note(et, dual_flat_fifth, i) .. '\n'
result = result .. string.format("| %s\n", ud_note(et, dual_flat_fifth, i))
wikitext = wikitext .. '|' .. ud_note(et, dual_sharp_fifth, i) .. '\n'
result = result .. string.format("| %s\n", ud_note(et, dual_sharp_fifth, i))
else
else
wikitext = wikitext .. '|' .. ud_note(et, fifth, i) .. '\n'
result = result .. string.format("| %s\n", ud_note(et, fifth, i))
end
end
wikitext = wikitext .. '|'
for j=1,#ratios_list do
-- Go through all ratios, if they are mapped to the current step in the ET's patent val with <1/3 relative error, add to approximate ratios column
t = mapping(u.get_monzo(ratios_list[j][1], ratios_list[j][2]), patent_val)
relative_error = math.abs(t - math.log2(ratios_list[j][1]/ratios_list[j][2]))
if t == i and relative_error < 1/3 then
wikitext = wikitext .. '[[' .. ratios_list[j][1] .. '/' .. ratios_list[j][2] .. ']], '
end
end
end
end
if wikitext:sub(-1, -1) == '|' then
wikitext = wikitext .. '\n'
else
wikitext = wikitext:sub(0, -3) .. '\n'
end
if (i + 2) <= #additional_split then
if (i + 2) <= #additional_split then
for j=1,#additional_split[i+2] do
for j = 1, #additional_split[i + 2] do
wikitext = wikitext .. '|' .. additional_split[i+2][j] .. '\n'
result = result .. string.format("| %s\n", additional_split[i + 2][j])
end
end
end
end
end
end


wikitext = wikitext .. '|}'
result = result .. "|}"
if wikitext:len() > 5000 then
if debugg == true then
wikitext = '{| class="wikitable mw-collapsible mw-collapsed"\n' .. wikitext:sub(20)
result = "<syntaxhighlight lang=\"wikitext\">" .. result .. "</syntaxhighlight>"
end
end
 
return wikitext
return frame:preprocess(result)
end
end


return p
return p
Retrieved from "https://en.xen.wiki/w/Module:Interval_table"