在 Dymola 中编译大型数组

Question

请考虑以下小型Modelica模型和函数：

model VectorizeDemo
  parameter Integer na=5;
  final parameter Integer nb=2*na;
  final parameter Real a[na] = {2*i for i in 1:na};
  final parameter Real b[nb] = {3*i for i in 1:nb};
  Real c[na];
  Real d[na,nb];

protected 
  function myFun
    input Real A;
    input Real B;
    output Real C;
  algorithm 
    C:=tanh(A)*sin(B);
  end myFun;

equation 
  c = sin(a);
  //d = myFun(a,b);
  // inner loop first
  d = {myFun(a[i], b[j]) for j in 1:nb, i in 1:na};
end VectorizeDemo;

这将在 Dymola 中编译和模拟，但是查看 dsmodel.c 中的 C 代码，每个数组元素都被声明为一个新变量：

...
DeclareVariable("d[4, 10]", "", 38.0, 0.0,0.0,0.0,0,513)
DeclareVariable("d[5, 1]", "", 13.0, 0.0,0.0,0.0,0,513)
DeclareVariable("d[5, 2]", "", 16.0, 0.0,0.0,0.0,0,513)
DeclareVariable("d[5, 3]", "", 19.0, 0.0,0.0,0.0,0,513)
...

因此，如果我通过设置 na=1000 来增加数组大小，我将声明 1000*2000 个变量。显示的示例仍然会编译，即使它需要很长时间，但我更复杂的用例因编译器 warning C4049: compiler limit, terminating line number emission 或 C1002 compiler is out of heap space.

而失败

旁注：较大的示例也需要几分钟的时间来检查，模拟后，在变量浏览器中展开变量时，GUI 将被阻塞很长时间。

是否有任何解决方法，例如重写我的代码或设置一些标志？临时增加堆space？我只需要运行模型一次。对正在发生的事情的任何见解也将不胜感激。使用 Dymola 2020 和 VisualStudio 2017。

Answer 1

默认情况下，所有 Modelica 编译器都会将所有数组（有一些例外，如函数中的记录中的数组）扩展为标量。 OpenModelica 在前端和后端开始了一些关于非扩展数组的工作。看： http://www.ep.liu.se/ecp/157/071/ecp19157071.pdf

Answer 2

IIRC，Modelica 规范的某些部分是以需要对所有变量进行标量化的方式编写的。在许多情况下这不是必需的，但对于一些常见的索引缩减情况来说更简单。对于某些离散化 PDE 的情况，它似乎特别没有必要。

Answer 3

是的，至少在Dymola 2020（可能还有更早的版本）中可以如下初步避免编译问题：

Real d[na,nb] annotation(__Dymola_HideArray=true);

但是，也可能存在其他可能性 - 该示例并未完全阐明它的用途。特别是我注意到 'd' 可以被评估并且根本没有被使用。

在 Dymola 中编译大型数组

Compile large array in Dymola

arrays

compiler-errors

modelica

dymola