Added missing header file

Samo Penic

2014-12-16 6bc95bfe19012a33f185adf1476e623d98d33465

 src/sh.c

@@ -2,9 +2,10 @@
#include<stdlib.h>
#include "general.h"
#include "sh.h"
#include "io.h"
#include <string.h>



/* do not use!!! */
ts_spharm *sph_init(ts_vertex_list *vlist, ts_uint l){
    ts_uint j,i;
    ts_spharm *sph=(ts_spharm *)malloc(sizeof(ts_spharm));
@@ -391,7 +392,7 @@
for(i=0;i<sph->l;i++){
    for(j=0;j<2*i+1;j++){
   /* DEBUG fprintf(stderr,"sph->sumUlm2[%d][%d]=%e\n",i,j,sph->ulm[i][j]* sph->ulm[i][j]); */
        sph->sumUlm2[i][j]+=sph->ulm[i][j]* sph->ulm[i][j];
        sph->sumUlm2Old[i][j]+=sph->ulm[i][j]* sph->ulm[i][j];
    }
}
   sph->N++;
@@ -402,8 +403,10 @@
ts_bool saveAvgUlm2(ts_vesicle *vesicle){

   FILE *fh;
	
   fh=fopen("sph2out.dat", "w");
    char filename[10000];
    strcpy(filename, command_line_args.path);
    strcat(filename, "sph2outOld.dat");
   fh=fopen(filename, "w");
   if(fh==NULL){
      err("Cannot open file %s for writing");
      return TS_FAIL;
@@ -414,7 +417,7 @@
   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,"%d,\t%d,\t%e\n", i, j-i, sph->sumUlm2Old[i][j]/(ts_double)sph->N);

          }
    fprintf(fh,"\n");