[Lme4-commits] r1602 - pkg/lme4Eigen/src

[noreply at r-forge.r-project.org]

Author: dmbates
Date: 2012-02-15 22:54:04 +0100 (Wed, 15 Feb 2012)
New Revision: 1602

Modified:
   pkg/lme4Eigen/src/external.cpp
Log:
Remove comments used in tracking down memory errors.


Modified: pkg/lme4Eigen/src/external.cpp
===================================================================
--- pkg/lme4Eigen/src/external.cpp	2012-02-15 21:53:19 UTC (rev 1601)
+++ pkg/lme4Eigen/src/external.cpp	2012-02-15 21:54:04 UTC (rev 1602)
@@ -314,16 +314,25 @@
 
     SEXP glmerWrkIter(SEXP pp_, SEXP rp_) {
 	BEGIN_RCPP;
-
+//	Rcpp::Rcout << "Begin glmerWrkIter" << std::endl;
 	XPtr<glmResp>    rp(rp_);
+//	Rcpp::Rcout << "Extracted response pointer" << std::endl;
+	const Eigen::VectorXd   wt(rp->sqrtWrkWt());
+//	Rcpp::Rcout << "Working weights: " << wt.head(7).adjoint() << std::endl;
 	XPtr<merPredD>   pp(pp_);
-	const Eigen::VectorXd   wt(rp->sqrtWrkWt());
+//	Rcpp::Rcout << "Extracted predictor pointer" << std::endl;
 	pp->updateXwts(wt);
+//	Rcpp::Rcout << "Updated predictor Xwts" << std::endl;
 	pp->updateDecomp();
+//	Rcpp::Rcout << "Updated predictor decomp" << std::endl;
 	pp->updateRes(rp->wrkResp());
+//	Rcpp::Rcout << "Updated predictor residuals" << std::endl;
 	pp->solve();
+//	Rcpp::Rcout << "Predictor solve" << std::endl;
 	rp->updateMu(pp->linPred(1.));
+//	Rcpp::Rcout << "Update mu" << std::endl;
 	pp->installPars(1.);
+//	Rcpp::Rcout << "Install pars" << std::endl;
 
 	return ::Rf_ScalarReal(rp->updateWts());