Trisurf Monte Carlo simulator
Samo Penic
2023-01-24 a8e354c7fad70eb7fdfda62ec83faf1be6c4ed44
src/tsmeasure.c
@@ -26,7 +26,7 @@
#include <snapshot.h>
#include<gsl/gsl_complex.h>
#include<gsl/gsl_complex_math.h>
#include "globals.h"
ts_vesicle *restoreVesicle(char *filename){
   ts_vesicle *vesicle = parseDump(filename);
@@ -34,7 +34,9 @@
}
void vesicle_calculate_ulm2(ts_vesicle *vesicle){
   vesicle->sphHarmonics=complex_sph_init(vesicle->vlist,21);
   //complex_sph_free(vesicle->sphHarmonics);
   //vesicle->sphHarmonics=complex_sph_init(vesicle->vlist,21);
   vesicle_volume(vesicle);
   preparationSh(vesicle,getR0(vesicle));
   calculateUlmComplex(vesicle);
@@ -65,7 +67,9 @@
      struct dirent *ent;
      ent=list[n];   
               i=rindex(ent->d_name,'.');
               if(i==NULL) continue;
               if(i==NULL) {
            continue;
      }
               if(strcmp(i+1,"vtu")==0){
                    j=rindex(ent->d_name,'_');
                    if(j==NULL) continue;
@@ -76,11 +80,16 @@
         vesicle=restoreVesicle(ent->d_name);
         vesicle_calculate_ulm2(vesicle);
                       tstep++;
         //vesicle_free(vesicle);
                    free(number);
         tape_free(vesicle->tape);
         vesicle_free(vesicle);
               }
      free(ent);
      }
   for (n = 0; n < count; n++)
     {
        free(list[n]);
     }
   free(list);
   return 0;
}