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

希望你能为你的问题检查同样的事情。