ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/group/trunk/mdtools/md_code/Bend.cpp
(Generate patch)

Comparing trunk/mdtools/md_code/Bend.cpp (file contents):
Revision 11 by mmeineke, Tue Jul 9 18:40:59 2002 UTC vs.
Revision 134 by chuckv, Fri Oct 11 15:09:09 2002 UTC

# Line 24 | Line 24 | void Bend::calc_forces(){
24  
25    double comf2, comf3, comf4;
26    double dcsidx, dcsidy, dcsidz, dcskdx, dcskdy, dcskdz;
27 <  double dcsjdx, dcsjdy, dcsjdz;
27 >  // double dcsjdx, dcsjdy, dcsjdz;
28    double dadxi, dadyi, dadzi;
29 <  double dadxk, dadyk, dadzk, dadxj, dadyj, dadzj;
29 >  double dadxk, dadyk, dadzk;//, dadxj, dadyj, dadzj;
30    double daxi, dayi, dazi, daxk, dayk, dazk, daxj, dayj, dazj;
31    
32  
# Line 80 | Line 80 | void Bend::calc_forces(){
80    dcskdy = dy*denom - comf3*gy;
81    dcskdz = dz*denom - comf3*gz;
82    
83 <  dcsjdx = -dcsidx - dcskdx;
84 <  dcsjdy = -dcsidy - dcskdy;
85 <  dcsjdz = -dcsidz - dcskdz;
83 > //   dcsjdx = -dcsidx - dcskdx;
84 > //   dcsjdy = -dcsidy - dcskdy;
85 > //   dcsjdz = -dcsidz - dcskdz;
86  
87    dadxi = -sinai*dcsidx;
88    dadyi = -sinai*dcsidy;
# Line 92 | Line 92 | void Bend::calc_forces(){
92    dadyk = -sinai*dcskdy;
93    dadzk = -sinai*dcskdz;
94  
95 <  dadxj = -dadxi - dadxk;
96 <  dadyj = -dadyi - dadyk;
97 <  dadzj = -dadzi - dadzk;
95 > //   dadxj = -dadxi - dadxk;
96 > //   dadyj = -dadyi - dadyk;
97 > //   dadzj = -dadzi - dadzk;
98  
99    daxi = comf4*dadxi;
100    dayi = comf4*dadyi;

Diff Legend

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