kalmar/SmallHyperbolic
1
0
Fork 0
mirror of https://github.com/kalmarek/SmallHyperbolic synced 2024-11-30 17:15:28 +01:00
Code Issues Releases Wiki Activity

small tweaks to the script

Browse Source
This commit is contained in:
kalmarek 2020-05-13 01:45:19 +02:00
parent ed6a987f70
commit a23f0ba464
No known key found for this signature in database
GPG Key ID: 8BF1A3855328FC15
1 changed files with 27 additions and 22 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

49
PSL.jl
View File

@ -21,8 +21,8 @@ function parse_evalzz(arg, expr_str)
end end
end end
function load_discrete_repr(i, q=109; CC=AcbField(128)) function load_discrete_repr(i, q=109; CC=AcbField(512))
ζ = root_of_unity(CC, (q+1)÷2) ζ = root_of_unity(CC, (q-1)÷2)
degree = q-1 degree = q-1
ra = readdlm("data/Discrete reps PSL(2, $q)/discrete_rep_$(i)_a.txt", ',', String) ra = readdlm("data/Discrete reps PSL(2, $q)/discrete_rep_$(i)_a.txt", ',', String)
@ -35,8 +35,8 @@ function load_discrete_repr(i, q=109; CC=AcbField(128))
end end
function load_principal_repr(i, q=109; CC=AcbField(128)) function load_principal_repr(i, q=109; CC=AcbField(512))
ζ = root_of_unity(CC, (q+1)÷2) ζ = root_of_unity(CC, (q-1)÷2)
degree = q+1 degree = q+1
ra = readdlm("data/Principal reps PSL(2, $q)/principal_rep_$(i)_a.txt", ',', String) ra = readdlm("data/Principal reps PSL(2, $q)/principal_rep_$(i)_a.txt", ',', String)
@ -48,29 +48,34 @@ function load_principal_repr(i, q=109; CC=AcbField(128))
return a,b return a,b
end end
# for i in 0:27 for i in 0:27
# try try
# a,b = load_principal_repr(i) a,b = load_principal_repr(i)
# adjacency = sum([[a^i for i in 1:4]; [b^i for i in 1:4]]) adjacency = sum([[a^i for i in 1:4]; [b^i for i in 1:4]])
# M = parent(adjacency) # @time evc = eigvals(adjacency)
#
# # X = M(rand(base_ring, size(adjacency))) CC = base_ring(adjacency)
# X = matrix(CC, rand(CC, size(adjacency)))
# # @time ev = eigvals(X*adjacency*inv(X)) @time evc = eigvals(X*adjacency*inv(X))
# @time evc = eigvals(adjacency)
# ev = sort(real.(first.(evc)), lt=<, rev=true) ev = sort(real.(first.(evc)), lt=<, rev=true)
# @info "Principal Series Representation $i" ev[1:4] @info "Principal Series Representation $i" ev[1:4]
# catch ex catch ex
# @error "Principal Series Representation $i : failed" @error "Principal Series Representation $i failed"
# ex isa InterruptException && throw(ex) ex isa InterruptException && throw(ex)
# end end
# end end
for i in 1:27 for i in 1:27
try try
a,b = load_discrete_repr(i) a,b = load_discrete_repr(i)
adjacency = sum([[a^i for i in 1:4]; [b^i for i in 1:4]]) adjacency = sum([[a^i for i in 1:4]; [b^i for i in 1:4]])
@time evc = eigvals(adjacency) # @time evc = eigvals(adjacency)
CC = base_ring(adjacency)
X = matrix(CC, rand(CC, size(adjacency)))
@time evc = eigvals(X*adjacency*inv(X))
ev = sort(real.(first.(evc)), lt=<, rev=true) ev = sort(real.(first.(evc)), lt=<, rev=true)
@info "Discrete Series Representation $i" ev[1:4] @info "Discrete Series Representation $i" ev[1:4]
catch ex catch ex