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

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