diff --git a/.github/workflows/CI.yaml b/.github/workflows/CI.yaml
index 9ddce1b39..fb472336a 100644
--- a/.github/workflows/CI.yaml
+++ b/.github/workflows/CI.yaml
@@ -18,6 +18,10 @@ jobs:
       fail-fast: false
       matrix:
         version:
+          - '1.10'
+          - '1.9'
+          - '1.8'
+          - '1.7'
           - '1.6'
           - '1'
         os:
diff --git a/Project.toml b/Project.toml
index d7f9cd4a0..c589fdff1 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,7 +1,7 @@
 name = "GXBeam"
 uuid = "974624c9-1acb-4ad6-a627-8ac40fc27a3e"
 authors = ["Taylor McDonnell <taylor.golden.mcdonnell@gmail.com> and Andrew Ning"]
-version = "0.5.1"
+version = "0.5.2"
 
 [deps]
 ArnoldiMethod = "ec485272-7323-5ecc-a04f-4719b315124d"
@@ -26,21 +26,22 @@ UnPack = "3a884ed6-31ef-47d7-9d2a-63182c4928ed"
 WriteVTK = "64499a7a-5c06-52f2-abe2-ccb03c286192"
 
 [compat]
-ArnoldiMethod = "0.0.4, 0.2"
+julia = "1.6"
+ArnoldiMethod = "0.0.4 - 0.5"
 ChainRulesCore = "1"
-FLOWMath = "0.3"
-FillArrays = "0.13, 1"
-ForwardDiff = "0.10"
-ImplicitAD = "0.3.1"
-IterativeSolvers = "0.9"
-LinearMaps = "2, 3"
+FLOWMath = "0.3 - 0.5"
+FillArrays = "0.13 - 2"
+ForwardDiff = "0.10 - 1"
+ImplicitAD = "0.3.1 - 0.4"
+IterativeSolvers = "0.9 - 0.10"
+LinearMaps = "2 - 4"
 NLsolve = "4"
 OffsetArrays = "1"
 ReverseDiff = "1"
-Roots = "1, 2"
+Roots = "1 - 3"
 SciMLBase = "1"
-SparseDiffTools = "1, 2"
-StaticArrays = "0.12, 1.0"
+SparseArrays = "1.6 - 2"
+SparseDiffTools = "2"
+StaticArrays = "1.9 - 2"
 UnPack = "1"
-WriteVTK = "1"
-julia = "1.6"
+WriteVTK = "1.21 - 2"
diff --git a/src/analyses.jl b/src/analyses.jl
index c99d0bb05..7e1e3c662 100644
--- a/src/analyses.jl
+++ b/src/analyses.jl
@@ -952,7 +952,7 @@ function solve_eigensystem(x, K, M, nev)
     A = LinearMap{T}(f!, fc!, nx, nx; ismutating=true)
 
     # compute eigenvalues and eigenvectors
-    λ, V = partialeigen(partialschur(A; nev=min(nx, nev), which=LM(), tol=1e-9)[1])
+    λ, V = partialeigen(partialschur(A; nev=min(nx, nev), which=ArnoldiMethod.LM(), tol=1e-9)[1])
 
     # sort eigenvalues by magnitude
     perm = sortperm(λ, by=(λ) -> (abs(λ), imag(λ)), rev=true)