Updating a mesh over a function

Question

I got a mesh and a solution of a PDE in FEniCS

omega = #some code for to get a mesh
V = FunctionSpace(omega)
u = TrialFunction(V)
v = TestFunction(V)
a = # ... some bilinear form
L = # ... some linear form
bc = # ... some Dirichlet bound conditions
solve(a = L, u , bc)

In this part we have the PDE solution, but, with this information we would like to redefine the omega mesh for catch better the information of solution u

¿How to redefine the mesh for this goal?

Answer 1

Hi

You may want to take a look at this:
http://fenicsproject.org/featured/2011/automated_error_control.html