I am working on writing out some boundary conditions which have a few symmetric trace-free components of tensors in them, for example the symmetric trace-free component of xi xj fk (where xi and fk are both vectors), the component of fi xj and also eye(3)*fk.
At the moment in the codes I represent these with
Normal[Symmetrize[Outer[Times, xx , f]]]
Normal[Symmetrize[Outer[Times, eye , f]]];
Normal[Symmetrize[Outer[Times, f , x]]]
However, I need to have these symmetric trace-free components written out explicitly so that I have an expression where the vector f_k is being multiplied by everything else remaining in the expression. Is there some piece of code that would enable me t