Solving system of PDEs/ODEs

[ni-sha-c]

Hello, I was wondering if the code currently supports solving a BVP of the form $L(u) = v,$ where $L$ is some differential operator, $u(x,y), v(x,y) \in \mathbb{R}^2$ and $x$ and $y$ are scalar coordinates.
I tried to modify the example BVP to be trained with a neural network with a vector output, but the functions in conditions.py and solvers.py assume that the output of the PDE is scalar.
This package is very nicely written, thank you!

[sathvikbhagavan]

Can't you substitute x, y into u, v and create a system of two equations?
Also, it would be great if you can post an MWE for this case.