OpenFOAM 如何从列表中删除一些元素？

Question

在 OpenFOAM 中，我可以访问我的模拟时间列表，如下所示：

const auto& tlist = mesh.time().times(); //or runTime.times();

只要在自定义函数对象的上下文中考虑这一点，您就可以在其中访问时间列表。

当我打印该列表时：

Foam::Info << tlist << Foam::endl;

然后通过postProcess命令运行函数对象，我得到：

9
(
0   constant
0   0
0.001   0.001
0.002   0.002
0.003   0.003
0.004   0.004
0.005   0.005
0.006   0.006
0.007   0.007
)

End

我想删除该列表中的前两个元素，即(0 constant) 和(0 0)。但是我找不到任何方法可以做到这一点，我发现的只是如果对象是一个 HashTable，那么有一个 erase 方法可以通过它的键删除元素。

任何想法如何从我的时间列表中删除前两个元素tlist？或者至少我怎样才能将该列表转换为另一个允许我这样做的数据结构？

谢谢

编辑：

这里是OpenFOAM上下文中List的定义链接：https://cpp.openfoam.org/v9/classFoam_1_1List.html

Answer 1

免责声明，我从未使用过 OpenFOAM。 看起来 List 有迭代器。 https://cpp.openfoam.org/v9/classFoam_1_1UList.html。 所以你可以尝试这样的事情（假设迭代器像我在其他库中习惯的那样工作）：

  const auto& tlist = mesh.time().times(); //or runTime.times();

  // assuming operator+ available on iterator. this will copy data
  auto sublist = List<instant>(tlist.begin()+2,tlist.end());

  Foam::Info << sublist << Foam::endl;

  // or you could try to loop over the elements manually, this won't copy data
  for (auto it = tlist.begin()+2; it != tlist.end(); ++it )
  {
      std::cout << *it << "\n";
  }
  std::cout << std::endl;

Answer 2

与使用迭代器相比，使用特定于 OpenFOAM 的习语更有意义。获取第一个列表，您可以使用 SubList 类获取原始列表的视图。从那里，您可以选择是要使用 SubList，还是要深拷贝到另一个列表中。

编辑（第一部分是从我的手机上潦草写的）：

这是一个简单的示例，说明如何使用 SubList:

instantList allTimes(runTime.times());

Info<< "I have " << allTimes.size() << " times:" << nl
    << allTimes << nl;

if (allTimes.size() > 2)
{
    SubList<instant> someTimes(allTimes, allTimes.size()-2, 2);

    Info<< "Subset of times: " << someTimes << nl;
}

SubList 构造函数的参数顺序如果你想到 string::substr() 可能会有点混乱，但本质上是：

• 原始列表
• 子集大小
• 偏移量（缺失为0）

请注意，这个 SubList 是一个非常轻量级的 视图 到现有列表中，您可以使用它来修改列表的子部分等。如果您有分支情况，您可能想要处理一些事情统一地，您可以简单地放入三元组。例如，

const SubList<instant>& myTimes =
(
    (allTimes.size() > 2)
  ? SubList<instant>(allTimes, allTimes.size()-2, 2)
  : SubList<instant>(allTimes)
);

您也可以使用 labelRange 来选择事物。例如，

labelRange myslice(0, allTimes.size());
if (allTimes.size() > 2)
{
    myslice.reset(2, allTimes.size()-2);
}

// Now
Info<< "Selected times: "
    << SubList<instant>(allTimes, myslice) << nl;

由于 SubList 也是列表的一种类型，因此您可以在其上使用 range-for 等。

请注意，在下一个版本（2021 年 12 月）中，您还可以直接获取视图。例如，

#if (OPENFOAM >= 2112)

Info<< "Selected times: " << allTimes.slice(2) << nl;

#endif