From aafd67f2d592799167864698ea07eb8aadcc86a7 Mon Sep 17 00:00:00 2001
From: Samo Penic <samo.penic@gmail.com>
Date: Wed, 25 May 2016 18:37:21 +0000
Subject: [PATCH] Merge branch 'master' of bitbucket.org:samop/trisurf-ng with work done on removing the memory leaks.

---
 src/restore.c |   19 +++++++++++++++----
 1 files changed, 15 insertions(+), 4 deletions(-)

diff --git a/src/restore.c b/src/restore.c
index 2cc83bb..14586c6 100644
--- a/src/restore.c
+++ b/src/restore.c
@@ -1,3 +1,4 @@
+/* vim: set ts=4 sts=4 sw=4 noet : */
 #include <stdio.h>
 #include <string.h>
 #include <stdlib.h>
@@ -56,6 +57,7 @@
 						if ((!xmlStrcmp(cur2->name, (const xmlChar *)"Points"))){
 							//fprintf(stderr,"Found point data\n");
 							if(vesicle!=NULL)
+								//fprintf(stderr,"Fils: %ld, Nfono: %ld\n", vesicle->tape->nfil, vesicle->tape->nfono);
 								parseXMLVertexPosition(vesicle, doc, cur2);
 						}
 						if ((!xmlStrcmp(cur2->name, (const xmlChar *)"Cells"))){
@@ -82,7 +84,7 @@
 
 /* TODO: filaments */
 
-	fprintf(stderr,"Restoration completed\n");
+//	ts_fprintf(stdout,"Restoration completed\n");
 //	write_vertex_xml_file(vesicle,999);
 //	vesicle_free(vesicle);
 //	exit(0);
@@ -285,7 +287,7 @@
 	xmlNodePtr child = cur->xmlChildrenNode;
 	xmlChar *points;
 	char *pts;
-	int i, idx, polyidx, monoidx;
+	int i, idx, polyidx, monoidx, filidx, fonoidx;
 	char *token[3];
 	while (child != NULL) {
 		if ((!xmlStrcmp(child->name, (const xmlChar *)"DataArray"))){
@@ -300,12 +302,19 @@
 					vesicle->vlist->vtx[idx]->x=atof(token[0]);
 					vesicle->vlist->vtx[idx]->y=atof(token[1]);
 					vesicle->vlist->vtx[idx]->z=atof(token[2]);
-				} else {
+				} else if(vesicle->tape->nmono && vesicle->tape->npoly && idx<vesicle->vlist->n+vesicle->tape->nmono*vesicle->tape->npoly) {
 					polyidx=(idx-vesicle->vlist->n)/vesicle->tape->nmono;
 					monoidx=(idx-vesicle->vlist->n)%vesicle->tape->nmono;
 					vesicle->poly_list->poly[polyidx]->vlist->vtx[monoidx]->x=atof(token[0]);
 					vesicle->poly_list->poly[polyidx]->vlist->vtx[monoidx]->y=atof(token[1]);
 					vesicle->poly_list->poly[polyidx]->vlist->vtx[monoidx]->z=atof(token[2]);
+				} else {
+					filidx=(idx-vesicle->vlist->n-vesicle->tape->nmono*vesicle->tape->npoly)/vesicle->tape->nfono;
+					fonoidx=(idx-vesicle->vlist->n-vesicle->tape->nmono*vesicle->tape->npoly)%vesicle->tape->nfono;
+					//fprintf(stderr,"filidx=%d, fonoidx=%d, coord=%s,%s,%s\n",filidx,fonoidx,token[0],token[1],token[2]);
+					vesicle->filament_list->poly[filidx]->vlist->vtx[fonoidx]->x=atof(token[0]);
+					vesicle->filament_list->poly[filidx]->vlist->vtx[fonoidx]->y=atof(token[1]);
+					vesicle->filament_list->poly[filidx]->vlist->vtx[fonoidx]->z=atof(token[2]);
 				}
 				for(i=0;i<2;i++)	token[i]=strtok(NULL," ");	
 				token[2]=strtok(NULL,"\n");
@@ -315,6 +324,7 @@
 		}
 		child=child->next;
 	}
+	//fprintf(stderr,"Came here\n");
 	//fprintf(stderr,"Vertices position j=%d\n",idx);	
 
 	return TS_SUCCESS;
@@ -339,7 +349,8 @@
 				}
 				else {
 				//find grafted vtx
-					if((vesicle->tape->nmono-1)==(idx-3*(vesicle->vlist->n-2))%(vesicle->tape->nmono)){
+					if(vesicle->tape->npoly && vesicle->tape->nmono && (vesicle->tape->nmono-1)==(idx-3*(vesicle->vlist->n-2))%(vesicle->tape->nmono)
+						&& idx<(3*vesicle->vlist->n-2+vesicle->tape->nmono*vesicle->tape->npoly+vesicle->tape->npoly)){
 						polyidx=(idx-3*(vesicle->vlist->n-2))/(vesicle->tape->nmono);
 						//fprintf(stderr,"poly=%d, vertex=%d\n",polyidx,atoi(token[0]));
 						vesicle->poly_list->poly[polyidx]->grafted_vtx=vesicle->vlist->vtx[atoi(token[0])];

--
Gitblit v1.9.3