Ipopt添加避障约束触发malloc错误：2D点集优化问题求助

你的问题核心是添加避障约束后触发内存越界写入，导致malloc(): invalid size (unsorted)错误，结合Valgrind的检测结果，从以下几个方向逐一排查：

1. 约束数量与数组大小不匹配 #

避障约束针对的是中间N-2个点（首尾点固定无需避障），每个点对应1个约束：(x_i-4)^2 + (y_i-4)^2 >= 1。

• 检查get_n_constraints函数的返回值是否为N-2，如果返回值错误（比如写成N或N-1），会导致Ipopt分配的约束上下限数组大小和实际需要的不匹配，触发越界写入。
• 确认存储约束上下限的数组长度严格等于N-2：比如在get_constraints_lower_bounds和get_constraints_upper_bounds中，手动分配的数组（如new double[n_c]）的n_c必须是N-2，循环赋值时索引不能超出数组范围（比如循环从0到N-3，而非N-2）。

2. 变量索引计算错误 #

轨迹的总变量数是2*N（每个点的x、y坐标各占1个变量），计算约束值或雅可比矩阵时，极易出现索引偏移错误：

• 第i个中间点（i从1到N-2）对应的x变量索引是2*i，y变量索引是2*i+1，不要误写成i或2*i+2，否则会访问超出变量数组的内存区域。
• 对照Valgrind报告的越界位置，比如某行代码访问了x[xxx]，计算xxx的数值是否超过2*N-1（变量数组的最大合法索引）。

3. Ipopt接口函数的内存管理问题 #

实现Ipopt的核心接口时，需严格遵循内存规则：

• 在eval_g（计算约束值）和eval_jac_g（计算雅可比矩阵）中，确保传给Ipopt的数组（如g、values）的大小完全匹配约束数量和雅可比非零元素数量。比如雅可比矩阵的非零元素数，每个避障约束对应2个非零元素（对x_i和y_i的偏导），所以总非零数是2*(N-2)，如果这个数值算错，会导致Ipopt写入越界。
• 避免使用栈上的固定大小数组（比如double g[100]），如果N是动态输入的，栈数组大小可能不足，改用动态分配的堆数组（new double[n_c]）并确保大小正确。

4. 避障约束的上下限设置验证 #

虽然你调整上下限无效，但仍需确认约束形式正确：

• 避障约束是下限约束：g_i = (x_i-4)^2 + (y_i-4)^2，所以get_constraints_lower_bounds中对应位置设为1.0，get_constraints_upper_bounds设为Ipopt提供的Ipopt::infinity()（或一个极大值如1e20）。不要误将下限设为0或其他值，虽然这不会直接导致内存越界，但如果约束数量错误，上下限数组的赋值会覆盖无关内存。

快速定位技巧 #

直接查看Valgrind的详细报告，找到具体的越界写入位置（比如函数名、行号）：

• 如果报告显示在get_constraints_lower_bounds中写入了超出数组的位置，说明约束数量或数组大小计算错误；
• 如果显示在eval_g中访问了x数组的非法索引，说明变量索引计算有误。

内容的提问来源于stack exchange，提问作者Vladimir Petrov