ensure that multi_normal_suff returns an allowable value

ecmerkle · Aug 26, 2023 · 9e7a2da · 9e7a2da
1 parent 4307b89
commit 9e7a2da
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 1 deletion.
diff --git a/DESCRIPTION b/DESCRIPTION
@@ -1,6 +1,6 @@
 Package: blavaan
 Title: Bayesian Latent Variable Analysis
-Version: 0.4-9.1155
+Version: 0.4-9.1156
 Authors@R: c(person(given = "Edgar", family = "Merkle",
                     role = c("aut", "cre"),
                     email = "merklee@missouri.edu",

diff --git a/inst/stan/stanmarg.stan b/inst/stan/stanmarg.stan
@@ -450,6 +450,8 @@ functions { // you can use these in R following `rstan::expose_stan_functions("f
     // using elementwise multiplication + sum here for efficiency
     out = -.5 * N * ( sum(Supdate[1:Nobs, 1:Nobs] .* (S + (xbar - Mu) * (xbar - Mu)')) + Supdate[Nobs + 1, Nobs + 1] + Nobs * log(2 * pi()) );
 
+    if(is_nan(out) || out == positive_infinity()) out = negative_infinity();
+
     return out;
   }