49 std::vector<double> Pk(N+1);
50 gsl_sf_legendre_Pl_array(N, x, &Pk[0]);
52 std::complex<double> P=0.0;
53 std::complex<double> I(0,1.0);
67 if (c->
type==
MAG)
return abs(P);
70 if (!finite(retVal)) {
71 throw std::runtime_error(
"Non-finite return value in LegendreExpansion");