If one of the operand is a MeshFucntion convert it to a scalar
potential.
this is not completely working when the first opernad is a
meshfunction.. I guess due to some numpy.array behavior
* added submesh class
* added library of geometrcal shapes functions
* meshFuncion can be integrated on a submesh
* meshFuncition norm and norm2
Calcualte the integral in the interaction region