CaseStudyForm  

Title  The Cosmic Computers Problem 
DateSubmitted  21 Sep 2012 
CaseStudyType  TeachingCaseStudy 
OperationsResearchTopics  IntegerProgramming, FacilityLocationProblem 
ApplicationAreas  Industrial Planning 
ProblemDescription 
The Cosmic Computer company wish to set up production in America. They are contemplating building production plants in up to 4 locations. Each location has planning restrictions that effectively determine the monthly production possible in each location. The monthly fixed costs of the different locations have been calculated by the accounting department, and are listed below:
 Location  Monthly Fixed Cost   San Francisco  $70,000   Los Angeles  $70,000   Phoenix  $65,000   Denver  $70,000  These plants will supply stores at 4 locations. Predicted demand and perunit transport costs for supplying the 4 production plants to the 4 stores are shown below, along with the quantity each factory can supply:
Where should the company set up their plants to minimise their total costs (fixed plus transportation)? Like most businesses, Cosmic Computers are big Excel users, but don't know much about optimisation. They would like to develop an AMPL model that uses SolverStudio to allow them to solve this problem using Excel. 
ProblemFormulation 
This problem is a facility location problem where the masterslave constraints control the existence of supply nodes for a transportation problem.
Traditional transportation problem formulations are as follows: Let's go through the mathematical programming formulation steps and consider what changes we might have to make. 1. Identify the Decision Variables
The decisions are twofold:
The flow variables determine where to ship the computers, so we need to add binary variables that govern which of the supply nodes exist, i.e., that are 1 if exists and 0 otherwise. 2. Formulate the Objective Function
The cost of transporting good through the network is already incorporated into the transportation problem objective, now we have to add the fixed (construction + maintenance) cost :
3. Formulate the Constraints
The constraints are the same as the transportation problem formulation, except that the supply is determined by the construction (or not) of a plant:
and we need to (at least) meet demand, but the formulation no longer needs to be balanced (otherwise we would always build all our construction plants):
4. Identify the data
The data in this problem is the same as that required for a transportation problem, i.e., , , , etc except that the fixed costs also needs to be specified.

