KPM fix scaling of Bose factor and regularization.

src/SpinWaveTheory/KPM.jl

@@ -45,7 +45,7 @@ interval [-1,1]. Regularization is treated using a cubic cutoff function and the
 
 """
 function get_all_coefficients(M, ωs, broadening, σ, kT,γ;η=0.05)
-    f(ω, x) = regularization_function(x,η*σ) * broadening(ω, x*γ, σ) * (1 + bose_function(kT, x))
+    f(ω, x) = regularization_function(x*γ,η*σ) * broadening(ω, x*γ, σ) * (1 + bose_function(kT, x*γ))
     output = OffsetArray(zeros(M, length(ωs)), 0:M-1, 1:length(ωs))
     for i in eachindex(ωs)
         output[:, i] = cheb_coefs(M, 2M, x -> f(ωs[i], x), (-1, 1))