# Definition de la zone interdite [-100;9900]x[-100;9900] avec la centrale nucléaire au centre
def zone_interdite(x,y):
    # x et y doivent etre de type numpy.ndarray
    if any((x >= -100) & (x <= 9.9e3) & (-100 <= y) & (y <= 9.9e3)):
        return True
    else:
        return False


def zone_interdite_margin(x,y):
    x_proj = np.maximum(-100, np.minimum(x, 2e3))
    y_proj = np.maximum(-100, np.minimum(y, 1e3))
    distance = np.sqrt(np.min((x_proj - x) ** 2 + (y_proj - y)**2))
    if distance > 0:
        return distance
    else:  # we are inside in at least one point
        x_x_proj2 = np.minimum(np.abs(np.minimum(-100, x) - x), np.abs(np.maximum(x, 2e3) - x))
        y_y_proj2 = np.minimum(np.abs(np.minimum(-100, y) - y), np.abs(np.maximum(y, 1e3) - y))
        return - np.sqrt(np.max((x_x_proj2) ** 2 + (y_y_proj2)**2))