For a matrix A, is there a function f(A) such that $$\frac{df_{ij}}{dA_{kl}}=f_i^k (A) \delta^l_j$$ ?
I am specifically interested in the case of a symmetric matrix $A$ and function f of rank 2.
Prior research: I tried $e^A$ but it does not have this property. The ij element of $e^A$ is $$e^A_{ij} = \sum_{n=0}^\infty \frac{1}{n!} A_i^{k_1} ... A_{k_{n-1}j} ~.$$ Each matrix element can be seen as an independent variable, so the derivative toward $A^{kl}$ is $$e^A_{ij} = \sum_{n=0}^\infty \frac{1}{n!} \sum_{p=0}^{n-1} A^p_{ik} A^{n-1-p}_{lj} ~.$$ This does not have the desired form.