[NotOverUnderated]

Hello,

OpenFOAM makes it very easy to assign values to fields thanks to operator overloading which makes it more natural and avoids explicit looping over all mesh cells, e.g:

// assuming p is volScalarField, we can use
p = 5.0
// Now every mesh cell has a value p=5

Now comes the question that arise naturally. Does this flexibility extends to conditional assignment, i.e. is it possible to assign values to fields conditionally without using an explicit looping (which in Python with numpy can be achieved using the where function), e.g:

# Assuming t and p are two 1D arrays and have the same size
p = numpy.where(t<2.0, 7.0, 0.0)

# That means if t <2.0 assign the value 7.0 to p, otherwise 0.0

My question: is there a similar functionality in OpenFOAM that allows that without explicit looping?

Thank you