I'm trying to wrap my head around an equation that involves the derivative of a second order tensor valued function of a second order tensor, then double-dot producted with another second order tensor: \begin{equation} \frac{\partial\mathbf{E}}{\partial\mathbf{F}}:\delta\mathbf{F} \end{equation} The tensor $\mathbf{F}$ is a gradient (so contra- vs. co-variant comes into this somehow, probably), and $\delta\mathbf{F}$ is an infinitesimal change in $\mathbf{F}$. I'm trying to reconcile this with what I've been learning about index notation, especially with respect to upper vs. lower indices, and which indices are summed over. E.g. an almost-definitely-wrong guess: \begin{equation} \frac{\partial E_{i}^{.j}}{\partial F_{k}^{.l}}\delta F_{j}^{.l} \end{equation}
For some background, this equation is from Finite Element Method Simulation of 3D Deformable Solids, Sifakis & Barbic, 2015. $\mathbf{F}$ is the deformation gradient tensor of continuum mechanics, and the full equation reads $\mathbf{E}=\frac{1}{2}(\mathbf{F}^T\mathbf{F} - \mathbf{I})$, so that $\frac{\partial\mathbf{E}}{\partial\mathbf{F}}:\delta\mathbf{F} = \frac{1}{2}(\delta\mathbf{F}^T\mathbf{F} + \mathbf{F}^T\delta\mathbf{F})$. I figure I may understand how one gets from the lhs to the rhs if I can play with the equation in index form.