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

better parsing

Browse Source
This commit is contained in:
kalmarek 2020-05-13 15:18:46 +02:00
parent ea407d2e98
commit 9698901921
No known key found for this signature in database
GPG Key ID: 8BF1A3855328FC15
1 changed files with 40 additions and 29 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

47
PSL.jl
View File

@ -3,19 +3,12 @@ using DelimitedFiles
include("src/nemo_utils.jl") include("src/nemo_utils.jl")
function parse_evalZ(arg, expr_str)
ex = Meta.parse(expr_str)
return @eval begin
let Z=$arg
$ex
end
end
end
function parse_evalzz(arg, expr_str) function parse_eval(arg, expr_str, var)
ex = Meta.parse(expr_str) ex = Meta.parse(expr_str)
svar = :($var)
return @eval begin return @eval begin
let zz=$arg let $svar = $arg
$ex $ex
end end
end end
@ -25,11 +18,19 @@ 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(
a = matrix(CC, [CC(parse_evalZ(ζ, s)) for s in ra[1:degree, 1:degree]]) "data/Discrete reps PSL(2, $q)/discrete_rep_$(i)_a.txt",
',',
String,
)
a = matrix(CC, [CC(parse_eval(ζ, s, :Z)) for s in ra[1:degree, 1:degree]])
rb = readdlm("data/Discrete reps PSL(2, $q)/discrete_rep_$(i)_b.txt", ',', String) rb = readdlm(
b = matrix(CC, [CC(parse_evalZ(ζ, s)) for s in rb[1:degree, 1:degree]]) "data/Discrete reps PSL(2, $q)/discrete_rep_$(i)_b.txt",
',',
String,
)
b = matrix(CC, [CC(parse_eval(ζ, s, :Z)) for s in rb[1:degree, 1:degree]])
return a, b return a, b
end end
@ -38,15 +39,25 @@ 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(
a = matrix(CC, [CC(parse_evalzz(ζ, s)) for s in ra[1:degree, 1:degree]]) "data/Principal reps PSL(2, $q)/principal_rep_$(i)_a.txt",
',',
String,
)
a = matrix(CC, [CC(parse_eval(ζ, s, :zz)) for s in ra[1:degree, 1:degree]])
rb = readdlm("data/Principal reps PSL(2, $q)/principal_rep_$(i)_b.txt", ',', String) rb = readdlm(
b = matrix(CC, [CC(parse_evalzz(ζ, s)) for s in rb[1:degree, 1:degree]]) "data/Principal reps PSL(2, $q)/principal_rep_$(i)_b.txt",
',',
String,
)
b = matrix(CC, [CC(parse_eval(ζ, s, :zz)) for s in rb[1:degree, 1:degree]])
return a, b return a, b
end end
function safe_eigvals(m::acb_mat) function safe_eigvals(m::acb_mat)
CC = base_ring(m) CC = base_ring(m)
X = matrix(CC, rand(CC, size(m))) X = matrix(CC, rand(CC, size(m)))