Error message for Blume-Maleev with single-q

SunnySuite · Sep 4, 2024 · 0c91dbf · 0c91dbf
1 parent 8acce27
commit 0c91dbf
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/src/Measurements/MeasureSpec.jl b/src/Measurements/MeasureSpec.jl
@@ -190,6 +190,9 @@ function ssf_custom_bm(f, sys::System; u, v, apply_g=true, formfactors=nothing)
     e3 = normalize(u × v)
 
     return ssf_custom(sys::System; apply_g, formfactors) do q, ssf
+        if iszero(q)
+            error("Blume-Maleev axis system not defined at zero q")
+        end
         if abs(q ⋅ e3) > 1e-12
             error("Momentum transfer q not in scattering plane")
         end