在Gekko上使用scipy.integrate.trapz时出错
我正在尝试使用Gekko库解决PDE，该PDE涉及通过在一定时间与空间集成速度来计算体积流量。
Qgap==scipy.integrate.trapz(2*np.pi*xpos*np.array(u1))
在每个XPOS 1000点处离散U1。 我得到这个错误：
APM model error: string > 15000 characters
Consider breaking up the line into multiple equations
The may also be due to only using newline character CR
instead of CR LF (for Windows) or LF (for MacOS/Linux)
To fix this problem, save APM file with appropriate newline characters
STOPPING...
解决这个问题我该怎么办？ M.Integral（U1，XPOS）选项是否起作用或M.Integral（）仅与时间不可或缺？
I am trying to solve a PDE with gekko library which involves calculating volume flow rate by integrating velocity with the space at a certain time.
Qgap==scipy.integrate.trapz(2*np.pi*xpos*np.array(u1))
where u1 is discretized at each xpos 1000 points.
I get this error:
APM model error: string > 15000 characters
Consider breaking up the line into multiple equations
The may also be due to only using newline character CR
instead of CR LF (for Windows) or LF (for MacOS/Linux)
To fix this problem, save APM file with appropriate newline characters
STOPPING...
What should I do to solve this problem?
Does m.integral(u1,xpos) option works or m.integral() only works with integral with time only?
如果你对这篇内容有疑问，欢迎到本站社区发帖提问 参与讨论，获取更多帮助，或者扫码二维码加入 Web 技术交流群。
绑定邮箱获取回复消息
由于您还没有绑定你的真实邮箱，如果其他用户或者作者回复了您的评论，将不能在第一时间通知您！
发布评论
评论（1）
scipy.integrate.trapz函数可与gekko一起起作用，因为它会产生符号结果，如下所示：问题有关解决1D Navier用可压缩的质量保护（液压阻尼器）的问题但是，这种符号形式与整合维度的复杂性增加。在关于早期问题的测试中，它仅限于150至200个空间离散步骤之间。在200或以上，我建议使用m.sum（）分解方程式编程您自己的梯形方程式The
scipy.integrate.trapzfunction works withgekkobecause it produces a symbolic result as shown here: Problem on Solving a 1D Navier Stokes with Compressible Mass Conservation (Hydraulic Damper) However, this symbolic form increases in complexity with the integrating dimension. In tests on the earlier problem, it is limited to somewhere between 150 and 200 spatial discretization steps. At 200 or beyond, I recommend programming your own trapezoidal equation that usesm.sum()to break up the equation