ComputationalModel 
Since a transportation problem is embedded within the Cosmic Computers Problem, we will start with the AMPL model file transportation.mod (see The Transportation Problem in AMPL for details):
set SUPPLY_NODES; set DEMAND_NODES; param Supply {SUPPLY_NODES} >= 0, integer; param Demand {DEMAND_NODES} >= 0, integer; param Cost {SUPPLY_NODES, DEMAND_NODES} default 0; param Lower {SUPPLY_NODES, DEMAND_NODES} integer default 0; param Upper {SUPPLY_NODES, DEMAND_NODES} integer default Infinity; var Flow {i in SUPPLY_NODES, j in DEMAND_NODES} >= Lower[i, j], <= Upper[i, j], integer; minimize TotalCost: sum {i in SUPPLY_NODES, j in DEMAND_NODES} Cost[i, j] * Flow[i, j]; subject to UseSupply {i in SUPPLY_NODES}: sum {j in DEMAND_NODES} Flow[i, j] = Supply[i]; subject to MeetDemand {j in DEMAND_NODES}: sum {i in SUPPLY_NODES} Flow[i, j] = Demand[j]; We add new decision variables that decide if we build the production plants or not: var Build {SUPPLY_NODES} binary; We define a new parameter for the fixed cost (of construction + maintenance): param FixedCost {SUPPLY_NODES};
We incorporate the fixed cost into the objective function and the minimize TotalCost: sum {i in SUPPLY_NODES, j in DEMAND_NODES} Cost[i, j] * Flow[i, j] + sum {i in SUPPLY_NODES} FixedCost[i] * Build[i]; subject to UseSupply {i in SUPPLY_NODES}: sum {j in DEMAND_NODES} Flow[i, j] <= Supply[i] * Build[i]; SolverStudio does not use a .run file, and so we simply add our data, solve and output commands below the model. data SheetData.dat; # Get data from the spreadsheet check : sum {s in SUPPLY_NODES} Supply[s] >= sum {d in DEMAND_NODES} Demand[d]; option solver cplex; # Use this for the student version of AMPL solve; option display_1col 9999999; # a magic line SolverStudio needs display Build > Sheet; display Flow > Sheet; Note that we have included a check command to ensure there is enough supply. The spreadsheet file with this model (but still with some more work to do, as we discuss next) is SolverStudioCosmicComputersAMPLIncomplete.xlsx. For this problem, our data includes supplies and demands, lower and upper bounds on the flows, and costs. Note that only some routes have lower and upper bounds; a blank entry in one of these tables means AMPL uses the default values (0 and infinity) given in the model.
We need to define (i.e. give names to) these 'data items' on the spreadsheet, as shown below. These names must match those given in the AMPL model for the parameters and decision variables. The parameter data items include the supplies (cells B2:B6), demands (cells I4:L4), the flow costs, and the bounds. The decision variables are the flows and build decisions. One possible way of defining these data items is shown below. Notice, for example, that the Supply (C3:C6) and FixedCosts (D3:D6) are indexed over the SUPPLY_NODES (B3:B6) given in first column of the same table. The Demand is similarly defined over DEMAND_NODES listed in the first row of the Demand table. The Cost (C9:F12) is similarly defined, but note now that that indices are not the rows and column titles for the table, but actually refer to the SUPPLY_NODES and DEMAND_NODES defined not in the cost table, but in the supply and demand tables.
Another way of defining these items, in which we always index over the rows and columns of the actual table defining the data, is shown below. Either of these approaches works fine.
The spreadsheet file SolverStudioCosmicComputersAMPLIncomplete.xlsx only has some of the data items defined. We will add the missing data items, and then solve this model. If you have not yet installed SolverStudio, then either please visit SolverStudio to download this software, or, in the Engineering Science labs, install it from the path FoE Shared drive (S:) : ENGSCI : ugcourses : Courses : EngSci355. Download and Open the Example Spreadsheet 1/ Download and save the spreadsheet file SolverStudioCosmicComputersAMPLIncomplete.xlsx to the hard disk. 2/ Open your downloaded spreadsheet. 3/ Switch to the Data tab in Excel, and use SolverStudio's "Show/Hide Data" and "Show Data in Color" buttons to see the current data items. Add the Missing Data Items 4/ Click SolverStudio's "Edit Data" button to open SolverStudio's Data Items editor. 5/ We wish to define DEMAND_NODES. Make sure that the "Add New Data Item" item is selected. Then, type in a 'Name' of DEMAND_NODES, click the 'Cell Range' box, and then drag over the Demand cells (I3:L3). (Excel may add dollars and the sheet name; that's ok.) This data item is an AMPL set, and so we leave the 'Index Range(s)' blank. Click 'Add Data Item'.
6/ Repeat the previous step to add the 'Demand' data item (cell range I4:L4), which is indexed over DEMAND_NODES. To create this index, click on the first (i.e. top) "Index Range(s)" entry box, and either type in the name of an already defined data item (DEMAND_NODES), or type in or drag over the appropriate cells (I3:L3 in this case). Then click 'Add Data Item'. 7/ Repeat the previous step to add the 'FixedCost' data item (cell range D3:D6), which is indexed over SUPPLY_NODES (or B3:B6 if you prefer). 8/ Create the 'Flow' data item (which has SUPPLY_NODES (B3:B6) as its first index, and DEMAND_NODES (I4:L4) as its second, or cells H9:H12 and I8:L8 if you prefer). 9/ Still staying in the Data Items editor, check each data item is correct by clicking on it in the data items list, and observing the cells SolverStudio highlights. Correct any mistakes you have made. 10/ Close the Data Items editor. Check and Solve using AMPL 11/ If you cannot see the AMPL model text, click SolverStudio's 'Show Model' button to view the SolverStudio model editor. Note that you can print this model using SolverStudio's File... Print menu. (This menu also lets you print the model output.) 12/ By carefully working through the AMPL model, check that you have now defined as data items all the parameters in the AMPL model that have data on the spreadsheet, and all the decision variables in the AMPL model that need to be written to the spreadsheet. 13/ Use SolverStudio's Language menu to check that AMPL is selected as the language. 14/ We now wish to solve this problem using AMPL. If you have downloaded SolverStudio yourself, then you will need to install the AMPL Student Version; the SolverStudio AMPL menu has a command to do this automatically. This will take a minute or two. 15/ Use the SolverStudio 'Solve Model' button (or File... Solve, which has the shortcut key F5) to run and solve this model, and check that the flows and build decisions appear on the spreadsheet. 16/ SolverStudio creates a data file from your spreadsheet data. You can view this using the AMPL menu item 'View Last Data File'. This can be helpful in debugging. 17/ Save your spreadsheet to keep the changes you have made. 18/ SolverStudio provides access to various sources of online AMPL documentation under the AMPL menu. Try opening one of these. Congratulations; you now have a working AMPL model running in SolverStudio, and have completed this SolverStudio tutorial.

Results 
The output shows that CPLEX has not required any branchandbound nodes again. Is this problem also naturally integer? We can check by solving the LP relaxation using Unfortunately, our problem is not naturally integer. CPLEX performs a lot of "clever tricks" to solve mixedinteger programmes quickly. Let's look at the effect of some of these techniques.
First, we can see what is happening in the branchandbound tree by setting the CPLEX options for displaying the branchandbound process:
This shows that
The RINS heuristic looks for a feasible integer solution starting from the LP relaxation solution and the repair heuristic looks for a feasible integer solution from the initial values (here all zero). Let's turn the heuristics off (CPLEX options
Now the LP relaxation
Here the optimal integer solution is found at the second node and then the rest of the tree (6 branchandbound nodes in all) is used to ensure this is the optimum. CPLEX is very careful with its selection of a branching variable. Here, it branches on
Notice that the branchandbound tree is larger (8 nodes as opposed to 6), so this selection strategy did not work as well as the CPLEX default. While solving the Cosmic Computers Problem we have explored many of the "behind the scenes" work that CPLEX does for you. See Integer Programming with AMPL for more examples. In the Extra for Experts we will see another technique called constraint branching. 
Conclusions 
As with the Brewery Problem, the solution to the Cosmic Computers Problem can be presented in many ways. Figure 1 gives a graphical solution to the Cosmic Computers Problem.
Figure 1 Solution to the Cosmic Computers Problem

ExtraForExperts 
Another form of branch and bound implements branches on constraints that should take integer values. For example, in the Cosmic Computers Problem the sum of any subset of Build variables should be integer. Consider the LP relaxation solution:
The sum of the three nonzero
Our solution is integer so we can fathom this node. Now let's drop the Up constraint and branch down:
This solution still has fractional

StudentTasks 

I  Attachment  History  Action  Size  Date  Who  Comment 

xlsx  SolverStudioCosmicComputersAMPLIncomplete.xlsx  r2 r1  manage  20.5 K  20120924  23:21  AndrewMason  CosmicComputers Spreadsheet model with data items still to be finished 