ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/group/trunk/OOPSE-3.0/src/UseTheForce/DarkSide/electrostatic.F90
(Generate patch)

Comparing trunk/OOPSE-3.0/src/UseTheForce/DarkSide/electrostatic.F90 (file contents):
Revision 2123 by gezelter, Sun Mar 13 07:39:01 2005 UTC vs.
Revision 2124 by gezelter, Tue Mar 15 14:20:50 2005 UTC

# Line 550 | Line 550 | contains
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
# Line 676 | Line 676 | contains
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))

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines