Adding ending energy to filaments. In progress. Continue in single_filament...

mihaf

2014-03-18 58230a2591414fb38b9ec8d3a76439b290cb0a6f

 src/initial_distribution.c

@@ -10,6 +10,7 @@
#include "initial_distribution.h"
#include "energy.h"
#include "poly.h"
#include "io.h"

ts_vesicle *initial_distribution_dipyramid(ts_uint nshell, ts_uint ncmax1, ts_uint ncmax2, ts_uint ncmax3, ts_double stepsize){
   ts_fprintf(stdout,"Starting initial_distribution on vesicle with %u shells!...\n",nshell);
@@ -37,9 +38,24 @@
ts_vesicle *create_vesicle_from_tape(ts_tape *tape){
   ts_vesicle *vesicle;
   vesicle=initial_distribution_dipyramid(tape->nshell,tape->ncxmax,tape->ncymax,tape->nczmax,tape->stepsize);
   vesicle->poly_list=init_poly_list(tape->npoly,tape->nmono, vesicle->vlist);
   // Nucleus:
   vesicle->R_nucleus=tape->R_nucleus;
   //Initialize grafted polymers (brush):
   vesicle->poly_list=init_poly_list(tape->npoly,tape->nmono, vesicle->vlist, vesicle);
   vesicle->spring_constant=tape->kspring;
   poly_assign_spring_const(vesicle);
   //Initialize filaments (polymers inside the vesicle):
   vesicle->filament_list=init_poly_list(tape->nfil,tape->nfono, NULL, vesicle);
ts_uint i,j;
   for(i=0;i<vesicle->filament_list->n;i++){
         for(j=0;j<vesicle->filament_list->poly[i]->vlist->n;j++){

   fprintf(stderr,"%e %e %e\n", vesicle->filament_list->poly[i]->vlist->vtx[j]->x,vesicle->filament_list->poly[i]->vlist->vtx[j]->y, vesicle->filament_list->poly[i]->vlist->vtx[j]->z );
   }
   }
//   vesicle->spring_constant=tape->kspring;
//   poly_assign_spring_const(vesicle);

   
   vesicle->nshell=tape->nshell;
   vesicle->dmax=tape->dmax*tape->dmax; /* dmax^2 in the vesicle dmax variable */