Dear Tomasz,


I would like to check the convergence of my solution as I refine the mesh.

Say, I have two solutions in CG1 space, f_coarse on mesh_coarse and f_fine on mesh_fine, which is a refined mesh_coarse. Can I map/interpolate/project f_fine onto mesh_coarse in Firedrake?


You can do this through the multigrid tools.


# make an empty coarse function (Here, Vc is the coarse FS)

f_coarse_ = Function(Vc)


# project fine onto coarse

inject(f_fine, f_coarse_)


# check error

norm(assemble(f_coarse - f_coarse_))


Many Thanks,

Alastair Gregory
-------------------------------------------------------------------
Research Assistant (Maths Foresees Grant)
Imperial College London
Office 759 Huxley Building, South Kensington
(Tel: 07794 243913) | (Email: a.gregory14@imperial.ac.uk)
------------------------------------------------------------------------------------
From: firedrake-bounces@imperial.ac.uk <firedrake-bounces@imperial.ac.uk> on behalf of Tomasz Salwa [RPG] <mmtjs@leeds.ac.uk>
Sent: 25 October 2016 11:25:45
To: firedrake
Subject: [firedrake] mapping between two meshes
 

Dear all,


I would like to check the convergence of my solution as I refine the mesh.

Say, I have two solutions in CG1 space, f_coarse on mesh_coarse and f_fine on mesh_fine, which is a refined mesh_coarse. Can I map/interpolate/project f_fine onto mesh_coarse in Firedrake?


Thank you,

Tomasz