From 0af0cf6ac748ba9c135c9ca0a23cdd6cbf50caf1 Mon Sep 17 00:00:00 2001 From: Samo Penic <samo.penic@gmail.com> Date: Fri, 13 Jun 2014 10:27:19 +0000 Subject: [PATCH] First attempt to calculate points of cross-section --- src/general.h | 98 ++++++++++++++++++++++++++++++++++++++++++------ 1 files changed, 85 insertions(+), 13 deletions(-) diff --git a/src/general.h b/src/general.h index 063e124..ef57b1e 100644 --- a/src/general.h +++ b/src/general.h @@ -3,7 +3,7 @@ #include<stdarg.h> #include<stdio.h> - +#include<gsl/gsl_complex.h> /* @brief This is a header file, defining general constants and structures. * @file header.h * @author Samo Penic @@ -48,6 +48,8 @@ #define TS_FAIL 1 /* CONSTANTS */ + +#define TS_ID_FILAMENT 1 /* DATA TYPES */ /** @brief Sets the default datatype for ts_double @@ -127,6 +129,10 @@ } ts_coord; +typedef struct { + ts_double N; + ts_coord *coord; +} ts_coord_list; /** @brief Data structure of all data connected to a vertex * @@ -155,6 +161,7 @@ ts_double projArea; ts_double relR; ts_double solAngle; + struct ts_poly *grafted_poly; }; typedef struct ts_vertex ts_vertex; @@ -165,12 +172,14 @@ } ts_vertex_list; struct ts_bond { - ts_uint idx; + ts_uint idx; ts_vertex *vtx1; ts_vertex *vtx2; - ts_double bond_length; - ts_double bond_length_dual; - ts_bool tainted; + ts_double bond_length; + ts_double bond_length_dual; + ts_bool tainted; //TODO: remove + ts_double energy; + ts_double x,y,z; }; typedef struct ts_bond ts_bond; @@ -213,12 +222,14 @@ ts_double dcell; ts_double shift; ts_double max_occupancy; + ts_double dmin_interspecies; } ts_cell_list; typedef struct { ts_uint l; ts_double **ulm; + gsl_complex **ulmComplex; ts_double **sumUlm2; ts_uint N; ts_double **co; @@ -227,20 +238,81 @@ +struct ts_poly { + ts_vertex_list *vlist; + ts_bond_list *blist; + ts_vertex *grafted_vtx; + ts_double k; +}; +typedef struct ts_poly ts_poly; + + +struct ts_poly_list { + ts_uint n; + ts_poly **poly; +}; +typedef struct ts_poly_list ts_poly_list; + + + +typedef struct { + long int nshell; + long int ncxmax; + long int ncymax; + long int nczmax; + long int npoly; + long int nmono; + long int nfil; + long int nfono; + long int R_nucleus; + long int pswitch; + long int constvolswitch; + ts_double constvolprecision; + char *multiprocessing; + long int brezveze0; + long int brezveze1; + long int brezveze2; + ts_double xk0; + ts_double dmax; + ts_double dmin_interspecies; + ts_double stepsize; + ts_double kspring; + ts_double xi; + ts_double pressure; + long int iterations; + long int inititer; + long int mcsweeps; + long int quiet; + long int shc; +} ts_tape; + + + + typedef struct { ts_vertex_list *vlist; ts_bond_list *blist; ts_triangle_list *tlist; - ts_cell_list *clist; + ts_cell_list *clist; ts_uint nshell; - ts_double bending_rigidity; - ts_double dmax; - ts_double stepsize; - ts_double cm[3]; - ts_double volume; - ts_spharm *sphHarmonics; -} ts_vesicle; + ts_double bending_rigidity; + ts_double dmax; + ts_double stepsize; + ts_double cm[3]; + ts_double volume; + ts_spharm *sphHarmonics; +// Polymers outside the vesicle and attached to the vesicle membrane (polymer brush): + ts_poly_list *poly_list; +// Filaments inside the vesicle (not attached to the vesicel membrane: + ts_poly_list *filament_list; + ts_double spring_constant; + ts_double pressure; + ts_int pswitch; + ts_tape *tape; + ts_double R_nucleus; + +} ts_vesicle; -- Gitblit v1.9.3