[Help-glpk] problem with glp_ios_heur_sol

[Le Dinh Danh]

Hi all,

I am now having a problem with using glp_ios_heur_sol() to provide an initial integer feasible solution by heuristic.
To do so I did as follows:

void callback(glp_tree *tree, void *info)

{  if (glp_ios_reason(tree) == GLP_IHEUR &&
       glp_ios_curr_node(tree) = 1)
   {   /* lp relaxation to the root subproblem has been just
          solved */
       glp_ios_heur_sol(tree, <array of x[j] found by heuristic>);
   }
   return;
}

And then in the main function:

int main(...)
{  glp_prob *mip;
   glp_iocp parm;
   . . .
   <apply your own heuristic to find all x[j], j = 1,...,n>
   . . .
   glp_init_iocp(&parm);
   parm.cb_func = callback;
   parm.cb_info = ... ;
   ret = glp_intopt(mip, &parm);
   . . .
}

Although the array x[] is set correctly (I am sure after carefully checked), the result is still incorrect. Is anyone has successfully applied this can share something?

--

Le Dinh Danh

LIRMM, 161 rue Ada, 34095 Montpellier Cedex 5

Email: address@hidden; address@hidden               

Tel: (+84) 912.449.666; (+33) 6.52.98.59.66