35 |
|
&trq, |
36 |
|
&Amat, |
37 |
|
&mu, |
38 |
< |
&ul); |
38 |
> |
&ul, |
39 |
> |
&quat); |
40 |
|
} |
41 |
|
else{ |
42 |
|
sprintf( painCave.errMsg, |
432 |
|
|
433 |
|
sprintf( painCave.errMsg, |
434 |
|
"Attempt to get Trq for atom %d before coords set.\n", |
435 |
+ |
index ); |
436 |
+ |
painCave.isFatal = 1; |
437 |
+ |
simError(); |
438 |
+ |
} |
439 |
+ |
} |
440 |
+ |
|
441 |
+ |
void DirectionalAtom::setTrq( double theT[3] ){ |
442 |
+ |
|
443 |
+ |
if( hasCoords ){ |
444 |
+ |
trq[offsetX] = theT[0]; |
445 |
+ |
trq[offsetY] = theT[1]; |
446 |
+ |
trq[offsetZ] = theT[2]; |
447 |
+ |
} |
448 |
+ |
else{ |
449 |
+ |
|
450 |
+ |
sprintf( painCave.errMsg, |
451 |
+ |
"Attempt to add Trq for atom %d before coords set.\n", |
452 |
|
index ); |
453 |
|
painCave.isFatal = 1; |
454 |
|
simError(); |
557 |
|
|
558 |
|
|
559 |
|
double phi,theta,psi,eps; |
560 |
< |
double pi; |
543 |
< |
double cphi,ctheta,cpsi; |
544 |
< |
double sphi,stheta,spsi; |
545 |
< |
double b[3]; |
546 |
< |
int flip[3]; |
560 |
> |
double ctheta,stheta; |
561 |
|
|
562 |
|
// set the tolerance for Euler angles and rotation elements |
563 |
|
|