Added spherical harmonics as-is into rbc

Samo Penic

2014-04-09 dc77e8ef10378188ee0a9657cedd1b7991a0a29f

 src/initial_distribution.c

@@ -11,6 +11,7 @@
#include "energy.h"
#include "poly.h"
#include "io.h"
#include "sh.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);
@@ -41,7 +42,9 @@

   vesicle=initial_distribution_dipyramid(tape->nshell,tape->ncxmax,tape->ncymax,tape->nczmax,tape->stepsize);
   // Nucleus:
   vesicle->R_nucleus=tape->R_nucleus;
   vesicle->R_nucleus=tape->R_nucleus*tape->R_nucleus;

   vesicle->clist->dmin_interspecies = tape->dmin_interspecies*tape->dmin_interspecies;

   //Initialize grafted polymers (brush):
   vesicle->poly_list=init_poly_list(tape->npoly,tape->nmono, vesicle->vlist, vesicle);
@@ -69,6 +72,9 @@
      }
   }

   for(i=0;i<vesicle->filament_list->n;i++){
      vertex_list_assign_id(vesicle->filament_list->poly[i]->vlist,TS_ID_FILAMENT);
   }

//   vesicle->spring_constant=tape->kspring;
//   poly_assign_spring_const(vesicle);
@@ -87,7 +93,7 @@

   vesicle->pressure= tape->pressure;
   vesicle->pswitch=tape->pswitch;

   vesicle->sphHarmonics=sph_init(vesicle->vlist,tape->shc);
    return vesicle;

}