21 Oct
2016
21 Oct
'16
12:35 p.m.
On 21 Oct 2016, at 12:32, Schwedes, Tobias <t.schwedes14@imperial.ac.uk> wrote:
from firedrake import *
mesh = UnitSquareMesh(3, 3)
# Make a viz mesh that refines twice from the solution mesh viz_mesh = MeshHierarchy(mesh, refinement_levels=2)[-1]
V = FunctionSpace(mesh, "CG", 2) Vviz = FunctionSpace(viz_mesh, "CG", 2)
# solve on mesh for "u" u = Function(V) u_viz = Function(Vviz)
prolong(u, u_viz)
This works in current Firedrake, please update. Cheers, Lawrence