The integral mesh and the mesh associated to test functions or unknown functions can be different in the case of varf form

Dear all,

I read from the FreeFem document (page 202, release 4.6) that " The integral mesh and the mesh associated to test functions or unknown functions can be different in the case of
varf form". I am wondering what this means? for example, if we have two overlapped meshes Th and Ths, and

fespace Vh(Th, P1);
Vh u, v;

fespace Vhs(Ths, P1);
Vhs us;

varf a(u, v) = int2d(Th)( u*v) + int2d(Ths)( u*v)

Would the above integral gives the following discretised form automatically:

(A + Pā€™ * B * P ) * u = f

where us = P* u with P is the interpolation matrix between mesh Th and Ths, A and B are discretised matrices from Th and Ths respectively?

Best,
Yongxing.

1 Like