-#ifdef CONSTRAINT_110_5832
- if(count==5833) {
- atom[5832].f.x=0.5*(atom[5832].f.x-atom[5832].f.y);
- atom[5832].f.y=-atom[5832].f.x;
- }
-#endif
-#ifdef CONSTRAINT_11X_5832
- if(count==5833) {
- // second trafo
- xt=atom[5832].f.x;
- yt=atom[5832].f.y*cos(-0.16935129)-atom[5832].f.z*sin(-0.16935129);
- zt=atom[5832].f.y*sin(-0.16935129)+atom[5832].f.z*cos(-0.16935129);
- // first trafo
- xtt=xt*cos(M_PI/4.0)-yt*sin(M_PI/4.0);
- ytt=xt*sin(M_PI/4.0)+yt*sin(M_PI/4.0);
- ztt=zt;
- // apply constraints
- ytt=0.0;
- // first trafo backwards
- xt=xtt*cos(M_PI/4.0)+ytt*sin(M_PI/4.0);
- yt=-xtt*sin(M_PI/4.0)+ytt*sin(M_PI/4.0);
- zt=ztt;
- // second trafo backwards
- atom[5832].f.x=xt;
- atom[5832].f.y=yt*cos(-0.16935129)+zt*sin(-0.16935129);
- atom[5832].f.z=-yt*sin(-0.16935129)+zt*cos(-0.16935129);
- }
-#endif