Output logs from pyscipopt

[sheepone12138]

I'm working on a project where I need to output the error during the iteration when the solver is running and plot it as a curve. How should I do this and how can I output the running log as a file?

[Joao-Dionisio]

Hello, @sheepone12138! When you say error, do you mean an error in SCIP, or some sort of error in your project?

If you know what is causing your error and when it's happening (in relation to the solving process), you can try implementing an event handler. For example, suppose your error shows up when solving the LP in some node. Then, in the event handler you need to catch the event SCIP_EVENTTYPE_LPSOLVED in eventinit. You would then need to implement eventexec and plot the error you want or whatever you want to do.

Here is a list of all events available in SCIP.

You can look here for an example on using event handlers.

If you're having problem with this, you can tell me more details and I can try to help out.

To ouput the log into a file, you can try using redirectOutput(). I've never used this, so I don't know how it works, sorry :(

[mmghannam]

@Joao-Dionisio Sure yeah, but more broadly I think it can also be implemented more generally as an attach_event_handler method that takes some event types and a callback (function) to run when the event is dispatched (with the type and data passed to it as inputs). Much more ergonomic than writing a full class, which we can definitely handle internally.

[Joao-Dionisio]

Oh, super great idea!

[sheepone12138]

I tried your method and because of some bugs I changed your code to:

gap_handler = GapEventhdlr()
m.includeEventhdlr(gap_handler, "gapEventHandler", "Event handler which collects primal and dual solution evolution")

But then I ran into some mistakes that I couldn't figure out

Traceback (most recent call last):
File "C:\Users\15572\PycharmProjects\pythonProject2(scip)\sa.py", line 6, in eventinit
self.model.catchEvent(SCIP_EVENTTYPE.BESTSOLFOUND)
File "src\pyscipopt\scip.pxi", line 5014, in pyscipopt.scip.Model.catchEvent
TypeError: catchEvent() takes exactly 2 positional arguments (1 given)
Exception ignored in: 'pyscipopt.scip.PyEventInit'
Traceback (most recent call last):
File "C:\Users\15572\PycharmProjects\pythonProject2(scip)\sa.py", line 6, in eventinit
self.model.catchEvent(SCIP_EVENTTYPE.BESTSOLFOUND)
File "src\pyscipopt\scip.pxi", line 5014, in pyscipopt.scip.Model.catchEvent
TypeError: catchEvent() takes exactly 2 positional arguments (1 given)
[event.c:205] ERROR: Error <0> in function call
[set.c:5329] ERROR: Error <0> in function call
[scip_solve.c:407] ERROR: Error <0> in function call
[event.c:226] ERROR: event handler not initialized
[set.c:5450] ERROR: Error <-8> in function call
[scip_solve.c:1853] ERROR: Error <-8> in function call
[scip_solve.c:3393] ERROR: Error <-8> in function call
[scip_solve.c:2358] ERROR: Error <-8> in function call
[scip_solve.c:2607] ERROR: Error <-8> in function call
Traceback (most recent call last):
File "C:\Users\15572\PycharmProjects\pythonProject2(scip)\sa.py", line 35, in
m.optimize()
File "src\pyscipopt\scip.pxi", line 3529, in pyscipopt.scip.Model.optimize
File "src\pyscipopt\scip.pxi", line 280, in pyscipopt.scip.PY_SCIP_CALL
Exception: SCIP: method cannot be called at this time in solution process!

[Joao-Dionisio]

Ah, I forgot that you need self as well. So for catching the events you need to do self.model.catchEvent(SCIP_EVENTTYPE.BESTSOLFOUND, self) and self.model.catchEvent(SCIP_EVENTTYPE.LPSOLVED, self)

I'll try to implement this myself and I'll get back to you when I get it working

[sheepone12138]

When I follow your method there is no more error, but what I don't quite understand is: when I use print(m.data)，The primal_solutions and dual_solutions obtained are empty sets, like this{'primal_solutions': [], 'dual_solutions': []}