double dot product

I have two tensors. How do I compute the double dot product of the two tensors in FreeFEM?