def METHOD_NAME(self) -> None:
    """Testing Bickelhaupt charges (unrestricted) against outputs from Multiwfn."""
    data, logfile = getdatafile(Gaussian, "basicGaussian09", ["dvb_un_sp.log"])
    bpa = Bickelhaupt(data)
    bpa.logger.setLevel(logging.ERROR)
    bpa.calculate()
    e_bpaalpha = numpy.loadtxt(f"{os.path.dirname(os.path.realpath(__file__))}/dvb_un_sp.bpa")
    e_bpaspin = numpy.loadtxt(f"{os.path.dirname(os.path.realpath(__file__))}/dvb_un_sp.bpaspin")
    assert numpy.all(bpa.fragcharges >= e_bpaalpha - 0.05)
    assert numpy.all(bpa.fragcharges <= e_bpaalpha + 0.05)
    assert numpy.all(bpa.fragspins >= e_bpaspin - 0.05)
    assert numpy.all(bpa.fragspins <= e_bpaspin + 0.05)