If $f:A\to \mathbb R^k$ and $g:B\to \mathbb R$ are two functions of class $C^2$ and their composition is well defined.
For $c \in A$ what is $$\frac{\partial^2(g\circ f)}{\partial x_i \partial x_j}(c)$$
Is it just $$\frac{\partial^2(g\circ f)}{\partial x_i \partial x_j}(c)=\frac{\partial^2g}{\partial x_i \partial x_j}(f(c))\frac{\partial^2f}{\partial x_i \partial x_j}(c) $$
If it is how to prove it. I get from $$\frac{\partial^2(g\circ f)}{\partial x_i \partial x_j}(c)=\frac{\partial}{\partial x_i}\frac{\partial(g\circ f)}{ \partial x_j}(c)=\frac{\partial}{\partial x_i}\left( \frac{\partial g}{ \partial x_j}(f(c))\frac{\partial f}{\partial x_j}(c)\right)$$
but not sure how to calculate this, is this a product rule?
Edit this is what I got from 1 day of thought :D
$f : A \to B, g: B \to \mathbb R^m, A\subseteq \mathbb R^n, B\subseteq \mathbb R^k$
$ \begin{aligned}\frac{\partial^2(g\circ f)}{\partial x_i\partial x_j}(x) &=\frac{\partial}{\partial x_i}\left(\frac{\partial(g\circ f)}{\partial x_j}(x)\right) \\&=\frac{\partial}{\partial x_i}\left(\nabla (g\circ f)(x)e_j \right) \\&=\frac{\partial}{\partial x_i}\left(\nabla g(f(x))\nabla f(x)e_j) \right) \\&=\frac{\partial}{\partial x_i}\left(\nabla g(f(x)) \frac{\partial f}{\partial x_j}(x) \right)\\& =\frac{\partial}{\partial x_i}\left(\sum_{p=1}^k \frac{\partial g}{\partial x_p}(f(x))\frac{\partial f_p}{\partial x_j}(x)\right) \\&=\sum_{p=1}^k\frac{\partial}{\partial x_i}\left(\frac{\partial g}{\partial x_p}(f(x))\frac{\partial f_p}{\partial x_j}(x) \right)\\ &=\sum_{p=1}^k\left[ \frac{\partial}{\partial x_i}\left(\frac{\partial g}{\partial x_p}(f(x)) \right)\frac{\partial f_p}{\partial x_j}(x)+\frac{\partial g}{\partial x_p}(f(x)) \frac{\partial^2f_p}{\partial x_i\partial x_j}(x) \right]\\ &= \sum_{p=1}^k \frac{\partial}{\partial x_i}\left(\frac{\partial g}{\partial x_p}(f(x)) \right)\frac{\partial f_p}{\partial x_j}(x)+ \sum_{p=1}^k \frac{\partial g}{\partial x_p}(f(x)) \frac{\partial^2f_p}{\partial x_i\partial x_j}(x)\\ &=\sum_{p=1}^k\sum_{l=1}^k\frac{\partial^2g}{\partial x_l\partial x_p}(f(x))\frac{\partial f_l}{\partial x_i}(x)\frac{\partial f_p}{\partial x_j}(x)+\sum_{p=1}^k \frac{\partial g}{\partial x_p}(f(x)) \frac{\partial^2f_p}{\partial x_i\partial x_j}(x)\end{aligned} $
Also from this we get $\begin{aligned}D^2(g\circ f)(x)(h,k)&=\sum_{i=1}^n\sum_{j=1}^n\frac{\partial^2(g\circ f)}{\partial x_i\partial x_j}(x) h_ik_j\\ &=\sum_{i=1}^n\sum_{j=1}^n\left(\sum_{p=1}^k\sum_{l=1}^k\frac{\partial^2g}{\partial x_l\partial x_p}(f(x))\frac{\partial f_l}{\partial x_i}(x)\frac{\partial f_p}{\partial x_j}(x)+\sum_{p=1}^k \frac{\partial g}{\partial x_p}(f(x)) \frac{\partial^2f_p}{\partial x_i\partial x_j}(x) \right)h_ik_j\\ \\&=\sum_{i=1}^n\sum_{j=1}^n\left(\sum_{p=1}^k\sum_{l=1}^k\frac{\partial^2g}{\partial x_l\partial x_p}(f(x))\frac{\partial f_l}{\partial x_i}(x)h_i\frac{\partial f_p}{\partial x_j}(x)k_j\right)\\+ \sum_{i=1}^n\sum_{j=1}^n\left(\sum_{p=1}^k \frac{\partial g}{\partial x_p}(f(x)) \frac{\partial^2f_p}{\partial x_i\partial x_j}(x)h_ik_j \right)\\ \\&=\sum_{p=1}^k\sum_{l=1}^k\frac{\partial^2g}{\partial x_l\partial x_p}(f(x))\left(\sum_{l=1}^n\frac{\partial f_l}{\partial x_i}(x)h_i\right)\left(\sum_{p=1}^n\frac{\partial f_p}{\partial x_j}(x)k_j \right)\\ +\sum_{p=1}^k \frac{\partial g}{\partial x_p}(f(x))\left(\sum_{i=1}^n\sum_{j=1}^n \frac{\partial^2f_p}{\partial x_i\partial x_j}(x)h_ik_j\right)\\ &=\sum_{p=1}^k\sum_{l=1}^k\frac{\partial^2g}{\partial x_l\partial x_p}(f(x))(Df_l(x)h)(Df_p(x)k)\\+\sum_{p=1}^k \frac{\partial g}{\partial x_p}(f(x))D^2f_p(x)(h,k)\\& =\sum_{p=1}^k\sum_{l=1}^k\frac{\partial^2g}{\partial x_l\partial x_p}(f(x))(Df(x)h)_l(Df(x)k)_p\\+\sum_{p=1}^k \frac{\partial g}{\partial x_p}(f(x))\left(D^2f(x)(h,k)\right)_p\\&=D^2g(f(x))(Df(x)h,Dg(x)k)+Dg(f(x))(D^2f(x)(h,k)) \end{aligned}$
That is we get $$D^2(g\circ f)(x)(h,k)=D^2g(f(x))(Df(x)h,Dg(x)k)+Dg(f(x))(D^2f(x)(h,k))$$