Trisurf Monte Carlo simulator
Samo Penic
2018-05-15 00251384c9c40368b04cfab4f4082666b5541a80
Fixes in dump ;)
1 files modified
10 ■■■■ changed files
src/restore.c 10 ●●●● patch | view | raw | blame | history
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)/2*pow((vesicle->tape->dmax+1.0),2);
        for(i=0;i<vesicle->tlist->n;i++){
            stretchenergy(vesicle, vesicle->tlist->tria[i]);
        }
    }
/* TODO: filaments */
//    ts_fprintf(stdout,"Restoration completed\n");