| | |
|---|
| | | #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); |
|---|
| | | } |
|---|
| | |
|---|
| | | 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"); |
|---|
