> I am wondering how we can do initialization in GLPK. > In my problems, I am optimizating several thousands of small LP > problems. I need their solutions to be similar to each other. However, > due to various noise in these problems, some of them are not feasible > or not easy to converge. In result, the output from GLPK for these > problems are not close to each other at all. One solution I am > thinking about is to use the solution from one problem as the initial > solution of the next problem. In this way, even if the next problem > will not converge, the result might still be close to the initial > solution.
> Could anybody tell me how to specify the initial solution for GLPK? Or > is it possible to do? You can do that by specifying the initial basis which the glpk simplex solver will start from; see api routines glp_set_row_stat and glp_set_col_stat in the reference manual. However, this will not help you, because there is no guarantee that the final basis obtained by the solver will be similar to the initial one. You should formalize what "to be similar" means in your case. Are solutions x = (0, 1, 0, .5) and x' = (.5, 0, 1, 0) similar to each other? If so/not, why? _______________________________________________ Help-glpk mailing list [email protected] http://lists.gnu.org/mailman/listinfo/help-glpk
