def METHOD_NAME(s0, x0, x1, points, out):
    S = points.shape[0]
    X = x0.shape[0]
    for n in range(X):
        out[n] = x0[n]
        for i in range(S):
            out[n] += x1[n, i] * (points[i] - s0[i])
    # return out