From dba66b4f8a25cb7eb18de5f7ffca97e179d7a5d9 Mon Sep 17 00:00:00 2001
From: Samo Penic <samo.penic@gmail.com>
Date: Sat, 19 Oct 2019 11:45:15 +0000
Subject: [PATCH] Moved compression/decompression routines to separate file
---
src/shcomplex.c | 24 +++++++++++++++++++++++-
1 files changed, 23 insertions(+), 1 deletions(-)
diff --git a/src/shcomplex.c b/src/shcomplex.c
index 5dc6a97..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>
@@ -69,7 +70,7 @@
if(sph->co != NULL) free(sph->co);
if(sph->ulm !=NULL) free(sph->ulm);
if(sph->ulmComplex !=NULL) free(sph->ulmComplex);
-
+ if(sph->sumUlm2 !=NULL) free(sph->sumUlm2);
if(sph->Ylmi!=NULL) {
for(i=0;i<sph->l;i++){
if(sph->Ylmi[i]!=NULL){
@@ -117,6 +118,27 @@
return TS_SUCCESS;
}
+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(tmpstrng,"%e ", gsl_complex_abs2(vesicle->sphHarmonics->ulmComplex[i][j]));
+ strcat(strng,tmpstrng);
+ }
+ }
+ //strcat(strng,"\n");
+
+ return strng;
+}
+
+ts_bool freeUlm2String(char *strng){
+ free(strng);
+ return TS_SUCCESS;
+}
+
+
ts_bool storeUlmComplex2(ts_vesicle *vesicle){
ts_spharm *sph=vesicle->sphHarmonics;
--
Gitblit v1.9.3