/** execution method of event handler */
static
SCIP_DECL_EVENTEXEC(eventExecSofttimelimit)
{ /*lint --e{715}*/
SCIP_EVENTHDLRDATA* eventhdlrdata;
assert(eventhdlr != NULL);
assert(strcmp(SCIPeventhdlrGetName(eventhdlr), EVENTHDLR_NAME) == 0);
assert(event != NULL);
assert(scip != NULL);
assert(SCIPeventGetType(event) == SCIP_EVENTTYPE_BESTSOLFOUND);
eventhdlrdata = SCIPeventhdlrGetData(eventhdlr);
assert(eventhdlrdata != NULL);
SCIPdebugMessage("exec method of event handler for soft time limit\n");
SCIP_CALL( SCIPsetRealParam(scip, "limits/time", eventhdlrdata->softtimelimit) );
/* notify SCIP that your event handler wants to drop the event type best solution found */
SCIP_CALL( SCIPdropEvent(scip, SCIP_EVENTTYPE_BESTSOLFOUND, eventhdlr, NULL, eventhdlrdata->filterpos) );
eventhdlrdata->filterpos = -1;
/* print best solution value */
SCIPinfoMessage(scip, NULL, "changed time limit to %.1f after first solution was found\n",
eventhdlrdata->softtimelimit);
return SCIP_OKAY;
}
/** execution method of objective change event handler */
static
SCIP_DECL_EVENTEXEC(eventExecIntobj)
{ /*lint --e{715}*/
SCIP_EVENTHDLRDATA* eventhdlrdata;
SCIP_SEPADATA* sepadata;
SCIP_VAR* var;
SCIP_Real objdelta;
eventhdlrdata = SCIPeventhdlrGetData(eventhdlr);
sepadata = (SCIP_SEPADATA*)eventhdlrdata;
assert(sepadata != NULL);
/* we don't have anything to do, if the objective value inequality doesn't yet exist */
if( sepadata->objrow == NULL )
return SCIP_OKAY;
var = SCIPeventGetVar(event);
switch( SCIPeventGetType(event) )
{
case SCIP_EVENTTYPE_VARADDED:
SCIPdebugMessage("variable <%s> with obj=%g was added to the problem\n", SCIPvarGetName(var), SCIPvarGetObj(var));
objdelta = SCIPvarGetObj(var);
if( !SCIPisZero(scip, objdelta) )
{
SCIP_CALL( SCIPaddVarToRow(scip, sepadata->objrow, var, SCIPvarGetObj(var)) );
}
break;
case SCIP_EVENTTYPE_OBJCHANGED:
SCIPdebugMessage("variable <%s> changed objective value from %g to %g\n",
SCIPvarGetName(var), SCIPeventGetOldobj(event), SCIPeventGetNewobj(event));
objdelta = SCIPeventGetNewobj(event) - SCIPeventGetOldobj(event);
SCIP_CALL( SCIPaddVarToRow(scip, sepadata->objrow, var, objdelta) );
break;
default:
SCIPerrorMessage("invalid event type %x\n", SCIPeventGetType(event));
return SCIP_INVALIDDATA;
}
return SCIP_OKAY;
}
/** execution method of event handler */
static
SCIP_DECL_EVENTEXEC(eventExecAddedVar)
{ /*lint --e{715}*/
assert(eventhdlr != NULL);
assert(strcmp(SCIPeventhdlrGetName(eventhdlr), EVENTHDLR_NAME) == 0);
assert(event != NULL);
assert(SCIPeventGetType(event) == SCIP_EVENTTYPE_VARADDED);
SCIPdebugMessage("exec method of event handler for added variable to probdata\n");
/* add new variable to probdata */
SCIP_CALL( SCIPprobdataAddVar(scip, SCIPgetProbData(scip), SCIPeventGetVar(event)) );
return SCIP_OKAY;
}
/* exec the event handler
*
* we interrupt the solution process
*/
static
SCIP_DECL_EVENTEXEC(eventExecZeroobj)
{
SCIP_HEURDATA* heurdata;
assert(eventhdlr != NULL);
assert(eventdata != NULL);
assert(strcmp(SCIPeventhdlrGetName(eventhdlr), EVENTHDLR_NAME) == 0);
assert(event != NULL);
assert(SCIPeventGetType(event) & SCIP_EVENTTYPE_NODESOLVED);
heurdata = (SCIP_HEURDATA*)eventdata;
assert(heurdata != NULL);
/* interrupt solution process of sub-SCIP */
if( SCIPgetLPSolstat(scip) == SCIP_LPSOLSTAT_ITERLIMIT || SCIPgetNLPIterations(scip) >= heurdata->maxlpiters )
{
SCIP_CALL( SCIPinterruptSolve(scip) );
}
return SCIP_OKAY;
}
/** execution method of event handler */
static
SCIP_DECL_EVENTEXEC(eventExecProbdatavardeleted)
{
SCIP_VAR* var;
SCIP_PROBDATA* probdata;
int idx;
assert(SCIPeventGetType(event) == SCIP_EVENTTYPE_VARDELETED);
var = SCIPeventGetVar(event);
probdata = (SCIP_PROBDATA*) eventdata;
assert(probdata != NULL);
assert(var != NULL);
/* get index of variable in stablesets array */
idx = (int)(size_t) SCIPvarGetData(var);
SCIPdebugMessage("remove variable %s [%d] from list of stable sets\n", SCIPvarGetName(var), idx);
assert(probdata->stablesetvars[idx] == var);
/* remove variable from stablesets array and release it */
SCIPfreeBlockMemoryArray(scip, &(probdata->stablesets[idx]), probdata->stablesetlengths[idx]); /*lint !e866*/
SCIP_CALL( SCIPreleaseVar(scip, &(probdata->stablesetvars[idx])) );
/* move all subsequent variables to the front */
for( ; idx < probdata->nstablesets - 1; idx++)
{
probdata->stablesets[idx] = probdata->stablesets[idx + 1];
probdata->stablesetlengths[idx] = probdata->stablesetlengths[idx + 1];
probdata->stablesetvars[idx] = probdata->stablesetvars[idx + 1];
SCIPvarSetData(probdata->stablesetvars[idx], (SCIP_VARDATA*) (size_t) idx); /*lint !e571*/
}
probdata->nstablesets--;
return SCIP_OKAY;
}/*lint !e715*/
/* exec the event handler
*
* we interrupt the solution process
*/
static
SCIP_DECL_EVENTEXEC(eventExecLocalbranching)
{
SCIP_HEURDATA* heurdata;
assert(eventhdlr != NULL);
assert(eventdata != NULL);
assert(strcmp(SCIPeventhdlrGetName(eventhdlr), EVENTHDLR_NAME) == 0);
assert(event != NULL);
assert(SCIPeventGetType(event) & SCIP_EVENTTYPE_LPSOLVED);
heurdata = (SCIP_HEURDATA*)eventdata;
assert(heurdata != NULL);
/* interrupt solution process of sub-SCIP */
if( SCIPgetNLPs(scip) > heurdata->lplimfac * heurdata->nodelimit )
{
SCIPdebugMessage("interrupt after %"SCIP_LONGINT_FORMAT" LPs\n",SCIPgetNLPs(scip));
SCIP_CALL( SCIPinterruptSolve(scip) );
}
return SCIP_OKAY;
}
/** execution method of event handler */
static
SCIP_DECL_EVENTEXEC(eventExecBestsol)
{ /*lint --e{715}*/
SCIP_SOL* bestsol;
SCIP_Real solvalue;
assert(eventhdlr != NULL);
assert(strcmp(SCIPeventhdlrGetName(eventhdlr), EVENTHDLR_NAME) == 0);
assert(event != NULL);
assert(scip != NULL);
assert(SCIPeventGetType(event) == SCIP_EVENTTYPE_BESTSOLFOUND);
SCIPdebugMessage("exec method of event handler for best solution found\n");
bestsol = SCIPgetBestSol(scip);
assert(bestsol != NULL);
solvalue = SCIPgetSolOrigObj(scip, bestsol);
/* print best solution value */
SCIPinfoMessage(scip, NULL, "found new best solution with solution value <%g> in SCIP <%s>\n",
solvalue, SCIPgetProbName(scip) );
return SCIP_OKAY;
}
/** execution method of event handler */
static
SCIP_DECL_EVENTEXEC(eventExecBoundwriting)
{ /*lint --e{715}*/
SCIP_EVENTHDLRDATA* eventhdlrdata;
assert(scip != NULL);
assert(eventhdlr != NULL);
assert(strcmp(SCIPeventhdlrGetName(eventhdlr), EVENTHDLR_NAME) == 0);
assert(event != NULL);
assert(((SCIPeventGetType(event) & SCIP_EVENTTYPE_NODESOLVED) == SCIP_EVENTTYPE_NODEFEASIBLE) || ((SCIPeventGetType(event) & SCIP_EVENTTYPE_NODESOLVED) == SCIP_EVENTTYPE_NODEINFEASIBLE) || ((SCIPeventGetType(event) & SCIP_EVENTTYPE_NODESOLVED) == SCIP_EVENTTYPE_NODEBRANCHED));
SCIPdebugMessage("exec method of event handler for writing primal- and dualbounds\n");
eventhdlrdata = SCIPeventhdlrGetData(eventhdlr);
assert(eventhdlrdata != NULL);
#ifdef ONEFILE
/* check if we need to open the file */
if( strlen(eventhdlrdata->filename) > 0 && !eventhdlrdata->isopen )
{
assert(eventhdlrdata->file == NULL);
assert(eventhdlrdata->oldfilename[0] == '\0');
eventhdlrdata->file = fopen(eventhdlrdata->filename, "w");
(void)strncpy(eventhdlrdata->oldfilename, eventhdlrdata->filename, SCIP_MAXSTRLEN);
if( eventhdlrdata->file == NULL )
{
SCIPerrorMessage("cannot create file <%s> for writing\n", eventhdlrdata->filename);
SCIPprintSysError(eventhdlrdata->filename);
return SCIP_FILECREATEERROR;
}
eventhdlrdata->isopen = TRUE;
#ifdef LONGSTATS
SCIPinfoMessage(scip, eventhdlrdata->file, "Problem: %s (%d Original Constraints, %d Original Variables)\n", SCIPgetProbName(scip), SCIPgetNOrigConss(scip), SCIPgetNOrigVars(scip) );
SCIPinfoMessage(scip, eventhdlrdata->file, "\t (%d Presolved Constraints, %d Presolved Variables, (%d binary, %d integer, %d implicit integer, %d continuous))\n", SCIPgetNConss(scip), SCIPgetNVars(scip), SCIPgetNBinVars(scip), SCIPgetNIntVars(scip), SCIPgetNImplVars(scip), SCIPgetNContVars(scip));
SCIPinfoMessage(scip, eventhdlrdata->file, "\n");
#endif
}
#endif
/* call writing only at right moment */
if( eventhdlrdata->freq == 0 || (SCIPgetNNodes(scip) % eventhdlrdata->freq) != 0 )
return SCIP_OKAY;
#ifndef ONEFILE
if( strlen(eventhdlrdata->filename) > 0 )
{
char name[SCIP_MAXSTRLEN];
char number[SCIP_MAXSTRLEN];
char* pch;
int n;
assert(eventhdlrdata->file == NULL);
assert(!eventhdlrdata->isopen);
if( eventhdlrdata->oldfilename[0] == '\0' )
(void)strncpy(eventhdlrdata->oldfilename, eventhdlrdata->filename, SCIP_MAXSTRLEN);
/* find last '.' to append filenumber */
pch=strrchr(eventhdlrdata->filename,'.');
assert(eventhdlrdata->filenumber > 0);
n=sprintf(number, "%"SCIP_LONGINT_FORMAT"", eventhdlrdata->filenumber * eventhdlrdata->freq);
assert(n > 0);
assert(n < SCIP_MAXSTRLEN);
/* if no point is found, extend directly */
if( pch == NULL )
{
(void)strncpy(name, eventhdlrdata->filename, (unsigned int)(SCIP_MAXSTRLEN - n));
strncat(name, number, (unsigned int)n);
}
else
{
int len;
if( (pch-(eventhdlrdata->filename)) > (SCIP_MAXSTRLEN - n) ) /*lint !e776*/
len = SCIP_MAXSTRLEN - n;
else
len = (int) (pch-(eventhdlrdata->filename));
(void)strncpy(name, eventhdlrdata->filename, (unsigned int)len);
name[len] = '\0';
strncat(name, number, (unsigned int)n);
assert(len+n < SCIP_MAXSTRLEN);
name[len+n] = '\0';
if( len + n + strlen(&(eventhdlrdata->filename[len])) < SCIP_MAXSTRLEN ) /*lint !e776*/
{
strncat(name, &(eventhdlrdata->filename[len]), strlen(&(eventhdlrdata->filename[len])));
name[strlen(eventhdlrdata->filename)+n] = '\0';
}
}
eventhdlrdata->file = fopen(name, "w");
if( eventhdlrdata->file == NULL )
{
SCIPerrorMessage("cannot create file <%s> for writing\n", eventhdlrdata->filename);
SCIPprintSysError(eventhdlrdata->filename);
return SCIP_FILECREATEERROR;
}
eventhdlrdata->isopen = TRUE;
#ifdef LONGSTATS
SCIPinfoMessage(scip, eventhdlrdata->file, "Problem: %s (%d Original Constraints, %d Original Variables)\n", SCIPgetProbName(scip), SCIPgetNOrigConss(scip), SCIPgetNOrigVars(scip) );
SCIPinfoMessage(scip, eventhdlrdata->file, "\t (%d Active Constraints, %d Active Variables, (%d binary, %d integer, %d implicit integer, %d continuous))\n", SCIPgetNConss(scip), SCIPgetNVars(scip), SCIPgetNBinVars(scip), SCIPgetNIntVars(scip), SCIPgetNImplVars(scip), SCIPgetNContVars(scip));
SCIPinfoMessage(scip, eventhdlrdata->file, "\n");
#endif
}
#endif
#ifndef NDEBUG
/* check the filename did not change during the solving */
if( strlen(eventhdlrdata->filename) > 0 && eventhdlrdata->isopen )
{
char tmp[SCIP_MAXSTRLEN];
(void)strncpy(tmp, eventhdlrdata->filename, SCIP_MAXSTRLEN);
/* the name should stay the same */
assert(strcmp(tmp, eventhdlrdata->oldfilename) == 0);
}
#endif
#ifdef FOCUSNODE
/* call writing method */
SCIP_CALL( writeBoundsFocusNode(scip, eventhdlrdata) );
#else
/* call writing method */
SCIP_CALL( writeBounds(scip, eventhdlrdata->file, eventhdlrdata->writesubmipdualbound) );
#endif
#ifndef ONEFILE
if( strlen(eventhdlrdata->filename) > 0 )
{
assert(eventhdlrdata->isopen);
(void) fclose(eventhdlrdata->file);
eventhdlrdata->isopen = FALSE;
eventhdlrdata->file = NULL;
++(eventhdlrdata->filenumber);
}
#endif
return SCIP_OKAY;
}
