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 2153 by chrisfen, Fri Apr 8 14:43:06 2005 UTC vs.
Revision 2156 by chrisfen, Fri Apr 8 22:13:41 2005 UTC

# Line 564 | Line 564 | contains
564            vpair = vpair + vterm
565            epot = epot + sw * vterm
566  
567 <          dudx = dudx - 5.0_dp*sw*vterm*riji*xhat - pref * sw * ri4 * ( &
567 >          dudx = dudx - 5.0_dp*sw*vterm*riji*xhat + pref * sw * ri4 * ( &
568                 qxx_j*(6.0_dp*cx_j*ux_j(1) - 2.0_dp*xhat) + &
569                 qyy_j*(6.0_dp*cy_j*uy_j(1) - 2.0_dp*xhat) + &
570                 qzz_j*(6.0_dp*cz_j*uz_j(1) - 2.0_dp*xhat) )
571 <          dudy = dudy - 5.0_dp*sw*vterm*riji*yhat - pref * sw * ri4 * ( &
571 >          dudy = dudy - 5.0_dp*sw*vterm*riji*yhat + pref * sw * ri4 * ( &
572                 qxx_j*(6.0_dp*cx_j*ux_j(2) - 2.0_dp*yhat) + &
573                 qyy_j*(6.0_dp*cy_j*uy_j(2) - 2.0_dp*yhat) + &
574                 qzz_j*(6.0_dp*cz_j*uz_j(2) - 2.0_dp*yhat) )
575 <          dudz = dudz - 5.0_dp*sw*vterm*riji*zhat - pref * sw * ri4 * ( &
575 >          dudz = dudz - 5.0_dp*sw*vterm*riji*zhat + pref * sw * ri4 * ( &
576                 qxx_j*(6.0_dp*cx_j*ux_j(3) - 2.0_dp*zhat) + &
577                 qyy_j*(6.0_dp*cy_j*uy_j(3) - 2.0_dp*zhat) + &
578                 qzz_j*(6.0_dp*cz_j*uz_j(3) - 2.0_dp*zhat) )
# Line 690 | Line 690 | contains
690            vpair = vpair + vterm
691            epot = epot + sw * vterm
692            
693 <          dudx = dudx - 5.0_dp*sw*vterm*riji*xhat - pref * sw * ri4 * ( &
693 >          dudx = dudx - 5.0_dp*sw*vterm*riji*xhat + pref * sw * ri4 * ( &
694                 qxx_i*(6.0_dp*cx_i*ux_i(1) - 2.0_dp*xhat) + &
695                 qyy_i*(6.0_dp*cy_i*uy_i(1) - 2.0_dp*xhat) + &
696                 qzz_i*(6.0_dp*cz_i*uz_i(1) - 2.0_dp*xhat) )
697 <          dudy = dudy - 5.0_dp*sw*vterm*riji*yhat - pref * sw * ri4 * ( &
697 >          dudy = dudy - 5.0_dp*sw*vterm*riji*yhat + pref * sw * ri4 * ( &
698                 qxx_i*(6.0_dp*cx_i*ux_i(2) - 2.0_dp*yhat) + &
699                 qyy_i*(6.0_dp*cy_i*uy_i(2) - 2.0_dp*yhat) + &
700                 qzz_i*(6.0_dp*cz_i*uz_i(2) - 2.0_dp*yhat) )
701 <          dudz = dudz - 5.0_dp*sw*vterm*riji*zhat - pref * sw * ri4 * ( &
701 >          dudz = dudz - 5.0_dp*sw*vterm*riji*zhat + pref * sw * ri4 * ( &
702                 qxx_i*(6.0_dp*cx_i*ux_i(3) - 2.0_dp*zhat) + &
703                 qyy_i*(6.0_dp*cy_i*uy_i(3) - 2.0_dp*zhat) + &
704                 qzz_i*(6.0_dp*cz_i*uz_i(3) - 2.0_dp*zhat) )

Diff Legend

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