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

move safe_eigvals

Browse Source
This commit is contained in:
kalmarek 2020-05-26 00:09:23 +02:00
parent 6cd237bb8e
commit d781aa46a9
No known key found for this signature in database
GPG Key ID: 8BF1A3855328FC15
2 changed files with 11 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
precomputed_reps_spectral_gap.jl
View File

@ -68,16 +68,6 @@ function load_principal_repr(i, q = 109; CC = AcbField(PRECISION))
return a, b return a, b
end end
function safe_eigvals(m::acb_mat)
evs = eigvals(m)
all(isfinite.(evs)) && return evs
CC = base_ring(m)
X = matrix(CC, rand(CC, size(m)))
evs = eigvals(X * m * inv(X))
all(isfinite.(evs)) && return evs
throw(ArgumentError("Could not compute eigenvalues"))
end
if !isinteractive() if !isinteractive()
for i = 0:27 for i = 0:27

11
src/nemo_utils.jl
View File

@ -163,3 +163,14 @@ function _count_multiplicites(evs)
rev = true, rev = true,
) )
end end
function safe_eigvals(m::acb_mat)
evs = eigvals(m)
all(isfinite.(evs)) && return evs
CC = base_ring(m)
X = matrix(CC, rand(CC, size(m)))
evs = eigvals(X * m * inv(X))
return evs
all(isfinite.(evs)) && return evs
throw(ArgumentError("Could not compute eigenvalues"))
end