fix

Author: wsmoses

src/periodograms.jl

@@ -190,11 +190,11 @@ function fft2pow2radial!(out::Array{T}, s_fft::Matrix{Complex{T}}, n1::Int, n2::
     kmax = length(out)   # the highest wavenumber
     wc = zeros(Int, kmax) # wave count for radial average
     if n1 == nmin        # scale the wavevector for non-square s_fft
-        c1 = 1
         c2 = div(n1, n2)
+        c1 = one(c2)
     else
         c1 = div(n2, n1)
-        c2 = 1
+        c2 = one(c1)
     end
 
     sqrt_muladd(a, k) = sqrt(muladd(a, a, k))