def crossings(y, value, x=None):
    """Finds the crossing points where a (poly-)line defined by a 1D y array has the given
    values and return the (linearly) interpolated index or x value if an x array is given 
    """
    raise NotImplementedError
