Bug fix for options montecarlo and cubature.

matlab/nonlinear-filters/nonlinear_kalman_filter.m

@@ -187,8 +187,8 @@ for t=1:sample_size
             lik(t)=-Inf;
             return
         end
+        PredictedObservedVarianceSquareRoot = chol(PredictedObservedVariance);
     end
-    PredictedObservedVarianceSquareRoot = chol(PredictedObservedVariance);
     PredictedObservedInverseVarianceSquareRoot = PredictedObservedVarianceSquareRoot\eye(n); % Inverse of the Cholesky -> inv(PredictedObservedVariance) = A*A'
     StandardPredictionError = PredictedObservedInverseVarianceSquareRoot'*PredictionError;
     lik(t) = -.5*n*log(2*pi) + sum(log(diag(PredictedObservedInverseVarianceSquareRoot))) - .5*dot(StandardPredictionError, StandardPredictionError);