Code for incremental adding of inclusions done

Samo Penic

2020-04-26 3006183b769f2e126b1a96e6bf697c2b7f657df7

 src/snapshot.c

@@ -23,8 +23,8 @@
/* outputs additional data into paraview xml file */
ts_bool xml_trisurf_data(FILE *fh, ts_vesicle *vesicle){

   ts_string *data=(ts_string *)malloc(sizeof(ts_sprintf));
   data->string=(char *)malloc(512000*sizeof(char)); /*TODO: warning, can break if the string is to long */
   ts_string *data=(ts_string *)malloc(sizeof(ts_string));
   data->string=(char *)malloc(5120000*sizeof(char)); /*TODO: warning, can break if the string is to long */
   data->beg=0;
   
   xml_trisurf_header(fh, vesicle);
@@ -32,6 +32,8 @@
   xml_trisurf_tria_neigh(data,vesicle->tlist);
   xml_trisurf_vtx_neigh(data,vesicle->vlist);   
   xml_trisurf_vtx_tristar(data,vesicle->vlist);
   xml_trisurf_nucleus(data,vesicle);
   xml_trisurf_constvolarea(data,V0,A0);
#ifdef COMPRESSION
   char *compressed;
   ts_uint nbytes=ts_compress_string64(data->string, data->beg-1, &compressed); //suppress null character at the end with by substracting 1
@@ -57,7 +59,7 @@
   fprintf(fh, "<trisurfversion>Trisurf (commit %s), compiled on %s %s</trisurfversion>\n",TS_VERSION, __DATE__,  __TIME__);
   fprintf(fh, "<dumpdate>%s</dumpdate>\n", c_time_string);

   fprintf(fh, "<tape>\n");
   fprintf(fh, "<tape>");
      fprintf(fh,"%s",tapetxt);   
   fprintf(fh, "</tape>\n");
   if(vesicle->poly_list!=NULL){
@@ -124,7 +126,18 @@
   return TS_SUCCESS;
}

ts_bool xml_trisurf_nucleus(ts_string *data, ts_vesicle* vesicle){
   if(vesicle->R_nucleus>0.0 || (vesicle->R_nucleusX>0.0 && vesicle->R_nucleusY>0.0 && vesicle->R_nucleusZ>0.0)){
      ts_sprintf(data,"<nucleus>%.17e %.17e %.17e</nucleus>",vesicle->nucleus_center[0], vesicle->nucleus_center[1], vesicle->nucleus_center[2]);
   }
   return TS_SUCCESS;
}
ts_bool xml_trisurf_constvolarea(ts_string *data, ts_double volume, ts_double area){
   ts_sprintf(data,"<constant_volume>%.17e</constant_volume>",volume);
   ts_sprintf(data,"<constant_area>%.17e</constant_area>",area);

   return TS_SUCCESS;
}


/* UTILITIES */