def METHOD_NAME(self):
    """Unit tests with a single cube"""
    # No isosurface
    cube_zero = numpy.zeros((2, 2, 2), dtype=numpy.float32)
    result = marchingcubes.MarchingCubes(cube_zero, 1.)
    self.assertEqual(result.shape, cube_zero.shape)
    self.assertEqual(result.isolevel, 1.)
    self.assertEqual(result.invert_normals, True)
    vertices, normals, indices = result
    self.assertEqual(len(vertices), 0)
    self.assertEqual(len(normals), 0)
    self.assertEqual(len(indices), 0)
    # Cube array dimensions: shape = (dim 0, dim 1, dim2)
    #
    #      dim 0 (Z)
    #        ^
    #        |
    #      4 +------+ 5
    #       /|     /|
    #      / |    / |
    #   6 +------+ 7|
    #     |  |   |  |
    #     |0 +---|--+ 1 -> dim 2 (X)
    #     | /    | /
    #     |/     |/
    #   2 +------+ 3
    #    /
    #   dim 1 (Y)
    # isosurface perpendicular to dim 0 (Z)
    cube = numpy.array(
        (((0., 0.), (0., 0.)),
         ((1., 1.), (1., 1.))), dtype=numpy.float32)
    level = 0.5
    vertices, normals, indices = marchingcubes.MarchingCubes(
        cube, level, invert_normals=False)
    self.assertAllClose(vertices[:, 0], level)
    self.assertAllClose(normals, (1., 0., 0.))
    self.assertEqual(len(indices), 2)
    # isosurface perpendicular to dim 1 (Y)
    cube = numpy.array(
        (((0., 0.), (1., 1.)),
         ((0., 0.), (1., 1.))), dtype=numpy.float32)
    level = 0.2
    vertices, normals, indices = marchingcubes.MarchingCubes(cube, level)
    self.assertAllClose(vertices[:, 1], level)
    self.assertAllClose(normals, (0., -1., 0.))
    self.assertEqual(len(indices), 2)
    # isosurface perpendicular to dim 2 (X)
    cube = numpy.array(
        (((0., 1.), (0., 1.)),
         ((0., 1.), (0., 1.))), dtype=numpy.float32)
    level = 0.9
    vertices, normals, indices = marchingcubes.MarchingCubes(
        cube, level, invert_normals=False)
    self.assertAllClose(vertices[:, 2], level)
    self.assertAllClose(normals, (0., 0., 1.))
    self.assertEqual(len(indices), 2)
    # isosurface normal in dim1, dim 0 (Y, Z) plane
    cube = numpy.array(
        (((0., 0.), (0., 0.)),
         ((0., 0.), (1., 1.))), dtype=numpy.float32)
    level = 0.5
    vertices, normals, indices = marchingcubes.MarchingCubes(cube, level)
    self.assertAllClose(normals[:, 2], 0.)
    self.assertEqual(len(indices), 2)