Im solving a coupled PDE-problem using NeuralPDE.jl, the boundary conditions are now set by the penalization method.
Instead I want to apply these as hard boundary constraints, or as an output transformation, as it is done in DeepXDE in Python.
The only option I’ve found for Julia, is an approach used in this blog post, however this example is for DAEs, and not PDEs.
Can hard boundary constraints be set when using NeuralPDE.jl?