I am looking for results that show smooth dependence of a solution to a parabolic equation, from the initial data.

More specifically I have the following problem:

CONSIDER spaces $P:=\mathbb{R}^k$ ("parameter space"), $\Omega:=\mathbb{R}^n$ ("actual space") and $T=[0,t_0]$ ("time").

Consider also $u_0,f_0\in C^{\infty}(P\times \Omega)$, and $L$ a linear elliptic operator on $\Omega$, whose components depend smoothly on $P$. For every $p\in P$, there is a solution $u_p\in C^{\infty}(\Omega\times T)$ of the PDE

${\partial\over \partial t}u_p+L_p u_p=f_p \qquad on\;\{p\}\times\Omega\times T$

$u_p=(u_0)_p\qquad on\; \{p\}\times\Omega\times\{0\}$

QUESTION: is the function $\overline{u}(p,x,t):=u_p(x,t)$ a function in $C^{\infty}(P\times\Omega\times T)$?

The problem comes from the following, geometric situation: I have a space $P\times \Omega$ with a metric on it, and $L$ is just the Laplacian $\Delta$ of $P\times \Omega$, restricted to the levels $\{p\}\times \Omega$.

Any reference will be very welcome. Thanks in advance.