My code creates protein models that need to be relaxed in order to lower their energy as these models may have some close contacts and other non-protein like features. According to the online documentation, the Rosetta command for relaxing non-idealized structures is
rosetta.exe aa xxx_ -relax -farlx -minimize -s xxx.start -fa_input -fa_output -use_input_bond
My questions are: 1) Is this the correct command for relaxing my models and make them "fall" in the nearest local minimum? 2) If so, how can I call these functions from inside my C++ code?
Thanks in advance,