[Genabel-commits] r1100 - pkg/MixABEL/src/MXlib

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

Author: wja
Date: 2013-02-12 21:20:19 +0100 (Tue, 12 Feb 2013)
New Revision: 1100

Modified:
   pkg/MixABEL/src/MXlib/twovarcomp.cpp
Log:
fixed overflow bug in twovarcomp.cpp

Modified: pkg/MixABEL/src/MXlib/twovarcomp.cpp
===================================================================
--- pkg/MixABEL/src/MXlib/twovarcomp.cpp	2013-02-12 17:53:53 UTC (rev 1099)
+++ pkg/MixABEL/src/MXlib/twovarcomp.cpp	2013-02-12 20:20:19 UTC (rev 1100)
@@ -854,7 +854,8 @@
 
   for(pdev_work->it=0;pdev_work->it<n;pdev_work->it++)
 	{
-		gsl_vector_view rowview;
+		gsl_vector_view rowview=gsl_matrix_row(pdev_work->pa_mat,pdev_work->it);
+		
     		gsl_vector_scale(&(rowview.vector), gsl_vector_get(pdev_work->pa2_vec,pdev_work->it));
 	}