Thanks mauro3 for your relpy.
Is i correct that isoutofdomain=(u,p,t) -> any(x -> x < 0, u) ensure that all variables to be within the positive domain?
If so, is there a way to specify that u[2] > 0 and is free u[2]?
Thanks mauro3 for your relpy.
Is i correct that isoutofdomain=(u,p,t) -> any(x -> x < 0, u) ensure that all variables to be within the positive domain?
If so, is there a way to specify that u[2] > 0 and is free u[2]?