inst/tests/test_sim2.R

In rcqls/EBSpatCGAL: EBSpat using CGAL

require(EBSpatCGAL)

new(Del2TermType2D) -> del2Term

#new(Del2TermType2D) -> del2Term2


del2 <- Delaunay()

insert(del2,matrix(runif(200,-350,350),ncol=2))

#### OLD CALL
# gd<-EBGibbs(~ 2 +Del2(th[1]*(l<=20)+th[2]*(20<l & l<=80),th=c(2,4)))

del2Term$set_struct(del2$rcpp())
#del2Term2$set_struct(del2$rcpp())

del2Term$infos <- c("l")
#del2Term2$infos <- c("x","l","l2")

del2Term$params <- list(th=c(2,4)) #,theta2=1)
#del2Term2$params <- list(theta=2)


del2Term$exprs<-list(substitute(th[1]*(l<=20)+th[2]*(20<l & l<=80))) #,a1=substitute(theta*c(l,sqrt(l2))))
del2Term$exprs.size <- c(1)
#del2Term2$exprs<-list(a3=substitute(theta*l))
#del2Term2$exprs.size <- c(1)

#new(Interaction,list(del2Term,del2Term2)) -> inter

#new(Interaction,list(del2Term)) -> inter


sim2 <- new(SimGibbsDel2D,list(del2Term),del2$rcpp(),c(-350,-350),c(350,350))

sim2$single <- 2

sim2$nb_runs <- 10000

(sc <- Scene()) %<<% window2d(c(-350,350),c(-350,350),xlab="",ylab="",main="what a beautiful plot!") 
sc %<<% points(del2,col="blue")  %<<% lines(del2,type="vor") 
sc %<<% lines(del2)

sim2$run( )

del2$save()

plot(sc)