auxiliary_functions_dealii Auxiliary functions for dealii such as special tensor operators outer_product_sym outer product/dyadic product/"x" D = outer_product_sym(A,B) = AxB + BxA; D = outer_product_sym(A) = AxA; ToDo check use of C++ "inline" and suchlike for efficiency