Skip to content

Instantly share code, notes, and snippets.

@vznvzn
Created January 9, 2021 05:28
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save vznvzn/f421a66aecc8890455ed055c264ed761 to your computer and use it in GitHub Desktop.
Save vznvzn/f421a66aecc8890455ed055c264ed761 to your computer and use it in GitHub Desktop.
def f2(n)
return n.odd? ? (n * 3 + 1) / 2 : n / 2
end
def adj(x, m, n, s2, p2, p, x3)
s2.replace(p[0...x] + (x3 % 2).to_s)
m = (s2[x, 1] == p[x, 1])
n = p2.reverse.to_i(2)
p2[x, 1] = m ? '11' : '01'
return m, n, x3
end
def terras121(p)
p2 = ['01', '11'][p[0, 1].to_i]
n = 1
n3 = 0
m = p2 == '11'
x3 = (p2.reverse.to_i(2)) >> 1
(1...p.length).each \
{
|x|
if (!m)
ns = n.to_s(2)
ns[0, 1] = ''
m1 = x3 - 3**n3
n1 = f2(m1)
end
n3 = p[0...x].split('').select { |z| z == '1' }.size
s2 = ''
m, n, x3 = adj(x, m, n, s2, p2, p, m ? (f2(x3) + 3**n3) : (n1 + 3**n3))
}
n = p2.reverse.to_i(2)
return n
end
def seq(n)
l = [n]
while (n != 1)
n = f2(n)
l << n
end
return l
end
def outd(f, l)
f.puts('$dat << eof')
k = l[0].keys
f.puts(k.join("\t"))
l.each { |x| f.puts(x.values.join("\t")) }
f.puts('eof')
return k
end
def outa(f, l, a = {}, t = '')
k = outd(f, l)
f.puts("set colors classic; set title '#{t}'; ")
f.puts("set ytics nomirror; set y2tics;")
f.puts("plot \\")
k, ct = [k - ['t'], "(column('t')):"] if (k.member?('t'))
k.each \
{
|x|
next if (a.member?(x) && a[x].nil?)
opt = a.fetch(x, '')
opt += ' lw 2 ' if (!opt.include?('lw'))
opt += ' with line ' if (!opt.include?('with') && !opt.include?('pt'))
f.puts("'$dat' using #{ct}(column('#{x}')) #{opt} title '#{x}',\\")
}
f.puts
# f.puts("reset; pause -1;")
end
def outafn(l, a = {}, fno = nil, t = '')
fn = "gnuplot#{fno}.cmd"
outa(f = File.open(fn, 'w'), l, a, t)
f.close
$stderr.puts([fn, t, l.size].inspect)
end
def dense(w, d)
w2 = w - 1
a = (0...w2).to_a
s = '0' * w2
(1..(d * w - 1)).map { a.delete_at(rand(a.size)) }.each { |x| s[x, 1] = '1' }
return ('1' + s)
end
def d(s)
c = s.split('').select { |x| x == '1' }.size
d = c.to_f / s.length
return d
end
def add(a, x)
a[x] = a.fetch(x, 0) + 1
end
def val(a)
a.fetch(1, 1).to_f / a.fetch(2, 1)
end
def add1(l, w, a, k)
x = l[w]
x = 0 if (x.nil?)
add(a, x % 3)
Hash[[["#{k}0", "#{k}1", "#{k}2"], a.sort.transpose[1]].transpose]
end
c = 500
a = {0 => 0, 1 => 0, 2 => 0}
a2 = a.dup
w = 100
l = []
c.times \
{
|i|
p = dense(w, i.to_f / (c - 1))
n = terras121(p)
l1 = seq(n)
l2 = seq(p.to_i(2))
d1 = add1(l1, w, a, 'x')
d2 = add1(l2, w, a2, 'y')
l << d1.merge(d2)
}
outafn(l)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment