From 59a0e085dc9824d31ad141ef51c9f69c1c575db2 Mon Sep 17 00:00:00 2001 From: Paul Tiede Date: Mon, 23 Sep 2024 22:15:34 -0400 Subject: [PATCH] Tweak hybrid example --- examples/advanced/HybridImaging/main.jl | 3 ++- src/mrf_image.jl | 2 +- 2 files changed, 3 insertions(+), 2 deletions(-) diff --git a/examples/advanced/HybridImaging/main.jl b/examples/advanced/HybridImaging/main.jl index 756f4597..9cbecc06 100644 --- a/examples/advanced/HybridImaging/main.jl +++ b/examples/advanced/HybridImaging/main.jl @@ -40,7 +40,7 @@ using Pyehtim # For reproducibility we use a stable random number genreator using StableRNGs -rng = StableRNG(44) +rng = StableRNG(49) # To download the data visit https://doi.org/10.25739/g85n-f134 @@ -246,6 +246,7 @@ fig |> DisplayAs.PNG |> DisplayAs.Text # Finally, let's take a look at some of the ring parameters + figd = CM.Figure(;resolution=(650, 400)); p1 = CM.density(figd[1,1], rad2μas(chain.sky.r)*2, axis=(xlabel="Ring Diameter (μas)",)) p2 = CM.density(figd[1,2], rad2μas(chain.sky.σ)*2*sqrt(2*log(2)), axis=(xlabel="Ring FWHM (μas)",)) diff --git a/src/mrf_image.jl b/src/mrf_image.jl index ad37080c..54fa5419 100644 --- a/src/mrf_image.jl +++ b/src/mrf_image.jl @@ -44,7 +44,7 @@ end function _apply_fluctuations(t::VLBIImagePriors.LogRatioTransform, mimg::AbstractArray, δ::AbstractArray) - @argcheck _checknorm(mimg) "Mean image must have unit flux when using log-ratio transformations in apply_fluctuations" + @argcheck _checknorm(mimg) "Mean image must have unit flux when using log-ratio transformations in apply_fluctuations while it seems to be $(sum(mimg))" r = to_simplex(t, baseimage(δ)) r .= r.*baseimage(mimg) r .= r./_fastsum(r)