In this paper we propose an original implementation for a large family of lowest-order methods to solve diffusive problems with a \lstinline|FreeFEM|-like domain specific language targeted at defining discrete linear and bilinear forms. We discuss how by using the Boost Proto framework we have developped the back-end and the front-end of the language.We validate the proposed DSEL design by the implementation of several academic problems. The overhead of the language is evaluated by comparing with a more traditional implementation.
