Dear colleagues, greetings!
I am trying to demonstrate that a Lagrangian is gauge invariant. The model I am working on is based on a fully symmetric rank-2 field that I have called "h". This theory is set in a curved spacetime. The final result in the attached notebook (please focus only on section 3) should be zero, or at least a total derivative, however, it is not trivial to see this from the xAct output. Any suggestions on what else I could try? Or should I proceed with the calculations by hand?
All the best!