SphericalHarmonics_Coefficients.h

Go to the documentation of this file.

/*******************************************************************
 *  Copyright (C) 2010 by Federico Marulli                         *
 *  federico.marulli3@unibo.it                                     *
 *                                                                 *
 *  This program is free software; you can redistribute it and/or  *
 *  modify it under the terms of the GNU General Public License as *
 *  published by the Free Software Foundation; either version 2 of *
 *  the License, or (at your option) any later version.            *
 *                                                                 *
 *  This program is distributed in the hope that it will be useful,*
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of *
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the  *
 *  GNU General Public License for more details.                   *
 *                                                                 *
 *  You should have received a copy of the GNU General Public      *
 *  License along with this program; if not, write to the Free     *
 *  Software Foundation, Inc.,                                     *
 *  59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.      *
 *******************************************************************/
 
#ifndef __SHCOEFF__
#define __SHCOEFF__
 
#include "Func.h"
 
namespace cbl {
 
  namespace glob {
 
    class SphericalHarmonics_Coefficients {
 
    protected:
 
      int m_nbins;
 
      int m_norder;
 
      int m_lmax;
 
      int m_n_sph;
 
      std::vector<int> m_n_sph_l;
 
      std::vector<std::vector<std::complex<double>>> m_alm;
 
      std::vector<double> m_normalization;
 
      std::vector<double> m_Plm;
    
      std::vector<std::complex<double>> m_sph;
 
    public:
 
      SphericalHarmonics_Coefficients () {}
 
      SphericalHarmonics_Coefficients (const int norder, const int nbins=1) { initialize(norder, nbins); }
 
      ~SphericalHarmonics_Coefficients () {}
 
 
      double real (const int n, const int bin=0) { return m_alm[bin][n].real();} 
 
      double imag (const int n, const int bin=0) { return m_alm[bin][n].imag();} 
 
      void initialize (const int norder, const int nbins=1);
 
      void reset ();
 
      std::vector<std::complex<double>> alm(const double xx, const double yy, const double zz);
 
      void add (const std::vector<std::complex<double>> alm, const double ww, const int bin=0);
 
      void add (const double xx, const double yy, const double zz, const double ww, const int bin=0);
 
      double power (const int l, const int bin1, const int bin2);
    };
 
  }
}
 
#endif