From 9166cbcd0e28d61a69646911af35bb7895ff9203 Mon Sep 17 00:00:00 2001
From: Samo Penic <samo.penic@fe.uni-lj.si>
Date: Mon, 28 Apr 2014 11:16:45 +0000
Subject: [PATCH] Did major change in vesicle datastructure. Added tape information into vesicle. This means the dump files will be broken each time tape structure changes. Beware.

---
 src/initial_distribution.c |   10 +++++++++-
 1 files changed, 9 insertions(+), 1 deletions(-)

diff --git a/src/initial_distribution.c b/src/initial_distribution.c
index 0587d18..35d533d 100644
--- a/src/initial_distribution.c
+++ b/src/initial_distribution.c
@@ -11,6 +11,8 @@
 #include "energy.h"
 #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);
@@ -40,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;
 
@@ -92,7 +95,12 @@
 
 	vesicle->pressure= tape->pressure;
 	vesicle->pswitch=tape->pswitch;
-
+    if(tape->shc>0){
+	    vesicle->sphHarmonics=complex_sph_init(vesicle->vlist,tape->shc);
+    }
+    else {
+        vesicle->sphHarmonics=NULL;
+    }
     return vesicle;
 
 }

--
Gitblit v1.9.3