Skip to content

GitLab

Commit 48c4e4da authored by BRAMAS Berenger's avatar BRAMAS Berenger
Browse files

add a cpp file for cheb roots, in case of multiple inclusion of the header

parent 5684b4b4
Hide whitespace changes
Inline Side-by-side
Src/Kernels/Chebyshev/FChebRoots.cpp 0 → 100644
View file @ 48c4e4da
#include "FChebInterpolator.hpp"
// In case of multiple include of FChebInterpolator.hpp, this has to be defined only once!!
// order 2
template<> const double FChebRoots<2>::roots[] = {-0.707106781186548,
0.707106781186547};
// order 3
template<> const double FChebRoots<3>::roots[] = {-8.66025403784439e-01,
0.0,
8.66025403784438e-01};
// order 4
template<> const double FChebRoots<4>::roots[] = {-0.923879532511287,
-0.382683432365090,
0.382683432365090,
0.923879532511287};
// order 5
template<> const double FChebRoots<5>::roots[] = {-9.51056516295154e-01,
-5.87785252292473e-01,
0.0,
5.87785252292473e-01,
9.51056516295154e-01};
// order 6
template<> const double FChebRoots<6>::roots[] = {-0.965925826289068,
-0.707106781186548,
-0.258819045102521,
0.258819045102521,
0.707106781186547,
0.965925826289068};
// order 7
template<> const double FChebRoots<7>::roots[] = {-9.74927912181824e-01,
-7.81831482468030e-01,
-4.33883739117558e-01,
0.0,
4.33883739117558e-01,
7.81831482468030e-01,
9.74927912181824e-01};
// order 8
template<> const double FChebRoots<8>::roots[] = {-0.980785280403230,
-0.831469612302545,
-0.555570233019602,
-0.195090322016128,
0.195090322016128,
0.555570233019602,
0.831469612302545,
0.980785280403230};
// order 9
template<> const double FChebRoots<9>::roots[] = {-9.84807753012208e-01,
-8.66025403784439e-01,
-6.42787609686539e-01,
-3.42020143325669e-01,
0.0,
3.42020143325669e-01,
6.42787609686539e-01,
8.66025403784438e-01,
9.84807753012208e-01,};
// order 10
template<> const double FChebRoots<10>::roots[] = {-0.987688340595138,
-0.891006524188368,
-0.707106781186548,
-0.453990499739547,
-0.156434465040231,
0.156434465040231,
0.453990499739547,
0.707106781186547,
0.891006524188368,
0.987688340595138};
// order 11
template<> const double FChebRoots<11>::roots[] = {-9.89821441880933e-01,
-9.09631995354518e-01,
-7.55749574354258e-01,
-5.40640817455598e-01,
-2.81732556841430e-01,
0.0,
2.81732556841430e-01,
5.40640817455597e-01,
7.55749574354258e-01,
9.09631995354518e-01,
9.89821441880933e-01};
// order 12
template<> const double FChebRoots<12>::roots[] = {-0.991444861373810,
-0.923879532511287,
-0.793353340291235,
-0.608761429008721,
-0.382683432365090,
-0.130526192220052,
0.130526192220052,
0.382683432365090,
0.608761429008721,
0.793353340291235,
0.923879532511287,
0.991444861373810};
// order 13
template<> const double FChebRoots<13>::roots[] = {-9.92708874098054e-01,
-9.35016242685415e-01,
-8.22983865893656e-01,
-6.63122658240795e-01,
-4.64723172043769e-01,
-2.39315664287558e-01,
0.0,
2.39315664287557e-01,
4.64723172043769e-01,
6.63122658240795e-01,
8.22983865893656e-01,
9.35016242685415e-01,
9.92708874098054e-01};
// order 14
template<> const double FChebRoots<14>::roots[] = {
-0.99371220989324258353,
-0.94388333030836756290,
-0.84672419922828416835,
-0.70710678118654752440,
-0.53203207651533656356,
-0.33027906195516708177,
-0.11196447610330785847,
0.11196447610330785847,
0.33027906195516708177,
0.53203207651533656356,
0.70710678118654752440,
0.84672419922828416835,
0.94388333030836756290,
0.99371220989324258353};
Src/Kernels/Chebyshev/FChebRoots.hpp
View file @ 48c4e4da
......@@ -91,141 +91,46 @@ struct FChebRoots : FNoCopyable
}
};
// We declare the roots here only once Please look to .cpp for definitions
// order 2
template<> const double FChebRoots<2>::roots[] = {-0.707106781186548,
0.707106781186547};
template<> const double FChebRoots<2>::roots[];
// order 3
template<> const double FChebRoots<3>::roots[] = {-8.66025403784439e-01,
0.0,
8.66025403784438e-01};
template<> const double FChebRoots<3>::roots[];
// order 4
template<> const double FChebRoots<4>::roots[] = {-0.923879532511287,
-0.382683432365090,
0.382683432365090,
0.923879532511287};
template<> const double FChebRoots<4>::roots[];
// order 5
template<> const double FChebRoots<5>::roots[] = {-9.51056516295154e-01,
-5.87785252292473e-01,
0.0,
5.87785252292473e-01,
9.51056516295154e-01};
template<> const double FChebRoots<5>::roots[];
// order 6
template<> const double FChebRoots<6>::roots[] = {-0.965925826289068,
-0.707106781186548,
-0.258819045102521,
0.258819045102521,
0.707106781186547,
0.965925826289068};
template<> const double FChebRoots<6>::roots[];
// order 7
template<> const double FChebRoots<7>::roots[] = {-9.74927912181824e-01,
-7.81831482468030e-01,
-4.33883739117558e-01,
0.0,
4.33883739117558e-01,
7.81831482468030e-01,
9.74927912181824e-01};
template<> const double FChebRoots<7>::roots[];
// order 8
template<> const double FChebRoots<8>::roots[] = {-0.980785280403230,
-0.831469612302545,
-0.555570233019602,
-0.195090322016128,
0.195090322016128,
0.555570233019602,
0.831469612302545,
0.980785280403230};
template<> const double FChebRoots<8>::roots[];
// order 9
template<> const double FChebRoots<9>::roots[] = {-9.84807753012208e-01,
-8.66025403784439e-01,
-6.42787609686539e-01,
-3.42020143325669e-01,
0.0,
3.42020143325669e-01,
6.42787609686539e-01,
8.66025403784438e-01,
9.84807753012208e-01,};
template<> const double FChebRoots<9>::roots[];
// order 10
template<> const double FChebRoots<10>::roots[] = {-0.987688340595138,
-0.891006524188368,
-0.707106781186548,
-0.453990499739547,
-0.156434465040231,
0.156434465040231,
0.453990499739547,
0.707106781186547,
0.891006524188368,
0.987688340595138};
template<> const double FChebRoots<10>::roots[];
// order 11
template<> const double FChebRoots<11>::roots[] = {-9.89821441880933e-01,
-9.09631995354518e-01,
-7.55749574354258e-01,
-5.40640817455598e-01,
-2.81732556841430e-01,
0.0,
2.81732556841430e-01,
5.40640817455597e-01,
7.55749574354258e-01,
9.09631995354518e-01,
9.89821441880933e-01};
template<> const double FChebRoots<11>::roots[];
// order 12
template<> const double FChebRoots<12>::roots[] = {-0.991444861373810,
-0.923879532511287,
-0.793353340291235,
-0.608761429008721,
-0.382683432365090,
-0.130526192220052,
0.130526192220052,
0.382683432365090,
0.608761429008721,
0.793353340291235,
0.923879532511287,
0.991444861373810};
template<> const double FChebRoots<12>::roots[];
// order 13
template<> const double FChebRoots<13>::roots[] = {-9.92708874098054e-01,
-9.35016242685415e-01,
-8.22983865893656e-01,
-6.63122658240795e-01,
-4.64723172043769e-01,
-2.39315664287558e-01,
0.0,
2.39315664287557e-01,
4.64723172043769e-01,
6.63122658240795e-01,
8.22983865893656e-01,
9.35016242685415e-01,
9.92708874098054e-01};
template<> const double FChebRoots<13>::roots[];
// order 14
template<> const double FChebRoots<14>::roots[] = {
-0.99371220989324258353,
-0.94388333030836756290,
-0.84672419922828416835,
-0.70710678118654752440,
-0.53203207651533656356,
-0.33027906195516708177,
-0.11196447610330785847,
0.11196447610330785847,
0.33027906195516708177,
0.53203207651533656356,
0.70710678118654752440,
0.84672419922828416835,
0.94388333030836756290,
0.99371220989324258353};
template<> const double FChebRoots<14>::roots[];
#endif
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment