def METHOD_NAME():
    # only uses the basis set
    fname = "rot_multi_2det_Be_STO.wfnoj.xml"
    basis, mo = read_qmcpack.parse_qmc_wf(fname, ["Be"])
    wf = Wavefunction_Be_STO(basis["Be"], mo)
    r = np.array([[0.7, 2.0, 3.0], [1.2, 1.5, 0.5], [1.5, 1.6, 1.5], [0.7, 1.0, 1.2]])
    VP = np.zeros(wf.rot_param_size + 1)
    VP[0] = 0.1
    p = wf.psi(r, VP)
    print("psi = ", p)
    print("log psi = ", np.log(abs(p)))
    dp = wf.dpsi(r, VP)
    print("dpsi = ", dp)
    print("dlogpsi = ", dp / p)
    dlogpsi = dp / p
    for i in range(dlogpsi.shape[0]):
        print("    CHECK(dlogpsi[{}] == ValueApprox({}));".format(i, dlogpsi[i]))
    en = wf.local_energy(r, VP)
    print("en = ", en)
    den = wf.dlocal_energy(r, VP)
    print("den = ", den)
    for i in range(den.shape[0]):
        print("   CHECK(dhpsioverpsi[{}] == ValueApprox({}));".format(i, den[i]))