OpenModelica中四维数组插值性能优化及技术问询

背景 #

建模过程中需获取依赖三个变量的值，因无法构建精度足够的多项式，采用三维数组线性插值方案（此前已在Wolfram Mathematica实现）。在OpenModelica中使用ExternData库加载CSV生成37×164×51的三维数组后，插值计算耗时长达20分钟（相同算法在Mathematica仅需数秒），针对以下问题解答：

1. 低性能是否由实现错误导致？若为是，该如何修复？ #

是，存在多处低效写法与潜在错误，修复点如下：

• 修正参数数组误用：将intarr函数中parameter Real A[:,:,:];改为Real A[:,:,:];。parameter是编译期常量，无法在algorithm段赋值，此写法会触发不必要的常量处理逻辑，大幅增加开销。
• 替换线性遍历为二分查找：当前用for循环线性遍历找索引（如找ii、jj、kk），对于y维度164个元素，每次遍历需162次判断。因你的x、y、z参考值均为有序数组，改用二分查找可将查找复杂度从O(n)降至O(logn)，大幅减少计算量。
• 修复循环边界错误：找jj的循环1:leny-2改为1:leny-1，找kk的循环1:lenz-2改为1:lenz-1，否则最后一个区间的变量无法匹配到索引。
• 移除冗余内存操作：loaddata函数中string[nz] := 0、list[ny,nz] := 0属于冗余赋值，每次循环都会覆盖这些值，直接删除即可减少内存操作开销。

2. 若实现无误，属于OpenModelica环境问题，该如何优化数组插值？ #

若代码修复后仍有性能问题，可从以下方向优化：

• 标记纯函数：在intarr函数前添加pure关键字（pure function intarr），告知编译器该函数无副作用，可进行编译期优化（如内联、常量折叠）。
• 分离参考坐标数组：提前将x、y、z的网格值从三维数组中分离为独立一维数组（如x_grid[lenx] = A[:, leny, lenz]），避免每次插值都从三维数组中提取参考值，减少数组访问开销。
• 启用编译优化：编译模型时添加-O3参数开启最高级别优化，或在函数上添加annotation(Inline=true)实现函数内联，减少函数调用的开销。
• 预生成静态常量数组：用loaddata生成数组后，将其导出为Modelica格式的常量数组文件，直接在模型中导入使用，避免运行时加载CSV与转换数组的开销。

3. 是否可不依赖ExternData库，仅用内置函数实现类似功能？ #

可以，使用OpenModelica内置的文件处理函数即可实现CSV加载，示例代码如下：

pure function loadCSVWithoutExternData
  input String fileName;
  input Integer nHeaderLines = 1;
  output Real data[:, :];
protected
  String fileContent = Modelica.Utilities.Files.readFile(fileName);
  String lines[:] = Modelica.Utilities.Strings.split(fileContent, "\n");
  // 过滤空行与表头行
  String dataLines[:] = {line for line in lines if not Modelica.Utilities.Strings.isBlank(line)}[nHeaderLines+1:end];
  Integer nRows = size(dataLines, 1);
  String firstLineCols[:] = Modelica.Utilities.Strings.split(dataLines[1], ",");
  Integer nCols = size(firstLineCols, 1);
  Real result[nRows, nCols];
algorithm
  for i in 1:nRows loop
    String cols[:] = Modelica.Utilities.Strings.split(dataLines[i], ",");
    for j in 1:nCols loop
      result[i,j] = String.toReal(cols[j]);
    end for;
  end for;
  data := result;
end loadCSVWithoutExternData;

注：需自行处理CSV格式不一致的情况（如分隔符、空值）。

4. 拥有Wolfram Mathematica授权，能否调用其API完成插值并获取结果？ #

可以，两种可行方案：

• 命令行调用脚本：编写Mathematica脚本，接受x、y、z参数与数据路径，计算插值后输出结果。在OpenModelica中用Modelica.Utilities.System.command调用该脚本，读取控制台输出转换为Real值。
• 预计算插值结果：若插值网格固定，用Mathematica提前计算所有可能的插值结果，保存为CSV或Modelica常量数组，直接在OpenModelica中使用，完全避免运行时跨环境调用的开销。
• C API调用：利用Mathematica LinkKit生成C语言接口，在OpenModelica中通过external关键字直接调用Mathematica的插值函数，适合需要动态计算的场景。

内容的提问来源于stack exchange，提问作者Volodymyr S