Dear all, How do I properly construct the node subset given function values? I'd like to get something like this: a = Function(V) # ... initialize a with some values ... b = Function(V) b.assign(1/a, subset=["where a != 0"] ) Thank you, Tomasz