code: def constrain(value, lower=-np.Inf, upper=np.Inf, allow_equal=False):
      
  
      ok = flib.constrain(value, lower, upper, allow_equal)
      if ok == 0:
          raise ZeroProbability
docs: Apply interval constraint on stochastic value.