In the process of time iteration, the finite element solution obtained at each step is actually the value on the grid node, but (3,0.5) is not the grid node. Why can it be obtained?
mesh Th = square(30, 5, [6.*x,y]);
…
ff << u(3., 0.5) << endl;
In the process of time iteration, the finite element solution obtained at each step is actually the value on the grid node, but (3,0.5) is not the grid node. Why can it be obtained?
mesh Th = square(30, 5, [6.*x,y]);
…
ff << u(3., 0.5) << endl;