From f977e8fe32b6a1a2ae85662ab35fd7c9f647cc3e Mon Sep 17 00:00:00 2001 From: Samo Penic <samo.penic@gmail.com> Date: Mon, 02 Jul 2018 10:50:58 +0000 Subject: [PATCH] Added a funciton that calculates everything needed for spherical harmonics --- src/general.h | 88 +++++++++++++++++++++++++++++++++++++++++-- 1 files changed, 83 insertions(+), 5 deletions(-) diff --git a/src/general.h b/src/general.h index e47ef1c..7676000 100644 --- a/src/general.h +++ b/src/general.h @@ -1,8 +1,10 @@ +/* vim: set ts=4 sts=4 sw=4 noet : */ #ifndef _GENERAL_H #define _GENERAL_H #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 @@ -154,7 +156,8 @@ ts_double projArea; ts_double relR; ts_double solAngle; - struct ts_poly *grafted_poly; + struct ts_poly *grafted_poly; + struct ts_cluster *cluster; }; typedef struct ts_vertex ts_vertex; @@ -192,11 +195,13 @@ ts_double znorm; ts_double area; // firstly needed for sh.c ts_double volume; // firstly needed for sh.c + ts_double energy; }; typedef struct ts_triangle ts_triangle; struct ts_triangle_list{ ts_uint n; + ts_double a0; ts_triangle **tria; }; typedef struct ts_triangle_list ts_triangle_list; @@ -222,6 +227,7 @@ typedef struct { ts_uint l; ts_double **ulm; + gsl_complex **ulmComplex; ts_double **sumUlm2; ts_uint N; ts_double **co; @@ -247,6 +253,50 @@ +typedef struct { + long int nshell; + long int ncxmax; + long int ncymax; + long int nczmax; + long int npoly; + long int nmono; + long int internal_poly; + long int nfil; + long int nfono; + long int R_nucleus; + ts_double R_nucleusX; + ts_double R_nucleusY; + ts_double R_nucleusZ; + long int pswitch; + long int constvolswitch; + long int constareaswitch; + long int stretchswitch; + ts_double xkA0; + 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; + long int number_of_vertices_with_c0; + ts_double c0; + ts_double w; + ts_double F; +} ts_tape; + + + typedef struct { ts_vertex_list *vlist; @@ -268,18 +318,38 @@ ts_double spring_constant; ts_double pressure; ts_int pswitch; - + ts_tape *tape; ts_double R_nucleus; - + ts_double R_nucleusX; + ts_double R_nucleusY; + ts_double R_nucleusZ; + ts_double nucleus_center[3]; + ts_double area; } ts_vesicle; + + +struct ts_cluster{ + ts_uint nvtx; + ts_uint idx; + ts_vertex **vtx; +}; + +typedef struct ts_cluster ts_cluster; + +typedef struct{ + ts_uint n; + ts_cluster **cluster; +} ts_cluster_list; /* GLOBAL VARIABLES */ int quiet; - - +ts_double V0; +ts_double A0; +ts_double epsvol; +ts_double epsarea; /* FUNCTIONS */ /** Non-fatal error function handler: @@ -300,4 +370,12 @@ #define VTX(n) &(vlist->vtx[n]) #define VTX_DATA(n) vlist->vtx[n].data + +/* FOR PID GENERATION ROUTINE */ +#define CPF_CLOEXEC 1 + +int createPidFile(const char *progName, const char *pidFile, int flags); + +int lockRegion(int fd, int type, int whence, int start, int len); +char *libVersion(); #endif -- Gitblit v1.9.3