def METHOD_NAME(self):
    eos = EOS(self.eos)
    eos_fit = eos.fit(self.volumes, self.energies)
    bulk_modulus = float(str(eos_fit.b0_GPa).split()[0])
    bulk_modulus_ans = float(str(self.qhda.bulk_modulus).split()[0])
    assert_allclose(bulk_modulus, bulk_modulus_ans, atol=1e-3)