Update deps

Project.toml

@@ -37,13 +37,13 @@ WGLMakiePrecompilesExt = "WGLMakie"
 
 [compat]
 CrystalInfoFramework = "0.6.0"
-DynamicPolynomials = "0.5.2"
+DynamicPolynomials = "0.6.0"
 ElasticArrays = "1.2.12"
 FFTW = "1.4.5"
 FilePathsBase = "0.9.18"
 GLMakie = "0.9, 0.10"
 HCubature = "1.7.0"
-JLD2 = "0.4.33"
+JLD2 = "0.5.4"
 LinearAlgebra = "1.10"
 Makie = "0.20, 0.21.1"
 OffsetArrays = "1.10.8"

src/Operators/Symbolic.jl

@@ -23,7 +23,7 @@ function stevens_as_spin_polynomials(k::Int)
     return map(𝒪s) do 𝒪
         # In the large-s limit only leading order terms contribute, yielding a
         # homogeneous polynomial of degree k
-        𝒪 = sum(t for t in 𝒪 if DP.degree(t) == k)
+        𝒪 = DP.filter_terms(DP.OfDegree(k), 𝒪)
         # Remaining coefficients must be real integers; make this explicit
         𝒪 = DP.map_coefficients(x -> Int(x), 𝒪)
         return 𝒪
@@ -35,7 +35,7 @@ end
 # expectation values
 function expand_as_spin_polynomial(p)
     𝒮 = spin_vector_symbol
-    return DP.subs(p, 
+    return DP.subs(p,
         spin_squared_symbol => 𝒮⋅𝒮,
         [stevens_symbols[k] => stevens_as_spin_polynomials(k) for k=0:6]...
     )