I have a 3D domain ( a box) , on which I'm trying to solve the following PDE
$
u_z = u_{xx} + u_{yy}
$
All FeniCS demos use a PDE, which has the full Laplace operator $\Delta u$ or the full gradient $\nabla u$ in it, and in my example I need to separate all variables, because I take fist derivative with respect to $z$ , and second derivatives with respect to the other variables $x$ and $y$.
How do I define in Fenics derivatives with respect to ${x,y,z}$ separately ?
I mean I would like to code something like
inner(xderivative(u), v) * dx
inner(yderivative(u), yderivative(v)) * dx