IpOpt拒绝解决不受约束的问题
IpOpt refuses to solve unconstrained problem
我对IpOpt很陌生,我正试图用它来解决简单的无约束优化问题。我的问题只是二次函数f(x) = (5x - 3)^2
。
我为这个问题创建了一个简单的类:
#include <cstdio>
#include "IpIpoptApplication.hpp"
#include "IpTNLP.hpp"
using namespace Ipopt;
class MyProblem : public Ipopt::TNLP
{
public:
const int nVars = 1;
virtual bool get_nlp_info(Index& n, Index& m, Index& nnz_jac_g,
Index& nnz_h_lag, IndexStyleEnum& index_style)
{
n = nVars;
m = 0;
nnz_jac_g = 0;
nnz_h_lag = 1;
index_style = IndexStyleEnum::C_STYLE;
return true;
}
virtual bool get_bounds_info(Index n, Number* x_l, Number* x_u,
Index m, Number* g_l, Number* g_u)
{
return true;
}
virtual bool get_starting_point(Index n, bool init_x, Number* x,
bool init_z, Number* z_L, Number* z_U,
Index m, bool init_lambda,
Number* lambda)
{
std::cout << "get_starting_point" << std::endl;
if(init_x){
x[0] = 0.0;
}
return true;
}
virtual bool eval_f(Index n, const Number* x, bool new_x,
Number& obj_value)
{
const Number residual = (5 * x[0] - 3);
obj_value = residual * residual;
std::cout << "obj_value " << obj_value << std::endl;
return true;
}
virtual bool eval_grad_f(Index n, const Number* x, bool new_x,
Number* grad_f)
{
const Number residual = (5 * x[0] - 3);
grad_f[0] = 2 * 10 * residual;
std::cout << "grad_f " << grad_f[0] << std::endl;
return true;
}
virtual bool eval_g(Index n, const Number* x, bool new_x,
Index m, Number* g)
{
std::cout << "eval_g was called m=" << m << " *g " << g << std::endl;
return true;
}
virtual bool eval_jac_g(Index n, const Number* x, bool new_x,
Index m, Index nele_jac, Index* iRow,
Index *jCol, Number* values)
{
std::cout << "eval_jac_g was called" << std::endl;
return true;
}
virtual void finalize_solution(SolverReturn status,
Index n, const Number* x, const Number* z_L, const Number* z_U,
Index m, const Number* g, const Number* lambda,
Number obj_value,
const IpoptData* ip_data,
IpoptCalculatedQuantities* ip_cq)
{
std::cout << "X final " << x[0] << std::endl;
}
};
int main()
{
SmartPtr<TNLP> mynlp = new MyProblem();
SmartPtr<IpoptApplication> app = new IpoptApplication();
app->Initialize();
ApplicationReturnStatus status = app->OptimizeTNLP(mynlp);
if (status == Solve_Succeeded) {
printf("nn*** The problem solved!n");
}
else {
printf("nn*** The problem FAILED!n");
}
return 0;
}
我设置了m = 0;
和nnz_jac_g = 0;
,以表明我没有限制
IpOpt给我以下输出:
******************************************************************************
This program contains Ipopt, a library for large-scale nonlinear optimization.
Ipopt is released as open source code under the Eclipse Public License (EPL).
For more information visit http://projects.coin-or.org/Ipopt
******************************************************************************
This is Ipopt version 3.12.11, running with linear solver mumps.
NOTE: Other linear solvers might be more efficient (see Ipopt documentation).
eval_g was called m=0 *g 0x5559e0661fd0
obj_value 9
X final 0
All variables are fixed and constraint violation 0.000000e+00
is below tolerance 1.000000e-08. Declaring success.
EXIT: Optimal Solution Found.
*** The problem solved!
但显然,问题的解决方案应该是类似x = 3/5
的东西。
我的问题课出了什么问题?
据我所知,IpOpt没有调用get_starting_point
方法,这看起来很奇怪。
我的IpOpt版本是3.12.11
,gcc是7.2.0-8ubuntu3.2
。我可以提供更多的信息,但我不知道我还能展示什么。源代码应该足够漂亮
PS抱歉我的英语不好
您应该定义var边界。
virtual bool get_bounds_info(Index n, Number* x_l, Number* x_u,
Index m, Number* g_l, Number* g_u)
{
*x_l = -1e19; //nlp_lower_bound_inf
*x_u = 1e19; //nlp_upper_bound_inf
return true;
}
变量被零边界随机冻结(因为其他编译器可能为未定义的边界提供不同的值(。
存在另一个错误。您应该定义eval_h或指定hessian近似方案。
相关文章:
- C++ 雷神库 - 使用资源加载器类时出现问题(不命名类型)
- C++问题:用户认为数字1-100,程序提出问题不超过6次即可得到答案。无法正确
- IpOpt拒绝解决不受约束的问题
- 受约束的成员函数和显式模板实例化
- 使用 NTAllocateVirtualMemory 和 GetProcAddress 的内存分配问题不起作用
- 循环仅对第一行正常工作.其他行不受 for 的影响
- 如何拥有受约束的运算符模板?
- 如何保护非托管应用程序中的字符串不受进程转储的影响
- 检查nullptr是否100%保护内存布局不受segfault影响
- 未使用的 asm() 在不受支持的体系结构上的行为
- OpenGL中使用Freetype进行字体渲染的问题不起作用
- 如何在不受其他文件影响的情况下"by itself" Visual Studio 项目中运行C++文件?
- IPOPT不遵守约束,但在使用CppAD时不记录违规行为
- 流的奇怪问题.不能存储具有特定值的成员变量的对象
- 可视化如何确定项目是由C 项目管理还是不受管理的C 项目
- 如何在QTextDocument中使用表格的框架属性,即使它不受支持?
- 为什么重载解析更喜欢不受约束的模板函数而不是更具体的模板函数?
- 将字节数组从c 不受管理的dll传递到c#unity
- 如何解决多线程绘图的问题不流畅
- 函数不受主内存约束的函数所需的复杂性是什么?