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