From f901ffe8aa9c458f0688d7deacc6d8708978f07d Mon Sep 17 00:00:00 2001
From: "Joachim Wuttke (h)" <j.wuttke@fz-juelich.de>
Date: Fri, 29 Jan 2016 10:36:44 +0100
Subject: [PATCH] ..
---
pub/lib/special.cpp | 4 ++--
pub/lib/special.hpp | 4 ++++
2 files changed, 6 insertions(+), 2 deletions(-)
diff --git a/pub/lib/special.cpp b/pub/lib/special.cpp
index 22e2fe7c..68bd6efc 100644
--- a/pub/lib/special.cpp
+++ b/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 );
diff --git a/pub/lib/special.hpp b/pub/lib/special.hpp
index 3cdda907..524e915e 100644
--- a/pub/lib/special.hpp
+++ b/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>
--
GitLab