From 7ff025cd5fe5bb4e51dfb92209f8221a3c8281ad Mon Sep 17 00:00:00 2001
From: Reinhold Willcox <33348824+reinhold-willcox@users.noreply.github.com>
Date: Tue, 17 Sep 2024 16:52:22 +0200
Subject: [PATCH] fixed typo with likelihood calculation of :\omega (#63)

* fixed typo with likelihood calculation of :\omega

* updated version number
---
 Project.toml        | 2 +-
 src/TuringModels.jl | 6 +++---
 2 files changed, 4 insertions(+), 4 deletions(-)

diff --git a/Project.toml b/Project.toml
index dec15ac..5a201cd 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,7 +1,7 @@
 name = "SideKicks"
 uuid = "ae7baef6-f8d7-4188-9988-8bf5c4952615"
 authors = ["Pablo Marchant ", "Reinhold Willcox "]
-version = "0.3.2"
+version = "0.3.3"
 
 [deps]
 CairoMakie = "13f3f980-e62b-5c42-98c6-ff1f3baf88f0"
diff --git a/src/TuringModels.jl b/src/TuringModels.jl
index 976460e..ebcc41d 100644
--- a/src/TuringModels.jl
+++ b/src/TuringModels.jl
@@ -343,7 +343,7 @@ function create_general_mcmc_model(;
             elseif obs_symbol == :ω_f
                 use_cauchy ?
                     obs_vals[ii] ~ WrappedCauchy(ω_f, obs_errs[ii]) :
-                    obs_vals[ii] ~ ModVonMises(Ω_f, 1/obs_errs[ii]^2)
+                    obs_vals[ii] ~ ModVonMises(ω_f, 1/obs_errs[ii]^2)
             elseif obs_symbol == :v_N
                 use_cauchy ?
                     obs_vals[ii] ~ Cauchy(v_N, obs_errs[ii]) :
@@ -360,9 +360,9 @@ function create_general_mcmc_model(;
         end
 
         dm2 = m2_i - m2_f
-        return     ( m1_i,  m2_i,  P_i,  e_i,  a_i,  ν_i,  vkick, θ, ϕ,  frac,  dm2,  i_f,  ω_f,  Ω_f,  m1_f, m2_f,  a_f,  P_f,  e_f,  K1,  K2,  v_N,  v_E,  v_r,  vsys)
+        return     ( m1_i,  m2_i,  P_i,  e_i,  a_i,  ν_i,  vkick,  θ,  ϕ,  frac,  dm2,  i_f,  ω_f,  Ω_f,  m1_f,  m2_f,  a_f,  P_f,  e_f,  K1,  K2,  v_N,  v_E,  v_r,  vsys)
     end
-    return_props = [:m1_i, :m2_i, :P_i, :e_i, :a_i, :ν_i, :vkick, :θ, :ϕ, :frac, :dm2, :i_f, :ω_f, :Ω_f, :m1_f,:m2_f, :a_f, :P_f, :e_f, :K1, :K2, :v_N, :v_E, :v_r, :vsys]
+    return_props = [:m1_i, :m2_i, :P_i, :e_i, :a_i, :ν_i, :vkick, :θ, :ϕ, :frac, :dm2, :i_f, :ω_f, :Ω_f, :m1_f, :m2_f, :a_f, :P_f, :e_f, :K1, :K2, :v_N, :v_E, :v_r, :vsys]
 
     # Need to combine some of the observations to compare against the predicted output
     obs_vals_cgs = observations.vals .* observations.units