A special finite element method is defined in
`getfem/getfem_interpolated_fem.h` which is not a real finite element
method, but a pseudo-fem which interpolates a finite element method defined on
another mesh. If you need to assemble a matrix with finite element methods
defined on different meshes, you may use the “interpolated fem” or “projected
fem” for that purpose:

```
// interpolation within a volume
getfem::new_interpolated_fem(getfem::mesh_fem mf, getfem::mesh_im mim);
// projection on a surface
getfem::new_projected_fem(getfem::mesh_fem mf, getfem::mesh_im mim);
```

Because each base function of the finite element method has to be interpolated, such a computation can be a heavy procedure. By default, the interpolated fem object store the interpolation data.

The interpolation is made on each Gauss point of the integration methods of
`mim`, so only this integration method can be used in assembly
procedures.

For instance if you need to compute the mass matrix between two different finite
element methods defined on two different meshes, this is an example of code which
interpolate the second FEM. on the mesh of the first FEM., assuming that `mf`
describes the finite element method and `mim` is the chosen integration method:

```
getfem::mesh_fem mf_interpole(mfu.linked_mesh());
pfem ifem = getfem::new_interpolated_fem(mf, mim);
dal::bit_vector nn = mfu.convex_index();
mf_interpole.set_finite_element(nn, ifem);
getfem::asm_mass_matrix(SM1, mim, mfu, mf_interpole);
del_interpolated_fem(ifem);
```

The object pointed by `ifem` contains all the information concerning the
interpolation. It could use a lot of memory. As pfem is a shared_ptr, the
interpolated fem will be automatically destroyed when the last pointer on it is
destroyed. To obtain a better accuracy, it is better to refine the integration
method (with `IM_STRUCTURED_COMPOSITE` for instance) rather than increase its
order.

Instead of using the previous tools (interpolated and projected fems), it is
possible to use a finite element variable defined on an another mesh than the one
on which an assembly is computed using the “interpolate transformation” tool of
the weak form language (see *Interpolate transformations* ), the finite element
variables will be interpolated on each Gauss point. There is no restriction
on the dimensions of the mesh used, which means in particular that a
two-dimensional fem variable can be interpolated on a one-dimensional mesh
(allowing the coupling of shell and beam elements, for instance).
It is also possible to use some transformations like polar coordinates to
euclidean ones.

Mortar methods are supported by *GetFEM++*. The coupling term between non matching
meshes can in particular be computed using the interpolate transformations of
the weak form language (see *Interpolate transformations*).