reduce allocs

Author: baggepinnen

src/AlphaStableDistributions.jl

@@ -277,7 +277,7 @@ function Distributions.fit(d::Type{<:AlphaSubGaussian}, x::AbstractVector{T}, m:
     α    = d1.α; scale=d1.scale
     cov  = zeros(T, m+1, m+1)
     xlen = length(x)
-    c    = ((sum(abs.(x).^p)/xlen)^(1/p))/scale
+    c    = ((sum(x->abs(x)^p, x)/xlen)^(1/p))/scale
     for i in 1:m
         tempxlen = xlen-mod(xlen, i)
         xtemp = reshape(x[1:end-mod(xlen, i)], i, tempxlen÷i)
@@ -286,7 +286,7 @@ function Distributions.fit(d::Type{<:AlphaSubGaussian}, x::AbstractVector{T}, m:
             tempxlen = size(xtemp, 1)*size(xtemp, 2)
         end
         xtemp = reshape(xtemp', 2, tempxlen÷2)
-        r = (2/(c^p))*(scale^(2-p))*(xtemp[1, :]'*((sign.(xtemp[2, :]).*(abs.(xtemp[2, :]).^(p-1)))))/(tempxlen/2)
+        @views r = (2/(c^p))*(scale^(2-p))*(xtemp[1, :]'*((sign.(xtemp[2, :]).*(abs.(xtemp[2, :]).^(p-1)))))/(tempxlen/2)
         cov[diagind(cov, i)] .+= r
     end
     cov = (cov+cov')+2*(scale^2)*I(m+1)