def METHOD_NAME():
    grid = make_pupil_grid(2)
    A = Field(np.random.randn(5, 10, grid.size), grid)
    U, S, Vh = field_svd(A)
    u, s, vh = field_svd(A, False)
    for i in range(grid.size):
        svd = np.linalg.svd(A[..., i])
        assert np.allclose(U[..., i], svd[0])
        assert np.allclose(S[..., i], svd[1])
        assert np.allclose(Vh[..., i], svd[2])
        svd2 = np.linalg.svd(A[..., i], full_matrices=False)
        assert np.allclose(u[..., i], svd2[0])
        assert np.allclose(s[..., i], svd2[1])
        assert np.allclose(vh[..., i], svd2[2])