$\DeclareMathOperator\SO{SO}$The following PDE defined on $\mathbb{R}^2$ $$\frac{\partial}{\partial x}\frac{\partial}{\partial y}f(x,y) = 0,$$ has solution $$f(x,y) = g(x) + h(y),$$ where $g,h : \mathbb{R} \to \mathbb{R}$ are arbitrary (nice enough) functions. I believe this to be a general solution to the equation, for any choice of boundary conditions. I don't have a proof of this and I'm happy to be shown otherwise if that's the case.
I want to solve a similar looking PDE for $(x,y) \in \mathbb{R}^{n}\times\mathbb{R}^{n},$ which is given by $$\frac{\partial}{\partial x}\cdot\frac{\partial}{\partial y}f(x,y) = 0,$$ where the dot product is with the Euclidean metric on $\mathbb{R}^{n}$. There is an $\SO(n)$ symmetry under $x \mapsto M x, y \mapsto M^T y$. I can see that there are still solutions of the form $f_1(x,y) = g(x) + h(y)$ where $g,h : \mathbb{R}^n \to \mathbb{R}$ , but that these are no longer the only solutions.
I've come up with something quite general that solves the equation;$$f_2(x,y) = \int d^na\,\left( g^x(x \wedge a) \,h^y(y\cdot a) + g^y(y \wedge a) \,h^x(x\cdot a)\right) + f_1(x,y),$$ where $g^x,g^y:\mathbb{R}^n\times\mathbb{R}^n\to \mathbb{R}$ and $h^x,h^y:\mathbb{R} \to \mathbb{R}$ are arbitrary (sufficiently nice) functions, which can depend on $a$. I added on $f_1(x,y)$ explicility, although I think it may be possible to fiddle the form of the integral and the choice of $g^x,g^y,h^x,h^y$ to include solutions of the form $f_1(x,y)$ without adding them on separately. One feature that I like about this solution is that it respects the $\SO(n)$ symmetry of the original PDE. The integral over the vector $a$ is also natural, in the context of the physical problem where the PDE comes from.
Is anyone able to find a general solution to this equation, in a form which respects the $\SO(n)$ symmetry, possibly with some (hopefully quite general) restriction on the boundary conditions?
Is my solution general for some choice of boundary conditions?