[x] Create wrapper type around Array{T,3} to make it behave like a convolution operator
[] Properly implement the convolution of two piecewise polynomial time basis functions
[] assemble and allocatestorage for tensor operators need to be more robust