TODO

   1 n = 100; m = 70; p = 5
   2 X = matrix(runif(n*p, -10, 10), nrow=n)
   3 Y = matrix(runif(n*m, -5, 15), nrow=n)
   4
   5 V1 = valse::valse(X, Y, fast=FALSE)
   6 Error in while (!pi2AllPositive) { :
   7   missing value where TRUE/FALSE needed
   8
   9 V2 = valse::valse(X, Y, fast=TRUE)
  10 list()
  11 Error in out[[ind_uniq[l]]] :
  12   attempt to select less than one element in get1index
  13
  14 ==> Error identified: line 61 in initSmallEM.R, division by 0
  15 It occurs also for smallers values of n and m, e.g.: n = 20; m = 20; p = 3