I would like to use the Python auto generated code for a MHE solver directly from C++. As it was stated in MHE using algebraic variables, I was able to solve my problem using the Python interface thanks to the help provided here. However, I couldn’t use the auto generated solver from C++. So, I started from scratch with minimal_example_mhe.py to see if I could achieve that task.
After digging around with the code, I was able to create a CMake project where the mhe_pendulum_ode_test.cpp file uses the auto generated code to solve the same problem that was solved using Python. One detail that is bugging me is that the errors between the true and the estimated states are a little different from the Python version.
The code is in this repository mhe_acados_cpp. You can try it by doing the following
Edit CMakeLists.txt to fit your acados installation.
Go to solver/mhe folder and run python3 minimal_example_mhe.py to get the auto generated code.
Go back to the project root, create a build directory and run cmake and make to get an executable
After runing the executable you should have a text file called simXest.txt inside the build directory that contains the values of the MHE estimated states obtained using C++.
Go into the src directory and run python3 plot_cpp_results.py to compare.
Could anyone point me if I’m missing something in the C++ code? It took me forever to realize how to set the parameters p (hint: it was the acados_update_params function).
and providing references for (x, u, x) then I get the result,
For native python:
Note that the result from cpp with
ocp_mhe.model.cost_y_expr = vertcat(x, u, x)
is exactly same as the one obtained with
ocp_mhe.model.cost_y_expr = vertcat(x, u)
I have been using NONLINEAR_LS cost in my formulation and face the same issue. In native python with NONLINEAR_LS cost my formulation works as expected but if I use generated code from native python in cpp then the solution is different in cpp from the native python.