R/optim.R

In llaniewski/ASEMOO:

#!/home/llaniewski/bin/RS

####################################################################
#                                                                  #
#  ASEMOO: ASynchronious Efficient Multi-Objctive Optimization     #
#                                                                  #
#  optim.R : R script to generate new designs                      #
#                                                                  #
####################################################################

#----------------------- Preprocessing ----------------------------#
lg = data.frame(matrix(NA,nrow=1,ncol=0))

optim.asemoo = function(case, fun, opt=list(), maxiter = 19, models = DiceKrigingModels, sub.optimize=nsga2_optim(), tab = data.frame(design=c())) {

	models$init()

	for (n in c("constraints", "objectives", "parameters"))
	{
		case[[n]]$name = as.character(case[[n]]$name)
		rownames(case[[n]]) = case[[n]]$name
	}



	#----------------------- newdesign generation ---------------------#
	comment = "Good"
	opt$count = 1
	while (nrow(tab) < maxiter) {
		doe = nrow(case$parameters)*2 + 2
		opt$sd_scale = smart_scale(nrow(tab), maxiter, doe)

		newdesign = try(gen.new(tab=tab, case=case, models=models, opt=opt, doe=doe, sub.optimize=sub.optimize))
		if ("try-error" %in% class(newdesign)) {
			warning("Encountered error in gen.new. quiting optimization")
			comment="Error"
			break;
		}
		newdesign[,c(case$objectives$name, case$constraints$name)] = as.numeric(NA)
		newdesign$state = "running"
		if (nrow(tab) > 0) {
			tab = rbind(tab,newdesign)
		} else {
			tab = newdesign
		}
		sel = apply( is.na(tab[,c(case$objectives$name, case$constraints$name),drop=F] ), 1, any)
		sel = sel & (tab$state != "done")
		tab[sel, case$objectives$name] = fun(tab[sel,case$parameters$name,drop=F])
		tab$state[sel] = "done"
		save(tab, file=paste(opt$actual,paste(opt$casename, "Save", "Rdata", sep="."), sep="/"))
#		Sys.sleep(5)
	}
	list(tab=tab, comment=comment)
}