I assume the solution is something along the lines of this:
mu = a + b * ( ((x-c) > 0) * (x-c) )
For those still looking for other options, I will also note that you can get similar behavior out of aesara.tensor.clip().
1 Like