From 86b69b3803b1b31c040cede6b47bf617561b3b81 Mon Sep 17 00:00:00 2001
From: Samo Penic <samo.penic@fe.uni-lj.si>
Date: Tue, 01 Mar 2016 08:36:37 +0000
Subject: [PATCH] Polymeres are restored
---
src/shcomplex.c | 10 +++++-----
1 files changed, 5 insertions(+), 5 deletions(-)
diff --git a/src/shcomplex.c b/src/shcomplex.c
index 7afe27e..7d8c9d8 100644
--- a/src/shcomplex.c
+++ b/src/shcomplex.c
@@ -130,9 +130,9 @@
}
-ts_double calculateKc(ts_vesicle *vesicle){
- ts_int min=4;
- ts_int max=vesicle->sphHarmonics->l-6;
+ts_double calculateKc(ts_vesicle *vesicle, ts_int lmin, ts_int lmax){
+ ts_int min=lmin;
+ ts_int max=lmax; //vesicle->sphHarmonics->l-3;
ts_long i,j;
ts_double retval, bval;
gsl_matrix *A=gsl_matrix_alloc(max-min,2);
@@ -146,7 +146,7 @@
for(i=min;i<max;i++){
gsl_matrix_set(A, i-min,0,(ts_double)((i-1)*(i+2)));
gsl_matrix_set(A, i-min,1,(ts_double)((i-1)*(i+2)*(i+1)*i));
- fprintf(stderr,"%e %e\n", gsl_matrix_get(A,i-min,0), gsl_matrix_get(A,i-min,1));
+// fprintf(stderr,"%e %e\n", gsl_matrix_get(A,i-min,0), gsl_matrix_get(A,i-min,1));
bval=0.0;
//average for m from 0..l (only positive m's)
for(j=0;j<=i;j++){
@@ -155,7 +155,7 @@
bval=bval/(ts_double)vesicle->sphHarmonics->N/(ts_double)(i+1);
gsl_vector_set(b,i-min,1.0/bval);
- fprintf(stderr,"%e\n", 1.0/gsl_vector_get(b,i-min));
+// fprintf(stderr,"%e\n", 1.0/gsl_vector_get(b,i-min));
}
// fprintf(stderr,"b[2]=%e\n",gsl_vector_get(b,1));
gsl_linalg_QR_decomp(A,tau);
--
Gitblit v1.9.3