Add createOrigSol and tests

Author: Joao-Dionisio

src/pyscipopt/scip.pxd

@@ -837,6 +837,7 @@ cdef extern from "scip/scip.h":
     SCIP_Real SCIPgetSolTransObj(SCIP* scip, SCIP_SOL* sol)
     SCIP_RETCODE SCIPcreateSol(SCIP* scip, SCIP_SOL** sol, SCIP_HEUR* heur)
     SCIP_RETCODE SCIPcreatePartialSol(SCIP* scip, SCIP_SOL** sol,SCIP_HEUR* heur)
+    SCIP_RETCODE SCIPcreateOrigSol(SCIP* scip, SCIP_SOL** sol, SCIP_HEUR* heur)    
     SCIP_RETCODE SCIPsetSolVal(SCIP* scip, SCIP_SOL* sol, SCIP_VAR* var, SCIP_Real val)
     SCIP_RETCODE SCIPtrySolFree(SCIP* scip, SCIP_SOL** sol, SCIP_Bool printreason, SCIP_Bool completely, SCIP_Bool checkbounds, SCIP_Bool checkintegrality, SCIP_Bool checklprows, SCIP_Bool* stored)
     SCIP_RETCODE SCIPtrySol(SCIP* scip, SCIP_SOL* sol, SCIP_Bool printreason, SCIP_Bool completely, SCIP_Bool checkbounds, SCIP_Bool checkintegrality, SCIP_Bool checklprows, SCIP_Bool* stored)

src/pyscipopt/scip.pxi

@@ -4429,6 +4429,25 @@ cdef class Model:
         partialsolution = Solution.create(self._scip, _sol)
         return partialsolution
 
+    def createOrigSol(self, Heur heur = None):
+        """Create a new primal solution.
+
+        :param Heur heur: heuristic that found the solution (Default value = None)
+
+        """
+        cdef SCIP_HEUR* _heur
+        cdef SCIP_SOL* _sol
+
+        if isinstance(heur, Heur):
+            n = str_conversion(heur.name)
+            _heur = SCIPfindHeur(self._scip, n)
+        else:
+            _heur = NULL
+            
+        PY_SCIP_CALL(SCIPcreateOrigSol(self._scip, &_sol, _heur))
+        solution = Solution.create(self._scip, _sol)
+        return solution
+
     def printBestSol(self, write_zeros=False):
         """Prints the best feasible primal solution."""
         user_locale = locale.getlocale(category=locale.LC_NUMERIC)

tests/test_solution.py

@@ -47,6 +47,33 @@ def test_solution_create():
     assert m.getSolObjVal(s1) == -1
     m.freeSol(s1)
 
+def test_createOrigSol():
+    m = Model()
+
+    x = m.addVar("x", lb=0, ub=2, obj=-1)
+    y = m.addVar("y", lb=1, ub=4, obj=1)
+    z = m.addVar("z", lb=1, ub=5, obj=10)
+    m.addCons(x * x <= y*z)
+    m.presolve()
+
+    s = m.createOrigSol()
+    s[x] = 2.0
+    s[y] = 5.0
+    s[z] = 10.0
+    assert not m.checkSol(s)
+    assert m.addSol(s, free=True)
+
+    s1 = m.createOrigSol()
+    m.setSolVal(s1, x, 1.0)
+    m.setSolVal(s1, y, 1.0)
+    m.setSolVal(s1, z, 1.0)
+    assert m.checkSol(s1)
+
+    m.optimize()
+
+    assert m.getSolObjVal(s1) == 10.0
+    m.freeSol(s1)
+
 
 def test_solution_evaluation():
     m = Model()