if (FALSE) {
A<-c(0,0); B<-c(1,0); C<-c(1/2,sqrt(3)/2);
Te<-rbind(A,B,C);
n<-10
set.seed(1)
Xp<-runif.std.tri(n)$gen.points
M<-as.numeric(runif.std.tri(1)$g) #try also M<-c(.6,.2)
inc.mat<-IncMatCSTe(Xp,t=1.25,M)
inc.mat
sum(inc.mat)-n
NumArcsCSTe(Xp,t=1.25)
dom.greedy(inc.mat) #try also dom.exact(inc.mat) #might take a long time for large n
IndUBdom(inc.mat,1)
}
Run the code above in your browser using DataLab