3
#include <libciomr/libciomr.h>
4
#include <libdpd/dpd.h>
14
if(params.ref == 0) { /*** RHF ***/
15
dpd_buf4_init(&E, CC_EINTS, 0, 11, 0, 11, 0, 0, "E <ai|jk>");
16
dpd_buf4_scmcopy(&E, CC_EINTS, "E 2<ai|jk> - <ai|kj>", 2);
17
dpd_buf4_sort(&E, CC_TMP0, pqsr, 11, 0, "E <ai|kj>");
20
dpd_buf4_init(&E, CC_EINTS, 0, 11, 0, 11, 0, 0, "E 2<ai|jk> - <ai|kj>");
21
dpd_buf4_init(&E1, CC_TMP0, 0, 11, 0, 11, 0, 0, "E <ai|kj>");
22
dpd_buf4_axpy(&E1, &E, -1);