Tweak hybrid example

ptiede · Sep 24, 2024 · 59a0e08 · 59a0e08
1 parent 12fc316
commit 59a0e08
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 2 deletions.
diff --git 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
@@ -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)