diff --git a/src/qmckl_ao.org b/src/qmckl_ao.org index 5511e74..ee2c413 100644 --- a/src/qmckl_ao.org +++ b/src/qmckl_ao.org @@ -299,9 +299,10 @@ integer function qmckl_ao_polynomial_vgl_f(context, X, R, lmax, n, L, ldl, VGL, vgl(3,n) = dble(b) * pows(a ,1) * pows(b-1,2) * pows(c ,3) vgl(4,n) = dble(c) * pows(a ,1) * pows(b ,2) * pows(c-1,3) - vgl(5,n) = dble(a) * dble(a-1) * pows(a-2,1) * pows(b ,2) * pows(c ,3) & - + dble(b) * dble(b-1) * pows(a ,1) * pows(b-2,2) * pows(c ,3) & - + dble(c) * dble(c-1) * pows(a ,1) * pows(b ,2) * pows(c-2,3) + vgl(5,n) = & + dble(a) * dble(a-1) * pows(a-2,1) * pows(b ,2) * pows(c ,3) + & + dble(b) * dble(b-1) * pows(a ,1) * pows(b-2,2) * pows(c ,3) + & + dble(c) * dble(c-1) * pows(a ,1) * pows(b ,2) * pows(c-2,3) exit end if end do