| 49 |
|
&Amat, |
| 50 |
|
&mu, |
| 51 |
|
&ul, |
| 52 |
< |
&rc, |
| 53 |
< |
&massRatio); |
| 52 |
> |
&quat); |
| 53 |
|
} |
| 54 |
|
else{ |
| 55 |
|
sprintf( painCave.errMsg, |
| 149 |
|
} |
| 150 |
|
} |
| 151 |
|
|
| 152 |
< |
void Atom::addFrc( double theF[3] ){ |
| 152 |
> |
void Atom::setFrc( double theF[3] ){ |
| 153 |
|
|
| 154 |
|
if( hasCoords ){ |
| 155 |
< |
frc[offsetX] += theF[0]; |
| 156 |
< |
frc[offsetY] += theF[1]; |
| 157 |
< |
frc[offsetZ] += theF[2]; |
| 155 |
> |
frc[offsetX] = theF[0]; |
| 156 |
> |
frc[offsetY] = theF[1]; |
| 157 |
> |
frc[offsetZ] = theF[2]; |
| 158 |
|
} |
| 159 |
|
else{ |
| 160 |
|
|
| 161 |
|
sprintf( painCave.errMsg, |
| 162 |
< |
"Attempt to add frc for atom %d before coords set.\n", |
| 162 |
> |
"Attempt to set frc for atom %d before coords set.\n", |
| 163 |
|
index ); |
| 164 |
|
painCave.isFatal = 1; |
| 165 |
|
simError(); |
| 166 |
|
} |
| 167 |
|
} |
| 168 |
|
|
| 169 |
< |
|
| 171 |
< |
void Atom::zeroForces( void ){ |
| 169 |
> |
void Atom::addFrc( double theF[3] ){ |
| 170 |
|
|
| 171 |
|
if( hasCoords ){ |
| 172 |
< |
frc[offsetX] = 0.0; |
| 173 |
< |
frc[offsetY] = 0.0; |
| 174 |
< |
frc[offsetZ] = 0.0; |
| 172 |
> |
frc[offsetX] += theF[0]; |
| 173 |
> |
frc[offsetY] += theF[1]; |
| 174 |
> |
frc[offsetZ] += theF[2]; |
| 175 |
|
} |
| 176 |
|
else{ |
| 177 |
|
|
| 178 |
|
sprintf( painCave.errMsg, |
| 179 |
< |
"Attempt to zero frc for atom %d before coords set.\n", |
| 179 |
> |
"Attempt to add frc for atom %d before coords set.\n", |
| 180 |
|
index ); |
| 181 |
|
painCave.isFatal = 1; |
| 182 |
|
simError(); |
| 183 |
|
} |
| 184 |
|
} |
| 185 |
|
|
| 188 |
– |
void Atom::getRc(double theRc[3]){ |
| 186 |
|
|
| 187 |
+ |
void Atom::zeroForces( void ){ |
| 188 |
+ |
|
| 189 |
|
if( hasCoords ){ |
| 190 |
< |
theRc[0] = rc[offsetX]; |
| 191 |
< |
theRc[1] = rc[offsetY]; |
| 192 |
< |
theRc[2] = rc[offsetZ]; |
| 190 |
> |
frc[offsetX] = 0.0; |
| 191 |
> |
frc[offsetY] = 0.0; |
| 192 |
> |
frc[offsetZ] = 0.0; |
| 193 |
|
} |
| 194 |
|
else{ |
| 195 |
|
|
| 196 |
|
sprintf( painCave.errMsg, |
| 197 |
< |
"Attempt to get rc for atom %d before coords set.\n", |
| 197 |
> |
"Attempt to zero frc for atom %d before coords set.\n", |
| 198 |
|
index ); |
| 199 |
|
painCave.isFatal = 1; |
| 200 |
|
simError(); |
| 201 |
|
} |
| 203 |
– |
|
| 202 |
|
} |
| 203 |
|
|
| 206 |
– |
void Atom::setRc(double theRc[3]){ |
| 207 |
– |
if( hasCoords ){ |
| 208 |
– |
rc[offsetX] = theRc[0] ; |
| 209 |
– |
rc[offsetY] = theRc[1]; |
| 210 |
– |
rc[offsetZ] = theRc[2]; |
| 211 |
– |
} |
| 212 |
– |
else{ |
| 213 |
– |
|
| 214 |
– |
sprintf( painCave.errMsg, |
| 215 |
– |
"Attempt to set rc for atom %d before coords set.\n", |
| 216 |
– |
index ); |
| 217 |
– |
painCave.isFatal = 1; |
| 218 |
– |
} |
| 219 |
– |
} |