From fda1ab6babed79842534b3a21a6ee96bc26f9d93 Mon Sep 17 00:00:00 2001
From: Samo Penic <samo.penic@fe.uni-lj.si>
Date: Tue, 16 Dec 2014 14:47:56 +0000
Subject: [PATCH] Real spherical harmonics are programmed. Inefficiently, but for testing purposes they works.
---
src/sh.c | 8 ++++----
1 files changed, 4 insertions(+), 4 deletions(-)
diff --git a/src/sh.c b/src/sh.c
index 9cb623e..add5a38 100644
--- a/src/sh.c
+++ b/src/sh.c
@@ -5,7 +5,7 @@
#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));
@@ -392,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++;
@@ -405,7 +405,7 @@
FILE *fh;
char filename[10000];
strcpy(filename, command_line_args.path);
- strcat(filename, "sph2out.dat");
+ strcat(filename, "sph2outOld.dat");
fh=fopen(filename, "w");
if(fh==NULL){
err("Cannot open file %s for writing");
@@ -417,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");
--
Gitblit v1.9.3