def METHOD_NAME(self):
    self.coords = np.array(
        [
            4.6067455554,
            8.8719311819,
            6.3886531197,
            4.0044515745,
            4.2449530507,
            7.7902855220,
            2.6453069446,
            0.8772647726,
            1.2804446790,
            1.1445332290,
            0.0067366438,
            1.8606485070,
            7.1002867706,
            5.0325506787,
            3.1805888348,
            4.5352891138,
            7.7389683929,
            9.4260970128,
            2.1833238914,
            9.0916071034,
            7.2299906064,
            4.1040157820,
            1.0496745045,
            5.4748315591,
        ],
        dtype=np.float64,
    )
    # 1.1445332290,    0.0067366438,    1.8606485070,
    # 2.1833238914,    9.0916071034,    7.2299906064,
    # 4.0044515745,    4.2449530507,    7.7902855220,
    # 7.1002867706,    5.0325506787,    3.1805888348,
    self.atype = np.array([0, 3, 2, 1, 3, 4, 1, 4], dtype=int)
    self.box = np.array([10.0, 0.0, 0.0, 0.0, 10.0, 0.0, 0.0, 0.0, 10.0])
    self.expected_e = np.array([3.671081837126222158e00])
    self.expected_f = np.array(
        [
            8.786854427753210128e-01,
            -1.590752486903602159e-01,
            -2.709225006303785932e-01,
            -4.449513960033193438e-01,
            -1.564291540964127813e-01,
            2.139031741772115178e-02,
            1.219699614140521193e00,
            -5.580358618499958734e-02,
            -3.878662478349682585e-01,
            -1.286685244990778854e00,
            1.886475802950296488e-01,
            3.904450515493615437e-01,
            1.605017382138404849e-02,
            2.138016869742287995e-01,
            -2.617514921203008965e-02,
            2.877081057057793712e-01,
            -3.846449683844421763e-01,
            3.048855616906603894e-02,
            -9.075632811311897807e-01,
            -6.509653472431625731e-03,
            2.302010972126376787e-01,
            2.370565856822822726e-01,
            3.600133435593881881e-01,
            1.243887532859055609e-02,
        ]
    )
    self.expected_v = np.array(
        [
            3.714071471995848417e-01,
            6.957130186032146613e-01,
            -1.158289779017217302e00,
            6.957130186032139951e-01,
            -1.400130091653774933e01,
            -3.631620234653316626e-01,
            -1.158289779017217302e00,
            -3.631620234653316626e-01,
            3.805077486043773050e00,
        ]
    )
    self.natoms = self.atype.size
    self.coords = self.coords.reshape([-1, self.natoms, 3])