Changeset 2634 in Sophya for trunk/SophyaLib/Samba

Timestamp:

Nov 19, 2004, 11:26:42 AM (21 years ago)

Author:

cmv

Message:

theta part of spherical harmonics, stand-alone routine cmv 19/11/04

File:

• trunk/SophyaLib/Samba/harmspher.cc (modified) (1 diff)

trunk/SophyaLib/Samba/harmspher.cc

@@ -79 +79 @@
  if(m>0) {          // Psph(0,0) == xlm[0]
    //double st = sqrt((1.-x)*(1.+x));
-   double st = sin(teta);
+   double st = fabs(sin(teta));
    for(int i=1;i<=m;i++) {
      xlm[0] *= -sqrt(((double)i+0.5)/(double)i) *st;