void glp_erase_prob(glp_prob *lp)
{ glp_tree *tree = lp->tree;
if (tree != NULL && tree->reason != 0)
xerror("glp_erase_prob: operation not allowed\n");
delete_prob(lp);
create_prob(lp);
return;
}
void BackProjection::prep_memory()
{
t_mem.start();
/* Create back projection problem */
create_prob();
int i, j;
guess = (unsigned char *)malloc(r * c * sizeof(unsigned char));
rscore = (F_TYPE *)malloc(r * sizeof(F_TYPE));
cscore = (F_TYPE *)malloc(c * sizeof(F_TYPE));
uscore = (F_TYPE *)malloc((r + c - 1) * sizeof(F_TYPE));
dscore = (F_TYPE *)malloc((r + c - 1) * sizeof(F_TYPE));
for (i = 0; i < r; i++)
{
for (j = 0; j < c; j++)
{
guess[ i * c + j ] = FALSE;
}
}
for (i = 0; i < r; i++) {
rscore[i] = (F_TYPE) rproj[i] / (F_TYPE) rband[i];
}
for (i = 0; i < c; i++) {
cscore[i] = (F_TYPE) cproj[i] / (F_TYPE) cband[i];
}
for (i = 0; i < r + c - 1; i++) {
uscore[i] = (F_TYPE) uproj[i] / (F_TYPE) uband[i];
dscore[i] = (F_TYPE) dproj[i] / (F_TYPE) dband[i];
}
t_mem.stop();
}
glp_prob *glp_create_prob(void)
{ glp_prob *lp;
lp = xmalloc(sizeof(glp_prob));
create_prob(lp);
return lp;
}
////////////////////////////////////////////////////////////////////////////////
//create, delete problem
void MtxLP::createProb(){
create_prob();
create_index();
}
void glp_erase_prob(glp_prob *lp)
{ delete_prob(lp);
create_prob(lp);
return;
}
