def METHOD_NAME(self):
    """
    test get_geometry() in array_analysis.py
    """
    ll = np.array([[24.5797167, 121.4842444, 385.106],
                   [24.5797611, 121.4842333, 384.893],
                   [24.5796694, 121.4842556, 385.106]])
    la = get_geometry(ll)
    np.testing.assert_almost_equal(la[:, 0].sum(), 0., decimal=8)
    np.testing.assert_almost_equal(la[:, 1].sum(), 0., decimal=8)
    np.testing.assert_almost_equal(la[:, 2].sum(), 0., decimal=8)
    ll = np.array([[10., 10., 10.],
                   [0., 5., 5.],
                   [0., 0., 0.]])
    la = get_geometry(ll, coordsys='xy')
    np.testing.assert_almost_equal(la[:, 0].sum(), 0., decimal=8)
    np.testing.assert_almost_equal(la[:, 1].sum(), 0., decimal=8)
    np.testing.assert_almost_equal(la[:, 2].sum(), 0., decimal=8)