def METHOD_NAME(egrid, axes):
    nx, ny, _ = egrid.dimensions
    assume(egrid.valid_mapaxes(egrid.mapaxes))
    coord = np.swapaxes(egrid.coord.reshape((ny + 1, nx + 1, 6)), 0, 1).astype(
        np.float64
    )
    assert_allclose(
        inverse_transform_xtgeo_coord_by_mapaxes(
            axes, transform_xtgeo_coord_by_mapaxes(axes, coord.copy())
        ),
        coord,
        atol=1e-7,
    )