3
\brief Enter brief description of file here
6
#include <libdpd/dpd.h>
12
namespace psi { namespace cclambda {
16
dpdfile2 R1, L1, LIA, Lia, RIA, Ria;
17
dpdbuf4 R2, L2, LIJAB, Lijab, LIjAb, RIJAB, Rijab, RIjAb;
20
if(params.ref == 0 || params.ref == 1) { /** RHF/ROHF **/
22
dpd_file2_init(&LIA, CC_LAMBDA, L_irr, 0, 1, "New LIA");
23
dpd_file2_init(&Lia, CC_LAMBDA, L_irr, 0, 1, "New Lia");
24
dpd_buf4_init(&LIJAB, CC_LAMBDA, L_irr, 2, 7, 2, 7, 0, "New LIJAB");
25
dpd_buf4_init(&Lijab, CC_LAMBDA, L_irr, 2, 7, 2, 7, 0, "New Lijab");
26
dpd_buf4_init(&LIjAb, CC_LAMBDA, L_irr, 0, 5, 0, 5, 0, "New LIjAb");
28
norm = dpd_file2_dot_self(&LIA);
29
norm += dpd_file2_dot_self(&Lia);
30
norm += dpd_buf4_dot_self(&LIJAB);
31
norm += dpd_buf4_dot_self(&Lijab);
32
norm += dpd_buf4_dot_self(&LIjAb);
33
fprintf(outfile,"size of L <L|L> %15.10lf\n",norm);
35
dpd_file2_close(&LIA);
36
dpd_file2_close(&Lia);
37
dpd_buf4_close(&LIJAB);
38
dpd_buf4_close(&Lijab);
39
dpd_buf4_close(&LIjAb);
43
}} // namespace psi::cclambda