Well, I actually need something like δBμν/δAρσ. At the moment I have to use 'vary' to compute δBμν and then set δAρσ→1 at the end to extract δBμν/δAρσ. But then if I have a more complicated expression such as δBμν/δAρσ+BμνBρσ, I cannot simply substitute whatever expression I have for Bμν into it as the derivatives aren't automatically calculated. So I calculate the derivative separately and substitute in the expression. The problem with that is that I still have to define $\delta{\mu\nu} (= \delta/\delta A^{\mu\nu})asaderivativesothatCadabragetstheindicesrightinthesubstitutions.Since\delta\{\mu\nu}$ is a second derivative, my calculation would clash if I perform an operation such as product_rule
, thus I don't use δμν for anything other than a dummy variable that gets replaced by the derivative I obtained with vary. I was just wondering if there is any easier/cleaner route to this and it seems that a derivative with multiple indices would solve my problem.