According to the priority of operators in Fortran it is necessary to use brackets for Tensor operators. For example the power operator ** which is overridden by the double dot product .ddot.has the highest priority. This is wrong, especially for a double contraction where one tensor is a result of a crossed dyadic product: A**B.cdya.C - the crossdyadic product has to be evaluated first. So the correct way of defining the equation is: A**(B.cdya.C)

Improve this behaviour...