code: def test_grid():
      # USGS Professional Paper 1395, p 89, Table 14
      globe = ccrs.Globe(semimajor_axis=1.0, ellipse=None)
      mill = ccrs.Miller(central_longitude=0.0, globe=globe)
      geodetic = mill.as_geodetic()
  
      other_args = {"a=1.0", "lon_0=0.0"}
      check_proj_params("mill", mill, other_args)
  
      assert_almost_equal(
          np.array(mill.x_limits), [-3.14159265, 3.14159265]
      )
      assert_almost_equal(
          np.array(mill.y_limits), [-2.3034125, 2.3034125]
      )
  
      lats, lons = np.mgrid[0:91:5, 0:91:10].reshape((2, -1))
      expected_x = np.deg2rad(lons)
      expected_y = np.array(
          [
              2.30341,
              2.04742,
              1.83239,
              1.64620,
              1.48131,
              1.33270,
              1.19683,
              1.07113,
              0.95364,
              0.84284,
              0.73754,
              0.63674,
              0.53962,
              0.44547,
              0.35369,
              0.26373,
              0.17510,
              0.08734,
              0.00000,
          ]
      )[::-1].repeat(10)
  
      result = mill.transform_points(geodetic, lons, lats)
      assert_almost_equal(result[:, 0], expected_x, decimal=5)
      assert_almost_equal(result[:, 1], expected_y, decimal=5)
  
file_path: cartopy/lib/cartopy/tests/crs/test_miller.py
imports: import numpy
  from numpy.testing import assert_almost_equal
  import pytest
  import cartopy.crs
  from .helpers import check_proj_params
label_desc: Use image node_id 6 for calling a global function with example usage: test_grid() without return types