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

remove old code based on precomputed reps

Browse Source
This commit is contained in:
Marek Kaluba 2021-09-20 14:18:46 +02:00
parent 048a35800f
commit 051290ed07
No known key found for this signature in database
GPG Key ID: 8BF1A3855328FC15
1 changed files with 0 additions and 102 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

102
precomputed_reps_spectral_gap.jl
View File

@ -1,102 +0,0 @@
using Nemo
using DelimitedFiles
using LinearAlgebra
include("src/nemo_utils.jl")
const PRECISION = 256
function parse_eval(expr_str, value, var_name)
ex = Meta.parse(expr_str)
svar = :($var_name)
return @eval begin
let $svar = $value
$ex
end
end
end
function read_eval(fname, var_name, value)
a = readdlm(fname, ',', String)
a .= replace.(a, '/' => "//")
return parse_eval.(a, value, var_name)
end
function load_discrete_repr(i, q = 109; CC = AcbField(PRECISION))
ζ = root_of_unity(CC, (q + 1) ÷ 2)
degree = q - 1
ra = read_eval(
"data/Discrete reps PSL(2, $q)/discrete_rep_$(i)_a.txt",
:Z,
ζ,
)
a = matrix(CC, [CC(s) for s in ra[1:degree, 1:degree]])
rb = read_eval(
"data/Discrete reps PSL(2, $q)/discrete_rep_$(i)_b.txt",
:Z,
ζ,
)
b = matrix(CC, [CC(s) for s in rb[1:degree, 1:degree]])
@assert contains(det(a), 1)
@assert contains(det(b), 1)
return a, b
end
function load_principal_repr(i, q = 109; CC = AcbField(PRECISION))
ζ = root_of_unity(CC, (q - 1) ÷ 2)
degree = q + 1
ra = read_eval(
"data/Principal reps PSL(2, $q)/principal_rep_$(i)_a.txt",
:zz,
ζ,
)
a = matrix(CC, [CC(z) for z in ra[1:degree, 1:degree]])
rb = read_eval(
"data/Principal reps PSL(2, $q)/principal_rep_$(i)_b.txt",
:zz,
ζ,
)
b = matrix(CC, [CC(z) for z in rb[1:degree, 1:degree]])
@assert contains(det(a), 1)
@assert contains(det(b), 1)
return a, b
end
if !isinteractive()
for i = 0:27
try
a, b = load_principal_repr(i)
adjacency = sum(a^i for i = 1:4) + sum(b^i for i = 1:4)
@time ev = let evs = safe_eigvals(adjacency)
_count_multiplicites(evs)
end
@info "Principal Series Representation $i" ev[1:2] ev[end]
catch ex
@error "Principal Series Representation $i failed" ex
ex isa InterruptException && throw(ex)
end
end
for i = 1:27
try
a, b = load_discrete_repr(i)
adjacency = sum(a^i for i = 1:4) + sum(b^i for i = 1:4)
@time ev = let evs = safe_eigvals(adjacency)
_count_multiplicites(evs)
end
@info "Discrete Series Representation $i" ev[1:2] ev[end]
catch ex
@error "Discrete Series Representation $i : failed" ex
ex isa InterruptException && rethrow(ex)
end
end
end