Trisurf Monte Carlo simulator
Samo Penic
2014-05-06 7d5e22d9eda3020e8b2b3fd8a70f25306561b401
src/initial_distribution.c
@@ -12,6 +12,7 @@
#include "poly.h"
#include "io.h"
#include "sh.h"
#include "shcomplex.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,6 +42,7 @@
   ts_vertex *vtx;
   vesicle=initial_distribution_dipyramid(tape->nshell,tape->ncxmax,tape->ncymax,tape->nczmax,tape->stepsize);
    vesicle->tape=tape;
   // Nucleus:
   vesicle->R_nucleus=tape->R_nucleus*tape->R_nucleus;
@@ -93,7 +95,12 @@
   vesicle->pressure= tape->pressure;
   vesicle->pswitch=tape->pswitch;
   vesicle->sphHarmonics=sph_init(vesicle->vlist,tape->shc);
    if(tape->shc>0){
       vesicle->sphHarmonics=complex_sph_init(vesicle->vlist,tape->shc);
    }
    else {
        vesicle->sphHarmonics=NULL;
    }
    return vesicle;
}