| | |
| | | #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); |
| | |
| | | } |
| | | |
| | | void vesicle_calculate_ulm2(ts_vesicle *vesicle){ |
| | | complex_sph_free(vesicle->sphHarmonics); |
| | | //complex_sph_free(vesicle->sphHarmonics); |
| | | |
| | | vesicle->sphHarmonics=complex_sph_init(vesicle->vlist,21); |
| | | //vesicle->sphHarmonics=complex_sph_init(vesicle->vlist,21); |
| | | vesicle_volume(vesicle); |
| | | preparationSh(vesicle,getR0(vesicle)); |
| | | calculateUlmComplex(vesicle); |
| | |
| | | ent=list[n]; |
| | | i=rindex(ent->d_name,'.'); |
| | | if(i==NULL) { |
| | | continue; |
| | | continue; |
| | | } |
| | | if(strcmp(i+1,"vtu")==0){ |
| | | j=rindex(ent->d_name,'_'); |
| | |
| | | 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; |
| | | } |