From snapshots, raw tape configuration is extracted and written into tape f...

Samo Penic

2016-05-24 cc95bd1ca5e69ce13766ca0acea11560531ed620

 src/initial_distribution.c

@@ -1,3 +1,4 @@
/* vim: set ts=4 sts=4 sw=4 noet : */
#include<stdlib.h>
#include<math.h>
#include<stdio.h>
@@ -51,7 +52,9 @@
   ts_vertex *vtx;
   ts_tape *tape=vesicle->tape;
   vesicle->R_nucleus=tape->R_nucleus*tape->R_nucleus;

   vesicle->R_nucleusX=tape->R_nucleusX*tape->R_nucleusX;
   vesicle->R_nucleusY=tape->R_nucleusY*tape->R_nucleusY;
   vesicle->R_nucleusZ=tape->R_nucleusZ*tape->R_nucleusZ;
   vesicle->clist->dmin_interspecies = tape->dmin_interspecies*tape->dmin_interspecies;

   //Initialize grafted polymers (brush):
@@ -92,7 +95,7 @@
   vesicle->dmax=tape->dmax*tape->dmax; /* dmax^2 in the vesicle dmax variable */
   vesicle->bending_rigidity=tape->xk0;
   vtx_set_global_values(vesicle); /* make xk0 default value for every vertex */ 
   ts_fprintf(stdout, "Tape setting: xk0=%e\n",tape->xk0);
//   ts_fprintf(stdout, "Tape setting: xk0=%e\n",tape->xk0);
   vesicle->stepsize=tape->stepsize;
   vesicle->clist->ncmax[0]=tape->ncxmax;
   vesicle->clist->ncmax[1]=tape->ncymax;