From 3197854e89eb9ffe5ff3137e274d16806c2028a5 Mon Sep 17 00:00:00 2001
From: Samo Penic <samo.penic@gmail.com>
Date: Tue, 27 Jul 2021 12:51:09 +0000
Subject: [PATCH] Found one bug in the code when calculating mprod. Fixed.
---
src/shcomplex.c | 7 +++++--
1 files changed, 5 insertions(+), 2 deletions(-)
diff --git a/src/shcomplex.c b/src/shcomplex.c
index 8492609..4090346 100644
--- a/src/shcomplex.c
+++ b/src/shcomplex.c
@@ -1,6 +1,7 @@
/* vim: set ts=4 sts=4 sw=4 noet : */
#include<math.h>
#include<stdlib.h>
+#include<string.h>
#include<gsl/gsl_complex.h>
#include<gsl/gsl_complex_math.h>
#include<gsl/gsl_sf_legendre.h>
@@ -120,12 +121,14 @@
char *Ulm2Complex2String(ts_vesicle *vesicle){
ts_int i,j;
char *strng=(char *)calloc(5000, sizeof(char));
+ char tmpstrng[255];
for(i=0;i<vesicle->sphHarmonics->l;i++){
for(j=i;j<2*i+1;j++){
- sprintf(strng,"%e ", gsl_complex_abs2(vesicle->sphHarmonics->ulmComplex[i][j]));
+ sprintf(tmpstrng,"%e ", gsl_complex_abs2(vesicle->sphHarmonics->ulmComplex[i][j]));
+ strcat(strng,tmpstrng);
}
}
- sprintf(strng,"\n");
+ //strcat(strng,"\n");
return strng;
}
--
Gitblit v1.9.3