phonebook of internal addresses

Author: jeffwack

interactive.jl

@@ -7,7 +7,7 @@ function colorax(ax)
     hidedecorations!(ax)
     hidespines!(ax)
     deactivate_interaction!(ax, :rectanglezoom)
-    #deactivate_interaction!(ax, :dragpan)
+    deactivate_interaction!(ax, :dragpan)
     #deactivate_interaction!(ax, :scrollzoom)
 
     sg = SliderGrid(fig[2,1],

sequences/AngleDoubling.jl

@@ -255,6 +255,25 @@ function denominators(S::Vector{Int})
     return q
 end
 
+function denominator(AIA::AngledInternalAddress,new::Int)
+    K = kneadingsequence(AIA.addr)
+    p = rhoSequence(K)
+
+    r = mod1(new,AIA.addr[end])
+    orb = p(r)
+
+    if AIA.addr[end] in orb
+        qk = Int((new - r)/AIA.addr[end] + 1) #why is this always an integer?
+
+    else
+        qk = Int((new - r)/AIA.addr[end] + 2)
+
+    end
+    
+    return qk
+end
+
+
 function numembeddings(I::Vector{Int})
     D = denominators(I)
     n = 1

sequences/Phonebook.jl

@@ -0,0 +1,69 @@
+using GLMakie
+include("AngleDoubling.jl")
+
+#for each internal address, calculate a screen which shows the future options
+
+function drawwedges!(ax,AIA)
+    N = AIA.addr[end]
+    n = 5
+    for ii in N+1:N+n 
+        k = denominator(AIA,ii)
+        angles = []
+        for jj in 1:k-1
+            if gcd(jj,k) == 1
+                push!(angles,(jj-0.5)*2*pi/(k-1))
+            end
+        end
+        r = (2^ii+2^(ii-1))/2
+        centers = [(r*cos(theta),r*sin(theta)) for theta in angles]
+        pie!(ax,ones(k-1),radius = 2^ii, inner_radius = 2^(ii-1))
+        text!(ax,centers,text = fill("$ii",length(centers)), align = (:center,:center))
+    end
+end
+
+function phonebook!(gl,AIA)
+    label = Label(gl[2,1],repr(AIA[]))
+    rowsize!(gl, 1, Aspect(1, 1.0))
+    ax1 = Axis(gl[1,1],aspect = 1)
+    drawwedges!(ax1,AIA[])
+
+    hidedecorations!(ax1)
+    hidespines!(ax1)
+    deactivate_interaction!(ax1, :rectanglezoom)
+    deactivate_interaction!(ax1, :dragpan)
+
+    register_interaction!(ax1, :select) do event::MouseEvent, ax1
+        if event.type === MouseEventTypes.leftdown
+            ang = atan(event.data[2],event.data[1])
+            rad = sqrt(sum([x*x for x in event.data]))
+            idx = Int(ceil(log2(rad)))
+            if idx > 0
+                k = denominator(AIA[],idx)
+                num = Int(mod1(ceil(ang*(k-1)/(2*pi)),k-1))
+                println("$idx and $num / $k")
+                if gcd(num,k) == 1
+                    newaddr = push!(copy(AIA[].addr),idx)
+                    newangles = push!(copy(AIA[].angles),Rational(num,k))
+                    AIA[] = AngledInternalAddress(newaddr,newangles)
+                    println(AIA[])
+                    empty!(ax1)
+                    drawwedges!(ax1,AIA[])
+                    label.text = repr(AIA[])
+                end
+            end
+        end
+    end
+
+    
+
+end
+
+
+fig = Figure()
+ga = fig[1, 1] = GridLayout()
+
+AIA = Observable(AngledInternalAddress([1],[]))
+
+phonebook!(ga,AIA)
+
+fig