mesh - FreeFem++ 中的自适应网格

Question

我正在尝试通过在某些时间步调整我的初始网格来改进我的 FreeFem++ 代码，但 FreeFem++ 出现了以下错误消息：

Exec error : Try to get unset x,y, ...

在网格确实适应之后。在计算新网格中的解决方案时似乎出现了问题，但在其他情况下我没有遇到这个问题。这就是我在 FreeFem++ 中调整网格并更新变量的方式：

Th = adaptmesh(Th,[u1,u2]);
plot(Th);
u1 = u1;
u2 = u2;
p  = p;

但现在我对如何在收到警告消息后继续操作感到困惑。

有谁知道如何解决这个问题？任何建议将不胜感激。

谢谢！

Answer 1

我知道这是一个非常古老的问题，但我下面的答案是为那些仍然想要的人准备的。

cout我在尝试输出 ( ) 一个 FE 变量时遇到了同样的错误uh，如果我们尝试一下，它就会解决uh[]。在这种情况下，它是一个数组。

Vh uh = x;
cout << uh << endl; // error
cout << uh[] << endl; // good

希望您可以针对您的问题检查相同的内容。