|
2 | 2 | #include <chrono> |
3 | 3 | #include <fmt/format.h> |
4 | 4 |
|
5 | | -using namespace dlr2d; |
| 5 | +using namespace cppdlr2d; |
6 | 6 |
|
7 | 7 | nda::vector<double> hubatom_allfuncs(double beta, double u, double lambda, |
8 | 8 | double eps, int niomtst, int nbos_tst, |
@@ -597,26 +597,20 @@ nda::vector<double> hubatom_allfuncs_3term(double beta, double u, double lambda, |
597 | 597 | // polarization_res(beta, ifops_fer, ifops_bos, gc, gc, lam_s_c, |
598 | 598 | // lam_s_csing); |
599 | 599 | // pol_s += polarization_const(beta, itops, ifops_bos, gc, gc); |
600 | | - auto pol_s = |
601 | | - dlr2d::polarization_3term(beta, lambda, eps, itops, ifops_fer, ifops_bos, |
602 | | - gc, gc, lam_s_c, lam_s_csing); |
| 600 | + auto pol_s = cppdlr2d::polarization_3term(beta, lambda, eps, itops, ifops_fer, ifops_bos, gc, gc, lam_s_c, lam_s_csing); |
603 | 601 | pol_s *= -1.0 / 2; |
604 | 602 |
|
605 | 603 | // auto pol_d = |
606 | 604 | // polarization_res(beta, ifops_fer, ifops_bos, grc, gc, lam_d_c, |
607 | 605 | // lam_d_csing); |
608 | 606 | // pol_d += polarization_const(beta, itops, ifops_bos, grc, gc); |
609 | | - auto pol_d = |
610 | | - dlr2d::polarization_3term(beta, lambda, eps, itops, ifops_fer, ifops_bos, |
611 | | - grc, gc, lam_d_c, lam_d_csing); |
| 607 | + auto pol_d = cppdlr2d::polarization_3term(beta, lambda, eps, itops, ifops_fer, ifops_bos, grc, gc, lam_d_c, lam_d_csing); |
612 | 608 |
|
613 | 609 | // auto pol_m = |
614 | 610 | // polarization_res(beta, ifops_fer, ifops_bos, grc, gc, lam_m_c, |
615 | 611 | // lam_m_csing); |
616 | 612 | // pol_m += polarization_const(beta, itops, ifops_bos, grc, gc); |
617 | | - auto pol_m = |
618 | | - dlr2d::polarization_3term(beta, lambda, eps, itops, ifops_fer, ifops_bos, |
619 | | - grc, gc, lam_m_c, lam_m_csing); |
| 613 | + auto pol_m = cppdlr2d::polarization_3term(beta, lambda, eps, itops, ifops_fer, ifops_bos, grc, gc, lam_m_c, lam_m_csing); |
620 | 614 |
|
621 | 615 | auto pol_s_c = ifops_bos.vals2coefs(beta, pol_s); // DLR expansion |
622 | 616 | auto pol_d_c = ifops_bos.vals2coefs(beta, pol_d); |
|
0 commit comments