Le 29 mars 2016 19:31, "Gregory, Alastair" <a.gregory14@imperial.ac.uk> a écrit :
>
> Hi George,
>
> Sorry about the Easter delay in getting back to you.
>
> You can only project 'w' to other Function Spaces that are on the same original mesh. Here, you are after the ability to prolong / inject between a hierarchy of meshes.
>
> For this you can use: MeshHierarchy(coarsest_mesh,number_of_levels)
>
> So, given your original mesh,
>
> mesh=UnitSquareMesh(size,size,quadrilateral=True)
>
> You can make a Mesh Hierarchy,
> mesh_hierarchy=MeshHierarchy(mesh,L)

That won't work on a quadrilateral mesh, will it?  Only meshes of intervals or triangles?  There's a test at the start of the constructor.