The idea would be to use a special boundary condition to allow setting the value, derivative, or directly the value of the virtual point using a python function (which must be jitable for numbs-compiled cases). This could build on ExpressionBC and the interface would probably look something like this:
def bc_value(value, dx, x, y, z, t):
return [some expression of the variables]
bc={'value_function': bc_value}
The idea would be to use a special boundary condition to allow setting the value, derivative, or directly the value of the virtual point using a python function (which must be jitable for numbs-compiled cases). This could build on
ExpressionBCand the interface would probably look something like this:A use case is demonstrated in #365.