Hello all,
I am not sure if this is a CVXPy, modelica, python, or maybe even a C problem (I believe CVXPy and/or the SCS solver calls functions written in C). I am simulating a microgrid using openmodelica that has a python module that calls a CVXPy function as a control algorithm. I tested the function using Python 3.8(modelica buildings does not use newer versions of python) and it runs without error, however when I run it in openmodelica I get a segmentation fault, and then it crashes. Is there a way to adjust memory in openmodelica, or what else may be causing the error? Note: For the optimization to work it needs to gather data from the modelica for the first 96 iterations, it is not until the 97th iteration when the modelica optimizer is called, and it runs fine. Both the modelica and CVXPy work independently, I only see the problem when Modelica tries to run CVXPy.
Python function test:
Input:
for i in range(100):
tou_cvx_one_day([abs(df.building_load[i]), abs(df.solar[i]), 0.5, i * 900])
Output:
===============================================================================
CVXPY
v1.3.1
===============================================================================
(CVXPY) Aug 01 02:03:02 PM: Your problem has 480 variables, 960 constraints, and 0 parameters.
(CVXPY) Aug 01 02:03:02 PM: It is compliant with the following grammars: DCP, DQCP
(CVXPY) Aug 01 02:03:02 PM: (If you need to solve this problem multiple times, but with different data, consider using parameters.)
(CVXPY) Aug 01 02:03:02 PM: CVXPY will first compile your problem; then, it will invoke a numerical solver to obtain a solution.
-------------------------------------------------------------------------------
Compilation
-------------------------------------------------------------------------------
(CVXPY) Aug 01 02:03:02 PM: Compiling problem (target solver=SCS).
(CVXPY) Aug 01 02:03:02 PM: Reduction chain: Dcp2Cone -> CvxAttr2Constr -> ConeMatrixStuffing -> SCS
(CVXPY) Aug 01 02:03:02 PM: Applying reduction Dcp2Cone
(CVXPY) Aug 01 02:03:02 PM: Applying reduction CvxAttr2Constr
(CVXPY) Aug 01 02:03:02 PM: Applying reduction ConeMatrixStuffing
(CVXPY) Aug 01 02:03:03 PM: Applying reduction SCS
(CVXPY) Aug 01 02:03:03 PM: Finished problem compilation (took 5.538e-01 seconds).
-------------------------------------------------------------------------------
Numerical solver
-------------------------------------------------------------------------------
(CVXPY) Aug 01 02:03:03 PM: Invoking solver SCS to obtain a solution.
------------------------------------------------------------------
SCS v3.2.3 - Splitting Conic Solver
(c) Brendan O Donoghue, Stanford University, 2012
------------------------------------------------------------------
problem: variables n: 483, constraints m: 1248
cones: z: primal zero / dual free vars: 288
l: linear vars: 960
settings: eps_abs: 1.0e-05, eps_rel: 1.0e-05, eps_infeas: 1.0e-07
alpha: 1.50, scale: 1.00e-01, adaptive_scale: 1
max_iters: 100000, normalize: 1, rho_x: 1.00e-06
acceleration_lookback: 10, acceleration_interval: 10
lin-sys: sparse-direct-amd-qdldl
nnz(A): 1917, nnz(P): 0
------------------------------------------------------------------
iter | pri res | dua res | gap | obj | scale | time (s)
------------------------------------------------------------------
0| 4.50e+02 1.96e+01 2.85e+03 1.43e+03 1.00e-01 1.59e-02
250| 4.05e-01 6.32e-03 1.90e+00 -1.81e+01 9.53e-03 2.55e-02
500| 2.89e-01 3.57e-03 1.79e+00 -1.95e+01 9.53e-03 3.37e-02
750| 1.07e-01 2.15e-03 2.78e+00 -2.09e+01 9.53e-03 4.31e-02
1000| 1.14e-01 1.38e-03 4.14e-01 -2.29e+01 9.53e-03 5.34e-02
1250| 8.00e-02 1.09e-03 1.10e+00 -2.40e+01 9.53e-03 6.46e-02
1500| 5.70e-02 1.43e-03 1.39e+00 -2.44e+01 9.53e-03 7.35e-02
1750| 3.43e+02 3.29e+00 1.20e+02 4.40e+01 9.53e-03 8.15e-02
2000| 4.36e-02 7.66e-04 8.12e-01 -2.36e+01 9.53e-03 8.88e-02
2250| 2.97e-02 5.15e-04 9.12e-01 -2.37e+01 9.53e-03 9.60e-02
2500| 2.24e-02 5.13e-04 9.24e-01 -2.38e+01 9.53e-03 1.03e-01
2750| 3.10e-02 4.82e-04 2.37e+00 -2.56e+01 9.53e-03 1.11e-01
3000| 3.38e-02 2.01e-04 1.73e+00 -2.53e+01 9.53e-03 1.18e-01
3250| 3.62e-02 3.74e-04 1.10e+00 -2.50e+01 9.53e-03 1.26e-01
3500| 1.48e-02 2.73e-04 1.64e-01 -2.45e+01 9.53e-03 1.33e-01
3750| 2.95e-02 2.56e-04 8.51e-01 -2.49e+01 9.53e-03 1.40e-01
4000| 2.20e-02 2.00e-04 4.30e-02 -2.45e+01 9.53e-03 1.47e-01
4250| 2.46e-02 2.71e-04 2.70e-01 -2.46e+01 9.53e-03 1.55e-01
4500| 2.12e-03 3.13e-05 2.26e-01 -2.44e+01 9.53e-03 1.62e-01
4750| 1.29e-03 3.72e-05 2.15e-01 -2.46e+01 9.53e-03 1.69e-01
5000| 6.15e-03 6.53e-05 2.88e-01 -2.46e+01 9.53e-03 1.76e-01
5250| 2.33e-03 4.50e-05 2.99e-01 -2.43e+01 9.53e-03 1.83e-01
5500| 4.44e-03 4.87e-05 2.29e-01 -2.46e+01 9.53e-03 1.90e-01
5750| 3.85e-03 3.49e-05 1.83e-01 -2.44e+01 9.53e-03 1.97e-01
6000| 8.10e-03 1.17e-04 5.12e-01 -2.47e+01 9.53e-03 2.05e-01
6250| 1.59e-02 5.62e-05 1.94e-01 -2.46e+01 9.53e-03 2.13e-01
6500| 1.24e-02 3.67e-05 9.37e-02 -2.44e+01 9.53e-03 2.25e-01
6750| 2.92e-03 2.45e-05 6.63e-02 -2.45e+01 9.53e-03 2.35e-01
7000| 2.63e-03 1.49e-05 9.56e-02 -2.44e+01 9.53e-03 2.42e-01
7250| 2.33e-03 1.58e-05 9.63e-02 -2.45e+01 9.53e-03 2.50e-01
7500| 4.15e+00 3.96e-02 2.37e+00 -2.33e+01 9.53e-03 2.58e-01
7750| 9.16e-04 1.74e-05 1.05e-01 -2.45e+01 9.53e-03 2.65e-01
8000| 1.63e-03 3.81e-05 1.82e-01 -2.44e+01 9.53e-03 2.72e-01
8250| 1.02e-03 6.16e-05 1.25e-02 -2.45e+01 9.53e-03 2.79e-01
8500| 3.78e-04 5.75e-06 3.63e-02 -2.45e+01 9.53e-03 2.86e-01
8750| 3.77e-04 3.33e-06 2.26e-04 -2.45e+01 9.53e-03 2.94e-01
------------------------------------------------------------------
status: solved
timings: total: 2.94e-01s = setup: 8.91e-03s + solve: 2.85e-01s
lin-sys: 2.00e-01s, cones: 2.23e-02s, accel: 1.19e-02s
------------------------------------------------------------------
objective = -24.469205
------------------------------------------------------------------
-------------------------------------------------------------------------------
Summary
-------------------------------------------------------------------------------
(CVXPY) Aug 01 02:03:03 PM: Problem status: optimal
(CVXPY) Aug 01 02:03:03 PM: Optimal value: -2.447e+01
(CVXPY) Aug 01 02:03:03 PM: Compilation took 5.538e-01 seconds
(CVXPY) Aug 01 02:03:03 PM: Solver (including time spent in interface) took 3.015e-01 seconds
OpenModelica Buildings Output:
Process crashed Simulation process failed. Exited with code 11. /tmp/OpenModelica_sigi-laptop/OMEdit/conf_paper_microgrid/conf_paper_microgrid -port=45477 -logFormat=xmltcp -override=startTime=0,stopTime=3.15351e+07,stepSize=899.974,tolerance=1e-06,solver=dassl,outputFormat=mat,variableFilter=.* -r=/tmp/OpenModelica_sigi-laptop/OMEdit/conf_paper_microgrid/conf_paper_microgrid_res.mat -w -lv=LOG_STATS -inputPath=/tmp/OpenModelica_sigi-laptop/OMEdit/conf_paper_microgrid -outputPath=/tmp/OpenModelica_sigi-laptop/OMEdit/conf_paper_microgrid The initialization finished successfully without homotopy method. Process crashed =============================================================================== CVXPY v1.3.2 =============================================================================== (CVXPY) Aug 04 06:37:08 PM: Your problem has 480 variables, 960 constraints, and 0 parameters. (CVXPY) Aug 04 06:37:08 PM: It is compliant with the following grammars: DCP, DQCP (CVXPY) Aug 04 06:37:08 PM: (If you need to solve this problem multiple times, but with different data, consider using parameters.) (CVXPY) Aug 04 06:37:08 PM: CVXPY will first compile your problem; then, it will invoke a numerical solver to obtain a solution. ------------------------------------------------------------------------------- Compilation ------------------------------------------------------------------------------- (CVXPY) Aug 04 06:37:08 PM: Compiling problem (target solver=SCS). (CVXPY) Aug 04 06:37:08 PM: Reduction chain: Dcp2Cone -> CvxAttr2Constr -> ConeMatrixStuffing -> SCS (CVXPY) Aug 04 06:37:08 PM: Applying reduction Dcp2Cone (CVXPY) Aug 04 06:37:08 PM: Applying reduction CvxAttr2Constr (CVXPY) Aug 04 06:37:08 PM: Applying reduction ConeMatrixStuffing (CVXPY) Aug 04 06:37:08 PM: Applying reduction SCS (CVXPY) Aug 04 06:37:08 PM: Finished problem compilation (took 5.609e-01 seconds). ------------------------------------------------------------------------------- Numerical solver ------------------------------------------------------------------------------- (CVXPY) Aug 04 06:37:08 PM: Invoking solver SCS to obtain a solution. Limited backtrace at point of segmentation fault /lib/x86_64-linux-gnu/libpthread.so.0(+0x14420)[0x7fa45551d420] /home/sigi-laptop/.local/lib/python3.8/site-packages/_scs_direct.cpython-38-x86_64-linux-gnu.so(QDLDL_factor+0x82)[0x7fa4203f1722] /home/sigi-laptop/.local/lib/python3.8/site-packages/_scs_direct.cpython-38-x86_64-linux-gnu.so(+0x4b0e)[0x7fa4203ecb0e] /home/sigi-laptop/.local/lib/python3.8/site-packages/_scs_direct.cpython-38-x86_64-linux-gnu.so(scs_init_lin_sys_work+0x580)[0x7fa4203ed280] /home/sigi-laptop/.local/lib/python3.8/site-packages/_scs_direct.cpython-38-x86_64-linux-gnu.so(scs_init+0x4c3)[0x7fa4203fde33] /home/sigi-laptop/.local/lib/python3.8/site-packages/_scs_direct.cpython-38-x86_64-linux-gnu.so(+0x17dc1)[0x7fa4203ffdc1] /lib/x86_64-linux-gnu/libpython3.8.so.1.0(+0x2491b9)[0x7fa4568491b9] /lib/x86_64-linux-gnu/libpython3.8.so.1.0(PyObject_Call+0x74)[0x7fa4568a8994] /lib/x86_64-linux-gnu/libpython3.8.so.1.0(_PyEval_EvalFrameDefault+0x590a)[0x7fa45667aa7a] /lib/x86_64-linux-gnu/libpython3.8.so.1.0(_PyEval_EvalCodeWithName+0x8fb)[0x7fa4567cae4b] /lib/x86_64-linux-gnu/libpython3.8.so.1.0(_PyFunction_Vectorcall+0x94)[0x7fa4568a8124] /lib/x86_64-linux-gnu/libpython3.8.so.1.0(_PyObject_FastCallDict+0x79)[0x7fa4568aa3a9] /lib/x86_64-linux-gnu/libpython3.8.so.1.0(_PyObject_Call_Prepend+0xcd)[0x7fa4568aa52d] /lib/x86_64-linux-gnu/libpython3.8.so.1.0(+0x243d47)[0x7fa456843d47] /lib/x86_64-linux-gnu/libpython3.8.so.1.0(+0x2491b9)[0x7fa4568491b9]
Thank you for the help
