R/hartman6.R

In DiceKriging: Kriging Methods for Computer Experiments

hartman6 <- function(x){ 
	x1 <- x[1]   
	x2 <- x[2]
	x3 <- x[3]
	x4 <- x[4]   
	x5 <- x[5]
	x6 <- x[6]

	a <- matrix(c(10.00,  0.05,  3.00, 17.00,
     3.00, 10.00,  3.50,  8.00,
     17.00, 17.00,  1.70,  0.05,
     3.50,  0.10, 10.00, 10.00,
     1.70,  8.00, 17.00,  0.10,
     8.00, 14.00,  8.00, 14.00), 6, 4, byrow=TRUE)
	p <-matrix(c(0.1312, 0.2329, 0.2348, 0.4047,
     0.1696, 0.4135, 0.1451, 0.8828,
     0.5569, 0.8307, 0.3522, 0.8732,
     0.0124, 0.3736, 0.2883, 0.5743,
     0.8283, 0.1004, 0.3047, 0.1091,
     0.5886, 0.9991, 0.6650, 0.0381), 6, 4, byrow=TRUE)
	c <- c(1.0, 1.2, 3.0, 3.2)
	d <- matrix(0,1,4)
	for (i in seq(1,4)){
		d[i]<-sum(a[,i]*(x - p[,i])^2)
	}
	f <- -sum(c*exp(-d))
	return(f)
}