bug in wpm

broadcasting dictionaries is not allowed

src/juobs_linalg.jl

@@ -126,7 +126,7 @@ function energies(evals::Union{Vector{Vector{uwreal}},Array{Array{uwreal}} }; wp
             ratio = evals[t][i] / evals[t+1][i]
             aux_en[t] = 0.5*log(ratio * ratio)
         end
-        isnothing(wpm) ? uwerr.(aux_en) : uwerr.(aux_en, wpm)
+        isnothing(wpm) ? uwerr.(aux_en) : [uwerr(a, wpm) for a in aux_en]
         eff_en[i] = copy(aux_en)
     end
     return eff_en

src/juobs_obs.jl

@@ -347,7 +347,7 @@ function dec_const_pcvc(corr::Vector{uwreal}, plat::Vector{Int64}, m::uwreal, mu
         else
             uwerr(f, wpm)
             uwerr(R_av, wpm)
-            uwerr.(R, wpm)
+            [uwerr(Raux, wpm) for Raux in R]
         end
     
         v = value(R_av)
@@ -464,7 +464,7 @@ function comp_t0(Y::YData, plat::Vector{Int64}; L::Int64, pl::Bool=false,
             uwerr.(t2E)
         else
             uwerr(t0, wpm)
-            uwerr.(t2E, wpm)
+            [uwerr(t2E_aux, wpm) for t2E_aux in t2E]
         end
 
         v = value.(t2E)
@@ -553,7 +553,7 @@ function comp_t0(Y::Vector{YData}, plat::Vector{Int64}; L::Int64, pl::Bool=false
             uwerr.(t2E)
         else
             uwerr(t0, wpm)
-            uwerr.(t2E, wpm)
+            [uwerr(t2E_aux, wpm) for t2E_aux in t2E]
         end
 
         v = value.(t2E)

src/juobs_tools.jl

@@ -300,7 +300,7 @@ function md_val(a::uwreal, obs::Corr, derm::Vector{Corr})
 end
 
 function plat_av(obs::Vector{uwreal}, plat::Vector{Int64}, wpm::Union{Dict{Int64,Vector{Float64}},Dict{String,Vector{Float64}}, Nothing}=nothing)
-    isnothing(wpm) ? uwerr.(obs) : uwerr.(obs, wpm)
+    isnothing(wpm) ? uwerr.(obs) : [uwerr(obs_aux, wpm) for obs_aux in obs]
     w = 1 ./ err.(obs)[plat[1]:plat[2]].^2
     av = sum(w .* obs[plat[1]:plat[2]]) / sum(w)
     return av 
@@ -329,7 +329,7 @@ m2_phys = fitp[1] + fitp[2] * phi2_phys
 ```
 """
 function lin_fit(x::Vector{<:Real}, y::Vector{uwreal}; wpm::Union{Dict{Int64,Vector{Float64}},Dict{String,Vector{Float64}}, Nothing}=nothing)
-    isnothing(wpm) ? uwerr.(y) : uwerr.(y, wpm)
+    isnothing(wpm) ? uwerr.(y) : [uwerr(yaux, wpm) for yaux in y]
     par = lin_fit(x, value.(y), err.(y))
     chisq(p, d) = sum((d .- p[1] .- p[2].*x).^2 ./ err.(y) .^2)
     (fitp, csqexp) = fit_error(chisq, par, y)
@@ -383,7 +383,7 @@ fit_routine(model, xdata, ydata, param=3, covar=true)
 ```
 """
 function fit_routine(model::Function, xdata::Array{<:Real}, ydata::Array{uwreal}, param::Int64=3; wpm::Union{Dict{Int64,Vector{Float64}},Dict{String,Vector{Float64}}, Nothing}=nothing)
-    #isnothing(wpm) ? uwerr.(ydata) : uwerr.(ydata, wpm)
+    isnothing(wpm) ? uwerr.(ydata) : [uwerr(yaux, wpm) for yaux in ydata]
     
     yval = value.(ydata)
     yer = err.(ydata)
@@ -411,8 +411,8 @@ function fit_routine(model::Function, xdata::Array{uwreal}, ydata::Array{uwreal}
         uwerr.(ydata)
         uwerr.(xdata)
     else
-        uwerr.(ydata, wpm)
-        uwerr.(xdata, wpm)
+        [uwerr(yaux, wpm) for yaux in ydata] 
+        [uwerr(xaux, wpm) for xaux in xdata]
     end
 
     yval = value.(ydata)