Bug fix in mpcac

src/juobs_obs.jl

@@ -129,6 +129,27 @@ function meff(corr::Vector{Corr},plat::Vector{Vector{Int64}}; kwargs...)
   return Tuple([r...] for r in eachrow(a))
 end
 
+function meff(corr::Vector{Vector{uwreal}},plat::Vector{Vector{Int64}}; kwargs...)
+  a = meff.(corr,plat; kwargs...)
+  if length(a[1])==1
+    return a
+  end
+
+  # Due to the broadcast, a is a Vector{Tuple{...}} where each element of the vector is the output of the function 
+  # The following lines of code reshape a into a Matrix such that each column 
+  # correspond to one Tuple and each row correspond to one component of the Tuple. The it will return each row 
+  # separately.
+  # 
+  # Example:
+  # V = [(a1,a2,a3),(b1,b2,b3)]
+  # V = collect.(V) #[[a1,a2,a3],[b1,b2,b3]]
+  # V stack(V) #[a1 b1; a2 b2; a3 b3]
+  # return Tuple([[r...] for e in eachrow(r)]) #([a1,b1],[a1,b2],[a3,b3])
+
+  a = stack(collect.(a));
+  return Tuple([r...] for r in eachrow(a))
+end
+
 @doc raw"""
     mpcac(a0p::Vector{uwreal}, pp::Vector{uwreal}, plat::Vector{Int64}; 
     ca::Float64=0.0, pl::Bool=false, data::Bool=false, kappa::Union{Vector{Float64}, Nothing}=nothing,
@@ -294,7 +315,25 @@ function mpcac(a0p::Vector{Vector{uwreal}},pp::Vector{Vector{uwreal}},plat::Vect
 end
 
 function mpcac(a0p::Vector{Corr},pp::Vector{Corr},plat::Vector{Vector{Int64}}; kwargs...)
-  return mpcac([a.obs for a in a0p],[p.obs for p in pp],plat;kwargs...)
+  a = mpcac.(a0p,pp,plat; kwargs...)
+  
+  if length(a[1])==1
+    return a
+  end
+
+  # Due to the broadcast, a is a Vector{Tuple{...}} where each element of the vector is the output of the function 
+  # with one set of a0p, pp and plat. The following lines of code reshape a into a Matrix such that each column 
+  # correspond to one Tuple and each row correspond to one component of the Tuple. The it will return each row 
+  # separately.
+  # 
+  # Example:
+  # V = [(a1,a2,a3),(b1,b2,b3)]
+  # V = collect.(V) #[[a1,a2,a3],[b1,b2,b3]]
+  # V stack(V) #[a1 b1; a2 b2; a3 b3]
+  # return Tuple([[r...] for e in eachrow(r)]) #([a1,b1],[a1,b2],[a3,b3])
+
+  a = stack(collect.(a));
+  return Tuple([r...] for r in eachrow(a))
 end
 
 ## Decay constants