From 4a1b277c5b9ad5369441d028e53f108adc8ef62f Mon Sep 17 00:00:00 2001
From: mihaf <miha.fosnaric@gmail.com>
Date: Fri, 21 Feb 2014 08:36:20 +0000
Subject: [PATCH] Nothing special.

---
 src/frame.c |   45 +++++++++++++++++++++++++++++++++------------
 1 files changed, 33 insertions(+), 12 deletions(-)

diff --git a/src/frame.c b/src/frame.c
index a4a9d32..1ee270b 100644
--- a/src/frame.c
+++ b/src/frame.c
@@ -9,34 +9,55 @@
     vesicle->cm[1]=0;
     vesicle->cm[2]=0;
     for(i=0;i<n;i++){
-        vesicle->cm[0]+=vtx[i]->data->x;
-        vesicle->cm[1]+=vtx[i]->data->y;
-        vesicle->cm[2]+=vtx[i]->data->z; 
+        vesicle->cm[0]+=vtx[i]->x;
+        vesicle->cm[1]+=vtx[i]->y;
+        vesicle->cm[2]+=vtx[i]->z; 
     } 
     vesicle->cm[0]/=(ts_float)n;
     vesicle->cm[1]/=(ts_float)n;
     vesicle->cm[2]/=(ts_float)n;
+
+    for(i=0;i<n;i++){
+        vtx[i]->x-=vesicle->cm[0];
+        vtx[i]->y-=vesicle->cm[1];
+        vtx[i]->z-=vesicle->cm[2]; 
+    } 
+
+    vesicle->cm[0]=0;
+    vesicle->cm[1]=0;
+    vesicle->cm[2]=0;
+
     return TS_SUCCESS;
 }
 
 ts_bool cell_occupation(ts_vesicle *vesicle){
-    ts_uint i,cellidx, n=vesicle->vlist->n;
-    ts_double shift;
-    ts_double dcell;
-    shift=(ts_double) vesicle->clist->ncmax[0]/2;
-    dcell=1.0/(1.0 + vesicle->stepsize);
+    ts_uint i,j,cellidx, n=vesicle->vlist->n;
+    //ts_double shift;
+    //ts_double dcell;
+    //shift=(ts_double) vesicle->clist->ncmax[0]/2;
+    //dcell=1.0/(1.0 + vesicle->stepsize);
     //`fprintf(stderr, "Bil sem tu\n"); 
 
     cell_list_cell_occupation_clear(vesicle->clist);
     for(i=0;i<n;i++){
     cellidx=vertex_self_avoidance(vesicle, vesicle->vlist->vtx[i]);
-    vesicle->vlist->vtx[i]->data->cell=vesicle->clist->cell[cellidx];
+//	already done in cell_add_vertex
+//    vesicle->vlist->vtx[i]->cell=vesicle->clist->cell[cellidx];
 
     cell_add_vertex(vesicle->clist->cell[cellidx],vesicle->vlist->vtx[i]);
-
     }
+
+    for(i=0;i<vesicle->poly_list->n;i++){
+	for(j=0;j<vesicle->poly_list->poly[i]->vlist->n;j++){
+    	cellidx=vertex_self_avoidance(vesicle, vesicle->poly_list->poly[i]->vlist->vtx[j]);
+    	cell_add_vertex(vesicle->clist->cell[cellidx],vesicle->poly_list->poly[i]->vlist->vtx[j]);
+	}
+    }
+
+    
+
     //fprintf(stderr, "Bil sem tu\n"); 
-	if(dcell);
-	if(shift);
+	//if(dcell);
+	//if(shift);
     return TS_SUCCESS;
 }

--
Gitblit v1.9.3