source("setup.R")
source("anneal.R")
prim7 <- read.table("prim7", header=T, sep="\t")


trace1 <- scagnostics.anneal(prim7,7,temp=0.25, stop.it="NRUNS")
#
summary(trace1)

#############################
# start up ggobi

g <- ggobi(prim7)
d <- displays(g)[[1]] 
pmode(d) <- "2D Tour"

# find best projection from previous run
pp <- trace1$projs[[which.max(trace1$index)]]
ggobi_display_set_tour_projection(d,pp)
summary(trace1$index)


# restart annealing process
trace2 <- scagnostics.anneal(prim7,7,temp=0.5, cool=0.6, start=pp, stop.it="CONV", show.index=T, gg=T)
#############################

# not run
trace3 <- scagnostics.anneal(prim7,1,temp=0.5 , stop.it="MAXRUNS")
prim7.scag <- scagnostics.trace(prim7, trace3)
#

g$trace <- as.data.frame(cbind(1:length(trace3$index), prim7.scag)) 
linker <- function(gg, gplot, idx, data) { 
	if (idx != -1) {
		idx <- idx + 1
		ggobi_display_set_tour_projection(d,trace3$proj[[idx]])
	}
}
ggsgnl <- gSignalConnect(g, "identify-point", linker)
####################
gSignalHandlerDisconnect(g,ggsgnl)