*** empty log message ***

miscprob.cc

@@ -205,7 +205,7 @@ ReturnMatrix mvnrnd(const RowVector& mu, const SymmetricMatrix& covar, int nsamp
 
  float mvnpdf(const RowVector& vals, const RowVector& mu, const SymmetricMatrix& covar)
  {
-   return std::exp(-0.5*((vals-mu).t()*covar.i()*(vals-mu)).AsScalar())/(std::pow(covar.Determinant(),0.5)*std::pow(2*M_PI,vals.Ncols()/2.0));
+   return std::exp(-0.5*((vals-mu)*covar.i()*(vals-mu).t()).AsScalar())/(std::pow(covar.Determinant(),0.5)*std::pow(2*M_PI,vals.Ncols()/2.0));
  }
 
 // ReturnMatrix gammarnd(const int dim1, const int dim2,