#\donttest{
nx<-20; ny<-5;
set.seed(1)
Xp<-cbind(runif(nx,0,1),runif(nx,0,1))
Yp<-cbind(runif(ny,0,.25),
runif(ny,0,.25))+cbind(c(0,0,0.5,1,1),c(0,1,.5,0,1))
M<-c(1,1,1)
t<-1.5
IM<-inci.mat.undCS(Xp,Yp,t,M)
IM
pcds::dom.num.greedy(IM)
#}
Run the code above in your browser using DataLab