Added a funciton that calculates everything needed for spherical harmonics

Samo Penic

2018-07-02 f977e8fe32b6a1a2ae85662ab35fd7c9f647cc3e

 src/shcomplex.c

@@ -1,3 +1,4 @@
/* vim: set ts=4 sts=4 sw=4 noet : */
#include<math.h>
#include<stdlib.h>
#include<gsl/gsl_complex.h>
@@ -68,7 +69,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){
@@ -170,3 +171,17 @@
    
    return retval;
}


ts_double get_absolute_ulm2(ts_vesicle *vesicle, int l, int m){
   return gsl_complex_abs2(vesicle->sphHarmonics->ulmComplex[l][m]);
}


void solve_for_ulm2(ts_vesicle *vesicle){
   ts_double r0=getR0(vesicle);
   preparationSh(vesicle,r0);
   calculateUlmComplex(vesicle);
   storeUlmComplex2(vesicle);

}