Index: tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/gibbs/_nonLinBlock.tol
===================================================================
--- tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/gibbs/_nonLinBlock.tol	(revision d29988b88b00db97b6849236f39973f136935eb9)
+++ tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/gibbs/_nonLinBlock.tol	(revision 7cac98fea236b2d206ab37792107fe7778db6157)
@@ -81,6 +81,6 @@
     0
   });
-  Code _.get.prior = If(FunctionExist("Matrix","_.filterEval::get.prior"),
-    _.filterEval::get.prior,
+  Code _.get.priorLogDens = If(FunctionExist("Real","_.filterEval::get.priorLogDens"),
+    _.filterEval::get.priorLogDens,
     Matrix (Real paramIdx) {Rand(0,1,0,0) } );
   //Names of affected parameters in linear block 
@@ -164,5 +164,5 @@
       Set bounds = _.filterEval::get.bounds(i,u);
       Code fltEv = _.filterEval::eval;
-      Code get.prior = _.get.prior;
+      Code get.priorLogDens = _.get.priorLogDens;
       Matrix Z = _.Z;
       Real x0 = bounds[1];
@@ -182,5 +182,5 @@
         Real Do PutMatDat(u,i,1,old);
         Matrix fn = VMat2Mat(si*Mat2VMat(Z-Fi*omega));
-        fn << get.prior(i)
+        fn
       };
       Real LogDens.WhiteNoise(Real x)
@@ -190,5 +190,5 @@
       //WriteLn("TRACE LogDens.WhiteNoise("+_.colNames[i]+","<<x+")="<<ld);
       //Real UpdateMaxLogLikelihood(x,ld);
-        ld
+        ld+get.priorLogDens(i)
       };
       Real LogDens.ARMA(Real x)
@@ -201,5 +201,5 @@
       //WriteLn("TRACE LogDens.ARMA("+_.colNames[i]+","<<x+")="<<ld);
       //Real UpdateMaxLogLikelihood(x,ld);
-        ld
+        ld+get.priorLogDens(i)
       };
       Real LogDens.ARIMA(Real x)
@@ -217,5 +217,5 @@
       //WriteLn("TRACE LogDens.ARIMA("+_.colNames[i]+","<<x+")="<<ld);
       //Real UpdateMaxLogLikelihood(x,ld);
-        ld
+        ld+get.priorLogDens(i)
       };
       Code LogDens = Case(