merged merging into alejandro, idm rework for several replicas

src/juobs_tools.jl

@@ -368,24 +368,8 @@ function corr_obs_TSM(cdata1::Array{CData, 1}, cdata2::Array{CData, 1}; real::Bo
     
     nms = isnothing(nms) ? sum(replica_sl) : nms # assuming vcfg_sl >= vcfg_corr
     
-    if isnothing(idm_sl)
-        delta_vcfg = vcfg_sl[1][2] - vcfg_sl[1][1]
-        idm_sl = collect(1:delta_vcfg:replica_sl[1])
-        for ii in eachindex(cdata1)[2:end]
-            delta_vcfg = vcfg_sl[ii][2] - vcfg_sl[ii][1]
-            aux = collect(1:delta_vcfg:replica_sl[ii]) .+ replica_sl[ii-1]
-            append!(idm_sl, aux)
-        end
-    end
-    if isnothing(idm_corr)
-        delta_vcfg = vcfg_corr[1][2] - vcfg_corr[1][1]
-        idm_corr = collect(1:delta_vcfg:replica_corr[1])
-        for ii in eachindex(cdata2)[2:end]
-            delta_vcfg = vcfg_corr[ii][2] - vcfg_corr[ii][1]
-            aux = collect(1:delta_vcfg:replica_corr[ii]) .+ replica_corr[ii-1]
-            append!(idm_corr, aux)
-        end
-    end
+    idm_sl = isnothing(idm_sl) ? Int64.([vcfg_sl[1]; vcfg_sl[1][end] .+ vcfg_sl[2]]) : idm_sl
+    idm_corr = isnothing(idm_corr) ? Int64.([vcfg_corr[1]; vcfg_corr[1][end] .+ vcfg_corr[2]]) : idm_corr
 
     data1 = real ? getfield.(cdata1, :re_data) ./ L^3 : getfield.(cdata1, :im_data) ./ L^3
     data2 = real ? getfield.(cdata2, :re_data) ./ L^3 : getfield.(cdata2, :im_data) ./ L^3
@@ -422,7 +406,7 @@ function corr_obs_TSM(cdata1::Array{CData, 1}, cdata2::Array{CData, 1}; real::Bo
     else
         return Corr(obs1 + obs2, cdata1)
     end
-end    
+end     
 @doc raw"""
     corr_sym(corrL::Corr, corrR::Corr, parity::Int64=1)