testing some distance computations

main.R

require(igraph)
require(RColorBrewer)
require(rgl)
G <- make_lattice(length=10,  dim=2)
r <- c()
for (i in 4:6) r <- c(r, (4:6) + 10 * i)
G <- delete_vertices(G, r)
n <- length(V(G))
u0 <- rep(0, n)
u0[1] <- 1
A <- as_adj(G)
                                        # laplacian_matrix(G)
t <- .00001
B <- (diag(n) - t * A)
u <- round(log(solve(B, u0)) / log(t), 0)
m <- length(unique(u))
redblue<-colorRampPalette(c("red","blue"))
plot(G, vertex.color=redblue(m)[u-min(u)+1])
rglplot(G, vertex.color=redblue(n)[u+1])
plot(distances(G)[,1] ~ u)