def METHOD_NAME(vector_array):
    vector_norms = np.linalg.norm(vector_array, ord=2, axis=1, keepdims=True)
    non_zero_norms = np.reshape(vector_norms != 0, newshape = len(vector_array))
    normalized_vectors = np.copy(vector_array)
    normalized_vectors[non_zero_norms] /= vector_norms[non_zero_norms]
    return normalized_vectors