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

希望你能检查同样的事情来解决你的问题。