I am wondering if there is an efficient way to define the convolution of two functions in FEniCS. That is, given functions $f$ and $g$, define the function $h=f*g$ by
$$h(x)=\int_\Omega f(x-y)g(y)dy.$$
It should be possible to do this by defining $h$ at each point in the mesh manually, but this method would be extremely time consuming if one has a large 2D or 3D mesh. I am hoping there is a more efficient way of doing it.
