I would like to see a detailed, step-by-step derivation of the following identity
$$\frac{\partial}{\partial A} \left( y^T A x \right) = y x^T$$
where $x, y \in \mathbb R^n$ and $A \in \mathbb R^{n \times n}$. I thought it would be easy to do using Einstein notation, but I am messing up with the reciprocal basis.
Since $y^\intercal A x = \sum_{i,j} a_{ij}y_i x_j$, where $A = (a_{ij})$, we can easily compute element-wise \begin{align*} \frac{\partial y^\intercal A x}{\partial a_{kl}} = y_k x_l \end{align*} So clearly \begin{align*} \frac{\partial y^\intercal A x}{\partial A} = \left(\frac{\partial y^\intercal A x}{\partial a_{kl}}\right)_{kl} = (y_k x_l)_{kl} = y x^\intercal \end{align*}