HNUORTools: Operations Research Tools

geo<-new("GeoSituation")
#example taken from
#Bloech Management Methoden und Optimalplanung S. 72
geo<-add(geo,new("Warehouse", id="L1", x=25,   y=70,   supply = 100   ))
geo<-add(geo,new("Warehouse", id="L2", x=150,  y=115,  supply = 130   ))
geo<-add(geo,new("Warehouse", id="L3", x=80,   y=140,  supply = 170   )) 

geo<-add(geo,new("Customer",  id="K1", x=15,   y=130,  demand = 150   ))
geo<-add(geo,new("Customer",  id="K2", x=60,   y=80,   demand = 120   ))
geo<-add(geo,new("Customer",  id="K3", x=175,  y=140,  demand = 80   ))
geo<-add(geo,new("Customer",  id="K4", x=50,   y=100,  demand = 50   ))

 
x <- matrix(rep(0, 12), ncol=4, byrow=TRUE) 
# setting values from example
x[1,1] <- 100
x[2,1] <- 50
x[2,2] <- 30
x[2,4] <- 50
x[3,2] <- 90
x[3,3] <- 80
geo$tpp$x <- x 
geo$tpp$cij <- matrix(c(3,5,7,11,1,4,6,3,5,8,12,7), ncol=4, byrow=TRUE)
cij <- geo$tpp$cij
 

totalcosts <- round(sum(x*cij))# 
totalcosts == 2300 # true 

#use Modi-Method
geo<- TPP.MODI(geo, log=TRUE) 	
x <- geo$tpp$x
totalcosts <- round(sum(x*cij))# 
totalcosts == 2040 # true 
# Proofe results
x[1,2] == 20 # true
x[2,1] == 80 # true
x[3,1] == 70 # true
x[1,3] == 80 # true
x[2,4] == 50 # true
x[3,2] == 100 # true