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

noreply at r-forge.r-project.org noreply at r-forge.r-project.org
Wed Dec 21 23:57:32 CET 2011


Author: dmbates
Date: 2011-12-21 23:57:32 +0100 (Wed, 21 Dec 2011)
New Revision: 1484

Modified:
   pkg/lme4Eigen/src/predModule.cpp
Log:
Delete some debugging code.


Modified: pkg/lme4Eigen/src/predModule.cpp
===================================================================
--- pkg/lme4Eigen/src/predModule.cpp	2011-12-21 22:56:56 UTC (rev 1483)
+++ pkg/lme4Eigen/src/predModule.cpp	2011-12-21 22:57:32 UTC (rev 1484)
@@ -70,19 +70,10 @@
 	    for(MSpMatrixd::InnerIterator rhsIt(d_Ut, j); rhsIt; ++rhsIt) {
 		Scalar                       y(rhsIt.value());
 		Index                        k(rhsIt.index());
-//		Rcpp::Rcout << "d_Ut row k = " << k << ", v = " << y;
-//		Rcpp::Rcout << ", Lambdat: col has " << d_Lambdat.innerNonZeros(k)
-//			    << " nonzeros" << std::endl;
 		MSpMatrixd::InnerIterator prdIt(d_LamtUt, j);
 		for (MSpMatrixd::InnerIterator lhsIt(d_Lambdat, k); lhsIt; ++lhsIt) {
 		    Index                    i = lhsIt.index();
-//		    Rcpp::Rcout << "lhsIt: i = " << i << ", v = " << lhsIt.value() << std::endl;
-//		    Rcpp::Rcout << "prdIt: i=" << prdIt.index() << std::endl;
-		    while (prdIt && prdIt.index() != i) {
-			++prdIt;
-//			if (!prdIt) Rcpp::Rcout << "End of product column" << std::endl;
-//			else Rcpp::Rcout << "i=" << prdIt.index() << std::endl;
-		    }
+		    while (prdIt && prdIt.index() != i) ++prdIt;
 		    if (!prdIt) throw runtime_error("logic error in updateLamtUt");
 		    prdIt.valueRef()          += lhsIt.value() * y;
 		}



More information about the Lme4-commits mailing list