45 #include "blockmemshell/memory.h"
47 #include "scip/pub_misc.h"
53 #define BMS_CALL(x) do \
57 SCIPerrorMessage("No memory in function call.\n"); \
58 return SCIP_NOMEMORY; \
64 struct SCIP_SDPiSolver
66 SCIP_MESSAGEHDLR* messagehdlr;
72 SCIP_Real sdpsolverfeastol;
87 SCIPerrorMessage(
"No SDP-solver available (SDP=none).\n");
88 SCIPerrorMessage(
"Ensure <relaxing/SDP/freq = -1>.\n");
98 SCIPerrorMessage(
"No SDP-solver available (SDP=none).\n");
99 SCIPerrorMessage(
"Ensure <relaxing/SDP/freq = -1>.\n");
123 return "no SDP-Solver linked currently";
136 assert( sdpisolver != NULL );
185 SCIP_MESSAGEHDLR* messagehdlr,
190 assert( sdpisolver != NULL );
191 assert( blkmem != NULL );
192 assert( bufmem != NULL );
193 SCIPdebugMessage(
"Calling SCIPsdpiCreate \n");
194 SCIPdebugMessage(
"Note that currently no SDP-Solver is linked to the binary. Ensure <relaxing/SDP/freq = -1>. \n");
196 BMS_CALL( BMSallocBlockMemory(blkmem, sdpisolver) );
198 (*sdpisolver)->messagehdlr = messagehdlr;
199 (*sdpisolver)->blkmem = blkmem;
200 (*sdpisolver)->bufmem = bufmem;
210 assert( sdpisolver != NULL );
211 assert( *sdpisolver != NULL );
212 SCIPdebugMessage(
"Freeing SDPISolver\n");
214 BMSfreeBlockMemory((*sdpisolver)->blkmem, sdpisolver);
224 SCIPdebugMessage(
"SDPs aren't counted as there is no SDP-solver.\n");
234 SCIPdebugMessage(
"SDPs aren't counted as there is no SDP-solver.\n");
274 int* sdpconstnblocknonz,
278 SCIP_Real** sdpconstval,
280 int** sdpnblockvarnonz,
290 int* blockindchanges,
302 int* startZnblocknonz,
308 SCIP_Real** startZval,
310 int* startXnblocknonz,
316 SCIP_Real** startXval,
352 SCIP_Real penaltyparam,
363 int* sdpconstnblocknonz,
367 SCIP_Real** sdpconstval,
369 int** sdpnblockvarnonz,
379 int* blockindchanges,
391 int* startZnblocknonz,
397 SCIP_Real** startZval,
399 int* startXnblocknonz,
405 SCIP_Real** startXval,
412 SCIP_Bool* penaltybound
460 SCIP_Bool* primalfeasible,
461 SCIP_Bool* dualfeasible
466 return SCIP_PLUGINNOTFOUND;
623 return SCIP_PLUGINNOTFOUND;
634 return SCIP_PLUGINNOTFOUND;
651 return SCIP_PLUGINNOTFOUND;
659 int* startXnblocknonz
665 return SCIP_PLUGINNOTFOUND;
682 int* startXnblocknonz,
686 SCIP_Real** startXval
691 return SCIP_PLUGINNOTFOUND;
711 return SCIP_PLUGINNOTFOUND;
718 int* startXnblocknonz
723 return SCIP_PLUGINNOTFOUND;
733 int* startXnblocknonz,
737 SCIP_Real** startXval
742 return SCIP_PLUGINNOTFOUND;
752 return SCIP_PLUGINNOTFOUND;
763 return SCIP_PLUGINNOTFOUND;
774 return SCIP_PLUGINNOTFOUND;
785 return SCIP_PLUGINNOTFOUND;
824 assert( sdpisolver != NULL );
825 assert( dval != NULL );
830 *dval = sdpisolver->epsilon;
833 *dval = sdpisolver->gaptol;
836 *dval = sdpisolver->feastol;
839 *dval = sdpisolver->sdpsolverfeastol;
842 *dval = sdpisolver->objlimit;
845 return SCIP_PARAMETERUNKNOWN;
858 assert( sdpisolver != NULL );
863 sdpisolver->epsilon = dval;
864 SCIPdebugMessage(
"Setting sdpisolver epsilon to %f.\n", dval);
867 sdpisolver->gaptol = dval;
868 SCIPdebugMessage(
"Setting sdpisolver gaptol to %f.\n", dval);
871 sdpisolver->feastol = dval;
872 SCIPdebugMessage(
"Setting sdpisolver feastol to %f.\n", dval);
875 sdpisolver->sdpsolverfeastol = dval;
876 SCIPdebugMessage(
"Setting sdpisolver sdpsolverfeastol to %f.\n", dval);
879 SCIPdebugMessage(
"Setting sdpisolver objlimit to %f.\n", dval);
880 sdpisolver->objlimit = dval;
883 return SCIP_PARAMETERUNKNOWN;
896 assert( sdpisolver != NULL );
901 *ival = (int) sdpisolver->sdpinfo;
902 SCIPdebugMessage(
"Getting sdpisolver information output (%d).\n", *ival);
905 return SCIP_PARAMETERUNKNOWN;
918 assert( sdpisolver != NULL );
923 sdpisolver->sdpinfo = (SCIP_Bool) ival;
924 SCIPdebugMessage(
"Setting sdpisolver information output (%d).\n", ival);
927 return SCIP_PARAMETERUNKNOWN;
939 SCIPdebugMessage(
"Lambdastar parameter only used by SDPA.");
948 SCIP_Real* penaltyparam
951 assert( penaltyparam != NULL );
953 *penaltyparam = 1E+10;
961 SCIP_Real penaltyparam,
962 SCIP_Real* maxpenaltyparam
965 assert( maxpenaltyparam != NULL );
967 *maxpenaltyparam = 1E+10;
992 return SCIP_PLUGINNOTFOUND;
1003 return SCIP_PLUGINNOTFOUND;
SCIP_Bool SCIPsdpiSolverDoesWarmstartNeedPrimal(void)
SCIP_RETCODE SCIPsdpiSolverIgnoreInstability(SCIP_SDPISOLVER *sdpisolver, SCIP_Bool *success)
SCIP_Bool SCIPsdpiSolverIsOptimal(SCIP_SDPISOLVER *sdpisolver)
SCIP_Real SCIPsdpiSolverGetMaxPrimalEntry(SCIP_SDPISOLVER *sdpisolver)
SCIP_RETCODE SCIPsdpiSolverGetPrimalBoundVars(SCIP_SDPISOLVER *sdpisolver, SCIP_Real *lbvars, SCIP_Real *ubvars, int *arraylength)
SCIP_Bool SCIPsdpiSolverIsAcceptable(SCIP_SDPISOLVER *sdpisolver)
SCIP_Bool SCIPsdpiSolverIsDualFeasible(SCIP_SDPISOLVER *sdpisolver)
SCIP_RETCODE SCIPsdpiSolverGetPrimalNonzeros(SCIP_SDPISOLVER *sdpisolver, int nblocks, int *startXnblocknonz)
const char * SCIPsdpiSolverGetSolverDesc(void)
SCIP_RETCODE SCIPsdpiSolverGetPrimalMatrix(SCIP_SDPISOLVER *sdpisolver, int nblocks, int *startXnblocknonz, int **startXrow, int **startXcol, SCIP_Real **startXval)
SCIP_Bool SCIPsdpiSolverIsPrimalInfeasible(SCIP_SDPISOLVER *sdpisolver)
enum SCIP_SDPSolverSetting SCIP_SDPSOLVERSETTING
SCIP_RETCODE SCIPsdpiSolverReadSDP(SCIP_SDPISOLVER *sdpisolver, const char *fname)
SCIP_Bool SCIPsdpiSolverFeasibilityKnown(SCIP_SDPISOLVER *sdpisolver)
SCIP_Real SCIPsdpiSolverInfinity(SCIP_SDPISOLVER *sdpisolver)
SCIP_Bool SCIPsdpiSolverIsPrimalUnbounded(SCIP_SDPISOLVER *sdpisolver)
int SCIPsdpiSolverGetInternalStatus(SCIP_SDPISOLVER *sdpisolver)
interface methods for specific SDP-solvers
SCIP_Real SCIPsdpiSolverGetDefaultSdpiSolverFeastol(void)
SCIP_Bool SCIPsdpiSolverIsDualInfeasible(SCIP_SDPISOLVER *sdpisolver)
SCIP_Bool SCIPsdpiSolverIsObjlimExc(SCIP_SDPISOLVER *sdpisolver)
SCIP_Bool SCIPsdpiSolverIsTimelimExc(SCIP_SDPISOLVER *sdpisolver)
SCIP_RETCODE SCIPsdpiSolverResetCounter(SCIP_SDPISOLVER *sdpisolver)
const char * SCIPsdpiSolverGetSolverName(void)
SCIP_RETCODE SCIPsdpiSolverLoadAndSolve(SCIP_SDPISOLVER *sdpisolver, int nvars, SCIP_Real *obj, SCIP_Real *lb, SCIP_Real *ub, int nsdpblocks, int *sdpblocksizes, int *sdpnblockvars, int sdpconstnnonz, int *sdpconstnblocknonz, int **sdpconstrow, int **sdpconstcol, SCIP_Real **sdpconstval, int sdpnnonz, int **sdpnblockvarnonz, int **sdpvar, int ***sdprow, int ***sdpcol, SCIP_Real ***sdpval, int **indchanges, int *nremovedinds, int *blockindchanges, int nremovedblocks, int nlpcons, int noldlpcons, SCIP_Real *lplhs, SCIP_Real *lprhs, int *rownactivevars, int lpnnonz, int *lprow, int *lpcol, SCIP_Real *lpval, SCIP_Real *starty, int *startZnblocknonz, int **startZrow, int **startZcol, SCIP_Real **startZval, int *startXnblocknonz, int **startXrow, int **startXcol, SCIP_Real **startXval, SCIP_SDPSOLVERSETTING startsettings, SCIP_Real timelimit)
int SCIPsdpiSolverGetDefaultSdpiSolverNpenaltyIncreases(void)
SCIP_RETCODE SCIPsdpiSolverComputeMaxPenaltyparam(SCIP_SDPISOLVER *sdpisolver, SCIP_Real penaltyparam, SCIP_Real *maxpenaltyparam)
SCIP_RETCODE SCIPsdpiSolverFree(SCIP_SDPISOLVER **sdpisolver)
SCIP_RETCODE SCIPsdpiSolverGetSolFeasibility(SCIP_SDPISOLVER *sdpisolver, SCIP_Bool *primalfeasible, SCIP_Bool *dualfeasible)
SCIP_RETCODE SCIPsdpiSolverGetIterations(SCIP_SDPISOLVER *sdpisolver, int *iterations)
SCIP_RETCODE SCIPsdpiSolverGetRealpar(SCIP_SDPISOLVER *sdpisolver, SCIP_SDPPARAM type, SCIP_Real *dval)
SCIP_RETCODE SCIPsdpiSolverSettingsUsed(SCIP_SDPISOLVER *sdpisolver, SCIP_SDPSOLVERSETTING *usedsetting)
SCIP_RETCODE SCIPsdpiSolverGetObjval(SCIP_SDPISOLVER *sdpisolver, SCIP_Real *objval)
SCIP_RETCODE SCIPsdpiSolverGetSdpCalls(SCIP_SDPISOLVER *sdpisolver, int *calls)
void * SCIPsdpiSolverGetSolverPointer(SCIP_SDPISOLVER *sdpisolver)
SCIP_RETCODE SCIPsdpiSolverLoadAndSolveWithPenalty(SCIP_SDPISOLVER *sdpisolver, SCIP_Real penaltyparam, SCIP_Bool withobj, SCIP_Bool rbound, int nvars, SCIP_Real *obj, SCIP_Real *lb, SCIP_Real *ub, int nsdpblocks, int *sdpblocksizes, int *sdpnblockvars, int sdpconstnnonz, int *sdpconstnblocknonz, int **sdpconstrow, int **sdpconstcol, SCIP_Real **sdpconstval, int sdpnnonz, int **sdpnblockvarnonz, int **sdpvar, int ***sdprow, int ***sdpcol, SCIP_Real ***sdpval, int **indchanges, int *nremovedinds, int *blockindchanges, int nremovedblocks, int nlpcons, int noldlpcons, SCIP_Real *lplhs, SCIP_Real *lprhs, int *rownactivevars, int lpnnonz, int *lprow, int *lpcol, SCIP_Real *lpval, SCIP_Real *starty, int *startZnblocknonz, int **startZrow, int **startZcol, SCIP_Real **startZval, int *startXnblocknonz, int **startXrow, int **startXcol, SCIP_Real **startXval, SCIP_SDPSOLVERSETTING startsettings, SCIP_Real timelimit, SCIP_Bool *feasorig, SCIP_Bool *penaltybound)
SCIP_Bool SCIPsdpiSolverIsPrimalFeasible(SCIP_SDPISOLVER *sdpisolver)
SCIP_RETCODE SCIPsdpiSolverSetIntpar(SCIP_SDPISOLVER *sdpisolver, SCIP_SDPPARAM type, int ival)
SCIP_RETCODE SCIPsdpiSolverCreate(SCIP_SDPISOLVER **sdpisolver, SCIP_MESSAGEHDLR *messagehdlr, BMS_BLKMEM *blkmem, BMS_BUFMEM *bufmem)
SCIP_Real SCIPsdpiSolverGetDefaultSdpiSolverGaptol(void)
SCIP_RETCODE SCIPsdpiSolverSetRealpar(SCIP_SDPISOLVER *sdpisolver, SCIP_SDPPARAM type, SCIP_Real dval)
SCIP_RETCODE SCIPsdpiSolverWriteSDP(SCIP_SDPISOLVER *sdpisolver, const char *fname)
SCIP_Bool SCIPsdpiSolverIsInfinity(SCIP_SDPISOLVER *sdpisolver, SCIP_Real val)
SCIP_RETCODE SCIPsdpiSolverComputeLambdastar(SCIP_SDPISOLVER *sdpisolver, SCIP_Real maxguess)
SCIP_RETCODE SCIPsdpiSolverGetPreoptimalPrimalNonzeros(SCIP_SDPISOLVER *sdpisolver, int nblocks, int *startXnblocknonz)
SCIP_RETCODE SCIPsdpiSolverGetIntpar(SCIP_SDPISOLVER *sdpisolver, SCIP_SDPPARAM type, int *ival)
SCIP_Bool SCIPsdpiSolverWasSolved(SCIP_SDPISOLVER *sdpisolver)
SCIP_RETCODE SCIPsdpiSolverIncreaseCounter(SCIP_SDPISOLVER *sdpisolver)
SCIP_Bool SCIPsdpiSolverIsDualUnbounded(SCIP_SDPISOLVER *sdpisolver)
struct SCIP_SDPiSolver SCIP_SDPISOLVER
SCIP_Bool SCIPsdpiSolverIsConverged(SCIP_SDPISOLVER *sdpisolver)
static void errorMessageAbort(void)
SCIP_RETCODE SCIPsdpiSolverGetSol(SCIP_SDPISOLVER *sdpisolver, SCIP_Real *objval, SCIP_Real *dualsol, int *dualsollength)
SCIP_RETCODE SCIPsdpiSolverGetPreoptimalSol(SCIP_SDPISOLVER *sdpisolver, SCIP_Bool *success, SCIP_Real *dualsol, int *dualsollength, int nblocks, int *startXnblocknonz, int **startXrow, int **startXcol, SCIP_Real **startXval)
enum SCIP_SDPParam SCIP_SDPPARAM
static void errorMessage(void)
SCIP_RETCODE SCIPsdpiSolverComputePenaltyparam(SCIP_SDPISOLVER *sdpisolver, SCIP_Real maxcoeff, SCIP_Real *penaltyparam)
SCIP_Bool SCIPsdpiSolverIsIterlimExc(SCIP_SDPISOLVER *sdpisolver)