/** constraint deactivation notification method of constraint handler */
static
SCIP_DECL_CONSDEACTIVE(consDeactiveSamediff)
{ /*lint --e{715}*/
SCIP_CONSDATA* consdata;
SCIP_PROBDATA* probdata;
assert(scip != NULL);
assert(strcmp(SCIPconshdlrGetName(conshdlr), CONSHDLR_NAME) == 0);
assert(cons != NULL);
consdata = SCIPconsGetData(cons);
assert(consdata != NULL);
assert(consdata->propagated || SCIPgetNChildren(scip) == 0);
probdata = SCIPgetProbData(scip);
assert(probdata != NULL);
/* check if all variables which are not fixed locally to zero are valid for this constraint/node */
assert( consdataCheck(scip, probdata, consdata) );
SCIPdebugMessage("deactivate constraint <%s> at node <%"SCIP_LONGINT_FORMAT"> in depth <%d>: ",
SCIPconsGetName(cons), SCIPnodeGetNumber(consdata->node), SCIPnodeGetDepth(consdata->node));
SCIPdebug( consdataPrint(scip, consdata, NULL) );
/* set the number of propagated variables to current number of variables is SCIP */
consdata->npropagatedvars = SCIPprobdataGetNVars(probdata);
/* check if all variables are valid for this constraint */
assert( consdataCheck(scip, probdata, consdata) );
return SCIP_OKAY;
}
/** node selection method of node selector */
static
SCIP_DECL_NODESELSELECT(nodeselSelectMaster)
{
SCIP_NODESELDATA* nodeseldata;
SCIP_NODE** nodes;
SCIP_CONS* origcons;
SCIP_CONS* parentorigcons;
SCIP_CONS* parentmastercons;
SCIP* origscip;
int nnodes;
SCIP_Longint orignodenumber;
assert(nodesel != NULL);
assert(strcmp(SCIPnodeselGetName(nodesel), NODESEL_NAME) == 0);
assert(scip != NULL);
assert(selnode != NULL);
nodeseldata = SCIPnodeselGetData(nodesel);
assert(nodeseldata != NULL);
origscip = GCGpricerGetOrigprob(scip);
*selnode = NULL;
orignodenumber = SCIPnodeGetNumber(SCIPgetCurrentNode(origscip));
if( orignodenumber != nodeseldata->lastorignodenumber )
{
nodeseldata->lastorignodenumber = orignodenumber;
origcons = GCGconsOrigbranchGetActiveCons(origscip);
parentorigcons = GCGconsOrigbranchGetParentcons(origcons);
/* check whether the current node is the root node and has no parent */
if( parentorigcons == NULL )
{
assert((GCGconsOrigbranchGetNode(origcons) == SCIPgetRootNode(origscip)) || ( GCGconsOrigbranchGetNode(origcons) == NULL) );
assert(GCGconsOrigbranchGetMastercons(origcons) != NULL);
assert((GCGconsMasterbranchGetNode(GCGconsOrigbranchGetMastercons(origcons)) == SCIPgetRootNode(scip)) || (GCGconsMasterbranchGetNode(GCGconsOrigbranchGetMastercons(origcons)) == NULL));
*selnode = SCIPgetRootNode(scip);
SCIPdebugMessage("selected root node in the master program\n");
}
else
{
parentmastercons = GCGconsOrigbranchGetMastercons(parentorigcons);
assert(parentmastercons != NULL);
assert( (GCGconsOrigbranchGetChild1cons(parentorigcons) == origcons)
!= (GCGconsOrigbranchGetChild2cons(parentorigcons) == origcons));
/* the original cons is the left child of its parentcons,
select the left child of the corresponding parentcons in the master*/
if( GCGconsOrigbranchGetChild1cons(parentorigcons) == origcons )
{
assert(GCGconsMasterbranchGetChild1cons(parentmastercons) != NULL);
assert(GCGconsMasterbranchGetNode(GCGconsMasterbranchGetChild1cons(parentmastercons)) != NULL);
*selnode = GCGconsMasterbranchGetNode(GCGconsMasterbranchGetChild1cons(parentmastercons));
SCIPdebugMessage("Master nodeselector selected node %"SCIP_LONGINT_FORMAT" corresponding to node %"SCIP_LONGINT_FORMAT" in the original program, since the parents (%"SCIP_LONGINT_FORMAT"/o, %"SCIP_LONGINT_FORMAT"/m) are linked\n",
SCIPnodeGetNumber(*selnode), SCIPnodeGetNumber(GCGconsOrigbranchGetNode(origcons)),
SCIPnodeGetNumber(GCGconsOrigbranchGetNode(parentorigcons)),
SCIPnodeGetNumber(GCGconsMasterbranchGetNode(parentmastercons)));
}
/* the original cons is the right child of its parentcons,
select the right child of the corresponding parentcons in the master */
else
{
assert(GCGconsOrigbranchGetChild2cons(parentorigcons) == origcons);
assert(GCGconsMasterbranchGetChild2cons(parentmastercons) != NULL);
assert(GCGconsMasterbranchGetNode(GCGconsMasterbranchGetChild2cons(parentmastercons)) != NULL);
*selnode = GCGconsMasterbranchGetNode(GCGconsMasterbranchGetChild2cons(parentmastercons));
SCIPdebugMessage("Master nodeselector selected node %"SCIP_LONGINT_FORMAT" corresponding to node %"SCIP_LONGINT_FORMAT" in the original program, since the parents (%"SCIP_LONGINT_FORMAT"/o, %"SCIP_LONGINT_FORMAT"/m) are linked\n",
SCIPnodeGetNumber(*selnode), SCIPnodeGetNumber(GCGconsOrigbranchGetNode(origcons)),
SCIPnodeGetNumber(GCGconsOrigbranchGetNode(parentorigcons)),
SCIPnodeGetNumber(GCGconsMasterbranchGetNode(parentmastercons)));
}
}
if( *selnode == NULL )
{
SCIPerrorMessage("nodesel_master could not find a node corresponding to the current original node!\n");
}
assert(*selnode != NULL);
/* set the dual bound to the lower bound of the corresponding original node */
SCIP_CALL( SCIPupdateNodeDualbound(scip, *selnode, SCIPgetNodeLowerbound(origscip, SCIPgetCurrentNode(origscip))) );
}
else
{
SCIPdebugMessage("select random node\n");
if( SCIPgetNChildren(scip) > 0 )
{
SCIP_CALL( SCIPgetChildren(scip, &nodes, &nnodes) );
*selnode = nodes[0];
}
else if( SCIPgetNSiblings(scip) > 0 )
{
SCIP_CALL( SCIPgetSiblings(scip, &nodes, &nnodes) );
*selnode = nodes[0];
}
else if( SCIPgetNLeaves(scip) > 0 )
{
SCIP_CALL( SCIPgetLeaves(scip, &nodes, &nnodes) );
*selnode = nodes[0];
}
}
#ifndef NDEBUG
GCGconsOrigbranchCheckConsistency(origscip);
GCGconsMasterbranchCheckConsistency(scip);
#endif
return SCIP_OKAY;
}
