- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我是 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("\n\n*** The problem solved!\n");
}
else {
printf("\n\n*** 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 抱歉我的英语不好
最佳答案
您应该定义变量范围。
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++ - IpOpt 拒绝解决无约束问题,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/53065342/
我只想允许一个国家/地区访问,但排除该国家/地区内的代理。 这就是我所拥有的(为了方便起见,缩短了版本) order deny,allow deny from all allow from 139.
这个问题在这里已经有了答案: What is an unhandled promise rejection? (9 个回答) 关闭 4 年前。 我目前正在尝试实现我自己的 Promise,以便在 A
我在使用 Gitolite 推送 git 时遇到问题。 当我尝试这个时: git push origin :refs/tags/deployment 我收到这个错误: remote: D NAME/i
我已经为我的 laravel 5.0-dev 项目配置了 mysql,如下所示: 'mysql' => [ 'driver' => 'mysql', 'host' =>
我对 Web 和 SOF 进行了一些研究,但发现对于该错误没有任何真正的帮助。 我使用 Windows 10 Ubuntu Bash 安装了 Node 和 Puppeteer,但未能使其工作,但我设法
在我的应用审核期间,我收到了以下信息: “17.2:要求用户共享个人信息(例如电子邮件地址和生日)才能正常运行的应用将被拒绝 具体来说,您的应用仅使用Facebook登录名进行身份验证,但不包括该网站
我正在开发 VeriFone VX 终端的接口(interface)。虽然,这确实是一个普遍的 EMV 问题。我们的处理器的下限为零,因此它将始终在线发送。但是,如果它发生变化,您如何知道(哪些标签)
我编写了一些宏代码,根据表单提交向经理发送电子邮件(用于费用/审批流程),这是我使用谷歌表单/电子表格的第一个项目,所以也许我可能会错过一些简单的东西,但我为此浏览了 2 个教程,我的代码与重要的部分
clang 3.4 接受以下代码;而 vc++ NOV 2013 CTP 拒绝它并出现错误: error C2668: 'AreEqual' : ambiguous call to overloade
使用 nginx,您可以允许和拒绝范围和 ips (https://www.nginx.com/resources/admin-guide/restricting-access/)。使用realip模
官方编辑: 非常感谢您的帮助,但我仍然遇到问题。 我的 ffserver.conf 文件是这样的: # Port on which the server is listening. You must
我有一个问题:我是 Ubuntu 系统的根。我想授予用户(比如用户名是 X)执行任何命令的权限,但同时我有一个文件夹,除了我的用户(当然不是 X,因为它是 Admin ) 或根。有什么建议么?谢谢!
我使用 Apache2.2 作为 tomcat 服务器的前端。我想限制对某个位置的访问,但允许对子位置的所有访问,但遇到了一些麻烦。 我目前拥有的是: AllowOverride None
就像 this person ,我一直在为浏览器缓存 SSL session 而苦苦挣扎。简而言之,如果选择了客户端证书,则无法以编程方式清除状态,除非在 IE 中使用 document.execCo
我的网站是在由 Apache 服务器提供服务的 Angular 上设置的。我通过 View 将内容动态加载到主页上。 现在以下是我的问题: 我建立这个网站的主要目的是通过 google adsense
我最近遇到了我的应用程序的问题,当它突然被 Google Play 拒绝时因为他们发现我使用的是背景位置 .但实际上我并没有使用这个功能。我只有 ACCESS_COARSE_LOCATION和 ACC
function sendPushNotification(subscription, urlEncodedData){ try { webpush.sendNotification(su
我包裹了一个 request-promise-native调用返回 promise 的函数。 import request from 'request-promise-native'; functio
我正在开发我的 meteor 项目,并开始设置我的第一个更复杂的允许/拒绝规则。我发现很难看出哪些允许触发,哪些不允许触发,以及这些函数中的某些变量包含什么。例如: List.allow({ u
我正在 AngularJS 中创建一个 Factory,它是这样的: if (href) { return $http({ method: method, url: item.href });
我是一名优秀的程序员,十分优秀!