550 |
|
if (j_is_Quadrupole) then |
551 |
|
ri2 = riji * riji |
552 |
|
ri3 = ri2 * riji |
553 |
< |
ri4 = ri2 * ri4 |
553 |
> |
ri4 = ri2 * ri2 |
554 |
|
cx2 = cx_j * cx_j |
555 |
|
cy2 = cy_j * cy_j |
556 |
|
cz2 = cz_j * cz_j |
676 |
|
|
677 |
|
ri2 = riji * riji |
678 |
|
ri3 = ri2 * riji |
679 |
< |
ri4 = ri2 * ri4 |
679 |
> |
ri4 = ri2 * ri2 |
680 |
|
cx2 = cx_i * cx_i |
681 |
|
cy2 = cy_i * cy_i |
682 |
|
cz2 = cz_i * cz_i |
683 |
|
|
684 |
< |
pref = pre14 * q_i / 6.0_dp |
684 |
> |
pref = pre14 * q_j / 6.0_dp |
685 |
|
vterm = pref * ri3 * (qxx_i * (3.0_dp*cx2 - 1.0_dp) + & |
686 |
|
qyy_i * (3.0_dp*cy2 - 1.0_dp) + & |
687 |
|
qzz_i * (3.0_dp*cz2 - 1.0_dp)) |