Trisurf Monte Carlo simulator
Samo Penic
2014-04-28 460c2a326b332e1b3e03621970a1b48dfb872050
src/shdiscover.c
@@ -42,17 +42,17 @@
/*we calculate new position of each vertex of vesicle */
for(i=0;i<n;i++){
   fi=atan2(vlist->vtx[i]->data->y, vlist->vtx[i]->data->x);
   fi=atan2(vlist->vtx[i]->y, vlist->vtx[i]->x);
/*   theta=atan2(
       sqrt(vlist->vtx[i]->data->x*vlist->vtx[i]->data->x + 
      vlist->vtx[i]->data->y*vlist->vtx[i]->data->y),
      vlist->vtx[i]->data->z 
       ); */
   theta=acos(
      vlist->vtx[i]->data->z /
       sqrt(vlist->vtx[i]->data->x*vlist->vtx[i]->data->x +
      vlist->vtx[i]->data->y*vlist->vtx[i]->data->y+
      vlist->vtx[i]->data->z*vlist->vtx[i]->data->z)
      vlist->vtx[i]->z /
       sqrt(vlist->vtx[i]->x*vlist->vtx[i]->x +
      vlist->vtx[i]->y*vlist->vtx[i]->y+
      vlist->vtx[i]->z*vlist->vtx[i]->z)
      );
@@ -67,9 +67,9 @@
      /*ts_fprintf(stdout,"l=%d, m=%d, u=%s\n",l,m,argv[j]);*/
   }
   vlist->vtx[i]->data->z=fabs(r)*cos(theta);
   vlist->vtx[i]->data->x=fabs(r)*sin(theta)*cos(fi);
   vlist->vtx[i]->data->y=fabs(r)*sin(theta)*sin(fi);
   vlist->vtx[i]->z=fabs(r)*cos(theta);
   vlist->vtx[i]->x=fabs(r)*sin(theta)*cos(fi);
   vlist->vtx[i]->y=fabs(r)*sin(theta)*sin(fi);
}
write_vertex_xml_file(vesicle,0);