From 0661f62a1a4b5bc87310421d198e7df8538c7363 Mon Sep 17 00:00:00 2001 From: Samo Penic <samo.penic@gmail.com> Date: Tue, 26 Jun 2018 10:33:59 +0000 Subject: [PATCH] Added function to calculate absolute value of ulm2 --- src/tape | 4 ++-- src/shcomplex.c | 6 ++++++ src/shcomplex.h | 1 + 3 files changed, 9 insertions(+), 2 deletions(-) diff --git a/src/shcomplex.c b/src/shcomplex.c index bd3a733..9002e7e 100644 --- a/src/shcomplex.c +++ b/src/shcomplex.c @@ -171,3 +171,9 @@ return retval; } + + +ts_double get_absolute_ulm2(ts_vesicle *vesicle, int l, int m){ + return gsl_complex_abs2(vesicle->sphHarmonics->ulmComplex[l][m]); +} + diff --git a/src/shcomplex.h b/src/shcomplex.h index 63c6af5..d298228 100644 --- a/src/shcomplex.h +++ b/src/shcomplex.h @@ -7,4 +7,5 @@ ts_bool complex_sph_free(ts_spharm *sph); ts_bool calculateUlmComplex(ts_vesicle *vesicle); ts_double calculateKc(ts_vesicle *vesicle, ts_int lmin, ts_int lmax); +ts_double get_absolute_ulm2(ts_vesicle *vesicle, int l, int m); #endif diff --git a/src/tape b/src/tape index bc70c2f..978e3c9 100644 --- a/src/tape +++ b/src/tape @@ -25,8 +25,8 @@ #Stretching -stretchswitch=0 -xkA0=100.0 +stretchswitch=1 +xkA0=1.0 ####### Polymer (brush) definitions ########### # npoly is a number of polymers attached to npoly distinct vertices on vesicle -- Gitblit v1.9.3