## Regular example with obvious solution (minimum distance 48)
m <- 10 # Generate some points in the plane
z <- cbind(c(rep(0,m), rep(2,m), rep(5,m), rep(7,m)), rep(seq(0,m-1),4))
n <- nrow(z)
d <- compute_distance_matrix(z)
b <- build_tour_2tree(d, n)
b$distance # Distance 57.868
bi <- improve_tour_2opt(d, n, b$tour)
bi$distance # Distance 48 (optimum)
plot_tour(z,b)
plot_tour(z,bi)
## Random points
set.seed(1)
n <- 25
z <- cbind(runif(n,min=1,max=10),runif(n,min=1,max=10))
d <- compute_distance_matrix(z)
b <- build_tour_2tree(d, n)
b$distance # Distance 48.639
bi <- improve_tour_2opt(d, n, b$tour)
bi$distance # Distance 37.351
plot_tour(z,b)
plot_tour(z,bi)
Run the code above in your browser using DataLab