Trisurf Monte Carlo simulator
Samo Penic
2019-02-08 4019ca8157136a68b6e43b31f166dc93acd81b4d
src/restore.c
@@ -16,13 +16,14 @@
#include "poly.h"
#include "initial_distribution.h"
#include "io.h"
#include <math.h>
ts_vesicle *parseDump(char *dumpfname) {
   xmlDocPtr doc;
   xmlNodePtr cur, cur1,cur2;
   ts_vesicle *vesicle=NULL;
   doc = xmlParseFile(dumpfname);
   int i;
   if (doc == NULL ) {
      fatal("Dump file could not be found or parsed. It is correct file?",1);
   }
@@ -86,7 +87,12 @@
   init_normal_vectors(vesicle->tlist);
   mean_curvature_and_energy(vesicle);
   sweep_attraction_bond_energy(vesicle);
   if(vesicle->tape->stretchswitch==1){
      vesicle->tlist->a0=sqrt(3)/4.0*pow((vesicle->tape->dmax+1.0)/2.0,2);
      for(i=0;i<vesicle->tlist->n;i++){
         stretchenergy(vesicle, vesicle->tlist->tria[i]);
      }
   }
/* TODO: filaments */
//   ts_fprintf(stdout,"Restoration completed\n");