..

pub/lib/special.cpp

@@ -84,9 +84,9 @@ void NSpecial::Test()
     static int Ni = 1000;
     
     for ( int j=0; j<Nj; ++j ){
-        double b = pow( 1.2, j-Nj );
+        double b = pow( 1.5, j-Nj );
         PSpec s( new CSpec );
-        s->z.push_back( PObjInt( new CObjInt (b) ) );
+        s->z.push_back( PObjDbl( new CObjDbl (b) ) );
         s->resize(Ni,false);
         for ( int i=0; i<Ni; ++i ){
             double a = b + pow( 1.2, i-Ni );

pub/lib/special.hpp

@@ -12,6 +12,8 @@
 #include <complex>
 using std::complex;
 
+#include "gsl/gsl_sf_trig.h"
+
 namespace NSpecial {
     void FourierCosine();
     void Test();
@@ -22,6 +24,8 @@ namespace NSpecial {
         complex<T> ci( 0.L, 1.L );
         complex<T> cone( 1.L, 0.L );
         return z==0 ? -ci : (cone-std::exp(ci*z))/z;
+//        T z2 = z / (T)(2.L);
+//        return -ci * std::sin(z2) / (z2) * std::exp(ci*z2);
     }
 
     template<typename T>