Index: tol/stdlib/math/stat/models/bayesian/arima/_arma_process.tol
===================================================================
--- tol/stdlib/math/stat/models/bayesian/arima/_arma_process.tol	(revision 188403563b622a77a0e8d4b5d98d12223af4e509)
+++ tol/stdlib/math/stat/models/bayesian/arima/_arma_process.tol	(revision f2ee099d0b2364ec6e0236aa07c3bf844a1a1eff)
@@ -175,7 +175,5 @@
     Real _.n = Max(_.p,_.q)+_.p+1;
   //Real _trc.6 = { WriteLn("TRACE FastCholeskiCovFactor _.n=" << _.n); 0 };
-    //Calculates de main minor of length _.n of joint ARMA covariance
-    VMatrix _.cov_arma_n = 
-      ARMAProcess::DoDrop(Pol2VMat(ARMATACov(_.ar,_.ma,_.n), _.n, _.n));
+
     //Calculates de main minor of length _.n of just pure AR covariance
     VMatrix _.cov_ar_n = 
@@ -193,4 +191,8 @@
       (ARMAProcess::ExtendBandedLowTriangToeplitz(_.Li_ar_n,_.m)+_.n^2)-_.n^2;
   //Real _trc.10 = { WriteLn("TRACE FastCholeskiCovFactor _.Li_ar=" << (Matrix VMat2Mat(_.Li_ar))); 0 };
+
+    //Calculates de main minor of length _.n of joint ARMA covariance
+    VMatrix _.cov_arma_n = 
+      ARMAProcess::DoDrop(Pol2VMat(ARMATACov(_.ar,_.ma,_.n), _.n, _.n));
     //Filters full ARMA covariance of AR factor 
     VMatrix _.cov_ma_n = (_.Li_ar_n*_.cov_arma_n*Tra(_.Li_ar_n)+_.n^2)-_.n^2;