def METHOD_NAME(units, dims):
    bl     = BravaisLattice(units)
    bz     = BrillouinZone(bl)
    k_mesh = MeshBrZone(bz, dims)
    k1, k2, k3 = bz.units
    assert np.array_equal(k_mesh.closest_index(k1), [0,0,0])
    assert np.array_equal(k_mesh.closest_index(k2), [0,0,0])
    assert np.array_equal(k_mesh.closest_index(k3), [0,0,0])
    assert linalg.norm(dot(bl.units, bz.units.T) - eye(3)*2*pi) < 1e-14