Fix division to multiplication of precision diagonal in minnesota

Author: ygeunkim

inst/include/bvhardraw.h

@@ -957,7 +957,7 @@ inline void minnesota_lambda(double& lambda, double& shape, double& rate, Eigen:
 														 boost::random::mt19937& rng) {
 	coef_prec.array() *= lambda;
 	// double gig_chi = (coef - coef_mean).squaredNorm();
-	double gig_chi = ((coef - coef_mean).array().square() / coef_prec.array()).sum();
+	double gig_chi = ((coef - coef_mean).array().square() * coef_prec.array()).sum();
 	lambda = sim_gig(1, shape - coef.size() / 2, 2 * rate, gig_chi, rng)[0];
 	cut_param(lambda);
 	coef_prec.array() /= lambda;