diff --git a/stim/optics/scalarmie.h b/stim/optics/scalarmie.h
index 4250ba3..48b6972 100644
--- a/stim/optics/scalarmie.h
+++ b/stim/optics/scalarmie.h
@@ -389,7 +389,7 @@ void cpu_scalar_mie_scatter(stim::complex<T>* E, size_t N, T* x, T* y, T* z, std
 		r = p.len();
 		if(r >= a){
 			for(size_t w = 0; w < W.size(); w++){
-				Ew = W[w].E() * exp(stim::complex<double>(0, W[w].kvec().dot(c)));
+				Ew = W[w].E() * exp(stim::complex<float>(0, W[w].kvec().dot(c)));
 				kr = p.len() * W[w].kmag();											//calculate k*r
 				stim::bessjyv_sph<double>(Nl, kr, vm, j_kr, y_kr, dj_kr, dy_kr);
 				cos_phi = p.norm().dot(W[w].kvec().norm());							//calculate the cosine of the angle from the propagating direction
@@ -649,7 +649,7 @@ void cpu_scalar_mie_internal(stim::complex<T>* E, size_t N, T* x, T* y, T* z, st
 		if(r < a){
 			E[i] = 0;
 			for(size_t w = 0; w < W.size(); w++){
-				Ew = W[w].E() * exp(stim::complex<double>(0, W[w].kvec().dot(c)));
+				Ew = W[w].E() * exp(stim::complex<float>(0, W[w].kvec().dot(c)));
 				knr = (stim::complex<double>)n * p.len() * W[w].kmag();							//calculate k*n*r
 
 				stim::cbessjyva_sph<double>(Nl, knr, vm, j_knr, y_knr, dj_knr, dy_knr);
--
libgit2 0.21.4