From 167f57c15bca276db4af347440837d8662d913df Mon Sep 17 00:00:00 2001
From: hlane33 <119619760+hlane33@users.noreply.github.com>
Date: Wed, 13 Dec 2023 10:01:17 -0500
Subject: [PATCH] Extend Lanczos buffer.

---
 src/SpinWaveTheory/KPM.jl | 5 +----
 1 file changed, 1 insertion(+), 4 deletions(-)

diff --git a/src/SpinWaveTheory/KPM.jl b/src/SpinWaveTheory/KPM.jl
index f599dfd61..f804a746e 100644
--- a/src/SpinWaveTheory/KPM.jl
+++ b/src/SpinWaveTheory/KPM.jl
@@ -272,9 +272,6 @@ function intensity_formula_kpm(f,swt::SpinWaveTheory,corr_ix::AbstractVector{Int
     sqrt_Nm_inv = 1.0 / √Nm #define prefactors
     S = (Ns-1) / 2
     sqrt_halfS  = √(S/2) #define prefactors   
-    sqrt_Nm_inv = 1.0 / √Nm #define prefactors
-    S = (Ns-1) / 2
-    sqrt_halfS  = √(S/2) #define prefactors   
     Ĩ = spdiagm([ones(nmodes); -ones(nmodes)]) 
     n_iters = 50
     Avec_pref = zeros(ComplexF64, Nm) # initialize array of some prefactors   
@@ -285,7 +282,7 @@ function intensity_formula_kpm(f,swt::SpinWaveTheory,corr_ix::AbstractVector{Int
         q_absolute = swt.sys.crystal.recipvecs * q_reshaped
         u = zeros(ComplexF64,3,2*nmodes)
         multiply_by_hamiltonian!(y, x) = (swt.sys.mode == :SUN) ? multiply_by_hamiltonian_SUN!(y, x, swt, q_reshaped) : multiply_by_hamiltonian_dipole!(y, x, swt, q_reshaped)
-        lo,hi = Sunny.eigbounds_MF(swt,q_reshaped,n_iters; extend=0.25) # calculate bounds
+        lo,hi = Sunny.eigbounds_MF(swt,q_reshaped,n_iters; extend=0.5) # calculate bounds
         γ=max(abs(lo),abs(hi)) # select upper bound (combine with the preceeding line later)
         # u(q) calculation
         for site = 1:Nm