diff options
author | Gael Guennebaud <g.gael@free.fr> | 2009-12-04 12:40:29 +0100 |
---|---|---|
committer | Gael Guennebaud <g.gael@free.fr> | 2009-12-04 12:40:29 +0100 |
commit | ea684af6b41cdcfd9a8d2edbda37a866827e5347 (patch) | |
tree | 0f3b9b03a1803327dab5b468e44152666304b58e | |
parent | 7aad4341609260356845b16057fd581302c5072f (diff) | |
parent | ffaea19a70d22657008645e97a0b9efe2fcc9451 (diff) |
merge
-rw-r--r-- | Eigen/src/Geometry/Quaternion.h | 2 | ||||
-rw-r--r-- | unsupported/Eigen/AutoDiff | 2 | ||||
-rw-r--r-- | unsupported/Eigen/src/FFT/ei_kissfft_impl.h | 8 | ||||
-rw-r--r-- | unsupported/Eigen/src/NumericalDiff/NumericalDiff.h | 5 | ||||
-rw-r--r-- | unsupported/test/FFT.cpp | 7 | ||||
-rw-r--r-- | unsupported/test/FFTW.cpp | 3 | ||||
-rw-r--r-- | unsupported/test/autodiff.cpp | 12 |
7 files changed, 20 insertions, 19 deletions
diff --git a/Eigen/src/Geometry/Quaternion.h b/Eigen/src/Geometry/Quaternion.h index 7f102c942..5578ce9f4 100644 --- a/Eigen/src/Geometry/Quaternion.h +++ b/Eigen/src/Geometry/Quaternion.h @@ -317,7 +317,7 @@ class Map<Quaternion<_Scalar>, PacketAccess > * The pointer \a coeffs must reference the four coeffecients of Quaternion in the following order: * \code *coeffs == {x, y, z, w} \endcode * - * If the template paramter PacketAccess is set to Aligned, then the pointer coeffs must be aligned. */ + * If the template parameter PacketAccess is set to Aligned, then the pointer coeffs must be aligned. */ EIGEN_STRONG_INLINE Map(const Scalar* coeffs) : m_coeffs(coeffs) {} inline Coefficients& coeffs() { return m_coeffs;} diff --git a/unsupported/Eigen/AutoDiff b/unsupported/Eigen/AutoDiff index 229c15e69..a8b3ea90a 100644 --- a/unsupported/Eigen/AutoDiff +++ b/unsupported/Eigen/AutoDiff @@ -25,7 +25,7 @@ #ifndef EIGEN_AUTODIFF_MODULE #define EIGEN_AUTODIFF_MODULE -#include <Eigen/Core> +#include <Eigen/Array> namespace Eigen { diff --git a/unsupported/Eigen/src/FFT/ei_kissfft_impl.h b/unsupported/Eigen/src/FFT/ei_kissfft_impl.h index 5c958d1ec..2dff2bd00 100644 --- a/unsupported/Eigen/src/FFT/ei_kissfft_impl.h +++ b/unsupported/Eigen/src/FFT/ei_kissfft_impl.h @@ -160,7 +160,7 @@ struct ei_kiss_cpx_fft scratch[0]=scratch[1]-scratch[2]; tw1 += fstride; tw2 += fstride*2; - Fout[m] = Complex( Fout->real() - .5*scratch[3].real() , Fout->imag() - .5*scratch[3].imag() ); + Fout[m] = Complex( Fout->real() - Scalar(.5)*scratch[3].real() , Fout->imag() - Scalar(.5)*scratch[3].imag() ); scratch[0] *= epi3.imag(); *Fout += scratch[3]; Fout[m2] = Complex( Fout[m].real() + scratch[0].imag() , Fout[m].imag() - scratch[0].real() ); @@ -377,10 +377,10 @@ struct ei_kissfft_impl std::vector<Complex> m_tmpBuf2; inline - int PlanKey(int nfft,bool isinverse) const { return (nfft<<1) | isinverse; } + int PlanKey(int nfft, bool isinverse) const { return (nfft<<1) | int(isinverse); } inline - PlanData & get_plan(int nfft,bool inverse) + PlanData & get_plan(int nfft, bool inverse) { // TODO look for PlanKey(nfft, ! inverse) and conjugate the twiddles PlanData & pd = m_plans[ PlanKey(nfft,inverse) ]; @@ -400,7 +400,7 @@ struct ei_kissfft_impl int ncfft= ncfft2<<1; Scalar pi = acos( Scalar(-1) ); for (int k=1;k<=ncfft2;++k) - twidref[k-1] = exp( Complex(0,-pi * ((double) (k) / ncfft + .5) ) ); + twidref[k-1] = exp( Complex(0,-pi * (Scalar(k) / ncfft + Scalar(.5)) ) ); } return &twidref[0]; } diff --git a/unsupported/Eigen/src/NumericalDiff/NumericalDiff.h b/unsupported/Eigen/src/NumericalDiff/NumericalDiff.h index 98872e0bc..db6f791df 100644 --- a/unsupported/Eigen/src/NumericalDiff/NumericalDiff.h +++ b/unsupported/Eigen/src/NumericalDiff/NumericalDiff.h @@ -45,10 +45,11 @@ enum NumericalDiffMode { * * Currently only "Forward" and "Central" scheme are implemented. */ -template<typename Functor, NumericalDiffMode mode=Forward> -class NumericalDiff : public Functor +template<typename _Functor, NumericalDiffMode mode=Forward> +class NumericalDiff : public _Functor { public: + typedef _Functor Functor; typedef typename Functor::Scalar Scalar; typedef typename Functor::InputType InputType; typedef typename Functor::ValueType ValueType; diff --git a/unsupported/test/FFT.cpp b/unsupported/test/FFT.cpp index 4a6d495cc..b029ba655 100644 --- a/unsupported/test/FFT.cpp +++ b/unsupported/test/FFT.cpp @@ -45,13 +45,10 @@ complex<long double> promote(long double x) { return complex<long double>( x); long double totalpower=0; long double difpower=0; cerr <<"idx\ttruth\t\tvalue\t|dif|=\n"; + long double pi = acos((long double)-1); for (size_t k0=0;k0<size_t(fftbuf.size());++k0) { complex<long double> acc = 0; -#ifdef _GNU_SOURCE - long double phinc = -2.*k0* M_PIl / timebuf.size(); -#else - long double phinc = -2.*k0* M_PI / timebuf.size(); -#endif + long double phinc = -2.*k0* pi / timebuf.size(); for (size_t k1=0;k1<size_t(timebuf.size());++k1) { acc += promote( timebuf[k1] ) * exp( complex<long double>(0,k1*phinc) ); } diff --git a/unsupported/test/FFTW.cpp b/unsupported/test/FFTW.cpp index cf7be75aa..c182cab9d 100644 --- a/unsupported/test/FFTW.cpp +++ b/unsupported/test/FFTW.cpp @@ -45,10 +45,11 @@ complex<long double> promote(long double x) { return complex<long double>( x); { long double totalpower=0; long double difpower=0; + long double pi = acos((long double)-1 ); cerr <<"idx\ttruth\t\tvalue\t|dif|=\n"; for (size_t k0=0;k0<fftbuf.size();++k0) { complex<long double> acc = 0; - long double phinc = -2.*k0* M_PIl / timebuf.size(); + long double phinc = -2.*k0* pi / timebuf.size(); for (size_t k1=0;k1<timebuf.size();++k1) { acc += promote( timebuf[k1] ) * exp( complex<long double>(0,k1*phinc) ); } diff --git a/unsupported/test/autodiff.cpp b/unsupported/test/autodiff.cpp index a96927b41..cac08cd4b 100644 --- a/unsupported/test/autodiff.cpp +++ b/unsupported/test/autodiff.cpp @@ -31,7 +31,7 @@ EIGEN_DONT_INLINE Scalar foo(const Scalar& x, const Scalar& y) // return x+std::sin(y); EIGEN_ASM_COMMENT("mybegin"); return static_cast<Scalar>(x*2 - std::pow(x,2) + 2*std::sqrt(y*y) - 4 * std::sin(x) + 2 * std::cos(y) - std::exp(-0.5*x*x)); -// return y/x;// - y*2; +// return x - y;//x*2 -std::pow(x,2);//(2*y/x);// - y*2; EIGEN_ASM_COMMENT("myend"); } @@ -137,10 +137,12 @@ template<typename Func> void forward_jacobian(const Func& f) void test_autodiff_scalar() { std::cerr << foo<float>(1,2) << "\n"; - AutoDiffScalar<Vector2f> ax(1,Vector2f::UnitX()); - AutoDiffScalar<Vector2f> ay(2,Vector2f::UnitY()); - std::cerr << foo<AutoDiffScalar<Vector2f> >(ax,ay).value() << " <> " - << foo<AutoDiffScalar<Vector2f> >(ax,ay).derivatives().transpose() << "\n\n"; + typedef AutoDiffScalar<Vector2f> AD; + AD ax(1,Vector2f::UnitX()); + AD ay(2,Vector2f::UnitY()); + foo<AD>(ax,ay); + std::cerr << foo<AD>(ax,ay).value() << " <> " + << foo<AD>(ax,ay).derivatives().transpose() << "\n\n"; } void test_autodiff_jacobian() |