tag:blogger.com,1999:blog-593563533834706486.post5716396360084887759..comments2015-02-10T07:23:43.000-05:00Comments on Yet Another Math Programming Consultant: Assignment ProblemErwin Kalvelagenhttps://plus.google.com/100547539949080099832noreply@blogger.comBlogger10125tag:blogger.com,1999:blog-593563533834706486.post-22734561711670015282009-11-19T05:34:25.274-05:002009-11-19T05:34:25.274-05:00I think GAMS has new link structure. As far as I u...I think GAMS has new link structure. As far as I understand our link should be rewritten soon and that will solve this problem.<br /><br />Btw the specialized network flow alg in MOSEK does much better than dual simplex on your assignment problem. We will most likely provide Mittleman with a few cases. They will be good for his network solver test.Erling D. Andersenhttps://www.blogger.com/profile/07306894197500659436noreply@blogger.comtag:blogger.com,1999:blog-593563533834706486.post-65237513042846482622009-11-18T12:04:45.110-05:002009-11-18T12:04:45.110-05:00These reformulations are explained in the document...These reformulations are explained in the documentation for the GAMS IO library. You may want to ask GAMS for a copy.Erwin Khttps://www.blogger.com/profile/09496091402502236997noreply@blogger.comtag:blogger.com,1999:blog-593563533834706486.post-33678482823676527982009-11-18T04:48:00.520-05:002009-11-18T04:48:00.520-05:00GAMS has an objective variable instead of an objec...GAMS has an objective variable instead of an objective function. If the objvar appears only once (linearly) in the model, and has no bounds, it can be substituted out. Most links do that automatically so the solver sees a proper objective function. This is especially important for nonlinear solvers.Erwin Khttps://www.blogger.com/profile/09496091402502236997noreply@blogger.comtag:blogger.com,1999:blog-593563533834706486.post-84447232744726109922009-11-18T03:45:25.278-05:002009-11-18T03:45:25.278-05:00Sorry, my last put should say if the objective is ...Sorry, my last put should say if the objective is in the constraints then the problem does not have network structure and hence we will detect the embedded network structure.<br /><br />So your examplepoint out a problem in the GAMS/MOSEK link.Erling D. Andersenhttps://www.blogger.com/profile/07306894197500659436noreply@blogger.comtag:blogger.com,1999:blog-593563533834706486.post-12268453398392279062009-11-18T03:42:19.381-05:002009-11-18T03:42:19.381-05:00Just a quick note. We downloaded your GAMS model. ...Just a quick note. We downloaded your GAMS model. Thanks. MOSEK has a mode where it detect the embedded network structure.<br />But GAMS or at least the GAMS/MOSEK link puts the objective in the constraints because in GAMS you do<br /><br />min z <br /><br />st c'x = z<br /> ...<br /><br />so we should change MOSEK so it take that into account.Erling D. Andersenhttps://www.blogger.com/profile/07306894197500659436noreply@blogger.comtag:blogger.com,1999:blog-593563533834706486.post-17310104819444282602009-11-18T01:43:26.927-05:002009-11-18T01:43:26.927-05:00I totally agree with both dual simplex is surprisi...I totally agree with both dual simplex is surprisingly fast. That was also the intended implication of my initial post. <br /><br />It is amazing it can compete with a specialized network simplex solver.Erling D. Andersenhttps://www.blogger.com/profile/07306894197500659436noreply@blogger.comtag:blogger.com,1999:blog-593563533834706486.post-46475339797492788552009-11-17T14:44:19.044-05:002009-11-17T14:44:19.044-05:00Erwin makes a good point. Sparse LP solvers are so...Erwin makes a good point. Sparse LP solvers are so fast that they are competitive with the performance of many special-purpose solvers. Furthermore, many problems with special structure - such as network models - have side constraints. And once you consider side constraints, you're usually better off with a general purpose solver.Greg Glocknerhttps://www.blogger.com/profile/14964879286706089383noreply@blogger.comtag:blogger.com,1999:blog-593563533834706486.post-17876176731766571792009-11-17T12:40:05.723-05:002009-11-17T12:40:05.723-05:00No doubt a better engineered network code will imp...No doubt a better engineered network code will improve this by a large margin (e.g. Cplex network algorithm solves these models in 0.25 and 1.58 seconds). Then we are comparing network codes. However, the point I wanted to make, somewhat inarticulate no doubt, was that a standard lp solver does not always have to be slow compared to specialized solvers. Anyway, the problems are in GAMS, and can be retrieved from http://amsterdamoptimization.com/models/assignmentR.zip.Erwin Khttps://www.blogger.com/profile/09496091402502236997noreply@blogger.comtag:blogger.com,1999:blog-593563533834706486.post-84254274791783437012009-11-17T06:53:39.631-05:002009-11-17T06:53:39.631-05:00If a GAMS program with your problem is available s...If a GAMS program with your problem is available somewhere then let me know.<br /><br />ErlingErling D. Andersenhttps://www.blogger.com/profile/07306894197500659436noreply@blogger.comtag:blogger.com,1999:blog-593563533834706486.post-46451277515692379752009-11-17T04:54:31.465-05:002009-11-17T04:54:31.465-05:00Dual simplex as I assume GuRoBi use is usually qui...Dual simplex as I assume GuRoBi use is usually quite fast. The specialized network flow optimizer in MOSEK (or in CPLEX) may be even faster.<br /><br /><br />ErlingErling D. Andersenhttps://www.blogger.com/profile/07306894197500659436noreply@blogger.com