Skip to content

J
juobs
Commit d85d97b7 authored by Javier's avatar Javier
Browse files
Options

minor bug md_sea 

fluc_obs computed twice
parent 30a86546
Hide whitespace changes
Inline Side-by-side
Showing with 1 addition and 2 deletions
src/juobs_tools.jl
View file @ d85d97b7
...@@ -227,7 +227,6 @@ function md_sea(a::uwreal, md::Vector{Matrix{Float64}}, ws::ADerrors.wspace=ADer ...@@ -227,7 +227,6 @@ function md_sea(a::uwreal, md::Vector{Matrix{Float64}}, ws::ADerrors.wspace=ADer
md_aux = cat(md_aux, md[k][:, 1:ivrep[k]], dims=2) md_aux = cat(md_aux, md[k][:, 1:ivrep[k]], dims=2)
end end
fluc_obs = getfield.(ws.fluc[p], :delta)
fluc_md = md_aux .- mean(md_aux, dims=2) fluc_md = md_aux .- mean(md_aux, dims=2)
uwerr(a) uwerr(a)
fluc_obs = mchist(a, id) fluc_obs = mchist(a, id)
...@@ -277,7 +276,7 @@ derm = [[corr_pp_d1[k], corr_pp_d2[k]] for k = 1:length(pp_d1)] ...@@ -277,7 +276,7 @@ derm = [[corr_pp_d1[k], corr_pp_d2[k]] for k = 1:length(pp_d1)]
m = meff(corr_pp[1], plat) m = meff(corr_pp[1], plat)
m_md1, m_md2 = md_val(m, corr_pp[1], derm[1]) m_md1, m_md2 = md_val(m, corr_pp[1], derm[1])
m_shifted = m + 2 * dm1 * m_md1 + dm2 * m_md2 m_shifted = m + dm1 * m_md1 + dm2 * m_md2
``` ```
""" """
function md_val(a::uwreal, obs::Corr, derm::Vector{Corr}) function md_val(a::uwreal, obs::Corr, derm::Vector{Corr})
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment