[adegenet-commits] r826 - pkg/R

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

Author: jombart
Date: 2011-02-28 19:33:39 +0100 (Mon, 28 Feb 2011)
New Revision: 826

Modified:
   pkg/R/glFunctions.R
Log:
Amazing bug in glDotProd: 
the C function works for 1000 SNPs and 10,000 SNPs, but fails with e.g 10, 1001, or 1010 SNPs...
fmio 


Modified: pkg/R/glFunctions.R
===================================================================
--- pkg/R/glFunctions.R	2011-02-28 18:02:34 UTC (rev 825)
+++ pkg/R/glFunctions.R	2011-02-28 18:33:39 UTC (rev 826)
@@ -556,7 +556,7 @@
 
 
 #### TESTING PCA ####
-## M <- matrix(sample(c(0,1), 20*1e3, replace=TRUE), nrow=20)
+## M <- matrix(sample(c(0,1), 20*1000, replace=TRUE), nrow=20)
 ## rownames(M) <- paste("ind", 1:20)
 
 ## x <- new("genlight",M)
@@ -600,6 +600,9 @@
 ## round(t(toto$loadings) %*% toto$loadings,10) # must be diag(1,4)
 
 
+
+
+
 ## ## SPEED TESTS ##
 ## ## perform glPca
 ## M <- matrix(sample(c(0,1), 100*1e5, replace=TRUE), nrow=100)