9- avoid hashmap key error in removal of doubletons and singletons in dual presolve of setppc constraints by skipping tripleton locks
10- when upgrading to xor constraint require parity variable to be enforced integral
11- provide the depth of the applied effective root with SCIPtreeGetEffectiveRootDepth() to avert invalid globalization of constraint set changes
12- disable objective limit during lexicographic dual minimization to preserve feasibility status
13- declare infinite bound inference infeasible in SCIPinferVarLbCons() and SCIPinferVarUbCons() to maintain correct loose LP values
14- fixed update of consssorted flags in a variable expressions data when adding an additional constraint using this variable
15- corrected computation of number of variables affected by symmetry
16- fixed computation of symmetry group size
17- correct comparison conditions in checkRedundancySide() of cons_varbound.c to ensure pair redundancy
18- also for the initial LP solve force acceptance of Farkas row in SCIPlpGetDualfarkas() if all other means for obtaining a reliable Farkas proof have failed
19- avoid overwriting the cutoff flag in tightenCoefs(), preprocessConstraintPairs(), and applyFixings() of cons_varbound.c to reject infeasible solutions
20- fix memory leak in CppAD interface w.r.t. user expressions
21- adjust bound before checking feasibility in SCIPvarAddVlb() and SCIPvarAddVub() to detect integrality cutoff
22- choose incumbent solution as reference in varFixingsRins() and varFixingsDins() of ALNS heuristic to ensure integrality
23- consider variable type in selectInitialVariableRandomly() of GINS heuristic to count integral variables
24- correct fractionality sign when standardizing bounds in simplifyInequalities() of cons_linear.c
25
26Build system
27------------
28
29### CMake
30
31- disable setting the boolean variable MACOSX_RPATH to a string value to actually enable @rpath usage
32- define SOPLEX_NEEDED and ZIMPL_NEEDED to avoid syntax error in configuration file
33
34Miscellaneous
35-------------
36
37- when writing a problem with non-generic names, warnings are printed if variable or constraint names are not unique
38
39@section RN922 SCIP 9.2.2
40*************************
41
42Fixed bugs
43----------
44
45- fixed definition of slack variables of >= constraints and ranged rows in LP interface to Gurobi
46- free memory leaking singleton arrays for ignored variables in SCIPapplyHeurDualval()
47- use epsilon tolerance for bound implications of dual substitutions in dualPresolve() of cons_linear
48- fixed timeouts due to race conditions when using TPI=omp
49- fixed issues with indexing of concurrent solvers when a solver finishes (because SCIP is too fast) before others have been created when using TPI=tny
50- respect fixing if variable is declared binary in bounds section of mps file
51- fixed reading mps files with short row and column names and some fixed-form mps files with spaces in row names
52- force local bounds to adjusted value when fixing variable
53- adds missing NULL check for Benders' subproblems in feasalt cut
54- apply objective scale to pseudoboolean non-linear terms
55- impose pseudoboolean cost bound without intermediate casting
56- flatten aggregation graph before finishing presolve to accelerate evaluating an empty solution
57- respect fixed resultants in binvarGetActiveProbindex() of branch_relpscost
58- drop fixed resultants when copying pseudoboolean constraint
59- correct matrix entry when resolving aggregated variables in checkSystemGF2() of cons_xor
60- corrected symmetry detection graph of cons_and, cons_or, cons_xor, and cons_indicator
61- correct redundancy filter in consdataCreateRedundant() of cons_bounddisjunction
62- replaced unreliable comparison of node pointers by node numbers in prop_genvbounds
63- allocate memory for all integral variables in applyCliqueFixings() of heur_clique
64- fix conflict resolution for one particular case in cons_indicator
65- make interface to nauty thread safe
66- use relative epsilon tolerance to check Farkas row in SCIPlpGetDualfarkas() to avoid invalid cutoff
67- fixed bug with SOS2 constraints of size that avoided a correct enforcing
68
69Build system
70------------
71
72### Cmake
73
74- allow to set HIGHS_DIR via -D flag as well
75- reenable building and testing of SCFLP example
76- fixed that compiler flags for sanitizers were not used for compiling objects of SCIP library
77- prioritize local includes for unittests, examples, and applications to link the intended SCIP version
78- fix detection of Gurobi 11 and 12 libs on Linux
79
80### Makefile
81
82- fix use of SANITIZE=full for debug builds with GCC on Linux: removed -fsanitize=leak because it cannot be used together with -fsanitize=address
83- fix that make install SHARED=true missed to install libscipbase
84
85@section RN921 SCIP 9.2.1
86*************************
87
88Fixed bugs
89----------
90
91- apply same feasibility definition in PaPILO as in SCIP to avoid numerical inconsistencies
92- check variable cancellation in SCIPvarAddVlb() and SCIPvarAddVub() to avert wrong infeasibility
93- SCIPfreeReoptSolve() now also clears partial solutions
94- fixed bug in calculation of "fractionality score" for spatial branching candidates in cons_nonlinear
95- fixed bug with installing symmetry dialog in sub-SCIPs
96- added copy callbacks for presolvers dualagg and redvub, branching rules lookahead and cloud, primal heuristics dualval and repair, propagator nlobbt, and separators gauge and convexproj
97- respect maximal separation rounds setting by running remaining delayed separators in last round
98- replace terms consistently in applyFixings() of cons_linear to handle infinite aggregations
99
100Build system
101------------
102
103### Cmake
104
105- fix SCIP Optimization Suite build with -DSYM=(s)bliss and bliss being downloaded and built during configuration by cmake
106- fix use of a bliss cmake installation (from https://github.com/scipopt/bliss) specified via BLISS_DIR if -DSYM=(s)bliss
107
108### Makefile
109
110- the check/solchecker build can now use also GMP installations by Homebrew on macOS/arm64
111
112Miscellaneous
113-------------
114
115- the output precision for writing CIP/OPB files has been increased at several places for nonlinear constraints
116- the solchecker tool now also supports SCIP solutions without unknown or infinite values
117- generate an error if the same linear constraint is used by several indicator constraints
118
119@section RN920 SCIP 9.2.0
120*************************
121
122Features
123--------
124
125- added check for absolute and relative gap limits at end of synchronization in concurrent solving mode, in order to terminate earlier;
126 note that if the concurrent solve is stopped due to a gap limit, the "winner" solver will have an interrupted solution status and
127 its primal and dual bounds may not be the best possible ones (use SCIPgetConcurrentPrimalbound() and SCIPgetConcurrentDualbound() instead)
128- parse pseudoboolean constraint from CIP format (and add linear-"and"-reformulation)
129
130Performance improvements
131------------------------
132
133- reoptimization now also stores propagations from propagators if reoptimization/saveconsprop is enabled;
134 the parameter will be renamed to reoptimization/saveprop in a next major release
135- imposed stricter limits on the size of disconnected components which may be solved separately during presolve
136- use individual slack variables also for constraints indicated by a common binary variable to use tighter formulation by default
137- when computing symmetries using Nauty, iteration limits are now available to terminate Nauty early
138
139Fixed bugs
140----------
141
142- Benders' decomposition subproblems that are always infeasible are correctly handled and the complete problem is declared infeasible
143- skip linear constraint propagation if the residual activity bound cancels the side precision
144- correct bound tracking to make the evaluation of primal-dual-integrals work
145- skip aggregations on fixed variables in milp presolver to avoid errors for earlier versions of PaPILO
146- use indices of negation counterparts and accept fixings when ordering and-resultants of pseudoboolean constraints
147- update locks on model variables before removing pseudoboolean constraints in presolve
148- reformulate soft pseudoboolean constraints with linear constraints if the indicator decomposition is disabled
149- add workaround for recent HiGHS versions resetting the model status when changing the presolve option after a solve
150- avoid hashmap key error in removal of doubletons and singletons in dual presolve of setppc constraints by updating constraints and corresponding hashmaps after each multi-aggregation
151
152Interface changes
153-----------------
154
155### New API functions
156
157- added SCIPtpiIsAvailable() to check whether a working task processing interface is available (TPI != none)
158- added SCIPtpiGetLibraryName() and SCIPtpiGetLibraryDesc()
159- SCIPdelCons() can now also be called in SCIP_STAGE_TRANSFORMED
160- added SCIPstrcasecmp() and SCIPstrncasecmp() for case-insensitive string comparison
161- added SCIPbendersSubproblemsAreInfeasible() to return if at least one subproblem has been identified as being
162 infeasible prior to performing any variable fixing
163
164### New parameters
165
166- presolving/milp/abortfacexhaustive to control the abort threshold for exhaustive presolving in PAPILO
167- presolving/milp/abortfacmedium to control the abort threshold for medium presolving in PAPILO
168- presolving/milp/abortfacfast to control the abort threshold for fast presolving in PAPILO
169- constraints/components/maxcompweight to determine the maximum weight for a disconnected component that is solved during presolve
170- constraints/components/contfactor counts the contributing factor of a single continuous variables with respect to the weight limit specified by constraints/components/maxcompweight
171- constraints/indicator/usesameslackvar to decide whether the same slack variable should be used for indicators constraints with common binary variable
172- propagating/symmetry/nautymaxncells and propagating/symmetry/nautymaxnnodes to set iteration limits in Nauty (only available if build with SYM=nauty or SYM=snauty)
173
174### Changed parameters
175
176- presolving/milp/threads is now only available if PaPILO is built with TBB
177- changed default of numerics/recomputefac to 1e+6 to aim at relative epsilon precision
178
179Build system
180------------
181
182### Cmake
183
184- attempted to fix detection of CPLEX library on macOS and Windows systems
185
186Testing
187-------
188
189- added parameter FILTER for tests/Makefile to run only tests with a specific pattern (ctest with -R FILTER)
190
191Miscellaneous
192-------------
193
194- adjusted Gurobi interface for Gurobi 12
195- reordered events: BESTSOLFOUND/NODE_FEASIBLE are now processed before the nodes are cut off and before NODE_DELETE events are processed
196- removed `#define` of `getcwd` (in case of Windows builds) in scip/def.h
197- the `#define` of `strcasecmp` and `strncasecmp` (in case of Windows builds) in scip/def.h will be removed with SCIP 10;
198 use `SCIPstr(n)casecmp()` (scip/pub_misc.h) instead
199
200@section RN911 SCIP 9.1.1
201*************************
202
203Fixed bugs
204----------
205
206- fix bug in separator aggregation that caused incorrect terminations when for an integral objective the objective limit was set to a fractional value
207- fix bug in memory management of doublelex matrices where indices were set beyond array bounds
208- store relevance flags in consdataTightenCoefs() of cons_linear.c to avoid dropping variables unintentionally due to numerical changes
209- fix possible segmentation fault when resetting local and global bounds of original variables after abnormal termination
210- compare coefficient scales to detect parallelism to avoid integral relaxations in detectRedundantConstraints() of cons_linear.c
211- remove coefficient information from linear constraint hashes in order to compare keys with tolerances
212- fixed NULL pointer dereference in SCIPtpiGetNumThreads() of TinyCThread interface if called before SCIPtpiInit(); now returns 0 in this case
213- calling SCIPsolveConcurrent() when SCIP was compiled without a TPI now results in a plugin-not-found error
214- fixed LPI status functions of lpi_cpx if the barrier is called
215- check cuts for redundancy after scaling in cutTightenCoefs() and cutTightenCoefsQuad() of cuts.c
216- maintain definition of artificial integral variables in xor constraints to avoid invalidation of conflict constraints
217- free childless nodes only after the deactivation of the old focus path in treeSwitchPath() of tree.c
218- only normalize rows for column norm evaluations in shift-and-propagate heuristic to avoid unscaled violations
219- expression handler for absolute value computes now also a cut when one of the bounds is infinite as long as the sign of the argument is known
220- check for more data to be freed before modifying nonlinear constraints in SCIPaddExprNonlinear()
221- allow for terms with zero coefficients when checking sum in SCIPexprCheckQuadratic()
222- do not abort writing of constraints to lp/pip files (and forget to free buffer memory) when encountering unsupported nonlinear constraint
223- fixed that strings with names and descriptions of the LPI's were not null-terminated
224- fixed decision on when to enforce nonlinear constraints on relaxation solutions if integer variables are present
225- fixed reading of very long lines (>78000 characters) in FlatZinc reader
226- fixed write of strings with at least 1024 characters when compiling with MSVS 2015 (14.0) or later (still an issue with earlier MSVS, which will not be fixed)
227- keep epsilon coefficients in applyFixings() of cons_linear.c until multiples are merged to avoid loss of relevant contributions
228- fixed that some variable type changes were not accounted correctly in presolve of linear constraints
229- fix constraint sorting in setpp-orbitope upgrade in symmetry code
230- fixed that SCIPnodeGetConsProps() also counted bound changes from propagators
231- use original and-constraints of an untransformed pseudoboolean constraint when copying original problems to avert segmentation faults
232- sort and-constraints of a pseudoboolean constraint with respect to the resultant problem indices to avoid shuffled copies
233- treat pseudoboolean constraints as nonlinear constraints in OPB reader to prevent duplication mistakes
234
235Examples and applications
236-------------------------
237
238- in example Binpacking implement Farkas pricing for packings because branching can actually lead to infeasible restricted master LPs
239
240Testing
241-------
242
243- allow unknown parameters to let tests pass if only optional dependencies are missing
244
245Build system
246------------
247
248### Cmake
249
250- fixed compilation and linking with HiGHS; to use HiGHS, cmake now requires to specify the directory containing the `highs-config.cmake` file via `-DHIGHS_DIR`
251- extended identification of Ipopt library on Windows to also look for ipopt.dll.lib, as this is used since Ipopt 3.13.0 (2019)
252
253
254@section RN910 SCIP 9.1.0
255*************************
256
257Features
258--------
259
260- added constraint handler that checks and enforces bounds on original variables which have been fixed or aggregated in the transformed problem
261- allow for implicit integer variables in xor constraints
262- implemented symmetry detection callbacks for pseudoboolean constraints
263- added option to consider enforcement of nonlinear constraints (including spatial branching) before integrality has been enforced
264
265Fixed bugs
266----------
267
268- make sure that, in dual proof analysis, upgraded constraints can report their number of variables
269- check knapsack constraints as linear constraints to avoid violations based on implicit integrality
270- make sure that symmetry detection callbacks report success
271- make arithmetics in cuts.c more robust to avoid invalid scg cuts due to numerical rounding sensitivity
272- allow to copy pseudoboolean constraints even if no AND constraints are present in order to avoid a warning
273- set value for variable that is introduced to reformulate nonlinear objective function when reading .nl files with initial solution
274- fixed that multi-aggregations from PaPILO were accepted without ensuring that their bounds hold
275
276Examples and applications
277-------------------------
278
279- TSP: removed `gzstream.h/hpp` and parse gzipped files differently
280
281Interface changes
282-----------------
283
284### New callbacks
285
286- SCIP_DECL_NLHDLRENFO takes additional parameter branchcandonly to indicate whether only branching candidates should be registered,
287 without looking for cuts or boundchanges
288
289### New API functions
290
291- new function SCIPgetExprActivityNonlinear() to get activity of nonlinear constraint
292- new function SCIPincludeConshdlrFixedvar() to include new constraint handler that checks bounds on fixed original variables
293
294### New parameters
295
296- limits/primal to terminate the solve as soon as the primal bound is at least as good as this value, same as limits/objectivestop (deprecated)
297- limits/dual to terminate the solve as soon as the dual bound is at least as good as this value
298- presolving/milp/internalmaxrounds to control the maximal rounds for each call of the milp presolver (PaPILO)
299- constraints/fixedvar/∗ to control the new constraint handler that checks bounds on fixed original variables
300- constraints/nonlinear/branching/mixfractional to specify minimal average pseudo cost count for discrete variables
301 at which to start considering spatial branching before branching on fractional integer variables
302 (default infinity means that integrality is always enforced before nonlinear constraints)
303- constraints/nonlinear/branching/fracweight to specify weight of fractionality of integer variables in spatial branching
304
305### New data structures
306
307- enum SCIP_Status has been extended by SCIP_STATUS_PRIMALLIMIT and SCIP_STATUS_DUALLIMIT to indicate termination at limits/primal and limits/dual respectively
308
309Build system
310------------
311
312- changed the default symmetry computation package from sassy+bliss to sassy+nauty
313- bliss source is no longer included in the source tree, but can be linked as external library (https://github.com/ds4dm/Bliss)
314- cmake now automatically downloads and installs the latest release version of SoPlex from GitHub if it is not found directly
315
316Miscellaneous
317-------------
318
319- updated nauty to version 2.8.8
320
321
322@section RN901 SCIP 9.0.1
323*************************
324
325Fixed bugs
326----------
327
328- fixed bug in detection of double lex matrices due to wrong memory allocation
329- fixed bug where parameter `limits/objectivestop` was copied to subscips in plugins which have different objectives
330- in count mode respect empty solution when problem is solved in presolving
331- fixed that SCIP may have tried to pass a negative random seed to an LP solver
332- fixed bug where conflict analysis in strong branching tried to analyze an unsuccessful LP solve
333- enable parameter limits/restarts by allowing to reach the restart limit in restartAllowed()
334- handle normalizing linear constraints in presolving with more caution to avoid numerical difficulties for integer variables with large coefficients
335- use quad precision to calculate residual activity bounds to reduce numerical errors for large variable bounds in cons_linear.c
336- in SCIPsolveProbingLP() when objective has been changed do not return cutoff due to exceeding the cutoff bound
337- set probing LP to be a relaxation only when objective has not been changed
338- declare contradicting infinities in getMinActivity() and getMaxActivity() as non-tight
339- reject farkas solution with large values to bound magnification of errors in SCIPgetFarkasProof()
340- if all variables are fixed, apply relative feasibility tolerance to avoid invalid infeasibility in applyFixings() of cons_linear.c
341- delete empty redundant knapsack and linear constraints after cleanup in SCIPmatrixCreate()
342- cleanup first constraints in SCIPcleanupConssVarbound(), SCIPcleanupConssSetppc(), and SCIPcleanupConssLogicor()
343- symmetry handling methods have not been added immediately after computing symmetries, which has been fixed
344- delay checking node limits after presolving to correctly continue up to the exceeding main node
345- simplified the objective sense assumption for the subproblems in the Lagromory separator
346- update effective root depth iteratively to avoid trapping it at a single child node
347- fixed that dualsparsify presolver did not always check do-not-multiaggregate variable flag
348- fixed bug in computing elements of orbits in orbital reduction
349
350Miscellaneous
351-------------
352
353- redefine isrelax to istight in cons_linear.c
354
355Build system
356------------
357
358### Cmake
359
360- A file cmake.log containing the changed options is written to the cmake build directory.
361- removed replacement of the settings of RelWithDebInfo
362- add compile option -ffp-contract=off or /fp:precise to enhance reproducibility across different systems
363
364### Makefile
365
366- add flags -ffp-contract=off and /fp:precise to enhance reproducibility across different systems
367
368Interface changes
369-----------------
370
371### New parameters
372
373- added parameter "propagating/symmetry/symtiming" to determine when symmetries are computed and handled
374
375### Changed parameters
376
377- Disabled "propagating/symmetry/addconsstiming" since it is not used anymore.
378- Disabled "propagating/symmetry/ofsymcomptiming" since it is not used anymore.
379
380@section RN900 SCIP 9.0.0
381*************************
382
383Features
384--------
385
386- added simplification of abs(x)^n to x^n when n is even (simplify rule POW12)
387- in expression simplification, added possibility to expand power of sum with exponents larger than 2 (simplify rule POW7)
388- in expression simplification, added possibility to expand products with sums as factor when there are more than two factors (simplify rule SP12b)
389- in expression simplification, added possibility to allow expanding power of products with fractional exponent (simplify rule POW5a)
390- can now build with Lapack (used for solving nonlinear problems) explicitly; if IPOPT is used, then Lapack is accessed through it;
391 an interface to Lapack is available in lapack_calls.h
392- added support for logical constraints in binary variables and basic logical operators (and, or, not, equal) to AMPL .nl reader
393- indicator constraints now use the same slack variable for the linear constraints if binary variables are used multiple times
394 for indicator constraints
395- new diving heuristic indicatordiving that focuses on indicator variables which control semicontinuous variables
396- updated statistics on NLP relaxation; information on convexity of NLP has been added
397- added new bandit method (EXP.3-IX) that is the new default of the ALNS heuristic
398- added dialog to display recognized permutation symmetries ("display symmetry")
399- The symmetry detection code has been completely restructured. Instead of completely encoding symmetry information in
400 prop_symmetry.c, there is a new callback for constraint handlers, which returns symmetry relevant information by
401 encoding a constraint as a node and edge colored graph. This way, also custom constraint handlers can provide
402 symmetry information.
403- added special presolving step for logicor and set covering/packing for constraints with two variables and one independent variable
404- increased the depth field of the tree structure from 16 to 30 bits and the max depth of the search tree from 65534 to 1073741822
405- added new cut selector plugin called cutsel_ensemble. Performs a superset of hybrid cutsel. Implements new filtering methods and scoring options.
406- added new branching rule based on Gomory mixed-integer cuts
407- added new branching score to the reliability pseudo-cost branching rule
408- added new cut selector plugin called cutsel_dynamic. Performs a superset of hybrid cutsel. Implements new filtering methods and scoring options.
409- added monoidal strengthening to intersection cut generation of nlhdlr_quadratic
410- extended heuristic dps to run after processing the node and to use strict limits for the reoptimization step
411- extended obbt propagator to handle variables of the linear constraints corresponding to indicator constraints
412- added new heuristic heur_scheduler that models the selection of LNS and diving heuristics as a multi-armed bandit problem
413- added nonlinear handler for signomial expressions: adds cuts based on a difference of convex functions representation
414- Delayed separators are now called if no other cuts have been found. Previously, such separators were only called near stalling or when the separation rounds limit was reached.
415- added a new separator sepa_lagromory to generate Lagromory cuts in a relax-and-cut framework, i.e., GMI cuts generated
416 at different bases that are explored in the process of solving the Lagrangian dual problem at a node. This separator
417 is OFF by default.
418- when parsing nonlinear constraints from CIP files, the * after the number in a term is now optional if followed by a variable,
419 i.e., instead of 3*<x>*<y>, now also 3<x>*<y> can be read, but 3<x><y> is not supported;
420 this allows to read some CIP files that were written with SCIP < 8
421- when running bliss, we now limit the number of actually created generators and not the ones considered by bliss
422
423Performance improvements
424------------------------
425
426- Variable conflict-based addition of Schreier-Sims cuts for symmetry handling in prop_symmetry.c does no longer use
427 an expensive explicit conflict graph construction.
428- Only run propagation of indicator constraints if relevant variables have changed.
429- Propagation of indicator constraints propagates maximal activity of linear constraint to upper bound of slack variable.
430- Added initialization of LP relaxation for SOC constraints.
431- Avoided some numerical issues when separating SOC constraints.
432- If differentation of convex expression fails in convex nlhdlr, then a perturbed point is tried.
433- Use sassy/bliss as default symmetry computation package.
434- No longer linearize indicator constraints in undercover heuristic by default.
435- Deactivate pre-root heuristic zeroobj in undercover heuristic.
436- Added new branching score tiebreaker for reliability pseudo-cost branching rule
437- Try nearest bound if zero solution value is out of bounds in heuristic trivial.
438
439Interface changes
440-----------------
441
442### New and changed callbacks
443
444- Added solution linearization callback for nonlinear handlers.
445 The nonlinear handler may use this callback to add tight cuts that use the given solution as reference point to the cutpool.
446 The callback is called if option constraints/nonlinear/linearizeheursol is enabled.
447- The callback for copying problem data SCIP_DECL_PROBCOPY() now has a parameter "original" indicating whether the original
448 or transformed problem should be copied. Similarly, SCIPprobCopy() in prob.h/c has such a parameter.
449- Added CONSGETPERMSYMGRAPH and CONSGETSIGNEDPERMSYMGRAPH callback for constraint handlers, which provides symmetry information about constraints to SCIP.
450
451### Deleted and changed API methods
452
453- SCIPtpiInitLock(), SCIPtpiDestroyLock(), SCIPtpiInitCondition(), SCIPtpiDestroyCondition() now take a double
454 pointer to SCIP_LOCK/SCIP_CONDITION, respectively.
455- Renamed SCIPtpiSumbitJob() to SCIPtpiSubmitJob().
456- SCIPcreateConsIndicatorLinCons(), SCIPcreateConsIndicatorGenericLinCons() do not allow binvar to be NULL anymore
457- Replaced SCIPnlrowSetCurvature() by SCIPsetNlRowCurvature().
458- SCIPsnprintf() now returns the required string length in case reallocation is necessary.
459- SCIPcomputeCoverUndercover() receives three more arguments which nonlinearities to cover
460- removed definition of SQRT
461- SCIPincludeObjConshdlr(), SCIPincludeConshdlr(), SCIPconshdlrCreate() receive two more arguments for the
462 conshdlr callbacks CONSGETPERMSYMGRAPH and CONSGETSIGNEDPERMSYMGRAPH
464- new parameter issigned for SCIPcomputeOrbitsSym()
465- new parameter symtype for SCIPcomputeComponentsSym()
466- new parameter symtype and permvardomaincenter for SCIPlexicographicReductionAddPermutation()
467- replaced parameters matrixdata and exprdata of SYMcomputeSymmetryGenerators() by graph
468
469### New API functions
470
471- new functions for interfacing with Lapack: SCIPlapackIsAvailable(), SCIPlapackIsAvailable(), SCIPlapackComputeEigenvalues(), and SCIPlapackSolveLinearEquations
472 These functions should be used instead of SCIPisIpoptAvailableIpopt(), SCIPcallLapackDsyevIpopt(), and
473 SCIPsolveLinearEquationsIpopt() in the context of computing eigenvalues or solving linear equations.
474 If IPOPT is available, the calls will be passed on to Ipopt's interface to Lapack.
497- SCIPincludeCutselEnsemble(), SCIPselectCutsEnsemble() to include cutsel_ensemble or use the selection algorithm
498- SCIPincludeBranchruleGomory() includes new gomory cut based branching rule
499- SCIPgetVarAvgGMIScore(), SCIPincVarGMISumScore(), SCIPgetVarLastGMIScore(), SCIPsetVarLastGMIScore() set historical values used by SCIP's default branching rule
500- SCIPincludeNlhdlrSignomial() includes the signomial nonlinear handler
501- SCIPgetExprMonomialData() returns a monomial representation of a product expression
502- SCIPestimateRoot() computes estimators for roots with exponent in [0,1]
503- SCIPincludeCutselDynamic(), SCIPselectCutsDynamic() to include cutsel_dynamic or use the selection algorithm
504- SCIPincludeSepaLagromory() to include the new Lagromory separator
505- SCIPconsCompCheck() to compare two constraints based on their check priority
506- SCIPgetConsPermsymGraph() and SCIPgetConsSignedPermsymGraph() to access the (signed) permutation symmetry detection graph of a constraint
507- SCIPsetConshdlrGetPermsymGraph() and SCIPsetConshdlrGetSignedPermsymGraph() to set the (signed) permutation symmetry detection graph of a constraint
508- SCIPconshdlrSupportsPermsymDetection() and SCIPconshdlrSupportsSignedPermsymDetection() to access whether a conshdlr supports detection of (signed) permutation symmetries
509- new API functions for creating and manipulating a symmetry detection graph: SCIPcreateSymgraph(), SCIPfreeSymgraph(), SCIPcopySymgraph(), SCIPextendPermsymDetectionGraphLinear(),
541 String values should not be quoted and spaces in string values are not supported.
542- sassy can be used now as preprocessor for nauty/traces
543- The solution status codes written to AMPL solution files now match those from https://mp.ampl.com/details.html#_CPPv4N2mp3sol6StatusE.
544
545### Changed parameters
546
547- change default value of parameter separating/filtercutpoolrel to FALSE
548- Shrunk range of parameter "propagating/symmetry/sstleaderrule" from [0,3] to [0,2], removed "MAXCONFLICTS"-option.
549- Removed reading/gmsreader/freeints. A default upper bound of +inf is now assumed for integer variables in GAMS (which is the case since GAMS 23.1).
550- Disabled "propagating/symmetry/symfixnonbinaryvars" since it is not used in the new symmetry code.
551- Replaced "propagating/symmetry/onlybinarysymmetry" with "propagating/symmetry/enforcecomputesymmetry" and set to FALSE by default.
552- Disabled "propagating/symmetry/performpresolving" since running orbital reduction during presolving is not effective.
553- Repurposed "propagating/symmetry/usedynamicprop" since it now also applies to lexicographic reduction.
554- Changed default value of parameter "propagating/symmetry/addsymresacks" from FALSE to TRUE.
555
556### New parameters
557
558- separating/maxcutsgenfactor and separating/maxcutsrootgenfactor as factors (w.r.t. separating/maxcuts and
559 separating/maxcutsroot) for limits on the number of cuts generated (which used to be constants '2')
560- expr/pow/expandmaxexponent to specify limit on exponent when to expand power of sum
561- expr/pow/distribfracexponent to enable expanding power of products with fractional exponents
562- expr/product/expandalways to enable expanding any product with sums as factor (that is, also when more than two factors)
563- nlhdlr/convex/maxperturb to specify maximal perturbation of non-differential reference points
564- benders/default/nlpiterlimit to adjust iteration limit for NLP solver for Benders' decomposition NLP subproblems
565- limits/objectivestop to interrupt the solve as soon as a primal solution is found
566 that is at least as good as the given value
567- heuristics/undercover/{coverand,coverind,covernl} to control which nonlinearities to consider (by default only "and" and nonlinear constraints)
568- cutselection/ensemble/∗ all parameters for new ensemble cut selector
569- cutselection/dynamic/∗ all parameters for new dynamic cut selector
570- branching/relpscost/gmiavgeffweight specifies the weight in the relpscost branching score of the normalized efficacy of all generated GMI cuts (for a variable)
571 branching/relpscost/gmilasteffweight specifies the weight in the elpscost branching score of the normalized efficacy of the last generated GMI cut (for a variable)
572- nlhdlr/quadratic/usemonoidal to enable monoidal strengthening when separating quadratics with intersection cuts
573- nlhdlr/quadratic/trackmore to collect more statistics for monoidal strengthening when separating quadratics with intersection cuts
574- nlhdlr/quadratic/useminrep to use the minimal representation to compute cut coefficients when separating quadratics with intersection cuts
575- nlhdlr/quadratic/sparsifycuts to enable sparsification of intersection cuts for quadratics
576- heuristics/dps/reoptlimits to set strict limits for the reoptimization step in dps
577- heuristics/dps/timing to run dps after processing the node
578- heuristics/scheduler/∗ all parameters for new scheduler heuristic
579- nlhdlr/signomial/enabled, nlhdlr/signomial/maxnundervars, nlhdlr/signomial/mincutscale: parameters for the signomial nonlinear handler
580- separating/lagromory/freq to denote the frequency for calling Lagromory separator
581- separating/lagromory/away to denote the minimal integrality violation of a basis variable to try separation
582- separating/lagromory/dynamiccuts to denote whether the generated cuts be removed from LP if they are no longer tight
583- separating/lagromory/maxroundsroot to denote the maximal number of separation rounds in the root node
584- separating/lagromory/maxrounds to denote the maximal number of separation rounds in a non-root node
585- separating/lagromory/nmaxcutsperlproot to denote, in the root node, the maximal number of Lagromory cuts that can be
586 separated per explored basis of the Lagrangian dual
587- separating/lagromory/nmaxcutsperlp to denote, in a non-root node, the maximal number of Lagromory cuts that can be
588 separated per explored basis of the Lagrangian dual
589- separating/lagromory/perroundnmaxlpiters to denote the maximal number of separating LP iterations that are allowed per
590 separation round
591- separating/lagromory/dualdegeneracyratethreshold to denote the minimum dual degeneracy rate for the separator execution
592- separating/lagromory/varconsratiothreshold to denote the minimum variable-constraint ratio on the optimal face for the
593 separator execution
594- many other advanced parameters for the Lagromory separator, which are also mentioned under separating/lagromory/∗
595
596### Data structures
597
598- new data structure SYM_GRAPH to encode symmetry information from a constraint
599- new data structure SYM_EXPRDATA to encode symmetry information from expressions
600- replaced enum SYM_SYMTYPE by a complete new one; removed SYM_RHSSENSE, added SYM_NODETYPE and SYM_CONSOPTYPE
601
602Deleted files
603-------------
604
605- deleted the now unused files: tpi_none.h, tpi_openmp.h, tpi_tnycthrd.h, type_tpi_none.h, type_tpi_openmp.h, type_tpi_tnycthrd.h
606
607Unit tests
608----------
609
610- new unittests for nlhdlr_quadratic that test monoidal strengthening and the use of the minimal representation to compute cut coefficients
611- unittests for nlhdlr_signomial
612
613Testing
614-------
615
616- changed `set lim obj` to `set lim objective` in ctest since `obj` is no longer unambiguous
617
618Build system
619------------
620
621### Cmake
622
623- New flag -DLAPACK=on for linking with Lapack (must be available in the system)
624- added flag option "SYM=sbliss" for using sassy/bliss as a graph automorphism package
625- use SYM=sbliss by default, since sassy and bliss are now shipped with SCIP.
626- added flag option "SYM=snauty" for using sassy/nauty as a graph automorphism package
627
628### Makefile
629
630- added flag option "nauty" for SYM variable to specify which graph automorphism package should be used
631- added flag option "sbliss" (sassy/bliss) for SYM variable to specify which graph automorphism package should be used
632- added flag option "snauty" (sassy/nauty) for SYM variable to specify which graph automorphism package should be used
633- use SYM=sbliss by default, since sassy and bliss are now shipped with SCIP
634- New flag LAPACK for linking with Lapack (must be available in the system)
635- A file scip/config.h is created in the build directory ($(OBJDIR)/include) now.
636 Defining NO_CONFIG_HEADER is no longer necessary and has no effect anymore.
637
638Fixed bugs
639----------
640
641- Number of runs in which separation is enabled will now be correct as per "separating/maxruns" parameter. Earlier,
642 separation was enabled for "maxruns + 1" number of runs, i.e., separation was enabled even when "maxruns = 0" was set.
643- add safeguard for number of selected cuts
644- use double double precision for delta in activity updates of cons_linear.c
645- sort constraints in checkSolOrig() to keep the logic of constraint handlers intact
646- fixed make install when SHARED=no and make uninstall
647- check absolute coefficients relative to extremal average to avoid invalid unifications in normalizeCons()
648- enforce and check symmetry-based constraints to take their locks (in dual reductions) into account
649- flush changes at the end of SCIPlpiSetBase() in lpi_grb.c to apply the newly set basis correctly
650
651Miscellaneous
652-------------
653
654- The TPI has been changed to become mostly independent from the actual parallel interface during compile time.
655 For this, the locks and conditions are now pointers to abstract structs.
656- fixed that scip/def.h, the header for C preprocessor definition (`#define`'s), included headers for some SCIP type and function definitions
668- fixed that scaling of rows in Gurobi interface (SCIPlpiScaleRow()) changed infinite sides to large finite ones (1e20) when scale factor is negative
669- check whether objective limit is truly reached in SCIPlpiIsStable()
670- if an objective limit is set and Mosek claims that it is exceeded, check that this is indeed the case in SCIPlpiIsStable()
671- the .mps file reader was extended to additionally accept .mps files with the objective sense section on one line
672- ensure that delayed cutpool separation does not happen if enough cuts have been found already
673- fixed that a Gurobi or Mosek interface (LPI) could not be freed from a thread different than the one it was created in, which lead to problems when, e.g., using the concurrent optimizer
674- in cons_varbound if both variables are fixed, instead of propagating the constraint check the feasibility of the constraint since propagating can lead to unnecessary numerical troubles
675- fixed memory leak in branching rule "distribution"
676- fixed that undercover did not enforce that sufficiently many variables in an and-constraints were fixed to have it linearized in the generated sub-scip
677- fixed two bugs in cons_varbound: avoid adding variable bounds if vbdcoef is infinite and unlock/lock variables if lhs/rhs/vbdcoef are changed
678- globalize constraint changes before bound changes during the effective root shift in treeSwitchPath() to correctly update variable locks by avoiding global bounds with outstanding bound change events when deleting constraints
679- fixed error when attempting to copy SCIP with activated pricers
680- remove integrality check for implicit integral variables in cons_integral.c, updated the definition for variable type implicit integer to be more explicit that these variables can take fractional values, and updated FAQ entry
681- updated and, or, and pseudoboolean constraint handlers to better handle fractional implicit integer variables in enforcement
682- avoid invalid strong Chvatal-Gomory cuts by ceiling with tolerance in the calculation of parameter k and flooring slack coefficients with standard epsilon tolerance as well as quad precision
683- add missing row links to columns in SCIPflushRowExtensions() required to maintain correct row norms
684- restrict separator intobj to solutions for which there can be no better feasible integral objective values to avoid cutting them off
685- delay removal of clique lists from SCIPvarRemoveCliquesImplicsVbs() into cliqueCleanup() to avoid segmentation faults by incomplete cleanups after variable conversions
686- reject sub-solutions with invalid objectives in SCIPtranslateSubSols() due to contradicting infinite contributions
687- correct column index in getMinMaxActivityResiduals() to make presolver dualinfer work
688- copy quiet flag from main scip's message handler to the one of solverscip in concurrent optimization
689- avoid bound tightening cycles towards huge redundancy bounds in fullDualPresolve() to avert fatal numerical trouble in representing vertices of linear relaxations with unbounded domains
690- fix problem where debug solution did not work if branching constraints are present
691- update index of dominating variable on coefficient deletion in dualPresolving() of constraint handler logicor
692- pass on flags from a redundant to the implying varbound constraint in preprocessConstraintPairs() to avoid suboptimal fixings by model relaxations
693- apply reductions of dualinfer only if a dual solution exists (primal is bounded)
694- before presolving is finished check solutions before retransforming them including tiny values to keep check flags consistent
695- disable heuristics before and during presolving in benders decomposition to avoid messing up the objective structure in the solution store by dropping auxiliary variables in the required retransformations
696- correct local flag of singleton conflict constraints in tightenSingleVar() against invalid bound globalization
697- respect unboundedness in the computation of activity bounds in conflict.c to avoid invalid huge bounds due to small coefficients on unbounded variables
698- ensure positive sides of a linear constraint when recognizing a set partition in rangedRowSimplify() to account for redundancy issues
699- relax numerical conditions for variable aggregations to avert invalid variable fixings
700- fixed harmless read of uninitialized data when creating parameters
701- fix call of SCIPnlpiOraclePrintProblem in nlpi_worhp
702
703Performance improvements
704------------------------
705
706- Use sassy/bliss as default symmetry computation package.
707
708Interface changes
709-----------------
710
711### Interfaces to external software
712
713- added interface to nauty/traces for symmetry computation
714- added interface to sassy, a preprocessor for symmetry computation
715- The directory src/sassy contains the source code of sassy.
716
717Build system
718------------
719
720### Cmake
721
722- added flag option "SYM=sbliss" for using sassy/bliss as a graph automorphism package
723
724### Makefile
725
726- added flag option "nauty" for SYM variable to specify which graph automorphism package should be used
727- added flag option "sbliss" for SYM variable to specify which graph automorphism package should be used
728- use SYM=sbliss by default, since sassy and bliss are now shipped with SCIP
729
730Miscellaneous
731-------------
732
733- the parameter change callback is no longer called at the moment a parameter is created
734
735@section RN804 SCIP 8.0.4
736*************************
737
738Features
739--------
740
741- use quad precision in cutsSubstituteMIR() for substituting slack variables into CMIR cuts
742- use quad precision for activity computations in cons_linear.c
743
744### Deleted and changed API methods
745
746- SCIPshrinkDisjunctiveVarSet() no longer attempts to cutoff the current node if it finds the
747 given disjunction to be infeasible; it is now in the responsibility of the caller to act accordingly
748
749### New API functions
750
751- added SCIPskipSpace
752- added XML parser routines from src/xml/xml.h
753
754### Interfaces to external software
755
756- updated XPRESS interface to work with XPRESS 41 as well
757- new interface to LP solver HiGHS (lpi_highs.cpp)
758
759### Changed parameters
760
761- changed integer parameter "propagating/symmetry/recomputrestart" from range {0, 1, 2} to being fixed at value 0
762
763Unit tests
764----------
765
766- added new checks for objective, bounds and sides inside the scaling tests in lpi/change.c
767
768Fixed bugs
769----------
770
771- relax condition of assert in lpSolve to hold only if lp/checkstability is true
772- OPB reader: dynamically increase the size of the line buffer for reading in getMaxAndConsDim()
773- Enabled the solution limit checks in SCIPsolveIsStopped in solve.c when in SCIP_STAGE_PRESOLVING
774- fix quad precision comparison for the calculation of knapsackcover constraints.
775- force copying of symmetry constraints for concurrent solve
776- increase the precision of double multiplications in quad precision environments
777- corrected handling of different variable types in symmetry
778- orbitope detection: refine computation of overlap structure and fixed bug in finding orbitope generators
779- ensure that winner is always determined in concurrent mode
780- when a row with invalid index range is added to the global cutpool, update the indices instead of failing with an assert
781- update of watched literals in cons_xor.c and ensure that fixed variables are counted in propacateCons()
782- avoid some numerical issue in check for redundant variables in linear constraints (introduced with version 8.0.2)
783- more rows are treated as being integral in cutsSubstituteMIR() for CMIR cuts
784- fixed several relaxations of the rhs of cuts to 0 in cuts.c
785- avoid aggregation into a single linear constraint if they should not be upgraded (needed for indicator constraints)
786- apply effective root changes before fork is repropagated during treeSwitchPath() in tree.c
787- avoid invalid objective value for solutions generated from an unbounded LP relaxation in case of a nonlinear objective function
788- disable recomputing symmetries after a restart (constraint removal is not permitted)
789- avoid cutting off optimal solutions for small feasibility tolerances in conflict analysis
790- initialize variable in presol_dualsparsify.c
791- fix in presolving of variable bound constraints to use absolute instead of relative difference
792- free variable histories in reoptimization mode
793- handle whitespaces in interactive shell history consistently
794- check and enforce equalities on which the objective is projected in cons_linear.c
795- ensure row feasibilities of solution shiftings in heur_oneopt.c and heur_twoopt.c
882- new parameters to turn on more statistics: misc/showdivingstats (diving heuristics), heuristics/alns/shownbstats (ALNS), estimation/showstats (tree size estimation)
883
884### Data structures
885
886- added nrootcalls and ncutsadded to SCIP_Cutpool
887- added fromcutpool flag to SCIP_Row and fixed 33 bit allocation
1038 Removed support for parameterized expressions.
1039- Added support for expression types sin, cos, and entropy. Added an expression type that stores a SCIP variable.
1040- Added handling of sin and cos expressions in OSiL reader.
1041- Added handling of sin, cos, and tan expressions in ZIMPL reader.
1042- Added expression handler plugin type, removed user expression type.
1043- Replaced expressions, expression trees, and expression graphs by a new implementation of expressions.
1044- Added differentiation methods (Gradients and Hessian directions) to expressions.
1045- Added integrality information to expressions.
1046- Added comparison and hashing of expressions.
1047- Improved simplification of expressions and identification of common subexpressions, in particular terms that appear
1048 equally in several quadratic functions are now identified.
1049- Expression simplification now expands powers also when the resulting exponents will be larger than 2.
1050- Added table to print statistics on expression handlers.
1051- Removed interval-arithmetic evaluation from expression interpreter (interfaces to automatic differentiation codes).
1052- Hessians for user-provided expression handlers are not available in the NLP solver interfaces for now.
1053- Sparse Hessians are now returned by expression interpreter. Removed dense Hessians.
1054
1055Constraint handlers for NLP:
1056- cons_nonlinear has been completely rewritten and now handles constraints of the form lhs <= expr <= rhs, where
1057 expr is an expression using the new expression type. A linear part of a nonlinear constraint is no longer distinguished.
1058- Methods that work on a particular nonlinear structure are now implemented via nonlinear handlers (nlhdlr).
1059 These handlers are plugins of cons_nonlinear and interact with cons_nonlinear by detecting structures in expressions and
1060 providing propagation and separation mechanisms.
1061- Added nlhdlr "default" as a mandatory nlhdlr: this is a nlhdlr that gets active on expressions that no other nlhdlr handles.
1062 It then uses the callbacks implemented in expression handlers to provide domain propagation and linear under/overestimates.
1063- Added nlhdlr for quadratic expressions: domain propagation that was previously part of cons_quadratic; separation via intersection cuts (currently disabled).
1064- Added nlhdlr for bilinear expressions: collects bilinear products and provides under/overestimation of a product over a 2D-projection of the LP.
1065- Added nlhdlr for convex and concave expressions: detection of convex/concave subexpressions and linear under/overestimation.
1066- Added nlhdlr for quotients: detects quotients, provides specialized linear under/overestimation and domain propagation for univariate quotients.
1067- Added nlhdlr for perspective reformulations: detects expressions in semi-continuous variables and tightens linear
1068 under/overestimates of these expressions by applying a perspective transformation.
1069- Added nlhdlr for second-order cones: detects second-order cone constraints (as ||x|| <= t or a quadratic that can written like that) and
1070 separates on a disaggregated reformulation (similar to cons_soc previously).
1071- Removed constraint handlers abspower, bivariate, quadratic, and soc.
1072- Removed simplification of duplicate abspower constraints for now.
1073 Removed simplification of abspower constraints with loose variables.
1074 Removed resolving of domain propagation (ConsResProp callback) and conflict analysis for abspower constraints.
1075 Removed propagation of varbounds for variables of abspower constraints.
1076 Removed reformulation of abs(x)^p * x into signpower expression for now.
1077- Removed reformulation of SOC constraints when adding to the NLP.
1078 Removed possibility to add linearizations a la Glineur or BenTal&Nemirovski for SOC constraints.
1079- Removed possibility to disaggregate quadratic constraints.
1080 Removed lifted tangent inequalities for bilinear terms in quadratic constraints.
1081 Removed option to project LP solution onto feasible region of a quadratic constraint via interior point or gauge function.
1082 Removed handling of unbounded LPs when enforcing quadratic constraints.
1083- Removed construction of cuts from interval gradients for nonlinear constraints.
1084 Removed possibility to pretend that all nonlinear constraints are convex.
1085 Removed possibility to create local or non-initial nonlinear constraints for now.
1086- Removed explicit reformulation of nonlinear constraints in presolve.
1087 An extended formulation is now generated implicitly and used to construct the LP relaxation,
1088 but the original formulation is used to check feasibility, propagate variable domains, etc.
1089 Variables that are added for the extended formulation are now marked as relaxation-only.
1090- Improved locking of nonlinear variables in nonlinear constraints by taking monotonicity into account.
1091- Improved mechanisms to detect and handle structures in nonlinear constraints, which allows for several algorithms
1092 to work simultaneously on equal or overlapping expressions.
1093- Improved linearization of (sums of) products of binary variables, e.g., clique information is now taken into account.
1094 Removed option to linearize products between a binary and a non-binary variable.
1095- Nonlinear constraints with the same function (but possibly different left- or right-hand-side) are now merged.
1096- When a variable appears in only one concave less-or-equal constraint, it will be fixed to one of its bounds in more cases now.
1097- New branching rule for variables in nonlinear constraints (implemented in cons_nonlinear) that scores candidates
1098 based on constraint violation, variable type, and pseudo-costs simultaneously.
1099- Added tables to print statistics of nonlinear constraint handler and nonlinear handlers.
1100
1101Performance improvements
1102------------------------
1103
1104- Use the objective cutoff row as base row for separation in sepa_aggregation.c.
1105- Separate lifted cover cuts based on newer lifting function of Letchford and Souli (2019) in sepa_aggregation.c.
1106- In-tree restarts due to tree size estimation have been made compatible with orbital fixing.
1107- Improved upgrade of orbisacks to packing/partitioning orbitopes.
1108- Improved running time of symresack cover separation.
1109- Improved propagation algorithm for orbisack and symresack.
1110- Improved handling for expressions of type signpower in CppAD interface by making use of CppADs conditional expressions.
1111- Faster reevaluation of Hessian for quadratic expressions in CppAD interface.
1112- Enabled optimization of taped expressions in the CppAD interface.
1113- Sparse Hessian evaluation for expressions with sparse Hessian in the CppAD interface.
1114- Disabled branching on variables that are added for the extended formulation of a nonlinear constraint.
1115- Prefer branching on variables in nonconvex terms of nonlinear constraints also if a cut could be generated but would yield a small reduction in the convexification gap only.
1116- When violations of nonlinear constraints are tiny, the feasibility tolerance for the LP relaxation may now be reduced temporarily to allow for cuts to have an effect.
1117- Improved under/overestimation of multidimensional vertex-polyhedral (e.g., multilinear or concave) functions by use of scaling and keeping the cut-generating LP around.
1118- Products in nonlinear constraints are no longer disaggregated for the extended formulation.
1119- No longer relax bounds of integral variables when computing bounds on activity in nonlinear constraints.
1120- Revised and improved Ipopt interface.
1121- Revised and improved subNLP primal heuristic, in particular choice of iteration limit and starting condition and reuse of subSCIP.
1122- Merge strongcg into the gomory separator to avoid computing the same base equality twice. The files
1123 sepa_strongcg.{h|c} have been removed.
1124- The RENS neighborhood in ALNS now fixes fractional variables if necessary for reaching its target fixing rate.
1125- Use cpsolver emphasis in sepa_rapidlearning.
1126- If decomposition/maxgraphedge = 0, the computation of the block graph is completely turned off.
1127- Fixed a performance bug in cons_integral.c where rounding information was passed incorrectly to the generic diving algorithm.
1128
1129Examples and applications
1130-------------------------
1131
1132- Add c++ Sudoku example
1133- Removed circle.c and renamed string.c to spring.c in CallableLibrary example.
1134- SCIP-Jack (Steiner tree application):
1135 Most files have been modified and many new ones have been added (code base more than doubled).
1136 Two new problem classes are supported: maximum-weight connected subgraph problem with budget constraints, partial-terminal node-weighted Steiner tree problem.
1137 Major performance improvements accross all problem classes that are supported.
1138
1139Interface changes
1140-----------------
1141
1142### New and changed callbacks
1143
1144- extend SEPAEXEC{LP,SOL} callbacks by parameter "depth", which should be used within separation instead of SCIPgetCurrentDepth()
1145- extend SCIPseparateSolCutpool() by parameter "pretendroot", which is passed to SCIPcutpoolSeparate() instead of SCIPtreeGetCurrentDepth(scip->tree) == 0
1163- added SCIP_NLPIPROBLEM* problem to SCIP_DECL_NLPIGETSOLVERPOINTER
1164
1165### Deleted and changed API methods
1166
1167- When calling SCIPdelVar() in the SCIP_STAGE_PROBLEM stage, it is now documented as the user's responsibility to ensure
1168 the variable is not present in any constraints. Otherwise, SCIPsolve() now returns SCIP_INVALIDDATA.
1169- removed SCIPselectCuts(), replaced by SCIPselectCutsHybrid() of the hybrid cut selector plugin
1170- If SCIPvarMarkRelaxationOnly() is called for a variable, then this now need to happen before the
1171 variable is added to the transformed problem (SCIPaddVar()).
1172 SCIPvarMarkRelaxationOnly() will automatically call SCIPvarMarkDeletable().
1173- SCIPcomputeOrbitsFilterSym() and SCIPcomputeComponentsSym(): Type of argument componentblocked
1174 changed from SCIP_Shortbool* to unsigned*.
1175- removed SCIPsetParam()
1176- SCIPcreateConsOrbitope() and SCIPcreateConsBasicOrbitope() have two new arguments "usedynamicprop"
1177 and "mayinteract" to encode whether a dynamic reordering of the orbitope rows shall be used in
1178 full orbitope propagation and whether the orbitope might interact with other symmetry handling
1179 methods, respectively.
1180- internal function SCIPgetPropertiesPerm() has been replaced by SCIPisInvolutionPerm()
1181- internal function SCIPgenerateOrbitopeVarsMatrix() has new arguments "storelexorder", "lexorder",
1182 "nvarsorder", "maxnvarsorder" to compute and store the lexicographic order defined by an orbitope
1183- renamed SCIPcomputeArraysIntersection() and SCIPcomputeArraysSetminus() to
1184 SCIPcomputeArraysIntersectionInt() and SCIPcomputeArraysSetminusInt(), respectively, and changed
1185 their return type from SCIP_RETCODE (which _always_ was SCIP_OKAY) to void
1186- removed SCIPincludeDialogDefault(), use SCIPincludeDialogDefaultBasic() instead and check whether you also need SCIPincludeDialogDefaultSet() and SCIPincludeDialogDefaultFix()
1187- removed exprint argument from SCIPaddNlRowGradientBenderscutOpt()
1188- added copyexprhdlrs argument to SCIPcopyPlugins() (19th position)
1189- added copycutselectors argument to SCIPcopyPlugins() (10th position)
1190- removed refpoint argument from SCIPaddSquareSecant()
1191- added argument exprdata to SYMcomputeSymmetryGenerators()
1206- renamed SCIPsetNlRowExprtree() to SCIPsetNlRowExpr() and changed type of expr(tree) argument to SCIP_EXPR*
1207- SCIPenableNLP() is no longer allowed in stages INITSOLVE and SOLVING
1208- removed SCIPaddLinearConsToNlpHeurSubNlp()
1209- SCIPhasNLPSolution() now returns true also if a infeasible solution to the NLP exists; SCIPcreateNLPSol() can now also be used when only an infeasible solution to the NLP exists
1210- renamed SCIPsolveNLP() to SCIPsolveNLPParam() and added parameter SCIP_NLPPARAM param
1211- removed SCIPsolveDiveNLP(), use SCIPsolveNLPParam() instead
1259- removed SCIPnlpiOracleGetVarDegree() and SCIPnlpiOracleGetVarDegrees(), use SCIPnlpiOracleIsVarNonlinear() and SCIPnlpiOracleGetVarCounts(), respectively, instead
1260- removed SCIPnlpiOracleGetConstraintDegree(), use SCIPnlpiOracleIsConstraintNonlinear() instead
1261- removed argument messagehdlr from SCIPnlpiOraclePrintProblem() and SCIPnlpiOraclePrintProblemGams()
1262- removed SCIPcreateNlpiProb(), use SCIPcreateNlpiProblemFromNlRows() instead
1263- renamed SCIPupdateNlpiProb() to SCIPupdateNlpiProblem() and moved into scip_nlpi.{h,c}
1264- renamed SCIPaddNlpiProbRows() to SCIPaddNlpiProblemRows() and moved into scip_nlpi.{h,c}
1266- removed SCIPgetIpoptApplicationPointerIpopt(), use SCIPgetNlpiSolverPointer() instead
1267- removed SCIPnlpStatisticsCreate(), SCIPnlpStatisticsFree(), SCIPnlpStatisticsGetNIterations(), SCIPnlpStatisticsGetTotalTime(), SCIPnlpStatisticsSetNIterations(), SCIPnlpStatisticsSetTotalTime(), the struct is public now
1268- renamed LapackDsyev() to SCIPcallLapackDsyevIpopt()
1269- renamed SCIPsolveLinearProb() to SCIPsolveLinearEquationsIpopt()
1270
1271Expression interpreter:
1272- removed functions that operated with SCIPInterval
1273- changed first argument of SCIPexprintCreate() to SCIP*
1274- added SCIP* as first argument to SCIPexprintFree()
1275- added SCIP* as first argument to SCIPexprintCompile(), changed third argument to SCIP_EXPR*, added fourth argument SCIP_EXPRINTDATA**
1276- added arguments SCIP*, SCIP_EXPRINT*, and SCIP_EXPR* as first arguments to SCIPexprintFreeData()
1277- renamed SCIPexprintGetExprtreeCapability() to SCIPexprintGetExprCapability();
1278 added SCIP* as first argument, changed third argument to SCIP_EXPR*, added fourth argument SCIP_EXPRINTDATA*
1279- removed SCIPexprintNewParametrization()
1280- SCIPexprintEval(): added SCIP* as first argument, changed third argument to SCIP_EXPR*, added fourth argument SCIP_EXPRINTDATA*
1281- removed SCIPexprintEvalInt()
1282- SCIPexprintGrad(): added SCIP* as first argument, changed third argument to SCIP_EXPR*, added fourth argument SCIP_EXPRINTDATA*
1283- removed SCIPexprintGradInt()
1284- removed SCIPexprintHessianSparsityDense() and SCIPexprintHessianDense, use SCIPexprintHessianSparsity() and SCIPexprintHessian() instead
1285
1286Expressions:
1287- removed SCIPexpropGetName(), use SCIPexprhdlrGetName() instead
1288- removed SCIPexpropGetNChildren()
1289- removed SCIPexprGetOperator(), use SCIPexprGetHdlr() instead
1290- removed SCIPexprGetOpIndex(), SCIPexprGetOpReal(), SCIPexprGetOpData(), use SCIPexprGetData() instead
1291- removed SCIPexprGetRealPowerExponent(), SCIPexprGetIntPowerExponent(), SCIPexprGetSignPowerExponent(), use SCIPgetExponentExprPow() instead
1292- removed SCIPexprGetLinearCoefs(), SCIPexprGetLinearConstant(), use SCIPgetCoefsExprSum() and SCIPgetConstantExprSum() instead
1412- renamed SCIPgetNAllBilinearTermsQuadratic() to SCIPgetNBilinTermsNonlinear() and changed argument from SCIP* to SCIP_CONSHDLR*
1413- renamed SCIPgetAllBilinearTermsQuadratic() to SCIPgetBilinTermsNonlinear() and replaced arguments by SCIP_CONSHDLR* and changed return type to SCIP_CONSNONLINEAR_BILINTERM*
1414- removed SCIPaddBilinearIneqQuadratic(), use SCIPinsertBilinearTermImplicitNonlinear() instead
1415- renamed SCIPcreateConsQuadratic() to SCIPcreateConsQuadraticNonlinear() and removed argument stickingatnode
1416- removed SCIPincludeConshdlrSOC()
1417- replaced SCIPcreateConsBasicSOC() by SCIPcreateConsBasicSOCNonlinear()
1418- deprecated SCIPcreateConsSOC(), use SCIPcreateConsBasicSOCNonlinear() instead
1419- deprecated SCIPgetNlRowSOC(), use SCIPgetNlRowNonlinear() instead
1639- Removed parameters for strongcg separator: separating/strongcg/..., i.e.,
1640 maxrounds, maxroundsroot, maxsepacutsroot, maxsepacuts, dynamiccuts. One can control the
1641 strongcg behavior throught the parameters of sepa_gomory.
1642- renamed separating/{minortho,minorthoroot,objparalfac,intsupportfac} to cutselection/hybrid/{minortho,minorthoroot,objparalweight,intsupportweight}
1643
1644### New parameters
1645
1646- It is now possible to add orbitope constraints to handle symmetries even if not all variables
1647 in a component of the symmetry group are binary. This feature can be enabled using the
1706- new parameters "branching/inference/conflictprio" and "branching/inference/cutoffprio" that allow
1707 using the different score criteria in inference branching hierarchically
1708- new parameters "heuristics/alns/nodesquotmin", "heuristics/alns/maxcallssamesol" and "heuristics/alns/initduringroot"
1709 to control neighborhoods in ALNS
1710- new parameter "benders/default/benderscut/optimality/mir" to control whether mixed-integer rounding procedure is applied to Benders optimality cuts
1711- new parameter "misc/avoidmemout" to try to avoid running into memory limit
1827- added test projection:transitivity in prob/obbt.c
1828- added reader/nl.c
1829- added test intervalarith:sincos in scip/intervalarith.c
1830- added sepa/minor.c
1831- added sepa/rlt/cuts.c
1832- added sepa/rlt/product_detection.c
1833- added sepa/rlt/separation.c
1834- added tests test_compute_symmetry:expr{1,2,3,4} in symmetry/compute.c
1835
1836Build system
1837------------
1838
1839- the separate NLPI library has been integrated into the main SCIP library;
1840 uses of NLPILIBFILE or NLPILIBSRC in Makefiles that use make/make.project
1841 should be removed
1842
1843### Cmake
1844
1845- replace flag PARASCIP by THREADSAFE
1846- option AMPL is now ON by default
1847- add -D_XOPEN_SOURCE=600 cflag for compilation on systems other than windows
1848- add compile options directly instead of modifying CMAKE_C_FLAGS and CMAKE_CXX_FLAGS
1849
1850### Makefile
1851
1852- prefix PKG_CONFIG_PATH with Ipopt's pkgconfig dir instead of postfixing
1853- replace flag PARASCIP by THREADSAFE
1854- flag AMPL is now true by default
1855- libscip has been renamed to libscipbase and libscipsolver has been renamed to libscip
1856 in order to have the same library names for makefiles and cmake. For downward
1857 compatibility libscipsolver still exists (as a link).
1858
1859Fixed bugs
1860----------
1861
1862- fixed handling of multi-aggregation in cons_symresack and cons_orbisack
1863- SCIPgetLowerbound() now consistently returns infinity when status is infeasible
1864- improved the likelihood that a solution that satisfies the nonlinear constraints in the transformed (presolved) problem is also feasible for the original problem
1865- worked around a problem with derivatives of pow at 0 with non-integer exponents in CppAD
1866- fixed that setting of feasiblity and optimiality tolerance for NLP relaxation was overwritten with defaults when solving NLP
1867- fixed variance computation in t-value computation in hypothesis testing of reliability branching
1868- fixed potential overflow in cuts.c:evaluateLiftingFunctionKnapsack()
1869- do not transfer dual reductions to main SCIP in sepa_rapidlearning if strong dual reductions are forbidden
1870- fixed a bug in cons_setppc.c:removeDoubleAndSingletonsAndPerformDualpresolve() which led to wrong variable aggregations
1871- added missing redundancy check to coefficient tightening routine in cons_varbound that led to wrong variable locks
1872- added a check of numerics to avoid wrong reductions in cons_linear.c:simplifyInequalities()
1873- fixed a bug in cons_logicor.c regarding fixed variables that were not removed
1874- presol_milp.c handles variable fixings, suggested by PaPILO, more carefully to cover all cases that can appear
1875- fixed several reproducability bugs that led to different solving paths across multiple runs with one SCIP environment
1876- fixed check of limits/absgap, which was incorrectly applied to the scaled objective of the transformed problem
1877
1878Miscellaneous
1879-------------
1880
1881- in LP, use absolute instead of relative tolerances to check whether solutions from LP solver are primal feasible
1882- when a restart occurs, SCIP now asserts (during EXITSOLVE) that all relaxation-only variables are unlocked and
1883 only captured by the transformed problem.
1884- Avoid redundant calls to LP solver when LP is already flushed and solved
1885- In SCIPdebugMsg(), SCIPsetDebugMsg(), SCIPstatDebugMsg() strip directory from filename.
1886- Recompute activity of rows when checking LP solutions instead of trusting the value given by the LP solver
1887- The define NPARASCIP has been replaced by SCIP_THREADSAFE
1888- Parameters are not reset to their default values before activating an emphasis setting; this now allows to combine
1889 different emphasis settings.
1890- If a lazy bound is set to a value that is tighter than the current global bound, then the global bound is now tightened accordingly.
1891
1892@page RN7 Release notes for SCIP 7
1893
1894@section RN703 SCIP 7.0.3
1895*************************
1896
1897Interface changes
1898-----------------
1899
1900### New API functions
1901
1902- added function SCIPwithDebugSol() to determine whether SCIP was compiled with WITH_DEBUG_SOLUTION
1903
1904- New function SCIPcreateConsBounddisjunctionRedundant() and SCIPcreateConsBasicBounddisjunctionRedundant() that allow
1905 the input to contain redundant literals.
1906
1907### Changed API functions
1908
1909- added __attribute__((format(printf,,,))) to variadic SCIP message printing functions;
1910 this enables warnings about mismatches between format string and variadic function arguments if using GCC.
1911
1912### Interfaces to external software
1913
1914- removed MATLAB interface from interfaces/matlab
1915
1916Build system
1917------------
1918
1919- Compilation will now also work in a raw tarball of the scip git repository
1920
1921### Cmake
1922
1923- Make find_package more verbose and require packages to be found. Now enabled packages that are not found result in an error.
1924
1925Fixed bugs
1926----------
1927
1928- Fixed numerical bug by creating and checking the empty solution when problem vanishes in presolving.
1929- Fixed several bugs in reoptimization.
1930- Fixed bug in cons_bounddisjunction in which a constraint leads to a cutoff.
1931- Fixed bug in reoptimization by avoiding redundant literals for bound disjunctions.
1932- Replace wrong usage of CMAKE_BINARY_DIR by PROJECT_BINARY_DIR in CMake system
1933- Fixed lpi_glop to use absl flag instead of GFlags
1934- Fixed invalid initialization of limits/softtime parameter in cons_components
1935- Fixed MIN macro not defined in some cases.
1936- Safeguard some API methods for querying solving statistics against segmentation faults after SCIPpresolve()
1937- Fixed unfreed solution in concurrent solving
1938- Make sure that logicor constraints are enforced/separated if genlogicor is true in cons_indicator.
1939 Thus, do not change constraints/logicor/sepafreq anymore, but output warning.
1940- Fixed assert that assumed infeasibility cannot occur in presol_boundshift.
1941- Fixed sign bug during creation of conflict constraints with almost zero left-hand side.
1942- Disabled call of conflict handler for conflict sets that contain relaxation-only variables.
1943- Fixed bug in cons_symresack by initializing all constraint data in constrans
1944- Avoid allowing the upgrade of linear constraints corresponding to indicators in propagation of nodes.
1945- Fixed bug in computing lex. min/max faces in dynamic orbitope propagation.
1946- Limit vbounds propagator to one sweep through all bounds per call in order to avoid an infinite loop.
1947- Fixed mismatches between formatter strings and arguments in calls to SCIP message printing functions.
1948
1949Miscellaneous
1950-------------
1951
1952- To avoid performance variablity during reoptimization, the active constraints are now stored in the same order as in the
1953 original problem.
1954- Added workaround for bug in MSVS 2019 when compiling CppAD source.
1955- The original order of variables is now only restored after presolve if `randomization/permutevars` is set to FALSE.
1956
1957@section RN702 SCIP 7.0.2
1958*************************
1959
1960Features
1961--------
1962
1963- Extended check for valid variable and constraint names in PIP writer.
1964- Indicator constraints now write the name of the linear constraint in CIP format. This avoids a complicated
1965 name matching during parsing. The old format is still supported.
1966- The objective scale and offset are now taken into account when reading/writing problems in lp/mps/fzn/pip format.
1967 This often changes the behavior when writing transformed problems.
1968
1969Performance improvements
1970------------------------
1971
1972Examples and applications
1973-------------------------
1974
1975- Add c++ Sudoku example
1976
1977Interface changes
1978-----------------
1979
1980### Deleted and changed API methods
1981
1982- deprecated SCIPsetParam(); use SCIPsetIntParam(), SCIPsetRealParam(), etc, instead
1983
1984### New API functions
1985
1986- new function SCIPcolGetVarProbindex() for accessing the probindex of a variable corresponding to a column
1987- added new sorting functions SCIPsortPtrRealRealBoolBool()
1988
1989Fixed bugs
1990----------
1991
1992- Fixed a sign bug occuring for some cases in the computation of the directed cutoff distance
1993- Fixed overflow in inference information of cumulative constraint handler
1994- Fixed problem with propfreq == 0 in cons_linear.c
1995- Fixed aggregation heuristic in sepa_aggregation.c to not use local rows when that is requested nor modifiable rows
1996- Fixed problem with incompatibility of symmetry handling in Benders decomposition
1997- Fixed parsing of small number in CIP format of nonlinear constraint
1998- Fixed problem in upgrading full orbitopes to partitioning orbitopes
1999- Fixed treatment of negated variables in parsing orbitope constraints
2000- Fixed potential indexing error in prop_vbounds.c
2001- Fixed bug in writing fzn format
2002- Revert to old parameter value when parameter change is rejected in PARAMCHGD callback
2003- Fixed order of linker flags when building AMPL interface with cmake system.
2004- Fixed bug with ofsymcomptiming == 0: now symmetry is actually computed
2005- Improved separation for SOC constraint with multiaggregated variable on right-hand-side.
2006- Fixed bug in presol_boundshift: do not shift non-active variables
2007- Fixed bug in domain propagation for SOC constraints with negative coefficient for variable on right-hand-side.
2008- Fixed bug with multi-aggregated integer variables in cons_xor.
2009- Fixed bug with changed number of variables in storing relaxation solution in probing/diving.
2010- Fixed issues with scip statistics
2011- Fixed issue with multi-aggregated variables
2012- Fixed issue in with buffer memory in cons_cumulative
2013
2014@section RN701 SCIP 7.0.1
2015*************************
2016
2017Build system
2018------------
2019
2020### Cmake
2021
2022- Fixed cmake module for Ipopt to support Ipopt 3.13
2023- Added cmake variable BLISS_LIBRARY_DIR to pass a bliss library
2024- Added warning to require minimal GCC version 4.8
2025- Renamed FindCriterion cmake findmodule to FindCRITERION
2026
2027### Makefile
2028
2029- Updated make/make.project for macOS
2030- Add SANITZERFLAGS (only OPT=dbg) to LDFLAGS also if SHARED=true (and rename SANITZERFLAGS to SANITIZERFLAGS)
2031- Fixed GCC version check in make/make.linux.gnu.dbg
2032
2033Fixed bugs
2034----------
2035
2036- Fixed handling of constant term in quadratic objective of pip-reader
2037- Fixed incorrect value string for parameter estimation/completiontype
2038- Fixed initialization of first LP which led to wrong output about LPs not resolved correctly
2039- Fixed bug in SCIPcutpoolClear() function not clearing internal datastructures properly
2040- Fixed inability to set a value in a solution if it was within epsilon to the current value
2041- Fixed problems for linking optimized code with debug library of SCIP (and vice versa) by moving debug-only members
2042 to end of SCIP structs
2043
2044Miscellaneous
2045-------------
2046
2047- Variables that are marked as relaxation-only are no longer checked against debug solutions if
2048 they are captured by the transformed problem only.
2049
2050@section RN700 SCIP 7.0.0
2051*************************
2052
2053Features
2054--------
2055
2056- Using the parameter "propagating/symmetry/recomputerestart" one can now decide to recompute symmetries after a
2057 restart or not. Previously one could just turn off symmetry computation after a restart. If orbital fixing
2058 has found a reduction before the restart occured, symmetries have to be updated to ensure correctness.
2059 To this end, the user can decide via "propagating/symmetry/disableofrestart" whether orbital fixing is disabled or
2060 whether symmetries are recomputed.
2061
2062- extended relaxators statistics in solve statistics about information on number of cutoffs, improved bounds, etc.
2063
2064- extends SMPS file reader for the stochastic information, the sto files, to read a lower bound for the discrete
2065 scenarios. The lower bound is used when creating the auxiliary variables for Benders' decomposition.
2066
2067- extended Benders framework to solve NLPs and generate optimality and feasibility cuts from their solution
2068 when the subproblem is convex nonlinear.
2069
2070- extended Benders framework to create copies of Benders decompositions that can be used in a multithreading environment
2071
2072- additional functionality has been added to enable the transfer of Benders' decomposition cuts between different SCIP
2073 instances, e.g., when used in UG
2074
2075- LP rows (SCIP_ROW) can now store pointer to a constraint from which the row originates
2076
2077- Trust region, a new LNS heuristic, has been added to SCIP as both a standalone primal heuristic heur_trustregion.c
2078 and as a neighborhood inside of Adaptive Large Neighborhood Search (heur_alns.c). This heuristic is designed to improve the
2079 heuristic performance of the Benders' decomposition algorithm. It builds upon the successful trust region approaches
2080 applied to Benders' decomposition.
2081
2082- Modularity of symmetry handling has been increased. It is now possible to use orbitopes (i.e., polyhedral symmetry
2083 handling techniques) and orbital fixing on the same instance.
2084
2085- cut strengthening enabled within the Benders' decomposition framework. This uses a mix of the Improved Magnanti-Wong
2086 method and Kelly's method. The cut strengthening is enabled by setting the paramemter
2087 "benders/<bendersname>/cutstrengthenenabled" to TRUE. The parameters "cutstrengthenmult", "noimprovelimit" and
2088 "corepointperturb" control the behavior of the cut strengthening method. Additionally, the parameter
2089 "cutstrengthenintpoint" allows the user to specify the solution that is used to initialize the core point. The options
2090 are the first LP solution, the first integer feasible solution, a relative interior point, a solution vector of all
2091 ones and a solution vector of all zeros. Also, the core point can be reinitialized after each update to the incumbent
2092 solution.
2093
2094- added option to adjust weights of different scores in relpscost (hybrid) branching rule based on degeneracy information
2095 and skip strong branching for very high degeneracy rates
2096
2097- added new SCIP_DECOMP* object to store user decompositions. The user can specify one or several decompositions by assigning variable
2098 and/or constraint labels either through the API or by reading a file in DEC format (which is one of the formats that GCG writes).
2099 This format specifies constraint labels, variable labels are inferred from that. The decomposition is transformed after presolving.
2100- statistics regarding the number of blocks, the largest and smallest blocks, the border, and the constraint graph are printed for the original decomposition,
2101 and for each decomposition after presolving.
2102- the decomposition can be used as initial decomposition for the Benders functionality of SCIP.
2103- new subsection "How to provide a problem decomposition" in the SCIP documentation
2104
2105- GINS heuristic can make use of a user-provided decomposition labels in two ways:
2106 1. by selecting a block of variables that maximizes the potential, and randomly selecting a start variable for the neighborhood and/or
2107 2. by selecting an interval of consecutive blocks as neighborhood, until fixing rate is reached. In this case, no variable is randomly selected.
2108
2109- extend potential parameter of GINS to allow computation based on local LP solution, as well
2110
2111- new primal heuristic Adaptive Diving, which registers all publicly available dive sets from other diving heuristics.
2112 At each call, it selects one of the available dive sets based on the user's score type choice (heuristics/adaptivediving/scoretype).
2113 During the solution process, the heuristics learns online which divesets reach the best score, and executes them more frequently.
2114 The statistic output for Diving heuristics has been extended to incorporate the statistics of each dive set within Adaptive diving.
2115
2116- Added new penalty alternating direction method (PADM) primal heuristic that splits the problem into several sub-SCIPs
2117 according to a user-provided decomposition. The sub-SCIPs are solved on an alternating basis until a feasible solution of
2118 the original problem is found.
2119
2120- Symmetry handling constraints (cons_orbisack, cons_orbitope, cons_symresack) now have an additional parameter to encode
2121 whether they are model constraints, i.e., define the problem, or whether they are only present to handle symmetries.
2122
2123- The symmetry code has been completely restructured. The presolvers presol_symbreak and presol_symmetry as well as the propagator
2124 prop_orbitalfixing have been merged into the single propagator prop_symmetry to avoid storing the same information
2125 multiple times. This propagator is now responsible for adding symmetry handling constraints as well as activating orbital
2126 fixing. Moreover, the new file symmetry.c contains general functions for symmetry computations like orbit computations.
2127
2128- Variables can now be marked as "relaxation-only". This flag should be used to introduce new variables that are required
2129 to define a relaxation, but that are not part of any checked constraints. Essentially, these variables should only be used
2130 in the current SCIP solve and disregarded after a restart or in SCIP copies. Hence, these variables are not copied by SCIPcopy and
2131 SCIPgetVarCopy, they are currently not used in conflict constraints, and cuts involving them are not upgraded to linear constraints.
2132 Relaxation-only variables cannot appear in the objective function.
2133
2134- The OSiL reader now supports nonlinear expressions of type "signpower".
2135
2136- Expressions of form abs(x)^p * x in a nonlinear constraint are now sometimes recognized and handled by abspower constraints.
2137
2138- If polyhedral symmetry handling methods are used (cons_orbisack, cons_orbitope, cons_symresack), it is now possible to
2139 recompute symmetries if a restart occured.
2140
2141- upgrade some more quadratic constraints to second-order cone constraints, that is,
2142 handle linear binary variables as if squared in simple upgrade and
2143 do not require bounds for variables that have a zero entry in the computed eigenvectors in the non-simple upgrade
2144
2145- new variable event when changing the variable type
2146
2147- It is no longer necessary to provide a SCIP pointer for the subproblem in SCIPaddBendersSubproblem if custom solving
2148 methods are defined. A NULL pointer can be supplied to SCIPaddBendersSubproblem. In this case, no internal Benders'
2149 cut generation methods can be used.
2150
2151- Using the parameter "constraints/symresack/checkmonotonicity" one can now decide to upgrade to packing/partitioning
2152 symresacks even if the underlying permutation is not monotone.
2153
2154- New branching rule `vanillafullstrong`, mostly for scientific purpose, with the following features:
2155 1) no cutoff or domain reduction: only branching;
2156 2) idempotent (optional): leave SCIP, as much as possible, in the same state before / after the strong branching calls
2157 - basically, do not update any statistic;
2158 3) donotbranch (optional): do no perform branching. So that the brancher can be called as an oracle only
2159 (on which variable would you branch ? But do not branch please);
2160 4) scoreall (optional): keep scoring variables, even if infeasibility is detected;
2161 5) collectscores (optional): store the candidate scores from the last call, which can then be retrieved by calling SCIPgetVanillafullstrongData();
2162 6) integralcands (optional): consider integral candidates for branching,
2163 i.e., get candidates from SCIPgetPseudoBranchCands() instead of SCIPgetLPBranchCands().
2164
2165- If a reference value (misc/referencevalue) is given, the primal-reference and reference-dual integrals are calculated automatically and printed within the SCIP statistics.
2166
2167- Locally valid cuts / rows are now considered for dual proof analysis when `conflict/uselocalrows` is set to TRUE.
2168
2169- Linking variables in the linking constraint handler (cons_linking.{ch}) can now be integer or continuous. The coefficients of the binary variables are therefore now stored as SCIP_Real.
2170
2171- To save memory, it is now possible to remove all variables from the internal symmetry data structures that are
2172 not affected by any symmetry.
2173
2174- Allow to filter first variables from orbits and transfer pseudo cost information to variables in orbit
2175
2176- Add integration of external MILP presolve library as a SCIP presolver plugin that runs on MILP problems
2177
2178- Parallelisation can be used when applying Benders' decomposition. There are two different forms of parallelisation
2179 available. The first is applying Benders' decomposition within a parallel branch-and-bound. This is achieved through
2180 the integration with the UG framework. The second is the use of shared memory parallelisation for solving the Benders'
2181 decomposition subproblems. A priority queue has been added to help with load balancing.
2182
2183- The Benders' decomposition framework can handle MINLPs. If a convex relaxation of the MINLP exists, then this is
2184 solved to generate optimality and feasibility cuts. The extensions to the framework are:
2185 - New generic solving methods to solve convex NLP subproblems.
2186 - Modification to benderscut_opt and benderscut_feas to enable the generation of cuts from convex NLPs.
2187 - Addition of benderscut_feasalt to generate feasibility cuts from an alternative subproblem that minimises the
2188 violation of infeasible problems.
2189 - Better handling of subproblem solution results
2190
2191- Adds a feasibility phase to the Benders' decomposition subproblem solving methods. The feasibility phase adds slack
2192 variables to each of the constraints to ensure feasibility of the subproblem. A large coefficient is given to these
2193 slack variables in the objective function to penalise constraint violations. The coefficients are progressively
2194 increased until all slack variables take the value 0.0.
2195
2196- Improved convexity check for Benders' decomposition subproblems. The constraints of the subproblem are now checked for
2197 convexity in the initialisation of the Benders' decomposition algorithm. This enables the solving of convex NLPs as
2198 Benders' decomposition subproblems.
2199
2200- Benders' decomposition can be applied using decomposition supplied in the DEC format. To apply Benders' decomposition
2201 the parameters decomposition/benderslabels and decomposition/applybenders must be set to TRUE.
2202
2203- new event handler event_estim.c/h that approximates search tree completion and estimates tree size
2204 to trigger restarts; many approximations of search tree completion and estimation, including WBE, SSG, and tree profile method
2205- new display column that reports approx. search tree completion during the search, and an overview in the statistics table
2206- added resources (script, tutorial, test data) to adapt tree size estimation to user instances.
2207
2208- Orbital Fixing uses a list of variables that have been fixed globally since the computation of symmetries to filter
2209 symmetries. Previously, some plugins were disabled, which is not necessary anymore.
2210
2211- A new presolver "dualsparsify" was added. It tries to combine columns (i.e. variables) to cancel nonzero coefficients in the constraint matrix.
2212
2213- The presolver "tworowbnd" was implemented with better performance.
2214
2215- To be able to calculate better bounds for the dual variables, the presolver "dualinfer" was extended by the ability to perform convex combinations of continuous columns.
2216
2217- allow disabling of pricers during solving process
2218
2219- added emphasis setting for numerically challenging instances
2220
2221Performance improvements
2222------------------------
2223
2224- Extended cut presolving by removing variables that been fixed at their bounds
2225
2226- Improved branching point selection when branching on externals branching candidates. Instead of using exactly the
2227 LP solution, a point closer to the middle of the variables domain is chosen.
2228
2229- Matrix presolvers that do not work on incomplete matrices now skip matrix creation if unsupported constraint types are detected.
2230
2231- consLockBenders callback implemented to add down locks on the Benders' decomposition auxiliary variables and
2232 up and down locks per subproblem for all master problem variables. This allows the use of presolving and propagation
2233 with Benders' decomposition.
2234
2235- improved performance of orbital fixing in several ways:
2236 store permutations in transposed form to improve cache efficency;
2237 reverse order to speed up filtering of permutations;
2238 handle variables globally fixed to 1 in list;
2239 use event handler to catch global fixings;
2240 speed up orbit computations;
2241 change handling of restarts;
2242 use only permutations that can contribute to a variable's orbit;
2243
2244- allow rapid learning at local nodes
2245
2246- allow to recompute cut without using fractional values for sepa_cgmip
2247
2248- restrict the number of the clique table nonzeros relative to the number of problem nonzeros, which could be a performance bottleneck.
2249
2250- variable fixings of LP face heuristic are now computed earlier; subproblem creation is skipped if not enough variables are fixed.
2251
2252- Improved domcol presolver to not require a complete representation of all constraints in the matrix
2253
2254- performance improvement of adaptive large neighborhood search heuristic on merely continuous problems. The heuristic
2255 stops if presolving in the sub-SCIP fixes less than 50 % of the current target fixing rate over all variables
2256 (including continuous).
2257
2258- reduce memory usage in symmetry detection by a staggered allocation with decreasing overhead for larger instances
2259
2260- improved full orbitope propagation using a static implementation or a dynamic reordering of orbitope rows by a global rank function
2261
2262- improved detection of packing/partitioning orbitopes
2263
2264- enable an in-tree restart if after a reasonable initialization, the estimated size of the remaining tree is large.
2265
2266
2267Examples and applications
2268-------------------------
2269
2270- added methods to set and get hmin and hmax for optcumulative constraints
2271
2272Interface changes
2273-----------------
2274
2275### New and changed callbacks
2276
2277- new optional callback `SCIP_DECL_DIVESETAVAILABLE` to check preconditions for this dive set, e.g., if an incumbent solution is available,
2278 which is passed as new argument to SCIPcreateDiveset(). SCIPcreateDiveset() has another new parameter "ispublic".
2279
2280- new callback `SCIP_DECL_CONSHDLRCOPY` and `SCIP_DECL_CONSCOPY` in cons_orbisack and cons_symresack
2281
2282- new `idempotent` argument to SCIPgetVarStrongbranchInt() and SCIPgetVarStrongbranchFrac(),
2283 so that statistics are not updated during the call. Likewise, new `updatecol` and `updatestat` arguments to SCIPcolGetStrongbranch().
2284
2285- callback `SCIP_DECL_CONSHDLRENFOLP` can now also return SCIP_SOLVELP as *result, which indicates to the SCIP core that the LP relaxation
2286 should be solved again because the primal feasibility tolerance of the LP has been tightened (using SCIPsetLPFeastol())
2287
2288- extension of SCIP_PQUEUE by a new callback SCIP_DECL_PQUEUEELEMCHGPOS to catch swaps as well as functionality to delete arbitrary elements from the priority queue.
2289
2290### Deleted and changed API methods
2291
2292- LPI:
2293 + now for all lp interfaces consistent requirements on SCIP_LPPAR:
2294 LPITLIM and BARRIERCONVTOL positive or zero; FEASTOL, DUALFEASTOL, LPTILIM strictly positive
2295 + now projecting SCIP_LPPAR values on feasible values for each lp interface
2296 + add interface to Glop
2297 + fixed mapping between scaling parameter values in Gurobi LPI lpi_grb
2298- Symmetry:
2299 + removed method SCIPseparateCoversOrbisack() in cons_orbisack.h since the orbitope constraint
2300 handler has its own implementation of this routine with advanced features now
2301 + renamed SCIPgetGeneratorsSymmetry() to SCIPgetSymmetry() and removed two arguments
2302 + extended method SCIPgetSymmetry(): It is possible to access both the original and
2303 transposed permutations matrix as well as the (independent symmetry) components of a permutation group now.
2304 + arguments of functions SCIPcreateConsOrbisack(), SCIPcreateConsBasicOrbisack(), SCIPcreateConsOrbitope(),
2305 SCIPcreateConsBasicOrbitope(), SCIPcreateConsSymresack(), SCIPcreateConsBasicSymresack(), and SCIPcreateSymbreakCons() extended
2306 by "ismodelcons" to encode whether the constraints are model constraints or not
2307 + the function SCIPgetSymmetry() no longer accepts the parameter recompute, but has parameter permvarmap as new input
2311 + add function SCIPcomputeOrbitsComponentsSym() to compute orbits without filtering permutations and indices of orbits for each variable
2312- SCIPallowObjProp() and SCIPallowDualReds() are deprecated and replaced by SCIPallowWeakDualReds() and SCIPallowStrongDualReds(), respectively
2313- Benders' decomposition
2314 + changed SCIPstoreBenderscutCut() in scip_benders.c to SCIPstoreBendersCut(). Where this function used to take a
2315 SCIP_BENDERSCUT pointer, it now accepts a SCIP_BENDERS pointer.
2316 + the functions SCIPsolveBendersSubproblem() no longer accepts the parameter type. The type is not a necessary
2317 argument for the subproblem solving method.
2318 + arguments of functions SCIPbendersSolveSubproblemLP(), SCIPbendersSolveSubproblemCIP(), and SCIPbendersOnlyCheckConvexRelax() changed
2319 - removed SCIPbenderscutGetNAddedCuts() and SCIPbenderscutGetAddedCutData()
2320
2321
2322### New API functions
2323
2324- new argument "onlyifcomplete" in SCIPmatrixCreate() to skip matrix creation right after detecting unsupported constraint types
2325 and new arguments to count statistics when doing a clean-up of inactive variables in the constraints before building the matrix
2326
2327- new argument "threadsafe" in SCIPcopy(), SCIPcopyConsCompression(), SCIPcopyOrig(), SCIPcopyOrigConsCompression and
2328 SCIPcopyBenders(). This argument must only be set to TRUE if the source and target SCIP instances are to be solved in
2329 parallel. Setting this argument to TRUE has a performance cost.
2330
2331- new argument "append" in SCIPsetModifiedDefaultSettingsIpopt()
2332- functions SCIPclearRelaxSolVals(), SCIPsetRelaxSolVal(), SCIPsetRelaxSolVals(), SCIPsetRelaxSolValsSol(), and SCIPmarkRelaxSolValid()
2333 receive an additional argument "relax" to store the relaxation handler as creator of the relaxation solution.
2334
2335- LP:
2336 - SCIProwGetOriginCons() now returns a SCIP_CONS* instead of a SCIP_CONSHDLR*, use SCIProwGetOriginConshdlr() for the previous behavior
2337 - SCIPcreateRowCons() and SCIPcreateEmptyRowCons() now expect a SCIP_CONS* instead of a SCIP_CONSHDLR*,
2338 use SCIPcreateRowConshdlr() and SCIPcreateEmptyRowConshdlr(), respectively, for the previous behavior
2339 - deprecated SCIPlpfeastol() and SCIPchgLpfeastol(), use SCIPgetLPFeastol() and SCIPsetLPFeastol()
2340
2341- new parameter "divecontext" for every method that queries statistics for a diveset. The context can be used to distinguish between the dive set
2342 as single (standalone) heuristic or within Adaptive Diving.
2343- new parameters "divecontext" and "iterlim" to SCIPperformGenericDivingAlgorithm() to control in which context (single,adaptive) statistics are updated.
2344
2345- SCIPcopyVars, SCIPcopy, SCIPcopyConsCompression, and SCIPgetVarCopy do not copy variables that are marked as relaxation-only,
2346 thus it cannot be assumed anymore that each active variable from the master SCIP also has a counterpart in the copy.
2347 SCIPcopy, SCIPcopyConsCompression, and SCIPcopyConss can now return *valid=TRUE if some non-checked and non-enforced constraints were not copied, e.g.,
2348 because they involved relaxation-only variables. Thus, a copy is already regarded as valid if all checked or enforced constraints were copied successfully.
2349
2350- linking constraint handler:
2351 - changed type of vals argument from int* to SCIP_Real* in SCIPcreateConsLinking() and SCIPcreateConsBasicLinking()
2352 - SCIPgetIntvarLinking() has been renamed to SCIPgetLinkvarLinking().
2353 - changed return value of SCIPgetValsLinking() from int* to SCIP_Real*.
2354 - new method SCIPgetBinvarsDataLinking().
2355
2356- SCIPbendersCheckSubproblemOptimality() now returns a boolean indicating whether the subproblem is optimal or not.
2357 Previously this result was returned through a parameter. The change was required to facilitate the integration with
2358 the UG framework.
2359
2360- deleted SCIPcombineTwoInt(), SCIPcombineThreeInt(), SCIPcombineFourInt(); use the appropriate SCIPhashTwo(), ..., SCIPhashSeven() method instead
2361
2362- SCIPsetupBendersSubproblem takes a parameter of the enforcement type.
2363
2364- SCIPcreateNlpiProb takes a hashmap to store the map between the nlrows and the index in the nlrow array.
2365
2366### New API functions
2367
2368- SCIPallowWeakDualReds() and SCIPallowStrongDualReds() replace the deprecated SCIPallowObjProp() and SCIPallowDualReds(), respectively
2369- methods have been added to facilitate the transfer of Benders' decomposition cuts between solvers in UG. These include
2370 SCIPapplyBendersStoredCuts(), SCIPbendersGetNStoredCuts(), SCIPbendersGetStoredCutData() and
2371 SCIPbendersGetStoredCutOrigData().
2372- added SCIPisConvexAbspower()
2373- new functions SCIPsolGetType(), SCIPsolGetRelax(), SCIPsolSetRelax(), SCIPsolSetLPRelaxation(), SCIPsolSetStrongbranch(),
2374 SCIPsolSetPseudo to set or query the new type attribute of a primal solution. The type attribute gives information
2375 about the origin of the solution, ie, whether it was created by a relaxation handler,
2376 by the LP relaxation, by strong branching, by the current pseudo solution, or by a primal heuristic.
2377 The meaning of the argument 'heur' in all creation methods for primal solutions such as SCIPcreateSol() stays unchanged.
2379- new API functions SCIPsetCommonSubscipParams(), SCIPtranslateSubSol(), and SCIPtranslateSubSols() shared by several Large Neighborhood Search heuristics.
2380- new API function SCIPgetLPDualDegeneracy() to get two measures for the dual degeneracy of the current LP
2381- new API functions SCIPdivesetIsAvailable() to check preconditions of a dive set and SCIPdivesetIsPublic() to check if the dive set can be used
2382 by other primal heuristics.
2383- new API functions SCIPcomputeOrbitsSym(), SCIPcomputeOrbitsFilterSym(), SCIPgetPropertiesPerm(), SCIPdetermineBinvarAffectedSym(),
2384 SCIPdetermineNVarsAffectedSym(), SCIPcomputeComponentsSym(), and SCIPextendSubOrbitope(), SCIPgenerateOrbitopeVarsMatrix() for symmetry computations
2385- new API functions SCIPvarIsRelaxationOnly() and SCIPvarMarkRelaxationOnly() to query and set, resp., whether a variable is marked as relaxation-only
2386- new API functions SCIPconshdlrGetNUpdateConss() and SCIPconshdlrGetUpdateConss(), for expert users only
2387- new API function SCIPgetNConflictDualproofsApplied()
2388- new API functions SCIPeventGetOldtype() and SCIPeventGetNewtype() for the new event when changing the variable type
2389- new API function SCIPisConvexConsQuadratic() to check whether a quadratic constraint is convex when a given set of variables would be fixed
2390- new API functions SCIPgetLPFeastol(), SCIPsetLPFeastol(), and SCIPresetLPFeastol() to get, set, and reset (to the default), respectively, the primal
2391 feasibility tolerance for the LP relaxation
2392- new API functions SCIPcleanupConss{Linear,Varbound,Setppc,Logicor,Knapsack}() to clean up inactive variables from those types of linear constraints
2393- new API function SCIPsetBendersSubproblemComp() used to add a custom comparison method for ordering the Benders'
2394 decomposition subproblem solves. The comparison method is used to help with load balancing.
2395- new API function SCIPgetRowObjParallelism to get the objective parallelism of a row
2396- new API function SCIPcolGetAge to get the age of a column
2397- added SCIPhashThree(), SCIPhashFive(), SCIPhashSix(), and SCIPhashSeven() that complement SCIPhashTwo(), SCIPhashFour() to combine 32bit integers to
2398 a 32bit hash value
2399- new API function SCIPgenerateAndApplyBendersOptCut is used to generate a Benders' optimality cut using the dual
2400 solutions. This function can be supplied vectors for the primal and dual solution for generating an optimality cut.
2401 This avoids the need for a SCIP instance to solve the Benders' decomposition subproblem and generating cuts.
2402- new API function SCIPconsAddCoef used for adding a coefficient to a linear-type constraint.
2403- new API functions SCIPconsNonlinearGetRhs, SCIPconsNonlinearGetLhs and SCIPconsNonlinearAddLinearCoef for getting the
2404 RHS and LHS from a nonlinear-type constraint and adding a linear coefficient to the constraint.
2405- new API function SCIPbendersSolSlackVarsActive for checking whether any slack variables from the feasibility phase are
2406 active in the subproblem solution.
2407- new API functions SCIPbendersSetSubproblemType and SCIPbendersGetSubproblemType sets and gets the subproblem type.
2408 This is either:
2409 - Convex constraints with continuous variables
2410 - Convex constraints with discrete variables
2411 - Non-convex constraints with continuous variables
2412 - Non-convex constraints with discrete variables
2413- new API functions SCIPbendersSetSubproblemIsNonlinear() and SCIPbendersSubproblemIsNonlinear() for setting and
2414 identifying whether the Benders' decomposition subproblems contain nonlinear constraints. Similarly, the functions
2415 SCIPbendersSetMasterIsNonlinear() and SCIPbendersMasterIsNonlinear() sets and identifies whether the Benders'
2416 decomposition master problem contains nonlinear constraints.
2417- new API function SCIPapplyBendersDecomposition for applying Benders' decomposition given a decomposition in the DEC
2418 format
2419- new API function SCIPwasNodeLastBranchParent to query if a node has been the parent of the most recent branching in the tree
2420- new API functions SCIPtreemodelInit(), SCIPtreemodelFree(), SCIPtreemodelIsEnabled(), SCIPtreemodelSelectCandidate() related to the new
2421 treemodel way of comparing branching candidates. These functions are only currently used for reliability pscost branching, but they can be used
2422 in other parts of the code.
2423- New function SCIPcalcChildEstimateIncrease() to compute the increase in the child estimation
2424- new API functions SCIPisOrbitalfixingEnabled() and SCIPgetSymmetryNGenerators() to check whether orbital fixing is
2425 enabled and to get the number of generators of the current symmetry group, respectively
2426- new API function SCIPdelNlRow() to remove a row from the NLP
2427
2428### Event system
2429
2430- new event type SCIP_EVENTTYPE_NODEDELETE to react on nodes that are about to be deleted from the tree
2431
2432### Changed parameters
2433
2434- renamed parameter "propagating/orbitalfixing/enableafterrestart" to ".../symmetry/recomputerestart"
2435- Parameter "misc/allowdualreds" is now called "misc/allowstrongdualreds"
2436- Parameter "misc/allowobjprop" is now called "misc/allowweakdualreds"
2437- changed default values of propagation (new value: 1, old value: 5) and separation frequency (new value: -1, old value: 5) in cons_orbitope.c
2438- all primal heuristics that use sub-SCIPs are disabled within the heuristics fast emphasis setting
2439- deleted parameter heuristics/localbranching/useuct, use heuristics/useuctsubscip instead
2440- changed default value of "presolving/symbreak/detectorbitopes" (new value: TRUE, old value: FALSE)
2441- extended range of "misc/usesymmetry" (new range: [0,3], old range: [0,2])
2456- changed default value of "heuristics/coefdiving/freq" (old: 10, new: -1)
2457- changed default value of "heuristics/conflictdiving/freq" (old: -1, new: 10)
2458- changed default value of "heuristics/conflictdiving/lockweight" (old: 1.0, new: 0.75)
2459- replaced parameter "numerics/lpfeastol" by "numerics/lpfeastolfactor" to specify which factor should be applied to the SCIP feasibility
2460 tolerance to initialize the primal feasibility tolerance of the LP solver
2461- enabling aggressive presolving now activates all available presolving plugins,
2462 and decreases the presolving/restartfac parameter correctly with respect to default.
2463- changed default value of heuristics/rins/nodesquot to 0.3 (was 0.1), to compensate
2464 the removal of a hard coded factor of 3.0 in the code without affecting the default behavior
2465 of the RINS heuristic.
2466- changed default value of "constraints/quadratic/empathy4and" (old: 0, new: 2)
2467- changed default value of "propagating/redcost/useimplics" (old: TRUE, new: FALSE)
2468
2469### New parameters
2470
2471- the possibility to define the Benders' decomposition auxiliary variables as implicit integer is provided. This
2472 behavior is controlled with an additional parameter in the Benders' decomposition framework.
2473- added parameter benders/<bendersname>/cutcheck to enable the generation of Benders' decomposition cuts during solution
2474 checking.
2475- constraints/orbitope/usedynamicprop: the possibility to propagate orbitope constraints by reordering the rows based on the
2476 branching strategy is provided (only possible for non-model constraints)
2477
2478- new parameters heuristics/shiftandpropagate/minfixingratelp and heuristics/locks/minfixingratelp to stop the
2479 heuristics after propagating integer fixings if no sufficient fixing of the all variables (including continuous)
2480 could be achieved. These parameters help to avoid solving LP's that are comparable in hardness to the main root LP.
2481
2482- Added parameters branching/midpull and branching/midpullreldomtrig to control by how much to move
2483 the branching point for an external branching candidate closer to the middle of the candidates domain. The
2484 default of 0.75 and 0.5, respectively, uses a point that is 75*alpha% closer to the middle of the domain, where
2485 alpha is the relative width of the candidates domain (width of local domain divided by width of global domain),
2486 if the latter is below 0.5, and alpha=1.0 otherwise. That is, with the default settings, a branching point is
2487 chosen closer to the middle of the candidates domain if the variables local domain is still similar to its
2488 global domain, but is chosen closer to the LP solution if the local domain is much smaller than the global
2489 domain.
2490
2491- Added parameter lp/minmarkowitz to set the Markowitz stability threshold (range 0.0001 to 0.9999).
2492 High values sacrifice performance for stability.
2493
2494- Added parameters benders/<bendersname>/lnsmaxcalls and benders/<bendersname>/lnsmaxcallsroot to the Benders'
2495 decomposition core. These parameters limit the number of Benders' decomposition subproblem checks, for the full
2496 branch-and-bound tree and root node respective, when solving the auxiliary problem of LNS hueristics. These
2497 parameters only have effect if the lnscheck parameter is set to TRUE.
2498
2499- Added parameter cons/linear/maxmultaggrquot to limit the maximum coefficient dynamism of an equation on which
2500 multiaggregation is performed. This replaces a compiler define of the same name.
2501 Default value is 1000, smaller values make multiaggregations numerically more stable.
2502- new global parameter heuristics/useuctsubscip that affects all LNS heuristics using common sub-SCIP parameters
2503
2504- new parameter branching/relpscost/degeneracyaware to switch degeneracy-aware hybrid branching
2505
2506- new parameter separation/rapidlearning/checkexec to check whether rapid learning is allowed to run locally
2507- new parameters separation/rapidlearning/check{degeneracy,dualbound,leaves,nsols,obj} to enable checking the respective feature for local rapid learning
2508- new parameter separation/rapidlearning/maxcalls to limit the number of rapid learning executions
2509- new parameter separation/rapidlearning/nwaitingnodes to set the number of waiting nodes before the dual bound is checked
2510- new parameter separation/rapidlearning/mindegeneracy to set the minimal threshold of degenerate basic-variables
2511- new parameters separation/rapidlearning/minvarconsratio to set the minimal ratio of unfixed variables in relation to basis size
2512
2513- new parameters to control the Benders' decomposition two-phase method.
2514 - constraints/benderslp/depthfreq: after the maxdepth is reached, then the two-phase method will only be called at
2515 nodes at a depth divisible by depthfreq.
2516 - constraints/benderslp/stalllimit: after the maxdepth is reached, if there has been no improvement in the dual bound
2517 for stalllimit number of nodes, then the two-phase method is executed for the next fractional LP solution that is
2518 encountered.
2519 - constraints/benderslp/iterlimit: after the root node, only iterlimit fractional LP solutions are used at each node
2520 to generate Benders' decomposition cuts.
2521
2522- new parameters for symmetry handling
2523 - new parameter "propagating/symmetry/maxgenerators"
2524 - new parameter "propagating/symmetry/checksymmetries"
2525 - new parameter "propagating/symmetry/displaynorbitvars"
2526 - new parameter "propagating/symmetry/conssaddlp"
2527 - new parameter "propagating/symmetry/addsymresacks"
2528 - new parameter "propagating/symmetry/detectorbitopes"
2529 - new parameter "propagating/symmetry/addconsstiming"
2530 - new parameter "propagating/symmetry/ofsymcomptiming"
2531 - new parameter "propagating/symmetry/performpresolving"
2532 - new parameter "propagating/symmetry/recomputerestart"
2533 - new parameter "constraints/symresack/checkmonotonicity"
2534 - new parameter "propagating/symmetry/compresssymmetries"
2535 - new parameter "propagating/symmetry/compressthreshold"
2536 - new parameter "propagating/symmetry/disableofrestart"
2537 - new parameter "propagating/symmetry/symfixnonbinaryvars"
2538
2539- new parameter for enabling shared memory parallelisation for solving Benders' decomposition subproblems. The parameter
2540 benders/<bendersname>/numthreads sets the number of threads used for parallel subproblem solving.
2541
2542- new parameters to control enhancements for solving MINLPs by Benders' decomposition
2543 - benders/<bendersname>/execfeasphase: enables the feasibility phase for solving the Benders' decomposition
2544 subproblems
2545 - benders/<bendersname>/slackvarcoef: the initial coefficient of the slack variable for the feasibility phase
2546 - benders/<bendersname>/checkconsconvexity: should the constraints be checked for convexity. This can be set to FALSE
2547 if you are certain that the NLP subproblem is convex.
2548
2549- new parameter presolving/clqtablefac (default value 2.0) as limit on number of entries in clique table relative to number of problem nonzeros
2550
2551- new parameter conflict/uselocalrows (default: TRUE) to incorporate locally valid cuts / rows for dual proof analysis
2552
2553- new return code SCIP_NOTIMPLEMENTED for functions, e.g., in the LPI that have not been implemented (yet)
2554
2555- new parameter separating/cgmip/genprimalsols that allows to generate initial primal solutions from Gomory cuts
2556
2557- new parameter branching/relpscost/filtercandssym to allow filtering from orbits
2558
2559- new parameter branching/relpscost/transsympscost to transfer pseudo cost information to orbit
2560
2561- new parameters for tree size estimation and restarts:
2562 - estimation/restarts/restartpolicy (default value n)
2563 - estimation/method (default value c)
2564 - estimation/restarts/restartlimit (default value 1)
2565 - estimation/restarts/minnodes (default value 1000)
2566 - estimation/restarts/countonlyleaves (default value FALSE)
2567 - estimation/restarts/restartfactor (default value 2)
2568 - estimation/coefmonoprog (default value 0.3667)
2569 - estimation/coefmonossg (default value 0.6333)
2570 - estimation/restarts/hitcounterlim (default value 50)
2571 - estimation/reportfreq (default value -1)
2572 - estimation/regforestfilename (default value "-")
2573 - estimation/completiontype (default value a)
2574 - estimation/treeprofile/enabled (default value FALSE)
2575 - estimation/treeprofile/minnodesperdepth (default value 20)
2576 - estimation/useleafts (default value TRUE)
2577 - estimation/ssg/nmaxsubtrees (default value -1)
2578 - estimation/ssg/nminnodeslastsplit (default value 0)
2579
2580- new parameter constraints/linear/extractcliques to turn clique extraction off
2581
2582- new emphasis setting emphasis/numerics to increase numerical stability of (mostly) presolving operations such as (multi-)aggregations at the cost of performance.
2583
2584- new parameters for treemodel:
2585 - new parameter branching/treemodel/enable to enable the treemodel in reliability pscost branching and possible
2586 future parts of the code where it could be used.
2587 - new parameter branching/treemodel/highrule to specify which branching rule to use when treemodel thinks the node
2588 is high in the tree.
2589 - new parameter branching/treemodel/lowrule to specify which branching rule to use when treemodel thinks the node
2590 is low in the tree.
2591 - new parameter branching/treemodel/height to specify at which (estimated) height a node is high or low in the tree.
2592 - new parameter branching/treemodel/filterhigh to specify whether to filter dominated candidates in nodes which are
2593 high in the tree.
2594 - new parameter branching/treemodel/filterlow to specify whether to filter dominated candidates in nodes which are
2595 low in the tree.
2596 - new parameter branching/treemodel/maxfpiter to specify the maximum number of fixed-point iterations to use when
2597 computing the ratio of a variable using the fixed-point method.
2598 - new parameter branching/treemodel/maxsvtsheight to specify the maximum height to compute the SVTS score exactly
2599 before approximating it using the ratio.
2600 - new parameter branching/treemodel/fallbackinf defines the fallback strategy to use when the tree size estimates
2601 obtained by SVTS are infinite.
2602 - new parameter branching/treemodel/fallbacknoprim defines the fallback strategy to use when no primal bound is known
2603 and thus SVTS would not be able to compute a tree size (it would be infinite).
2604 - new parameter branching/treemodel/smallpscost defines the value under which pscosts are considered too small to be
2605 the deciding factor for branching, in which case it may be better not to use the treemodel.
2606
2607- new parameters for symmetry handling constraint handlers to enforce that also non-model constraint are copied:
2608 - new parameter "constraints/orbisack/forceconscopy"
2609 - new parameter "constraints/orbitope/forceconscopy"
2610 - new parameter "constraints/symresack/forceconscopy"
2611
2612
2613### Data structures
2614
2615- small changes in constants of hash functions
2616- added fast 2-universal hash functions for two to seven 32bit elements with 32bit output
2617- extended SCIPpqueueCreate() by additional callback argument SCIP_DECL_PQUEUEELEMCHGPOS to catch position changes
2618- new methods SCIPpqueueDelPos() to delete elements at a specific position in the priority queue and SCIPpqueueFind() to
2619 find a specific position. It is recommended to track position changes using the new callback SCIP_DECL_PQUEUEELEMCHGPOS.
2620 In contrast, using SCIPpqueueFind() can be slow because it needs to compare the element it searches for
2621 with each slot in the queue.
2622
2623### Build system
2624
2625- The default value for DFLAGS in the non-cmake buildsystem has changed from -MM to -MMD. This will break the
2626 generation of depend.* files if that was done by a compiler call that relied on -MM. The new preferred way
2627 to handle compilation dependencies is to additionally use $(DFLAGS) when compiling the object files (.o) and
2628 to include the generated .d files in the Makefile, see also "Build system / Makefile" below.
2629
2630Unit tests
2631----------
2632 - new unit test for treemodel.
2633
2634Testing
2635-------
2636
2637 - fixed an issue that may have lead to wrong status reports in the evaluation scripts
2638
2639Build system
2640------------
2641
2642### Cmake
2643
2644- avoid problem with doubly defined object together with CPLEX
2645
2646### Makefile
2647
2648- Removed static object compilation dependency files (depend.*). If using a GCC compatible compiler, then dependency
2649 files are now dynamically created and updated during build. The new dependency files (*.d) reside next to each object
2650 file (.o) in the corresponding obj subdirectory.
2651- added support for building against Ipopt >= 3.13
2652- unify compiler switches for Intel compiler and avoid problem with doubly defined object together with CPLEX
2653
2654Fixed bugs
2655----------
2656
2657- fix and improve memory handling in symmetry computation
2658- fix shown number of applied conflicts in solving statistics
2659- fix wrongly skipping strong branching call and using old information if LP was solved with 0 iterations
2660- fix minor bug in cut score calculation
2661- fixed several bugs related to rounding locks of variables not being updated correctly
2662- small fix in cons_varbound.c to skip changing bounds of multi-aggregated variables in separation callback
2663- fixed issue in SCIPtightenVar* and SCIPinferVar* that occurs for small bound changes
2664- fixed rejecting minimal boundchange that changed sign of variable, even though SCIPisLb/UbBetter approved it
2665- fixed issue in generateCutNonConvex() which is triggered when adding quadratic constraints during the solving process
2666- fixed bug in freeing the reoptimization data if no problem exists
2667- fixed bug in SCIPreoptReleaseData() when freeing all stored constraints
2668- fixed bug when freeing the transformed problem via interactive shell if reoptimization is enabled
2669- fixed two issues related to (near-)redundant logicor constraints in presolving
2670- fixed counting of aggregations in XOR constraint handler
2671- fixed handling of unbounded solutions
2672- fixed update of LP size information when an LP error occured during probing
2673- handle special case of variable bound constraints during aggregating variables
2674- tighten sides of linear constraints before trying to upgrade them to more specialized constraints (knapsack, logic-or etc.) when calling SCIPupgradeConsLinear()
2675- fixed an issue in repair heuristic in the case of loose (noncolumn) variables
2676- allow user to correctly set heuristics/alns/(un)fixtol
2677- fixed an issue in heur_completesol which is triggered during bound widening of unbounded continuous variables
2678- fixed bug in cons_indicator if addopposite is true
2679- fixed bug in sepa_disjunctive: treat case that conflictgraph is empty
2680- added safety check in conversion to rational number to avoid overflow
2681- fixed bug in interval evaluation with power-operator in certain situations
2682- fixed behavior of SCIPmatrixCreate() regarding memory management and column generation
2683- SCIPmatrixCreate() returns complete=FALSE when locks do not add up
2684- fixed bug in sepa_oddcylce when variables are fixed
2685- fixed numerical issues related to tighter constraint sides in varbound constraint handler
2686- fixed update of watchedvars in logicor constraint handler in case of a restart during the tree
2687- fixed treatment of multi-aggregated variables in logicor constraint handler
2688- handle special case of redundant implications
2689- fixed numerical issue related to almost-0-values in pseudosolution conflict analysis
2690- fixed numerical issue related to very large greatest common dividers in linear constraint handler
2691- avoid using implications on multiaggregated variables when propagating implications
2692- fixed creation of (Lagrangian) variable bounds in the OBBT propagator
2693- fixed sorting of primal solutions
2694- fixed cleaning of clean buffer in conflict analysis
2695- avoid probing on variables with huge bounds in shift and propagate heuristic
2696- fix issue in printing solutions for variables that have been added by the dual sparsify presolver
2697- fix issue related to fixing redundant logic-or constraints after presolving
2698- fixed bug when parsing logic-or and and-constraints
2699- fixed wrong assert in updateLazyBounds()
2700- fixed bug in pricestore, which resulted in too many problem variables being added
2701- fixed bug in cons_knapsack where weight of clique was not reset after an infeasibility was detected
2702- fixed bug in presol_inttobinary which did not take into account that the aggregation could be rejected due to numerics
2703- fixed bug in debug solution mechanism in connection to variables created by presol_inttobinary
2704- fixed wrong indexing while undoing the implications from a redundant variable in SCIPshrinkDisjunctiveVarSet
2705- redundancy checks in SCIPnodeAddBoundinfer now take a possible change to an active variable into account
2706- fixed adding already added quadratic rows to NLP relaxation during solve
2707- fixed issue related to variable locks in the varbound constraint handler
2708- fixed bug in the quadratic constraint handler when changing infinite constraint sides
2709- fixed sorting of variables in linear constraint handler
2710- added additional checks to ensure numerical stability of dual proofs
2711- fixed a case when activities of a linear constraint got unreliable but where still used for reductions
2712- ensure that lhs <= rhs for linear constraints (without tolerances)
2713- make handling of read errors in SCIPfread() consistent between version with and without ZLIB
2714- correctly drop variable events in cons_indicator in restart
2715- fixed bug in cons_orbitope with upgrading of orbitope constraints
2716- additional checks in some presolvers for time limit being exceeded
2717- fixed bug in presolving of cons_varbound with multi-aggregated variables
2718- improve numerics in conflict analysis by using double-double arithmetic
2719- fixed bound acceptance condition to avoid inconsistencies
2720- fixed numerics in pseudoobj propagator by using double-double arithmetic
2721
2722Miscellaneous
2723-------------
2724
2725- modified display column for memory usage ("mem"), which reports the memory usage most of the time, but shows the creator name
2726 (heuristic, relaxation handler, LP relaxation, strong branching, pseudo solution) of every new incumbent solution. Together with this change,
2727 heuristic display characters have been unified to represent the type of the heuristic
2728 (diving, Large neighborhood search, propagation, etc.), see also type_heur.h.
2729- added assert that ensures that the locks of a variable have been decreased to 0 when it is freed
2730- added more output for completing a partial solution
2731- checks in debug mode that clean buffer memory is really clean when being freed are now disabled by default
2732- don't compute symmetries if reoptimization is enabled
2733- prefer integral values when fixing an almost-fixed continuous variable in the trivial presolver
2734- changed the name of the variable that is added by the OSiL reader to represent the quadratic or nonlinear parts of the objective function
2735- SCIP_EXPORT is now defined as __attribute__((__visibility__("default"))) if GCC and no SCIP config header is used
2736
2737@page RN6 Release notes for SCIP 6
2738
2739@section RN602 SCIP 6.0.2
2740*************************
2741
2742Features
2743--------
2744
2745- The abspower constraint handler now passes more accurate convexity
2746 information to the NLP relaxation.
2747
2748Examples and applications
2749-------------------------
2750
2751- added parsing functionality for optcumulative constraints in CIP format
2752
2753Interface changes
2754-----------------
2755
2756### Interfaces to external software
2757
2758- Updated the Mosek LP solver interface to support Mosek 9.0.
2759
2760Build system
2761------------
2762
2763### Cmake
2764
2765- new target to 'doc' to build documentation
2766- ctests now fail if parameter file not found
2767- add flag STATIC_GMP and improve GMP find module
2768- remove non-API methods from library (API methods use new macro SCIP_EXPORT)
2769- increase minimal required CMake version to 3.3
2770- correct paths and dependency information when installing SCIP
2771
2772Fixed bugs
2773----------
2774- fixed SCIP-Jack presolving bug that could lead to wrong results for Steiner arborescence problems
2775- fixed wrong unboundedness result in case not all constraints were already in the LP and enforcement was
2776 skipped because an optimal solution was found
2777- fixed wrong enforcement of constraints in the disjunction constraint handler
2778- fixed wrong behavior of concurrent solve ignoring initial solutions
2779- fixed bug in concurrent solve when problem was already solved in presolving
2780- aggregate non-artificial integer variable for XOR constraints with two binary variables and delete constraint
2781- copy the objective offset when copying the original problem
2782- fixed bug in SCIPlpiGetBInvARow in lpi_cpx using wrong size of resulting vector
2783- fixed quadratic runtime behavior in sepa_aggregation
2784- fixed statistics of separators
2785- improve numerical stability in varbound constraint handler by using double-double arithmetic
2786- fixed bug in propagation of dual proofs
2787- fixed bugs that arise for multiaggregated indicator variables by disallowing multiaggregation for them
2788- improve numerical stability in SCIPcomputeBilinEnvelope* by using double-double arithmetic
2789- fixed bug related to releasing pending bound changes in tree.c
2790- set STD FENV_ACCESS pragma to on in code that changes floating-point rounding mode
2791- disable GCC optimizations in main interval arithmetic code to prevent wrong optimizations
2792- fixed wrong assert in cons_xor concerning the variable type
2793- fixed different behavior of SCIPisLbBetter and SCIPisUbBetter between having NDEBUG defined or not
2794- correctly handle bound disjunctions in symmetry detection
2795- fixed issue in reliability branching related to the LP error flag not being reset
2796- fixed treatment of near-infinite bounds in shiftandpropagate's problem transformation
2797- fixed handling of infinite values in SCIPcomputeHyperplaneThreePoints()
2798- fixed comparisons of infinite values in heur_intshifting.c and heur_shifting.c
2799- fixed bug related to updating unprocessed cuts in the cutpool
2800- fixed bug related to enabling quadratic constraints during `CONSINITLP`
2801- add missing SCIP_EXPORT for functions used by GCG
2802- fixed memory leak and wrong initialization for trival cases in cons_symresack.c
2803- fixed bug with upgrading to packing/partitioning orbitopes
2804- fixed bug with the status while upgrading in presol_symbreak.c
2805- fixed wrong stage while clearing the conflict store
2806- fixed behavior of SCIPfixVar() by setting infeasible pointer to TRUE if fixval lies outside variable domain
2807- allow tightenVar() in SCIP_STAGE_PROBLEM stage
2808- fixed bug in cumulative constraint handler when separating the LP solution
2809- fixed issues with integer overflow in cumulative constraint handler
2810- fixed bug where the convexity of Benders' decomposition subproblems was checked even when users defined subproblem
2811 solving methods. Now, as per the documentation, the user must explicitly state whether the subproblem is convex
2812- fixed wrong indexing in heur_dualval
2813- fixed issue with basis status in SoPlex LPi
2814
2815Miscellaneous
2816-------------
2817
2818- statistics now output primal/dual bounds if objective limit is reached
2819- memory check in debug mode is now disabled by default
2820- message is now provided to the user to inform that automatic Benders' auxiliary variable lower bound computations are
2821 not activated when user defined subproblem solving methods are present
2822- corrected documentation of the primalgap in SCIP; describe when it will be infinite
2823
2824@section RN601 SCIP 6.0.1
2825*************************
2826
2827Features
2828--------
2829
2830- when using a debug solution every (multi-)aggregation will be checked w.r.t. this solution
2831
2832Performance improvements
2833------------------------
2834
2835- try greedy solution first before solving knapsack exactly using dynamic programming in SCIPsolveKnapsackExactly,
2836 compute greedy solution by weighted median selection.
2837- don't consider implied redcost by default in the reduced cost propagator
2838
2839Interface changes
2840-----------------
2841
2842### Deleted and changed API methods and macros
2843
2844- The preprocessor macro NO_CONFIG_HEADER now needs to be defined when
2845 including SCIP header files from a SCIP build or installation that
2846 has been build via the Makefile-only build system.
2847
2848- The following preprocessor macros have been renamed:
2849 WITH_ZLIB to SCIP_WITH_ZLIB, WITH_GMP to SCIP_WITH_GMP, WITH_READLINE
2850 to SCIP_WITH_READLINE, NO_SIGACTION to SCIP_NO_SIGACTION, NO_STRTOK_R
2851 to SCIP_NO_STRTOK_R, ROUNDING_FE to SCIP_ROUNDING_FE, ROUNDING_FP to
2852 SCIP_ROUNDING_FP, ROUNDING_MS to SCIP_ROUNDING_MS. Note, however, that
2853 the names of macros NO_RAND_R and NO_STRERROR_R have not been changed
2854 so far.
2855
2856### New API functions
2857
2858- SCIPhashmapInsertInt(), SCIPhashmapSetImageInt(), and SCIPhashmapGetImageInt() to use integer values as images in hashmaps
2859
2860### Command line interface
2861
2862- warn about coefficients in MPS files with absolute value larger than SCIP's value for infinity
2863
2864### Changed parameters
2865
2866- default clock type for timing is now wallclock
2867
2868Unit tests
2869----------
2870
2871- added unit tests for exact knapsack solving and (weighted) median selection algorithms
2872
2873Build system
2874------------
2875
2876### Cmake
2877
2878- add missing GMP dependency when compiling with SYM=bliss
2879- add DL library when linking to CPLEX to avoid linker errors
2880- new config.h header defining the current build configuration, e.g. SCIP_WITH_GMP
2881
2882Fixed bugs
2883----------
2884
2885- fixed handling of weights in cons_sos1 and cons_sos2 (NULL pointer to weights)
2886- fixed handling of unbounded LPs in SCIP and in several LPIs; added heuristic method to guess solution
2887- the STO reader is capable of handling scenarios defined using lower case "rhs"
2888- fixed OPB reader for instances without explicit plus signs
2889- correct dual solution values for bound constraints
2890- fixed recognition of variable with only one lock in cons_bivariate, cons_quadratic, and cons_nonlinear
2891- fixed update of constraint violations in solution repair in cons_bivariate, cons_quadratic, and cons_nonlinear
2892- print error message and terminate if matrix entries of a column are not consecutive in mps format
2893- fixed incorrect handling of fixed variables when transfer of cuts from LNS heuristic for Benders' decomposition
2894- fix returning local infeasible status by Ipopt interface if Ipopt finds problem locally infeasible
2895- skip attempt to apply fixings in linear constraint handler during solving stage as LP rows cannot change anymore
2896- fixed bug when reading >= indicator constraints in MPS format
2897- fix issue with nodes without domain changes if we ran into solution limit in prop_orbitalfixing
2898- fixed unresolved reference to CppAD's microsoft_timer() function on builds with MS/Intel compilers on Windows
2899- ignore implications added through SCIPaddVarImplication() that are redundant to global bounds also in the
2900 special case of an implication between two binary variables; also, use implications instead of cliques in the case
2901 of a binary implied variable with nonbinary active representative
2902- fixed bug with aggregated variables that are aggregated in propagation of cons_sos1
2903- fixed some special cases in SCIPselect/SCIPselectWeighted methods
2904- relaxed too strict assertion in Zirounding heuristic
2905- fixed the upgrade routine to XOR constraints: aggregate integer variable if its coefficient has the wrong sign
2906- fixed handling of nonartificial parity variables when deleting redundant XOR constraints
2907- earlier deletion of trivial XOR constraints (at most 1 operator left)
2908- fixed wrong hashmap accesses and added sanity check for the correct hashmap type
2909- avoid copying of unbounded solutions from sub-SCIPs as those cannot be checked completely
2910- corrected the output of the first LP value in case of branch-and-price
2911- fixed possible integer overflow, which led to wrong conclusion of infeasibility, in energetic reasoning of cons_cumulative.c
2912
2913Miscellaneous
2914-------------
2915
2916- do not scale linear constraints to integral coefficients
2917
2918@section RN600 SCIP 6.0.0
2919*************************
2920
2921Features
2922--------
2923
2924- new diving heuristic farkasdiving that dives into the direction of the pseudosolution and tries to construct Farkas-proofs
2925- new diving heuristic conflictdiving that considers locks from conflict constraints
2926- restructuring of timing of symmetry computation that allows to add symmetry handling components within presolving
2927- `lp/checkstability` is properly implemented for SoPlex LPI (spx2)
2928- new branching rule lookahead that evaluates potential child and grandchild nodes to determine a branching decision
2929- limits on the number of presolving rounds a presolver (maxrounds) or propagator/constraint handler (maxprerounds)
2930 participates in are now compared to the number of calls of the particular presolving method, not the number of
2931 presolving rounds in general, anymore
2932- new miscellaneous methods for constraints that have a one-row linear representation in pub_misc_linear.h
2933- a Benders' decomposition framework has been added. This framework provides the functionality for a user to solve a
2934 decomposed problem using Benders' decomposition. The framework includes classical optimality and feasibility cuts,
2935 integer optimality cuts and no-good cuts.
2936- add statistic that presents the number of resolves for instable LPs
2937- new readers for stochastic programming problems in SMPS format (reader_sto.h, reader_smps.h)
2938
2939Performance improvements
2940------------------------
2941
2942- cuts generated from certain quadratic constraints with convex feasible region are now global
2943- performance improvements for Adaptive Large Neighborhood Search heur_alns.c
2944 + all neighborhoods now start conservatively from maximum fixing rate
2945 + new default parameter settings for bandit selection parameters
2946 + no adjustment of minimum improvement by default
2947- improved bound tightening for some quadratic equations
2948- constraint handler checking order for original solutions has been modified to check those with negative check priority
2949 that don't need constraints after all other constraint handlers and constraints have been checked
2950- deactivate gauge cuts
2951
2952Examples and applications
2953-------------------------
2954
2955- new example `brachistochrone` in CallableLibrary examples collection; this example implements a discretized model to
2956 obtain the trajectory associated with the shortest time to go from point A to B for a particle under gravity only
2957- new example `circlepacking` in CallableLibrary examples collection; this example models two problems about packing
2958 circles of given radii into a rectangle
2959- new price-and-branch application for the ringpacking problem
2960- new stochastic capacitated facility location example demonstrating the use of the Benders' decomposition framework
2961
2962Interface changes
2963-----------------
2964
2965### New and changed callbacks
2966
2967- added parameter locktype to `SCIP_DECL_CONSLOCK` callback to indicate the type of variable locks
2968
2969### Deleted and changed API methods
2970
2971- Symmetry:
2972 + removed function SCIPgetTimingSymmetry() in presol_symmetry.h since this presolver does not compute symmetries independent
2973 of other components anymore
2974 + additional argument `recompute` to SCIPgetGeneratorsSymmetry() to allow recomputation of symmetries
2975
2976- Random generators:
2977 + the seed of SCIPinitializeRandomSeed() is now an unsigned int
2978 + the seed of SCIPsetInitializeRandomSeed() is now an unsigned int and it returns an unsigned int
2979 + new parameter for SCIPcreateRandom() to specify whether the global random seed shift should be used in the creation of
2980 the random number generator
2981
2982- Miscellaneous:
2983 + additional arguments `preferrecent`, `decayfactor` and `avglim` to SCIPcreateBanditEpsgreedy() to choose between
2984 weights that are simple averages or higher weights for more recent observations (the previous default).
2985 The last two parameters are used for a finer control of the exponential decay.
2986 + functions SCIPintervalSolveUnivariateQuadExpression(), SCIPintervalSolveUnivariateQuadExpressionPositive(), and
2987 SCIPintervalSolveUnivariateQuadExpressionPositiveAllScalar() now take an additional argument to specify already
2988 existing bounds on x, providing an entire interval ([-infinity,infinity]) gives previous behavior
6083 - added writing for pip files (linear, quadratic, polynomial nonlinear, polynomial abspower, polynomial bivariate, and
6084 and constraints)
6085 - CIP format variable characters defined, e.g. `SCIP_VARTYPE_INTEGER_CHAR`
6086 - Improved support for wbo format for weighted PBO problems, IBM's xml-solution
6087 format and pip and zimpl format for polynomial mixed-integer programs
6088 - New reader for (standard) bounds on variables
6089 - Extended reader for CIP models to handle various new constraints, including all types of linear constraints
6090 - flatzinc reader is now capable to read cumulative constraints
6091 - changed opb(/wbo) reader which now creates pseudoboolean constraints instead of linear- and and-constraints, only a
6092 non-linear objective will create and-constraints inside the reader and while reading a wbo file the topcost constraint
6093 is created as well
6094 - added clock for determine the time for reading
6095 - added reader for variable bounds (reader_bnd.{c,h})
6096 - Removed method SCIPreadSol(); call solution reading via SCIPreadProb() which calls the solution reader for .sol files.
6097
6098- Nonlinear:
6099 - Major extensions for nonlinear CIP, new option for n-ary branching on nonlinear variables (within pseudocost branching rule)
6100 - added BETA version of constraint handler for nonlinear constraints (cons_nonlinear.{c,h}) to handle nonlinear
6101 equations given by algebraic expressions using operands like addition, multiplication, power, exp, log, bivariate
6102 nonlinear constraints; currently no trigonometric functions
6103 - added BETA version of constraint handler for bivariate nonlinear constraints (cons_bivariate.{c,h}) to compute tight
6104 estimators for 1-convex and convex-concave bivariate nonlinear functions (given as expression tree)
6105 - the gams writer can now write nonlinear, abspower and bivariate constraints
6106 - Extended writer for GAMS and pip format to write more types of nonlinear constraints
6107 - the pip and zimpl reader now create nonlinear constraints for polynomials of degree > 2
6108
6109- Presolving:
6110 - new dual presolving methods in cons_setppc and cons_logicor
6111 - new presolving step `removeConstraintsDueToNegCliques` in locigor constraint handler which updates logicor constraints
6112 to setppc constraints if a negated clique inside this constraint exist, by default is off
6113 - new presolving step in cons_knapsack (detectRedundantVars, deleteRedundantVars) which determines redundant variables
6114 in knapsack constraint with or without using clique information
6115 - cons_logicor is now able to replace all aggregated variables in presolving by there active or negation of an active
6116 variable counterpart
6117 - prop_pseudoobj is now working in presolving as well
6118 - implement presolving in exitpre() in cons_orbitope and cons_indicator
6119
6120- Propagators:
6121 - added counter for number calls and timing for resolve propagation calls for constraint handler and propagators
6122 - Propagators are now also called in node presolving
6123 - the probing presolver presol_probing.{c.h} is now a propagator prop_probing.{c,h}, all corresponding parameters moved as well
6124 - the redcost separator sepa_redcost.{c.h} is now a propagator prop_redcost.{c,h}, all corresponding parameters moved as well
6125 - outsourced propAndSolve() method in solve.c which calls domain propagation and solving of the lp and relaxation
6126
6127- Statistic:
6128 - solutions which are given by the user from the outside are now marked by `#` in the output
6129 - the `Solving Time` is now spitted into presolving, solving and reading time
6130 - Presolvers section has new column `AddCons` which states the number of added constraint
6131 - Constraints section has new column named \#ResProp which show the number of resolve propagation calls of certain
6132 constraint handler
6133 - Constraint Timing section has a new column \#ResProp which states the time spend in resolve propagation method of the
6134 constraint handler
6135 - improved output of propagators in display statistics
6136 - new section `Propagator Timing` which shows the time spend in different callbacks of the propagator
6137 - rearranged first two columns of Propagators section; \#Propagate and \#ResProp stating the number of call for
6138 propagation and resolve propagation; the Time column is moved into the new section Propagator Timings
6139 - Constraints section has new column named `MaxNumber` which the maximum number of active constraint of a certain
6140 constraint handler
6141 - added two columns `Time-0-It` and `Calls-0-It` in the LP section which states the number of LP call and time spend for
6142 solving LPs with zero iterations (only refactorization)
6143 - The display of statistics for presolvers, propagators, constraints and LP solving has changed.
6144
6145Performance improvements
6146------------------------
6147
6148- Reorganized filtering process of separation storage (allows adding cuts for different solutions)
6149- Improved presolving for various constraint handlers
6150- Improved propagation methods for variable bound constraints
6151- Improved performance for quadratic constraints
6152- performance improvements in prop_vbounds
6153- child selection rules now get also applied when the relaxation value is equal to the bound changed in branching
6154- added dual reduction to cons_cumulative.c
6155- for continuous variables, the pseudo costs update and the pscost branching rule now use the same strategies for
6156 updating the pseudo costs and estimating the improvement in the LP bound
6157- only perform probing if the variables are locked
6158- performance and memory consumption improvements in xmlparse.c
6159- Improved knapsack cover cuts
6160- avoid very long separation times of LEWIs in cons_knapsack for very large minimal covers
6161- used SCIPallocMemoryArray() instead of SCIPallocBlockMemoryArray() which leads to fewer memory consumption in
6162 getLiftingSequence() in cons_knapsack, also improved cache use bei using an extra array instead blockmemory chunks
6163- switched FASTMIP from 1 to 2 for CPLEX and changed default pricing rule back to steepest edge pricing instead of
6164 quickstart steepest edge pricing
6165- made sorting method more robust
6166- LNS heuristics now use SCIPcopy() by default
6167- considering inactive variables in undercover heuristic; limiting effort for solving covering problem
6168- if during probing mode the LP relaxation is solved from scratch, e.g., when calling the shiftandpropagate heuristic
6169 before root node solving, then we clear the resulting LP state, since it might be a bad starting basis for the next
6170 solve of the LP relaxation (controlled by new parameter `lp/clearinitialprobinglp`)
6171- included LP simplifier into SoPlex LP interface, applied when solving from scratch (lpi_spx.cpp)
6172- new presolving steps in varbound constraint handler, tightening bounds, coefficients, sides and pairwise presolving
6173
6174Interface changes
6175-----------------
6176
6177- Miscellaneous:
6178 - The emphasis setting types now distinguish between plugin-type specific parameter settings (default, aggressive, fast, off),
6179 which are changed by SCIPsetHeuristics/Presolving/Separating(), and global emphasis settings (default, cpsolver, easycip,
6180 feasibility, hardlp, optimality, counter), which can be set using SCIPsetEmphasis().
6181
6182### New and changed callbacks
6183
6184- added propagator timings `SCIP_PROPTIMING_BEFORELP`, `SCIP_PROPTIMING_DURINGLPLOOP` and `SCIP_PROPTIMING_AFTERLPLOOP` for
6185 all propagation callbacks (see propagators and constraint handlers) which lead to calling the propagation methods of a
6186 propagator before the lp is solved, during the lp loop and after the lp solving loop
6187
6188- Conflict Analysis:
6189 - Added parameter `separate` to conflict handler callback method SCIP_DECL_CONFLICTEXEC() that defines whether the conflict
6190 constraint should be separated or not.
6191
6192- Constraint Handler:
6193 - The new constraint handler callback SCIP_DECL_CONSDELVARS() is called after variables were marked for deletion.
6194 This method is optional and only of interest if you are using SCIP as a branch-and-price framework. That means,
6195 you are generating new variables during the search. If you are not doing that just define the function pointer
6196 to be `NULL`.
6197 If this method gets implemented you should iterate over all constraints of the constraint handler and delete all
6198 variables that were marked for deletion by SCIPdelVar().
6199
6200- NLP Solver Interface:
6201 - The callbacks SCIP_DECL_NLPIGETSOLUTION() and SCIP_DECL_NLPISETINITIALGUESS() got new parameters to get/set values of
6202 dual variables.
6203 - The callback SCIP_DECL_NLPICOPY() now passes the block memory of the target SCIP as an additional parameter.
6204
6205- Presolving:
6206 - New parameters `isunbounded` and `isinfeasible` for presolving initialization (SCIP_DECL_CONSINITPRE(),
6207 SCIP_DECL_PRESOLINITPRE(), SCIP_DECL_PROPINITPRE()) and presolving deinitialization (SCIP_DECL_CONSEXITPRE(),
6208 SCIP_DECL_PRESOLEXITPRE(), SCIP_DECL_PROPEXITPRE()) callbacks of presolvers,
6209 constraint handlers and propagators, telling the callback whether the problem was already declared to be
6210 unbounded or infeasible. This allows to avoid expensive steps in these methods in case the problem is already
6211 solved, anyway.
6212
6213 Note, that the C++ methods
6214 - scip::ObjConshdlr::scip_presol() corresponding to SCIP_DECL_CONSPRESOL()
6215 - scip::ObjConshdlr::scip_initpre() corresponding to SCIP_DECL_CONSINITPRE()
6216 - scip::ObjPresol::scip_initpre() corresponding to SCIP_DECL_PRESOLINITPRE()
6217 - scip::ObjProp::scip_initpre() corresponding to SCIP_DECL_PROPINITPRE()
6218 - scip::ObjConshdlr::scip_exitpre() corresponding to SCIP_DECL_CONSEXITPRE()
6219 - scip::ObjPresol::scip_exitpre() corresponding to SCIP_DECL_PRESOLEXITPRE()
6220 - scip::ObjProp::scip_exitpre() corresponding to and SCIP_DECL_PROPEXITPRE()
6221 are virtual functions. That means, if you are not adding the new parameters, your code will still compile, but these methods are not executed.
6222 - Propagators are now also called in during presolving, this is supported by the new callback methods SCIP_DECL_PROPINITPRE(),
6223 SCIP_DECL_PROPEXITPRE(), and SCIP_DECL_PROPPRESOL().
6224 - The new parameters `nnewaddconss` and `naddconss` were added to the constraint handler callback method SCIP_DECL_CONSPRESOL()
6225 and the presolver callback method SCIP_DECL_PRESOLEXEC(). These parameters were also added to corresponding C++ wrapper
6226 class methods (scip_presol() in objconshdlr.h and scip_exec() in objpresol.h)
6227
6228- Problem Data:
6229 - The callback SCIP_DECL_PROBCOPY() got a new parameter `global` to indicate whether the global problem or a local version is copied.
6230
6231### Deleted and changed API methods
6232
6233- implemented SCIPlpiGetPrimalRay() in SoPlex interface that has become available with SoPlex version 1.5.0.2
6234- allowed calling SCIPgetRowSolActivity() in `SCIP_STAGE_SOLVED`, since LP is still available
6235- various extensions and modifications for expressions and expression trees (too much to state here)
6236- The result value `SCIP_NEWROUND` has been added, it allows a separator/constraint handler to start a new separation round
6237 (without previous calls to other separators/conshdlrs).
6238- SCIPcalcNodeselPriority() got a new parameter `branchdir`, which defines the type of branching that was performed: upwards, downwards, or fixed.
6239
6240- Constraint Handlers:
6241 - Method SCIPincludeQuadconsUpgrade() of quadratic constraint handler got new parameter `active` to indicate whether the upgrading method is active by default.
6242 - Method SCIPseparateRelaxedKnapsack() in knapsack constraint handler got new parameter `cutoff`, which is a pointer to store whether a cutoff was found.
6243
6244- Nonlinear expressions, relaxation, and solver interface:
6245 - SCIPcreateNLPSol() now creates a `SCIP_SOL` that is linked to the solution of the current NLP relaxation
6246 - Various types and functions dealing with polynomial expressions have been renamed to use the proper terms `monomial` and
6247 `polynomial` in nonlinear expressions (nlpi/∗expr*); results in many renamings of types, structs and methods.
6248 - The methods SCIPnlpGetObjective(), SCIPnlpGetSolVals(), and SCIPnlpGetVarSolVal() have been removed, use SCIPgetNLPObjval(),
6249 SCIPvarGetNLPSol() and SCIPcreateNLPSol() to retrieve NLP solution values instead.
6250 SCIPcreateNLPSol() now returns an error if NLP or NLP solution is not available
6251 - Removed methods SCIPmarkRequireNLP() and SCIPisNLPRequired(), because the NLP is now always constructed if nonlinearities
6252 are present.
6253 - SCIPgetNLP() has been removed and NLP-methods from pub_nlp.h have been moved to scip.h, which resulted in some renamings, too.
6254 - renamed SCIPexprtreeEvalSol() to SCIPevalExprtreeSol() and now located in scip.h.
6255 - renamed SCIPexprtreeEvalIntLocalBounds() to SCIPevalExprtreeLocalBounds() and now located in scip.h.
6256 - renamed SCIPexprtreeEvalIntGlobalBounds() to SCIPevalExprtreeGlobalBounds() and now located in scip.h.
6257 - The functions SCIPnlpiGetSolution() and SCIPnlpiSetInitialGuess() got additional arguments to get/set dual values.
6258 - The method SCIPgetNLPI() got a new parameter `nlpiproblem`, which is a pointer to store the NLP solver interface problem.
6259
6260- Timing:
6261 - SCIPincludeProp() got additional parameters to set the timing mask of the propagator and the new callbacks and parameters
6262 related to calling the propagator in presolving.
6263 - SCIPincludeConshdlr() got additional parameters to set the variable deletion callback function and the timing mask for
6264 propagation.
6265 - removed parameters timelimit and memorylimit from SCIPapplyRens()
6266 - The parameters `timelimit` and `memorylimit` were removed from SCIPapplyRens().
6267
6268- Problem Data:
6269 - The method SCIPcopyProb() got a new parameter `global` to indicate whether the global problem or a local version is copied.
6270
6271- Writing and Parsing Constraints:
6272 - The methods SCIPwriteVarName(), SCIPwriteVarsList(), and SCIPwriteVarsLinearsum() got a new boolean parameter `type`
6273 that indicates whether the variable type should be written or not.
6274 - The methods SCIPparseVarName() and SCIPparseVarsList() got a new output parameter `endptr` that is filled with the position
6275 where the parsing stopped.
6276 - The method SCIPwriteVarsList() got additionally a new parameter `delimiter` that defines the character which is used for delimitation.
6277
6278- Variables:
6279 - SCIPmarkDoNotMultaggrVar()/SCIPvarMarkDoNotMultaggr() now allow to mark negated and aggregated variables
6280 - SCIPgetVarCopy() got a new parameter `success` that will be FALSE if method is called after problem creation stage and no hash map is
6281 given or no image for the given variable is contained in the given hash map.
6282 - SCIPchgVarType() got an extra boolean parameter to store if infeasibility is recognized while upgrading a variable from continuous
6283 type to an integer type.
6284 - SCIPdelVar() got a new parameter `deleted`, which stores whether the variable was successfully marked to be deleted.
6285
6286### New API functions
6287
6288- information about the quality of the solution of an LP (currently the condition number of the basis matrix) can now be:
6289 + requested from the LPI (currently only available for CPLEX): methods SCIPlpiGetRealSolQuality() and
6290 + SCIPprintLPSolutionQuality() command display lpsolquality in interactive shell display column lpcond to show
6291 + estimate on condition number, if available
6292- SCIPround() and SCIPfeasRound() to round to nearest integer
6293- SCIPsortRealRealIntInt() and corresponding sorting/inserting/deleting methods in pub_misc.h and necessary defines in misc.c
6294- SCIPsortRealIntLong(), SCIPsortPtrPtrRealInt() and corresponding sorting/inserting/deleting methods in
6295 pub_misc.h and necessary defines in misc.c
6296- SCIPcomputeLPRelIntPoint() to compute relative interior point of the current LP
6297- SCIPstartSolvingTime() and SCIPstopSolvingTime() which can be used to start or stop the solving time clock
6298- SCIPstrToRealValue() and SCIPstrCopySection() in pub_misc.h; these methods can be used to convert a string
6299 into a `SCIP_Real` value and to copy a substring.
6300- SCIPgetBinvarRepresentatives() which gets binary variables that are equal to some given binary variables,
6301 and which are either active, fixed, or multi-aggregated, or the negated variables of active, fixed, or multi-aggregated variables
6302- SCIPhasPrimalRay() and SCIPgetPrimalRayVal() that return whether a primal ray is stored and which value a
6303 given variable has in the primal ray, respectively
6304- SCIPsetConsModifiable()
6305- SCIPsetParam() which is a generic parameter setter method, independent of the parameter type
6306- SCIPpropInitpre(), SCIPpropExitpre(), SCIPpropPresol() which initializes, exists and executes the presolving phase
6307- SCIProwGetAge() to access the age of a row (pub_lp.h/lp.c)
6308- SCIPsolGetOrigObj() in pub_sol.h which returns for a solution in the original problem space the objective value
6309- SCIPretransformSol() in scip.h that allows to retransform a solution to the original space
6310- SCIPlpiClearState() to LP interfaces for clearing basis information in the LP solver
6311- SCIPgetSubscipDepth() to access the depth of the current SCIP as a copied subproblem
6312- SCIPdebugAddSolVal() and SCIPdebugGetSolVal() to add/get values to/from a debug solution
6313- SCIPsepastoreRemoveInefficaciousCuts() to remove non-efficious cuts from the separation storage
6314
6315- Nodes:
6316 - SCIPnodeGetParent() to get parent node of a node
6317 - SCIPnodesSharePath() in pub_tree.h that determines whether two nodes are on the same leaf-root path
6318 - SCIPnodesGetCommonAncestor() in pub_tree.h that finds the common ancestor node for two given nodes
6319
6320- Read and Write:
6321 - SCIPgetReadingTime() which returns the time for reading in seconds
6322 - SCIPparseVarsLinearsum(), SCIPparseVarsPolynomial() and SCIPwriteVarsPolynomial() and for writing and
6323 parsing polynomials in constraint handler writing/parsing methods
6324
6325- Memory:
6326 - SCIPcreateMesshdlrPThreads() and SCIPfreeMesshdlrPThreads() for allocating and deleting necessary memory
6327 for message handlers for parallel pthread version
6328 - SCIPallocClearMemoryArray() and BMSallocClearMemoryArray() for allocating cleared memory arrays in scip.h and memory.h
6329
6330- Intervals:
6331 - SCIPintervalPowerScalarInverse() to solve an equation y = x^p for given bounds on y and scalar exponent p
6332 - SCIPintervalQuadBivar() to compute tight bounds on a bivariate quadratic form
6333 - SCIPintervalSolveBivariateQuadExpressionAllScalar() to compute tight bounds on the solutions of a bivariate quadratic equation
6334
6335- Variables:
6336 - SCIPcomputeVarCurrent{L,U}b{Local,Global}() to compute local or global lower or upper bounds of a
6337 multiaggregated variable from the bounds of the aggregation variables
6338 - SCIPbranchVarValNary() for n-ary variable branching
6339 - SCIPgetNegatedVars() which returns all negated variables for a given array of variables, if the negated
6340 variables are not existing yet, they will be created
6341 - SCIPgetNTotalVars() that returns the total number of created vars, icluding variables that were deleted in the meantime
6342 - SCIPvarGetHashkey(), SCIPvarIsHashkeyEq(), SCIPvarGetHashkeyVal() in pub_var.h which can be used for `SCIP_HASHTABLE` of variables
6343 - SCIPvarGetNBdchgInfosLb() and SCIPvarGetNBdchgInfosUb() in pub_var.h returning the number of lower or upper bound changes on the active path
6344 - SCIPvarGetBdchgInfoLb() and SCIPvarGetBdchgInfoUb() returning the bound change information at the given position
6345 - SCIPvarMarkDeletable() to mark a variable to be deletable completely from the problem (for
6346 branch-and-price); can only be called before the variable is added to the problem
6347 - SCIPvarMarkNotDeletable() that marks a variable to be non-deleteable (used within SCIP for forbidding
6348 deletion of variables contained in solution, LP bases, (multi)aggregation, ...)
6349 - SCIPvarIsDeletable() that returns whether a variable is marked to be deletable (each variable is per default non-deletable)
6350
6351- NLP:
6352 - SCIPgetNLPVarsNonlinearity() to get for each variable in the NLP the number of NLP rows in which this variable appears in a nonlinear way
6353 - SCIPnlrowGetDualsol(), SCIPgetNLPVarsLbDualsol(), SCIPgetNLPVarsUbDualsol() to retrieve dual values from an NLP solution
6354 - SCIPgetNLPFracVars() to get the fractional variables in an NLP solution
6355
6356- Propagator:
6357 - SCIPpropSetPresolPriority() which changes the presolving priority of a given propagator