if (FALSE) {
A<-c(1,1); B<-c(2,0); C<-c(1.5,2);
Tr<-rbind(A,B,C);
n<-10
set.seed(1)
Xp<-runif.tri(n,Tr)$gen.points
M<-as.numeric(runif.tri(1,Tr)$g) #try also M<-c(1.6,1.2)
S<-rbind(Xp[1,],Xp[2,])
IndNAStri.domset(S,Xp,Tr,M)
S<-rbind(Xp[1,],Xp[2,],Xp[3,],Xp[5,])
IndNAStri.domset(S,Xp,Tr,M)
S<-rbind(c(.1,.1),c(.3,.4),c(.5,.3))
IndNAStri.domset(S,Xp,Tr,M)
IndNAStri.domset(c(.2,.5),Xp,Tr,M)
IndNAStri.domset(c(.2,.5),c(.2,.5),Tr,M)
IndNAStri.domset(Xp[5,],Xp[2,],Tr,M)
S<-rbind(Xp[1,],Xp[2,],Xp[3,],Xp[5,],c(.2,.5))
IndNAStri.domset(S,Xp[3,],Tr,M)
IndNAStri.domset(Xp,Xp,Tr,M)
P<-c(.4,.2)
S<-Xp[c(1,3,4),]
IndNAStri.domset(Xp,P,Tr,M)
IndNAStri.domset(S,P,Tr,M)
IndNAStri.domset(S,Xp,Tr,M)
IndNAStri.domset(rbind(S,S),Xp,Tr,M)
}
Run the code above in your browser using DataLab