In the jet bundle approach to differential equations
one identifies the equation with the set of a solution of the differential equations.
Now if I would formally want to capture a differential equation, I'd say it comes with a boundary condition as well. Like for example a Dirichlet/Neumann boundary conditions to a diffusion equation.
How are boundary conditions formally captured by the jet bundle approach to differential equations?