From 0634b97e436e480b203a22e1113c473d3ae4d281 Mon Sep 17 00:00:00 2001
From: Samo Penic <samo.penic@fe.uni-lj.si>
Date: Tue, 06 Dec 2016 13:35:11 +0000
Subject: [PATCH] Merge branch 'nirgov' into HEAD

---
 src/restore.c |   49 +++++++++++++++++++++++++++++++++++++++++++++++--
 1 files changed, 47 insertions(+), 2 deletions(-)

diff --git a/src/restore.c b/src/restore.c
index 086bf1c..c89ff6d 100644
--- a/src/restore.c
+++ b/src/restore.c
@@ -54,6 +54,10 @@
 				if ((!xmlStrcmp(cur1->name, (const xmlChar *)"Piece"))){
 					cur2=cur1->xmlChildrenNode;
 					while(cur2!=NULL){
+						if ((!xmlStrcmp(cur2->name, (const xmlChar *)"PointData"))){
+							if(vesicle!=NULL)
+								parseXMLPointData(vesicle,doc,cur2);
+						}
 						if ((!xmlStrcmp(cur2->name, (const xmlChar *)"Points"))){
 							//fprintf(stderr,"Found point data\n");
 							if(vesicle!=NULL)
@@ -81,6 +85,7 @@
 
 	init_normal_vectors(vesicle->tlist);
 	mean_curvature_and_energy(vesicle);
+	sweep_attraction_bond_energy(vesicle);
 
 /* TODO: filaments */
 
@@ -182,7 +187,6 @@
 
 	vesicle->tape=tape;
 	set_vesicle_values_from_tape(vesicle);
-
 	return vesicle;
 }
 
@@ -318,7 +322,48 @@
 	return TS_SUCCESS;
 }
 
-
+/* this parses the data for vertices (like spontaneous curvature, etc.) */
+ts_bool parseXMLPointData(ts_vesicle *vesicle,xmlDocPtr doc, xmlNodePtr cur){
+	xmlNodePtr child = cur->xmlChildrenNode;
+	xmlChar *property_name;
+	xmlChar *values;
+	char *vals;
+	char *token;
+	int idx, polyidx, monoidx, filidx, fonoidx;
+	while (child != NULL) {
+		if ((!xmlStrcmp(child->name, (const xmlChar *)"DataArray"))){
+			property_name=xmlGetProp(child, (xmlChar *)"Name");
+		//	fprintf(stderr,"Name: %s\n", property_name);
+			if(!xmlStrcmp(property_name,(const xmlChar *)"spontaneous_curvature")){
+				values=xmlNodeListGetString(doc,child->xmlChildrenNode,1);
+				vals=(char *)values;
+				token=strtok(vals," ");
+				idx=0;
+				while(token!=NULL){
+					if(idx<vesicle->vlist->n){
+						vesicle->vlist->vtx[idx]->c=atof(token);
+					} 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]->c=atof(token);
+					} 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]->c=atof(token);
+					}
+					idx++;
+					token=strtok(NULL," ");
+				}
+				xmlFree(values);		
+			}
+			xmlFree(property_name);
+		}
+		
+		child=child->next;
+	}
+	return TS_SUCCESS;
+}
 /* this is a parser of vertex positions and bonds from main xml data */
 ts_bool parseXMLVertexPosition(ts_vesicle *vesicle,xmlDocPtr doc, xmlNodePtr cur){
 	xmlNodePtr child = cur->xmlChildrenNode;

--
Gitblit v1.9.3