Trisurf Monte Carlo simulator
Samo Penic
2012-06-07 824cfd6fc38d1c583a33cb4d1a146d97b15cb543
src/sh.c
@@ -64,11 +64,11 @@
      fatal("Error using shY function!",1);
   fac1=1.0;
   for(i=1; i<=l-m;i++){
   for(i=1; i<=l-abs(m);i++){
      fac1 *= i;
   }
   fac2=1.0;
   for(i=1; i<=l+m;i++){
   for(i=1; i<=l+abs(m);i++){
      fac2 *= i;
   }
@@ -81,9 +81,9 @@
   else {
      //K=pow(-1.0,abs(m))*sqrt(1.0/(2.0*M_PI))*cos(m*fi);
      if(abs(m)%2==0)
      K=sqrt(1.0/(M_PI))*sin(m*fi);
      K=sqrt(1.0/(M_PI))*cos(m*fi);
      else
      K=-sqrt(1.0/(M_PI))*sin(m*fi);
      K=-sqrt(1.0/(M_PI))*cos(m*fi);
   }
   
   return K*sqrt((2.0*l+1.0)/2.0*fac1/fac2)*plgndr(l,abs(m),cos(theta));