| | |
|---|
| | | /* vim: set ts=4 sts=4 sw=4 noet : */ |
|---|
| | | #include<stdio.h> |
|---|
| | | #include<math.h> |
|---|
| | | #include "general.h" |
|---|
| | |
|---|
| | | |
|---|
| | | |
|---|
| | | |
|---|
| | | ts_bool saveAvgUlm2(ts_vesicle *vesicle){ |
|---|
| | | |
|---|
| | | FILE *fh; |
|---|
| | | |
|---|
| | | fh=fopen("sph2out.dat", "w"); |
|---|
| | | if(fh==NULL){ |
|---|
| | | err("Cannot open file %s for writing"); |
|---|
| | | return TS_FAIL; |
|---|
| | | } |
|---|
| | | |
|---|
| | | ts_spharm *sph=vesicle->sphHarmonics; |
|---|
| | | ts_int i,j; |
|---|
| | | fprintf(fh,"l,\tm,\tulm^2avg\n"); |
|---|
| | | for(i=0;i<sph->l;i++){ |
|---|
| | | for(j=0;j<2*i+1;j++){ |
|---|
| | | fprintf(fh,"%d,\t%d,\t%e\n", i, j-i, sph->sumUlm2[i][j]/(ts_double)sph->N); |
|---|
| | | |
|---|
| | | } |
|---|
| | | fprintf(fh,"\n"); |
|---|
| | | } |
|---|
| | | fclose(fh); |
|---|
| | | return TS_SUCCESS; |
|---|
| | | } |
|---|
