[Genabel-commits] r1637 - branches/ProbABEL-0.50/src

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

Author: maartenk
Date: 2014-03-09 22:47:13 +0100 (Sun, 09 Mar 2014)
New Revision: 1637

Modified:
   branches/ProbABEL-0.50/src/reg1.cpp
Log:
-r1636 introduced error causing wrong results: order of dimensions in a partial fixed matrix was swapped

Modified: branches/ProbABEL-0.50/src/reg1.cpp
===================================================================
--- branches/ProbABEL-0.50/src/reg1.cpp	2014-03-09 19:56:25 UTC (rev 1636)
+++ branches/ProbABEL-0.50/src/reg1.cpp	2014-03-09 21:47:13 UTC (rev 1637)
@@ -326,7 +326,7 @@
     }
     else if (X.data.cols() == 2)
     {
-        Matrix<double, 2, Dynamic> tXW =  W_masked.masked_data->data*X.data;
+        Matrix<double,  Dynamic,2> tXW =  W_masked.masked_data->data*X.data;
         Matrix2d xWx = tXW.transpose() * X.data;
         Ch = LDLT<MatrixXd> (xWx);
         Vector2d beta_2f = Ch.solve(tXW.transpose() * Y);