if (FALSE) {
A<-c(0,0); B<-c(1,0); C<-c(1/2,sqrt(3)/2)
Tr<-rbind(A,B,C)
n<-10 #try also n<-20
Xp<-runif.tri(n,Tr)$g
M<-as.numeric(runif.tri(1,Tr)$g) #try also M<-c(1,1,1)
r<-1.4
PEdomTri(Xp,Tr,r,M)
IM<-IncMatPEtri(Xp,Tr,r,M)
dom.greedy #try also dom.exact(IM)
gr.gam<-dom.greedy(IM)
gr.gam
Xp[gr.gam$i,]
PEdomTri(Xp,Tr,r,M=c(.4,.4))
}
Run the code above in your browser using DataLab