Gaunt.h

#ifndef gridripper_multipole_Gaunt_h
#define gridripper_multipole_Gaunt_h


#include <gridripper/math.h>
#include <gridripper/tvector.h>


namespace gridripper { namespace multipole {


extern GReal_t YYY_coupling(const int l1, const int l2, const int l3, 
                            const int m1, const int m2, const int m3);


extern bool YYY_coupling_zeros(const int l1, const int l2, const int l3, 
                               const int m1, const int m2, const int m3);


extern unsigned int YYY_coupling_size();


extern GReal_t YreYreYre_coupling(const int l1, const int l2, const int l3, 
                                  const int m1, const int m2, const int m3);


extern GComplex_t YYreY_coupling(const int l1, const int l2, const int l3, 
                                 const int m1, const int m2, const int m3);


extern GReal_t intYl1ExpMIPhi(const int l);


extern GReal_t ExpPmIPhiYY_coupling(const int l1, const int m1, 
                                    const int l2, const int sign);


extern unsigned int ExpPmIPhiYY_coupling_size();


extern GReal_t intYl0(const int l, const GReal_t theta);


extern GReal_t intBarYl1mYl2m(const int l1, const int l2, const int m, const GReal_t theta, tvector< tvector< tvector<GReal_t> > > *coeffstorage);


} } // namespace gridripper::multipole


#endif /* gridripper_multipole_Gaunt_h */