Add tensorProduct for DiscreteFunctionP0

The tensor product is performed cell by cell

Also give access to the tensor product in the language. In this
special case vectors have to be of the same dimension since non-square
matrices are not allowed there

Also add related tests