* Evaluate the Gamma function $ \Gamma(n) = (n-1)! $.
* @param n point of evaluation (integer).
*/
- long double gamma(const unsigned int n) const;
+ static long double gamma(const unsigned int n);
};
/**
* Compute the points of the quadrature formula.
*/
+ static
std::vector<double>
- get_quadrature_points(const unsigned int n) const;
+ get_quadrature_points(const unsigned int n);
/**
* Compute the weights of the quadrature formula.
*/
+ static
std::vector<double>
- get_quadrature_weights(const unsigned int n) const;
+ get_quadrature_weights(const unsigned int n);
};
JacobiP(const long double, const int, const int, const unsigned int) const;
template <>
long double
-QGaussLobatto<1>::gamma(const unsigned int n) const;
+QGaussLobatto<1>::gamma(const unsigned int n);
-template <> std::vector<double> QGaussLog<1>::get_quadrature_points(const unsigned int) const;
-template <> std::vector<double> QGaussLog<1>::get_quadrature_weights(const unsigned int) const;
+template <> std::vector<double> QGaussLog<1>::get_quadrature_points(const unsigned int);
+template <> std::vector<double> QGaussLog<1>::get_quadrature_weights(const unsigned int);
template <> QMidpoint<1>::QMidpoint ();
template <> QTrapez<1>::QTrapez ();
template <>
-long double QGaussLobatto<1>::gamma(const unsigned int n) const
+long double QGaussLobatto<1>::gamma(const unsigned int n)
{
long double result = n - 1;
for (int i=n-2; i>1; --i)
template <>
std::vector<double>
-QGaussLog<1>::get_quadrature_points(const unsigned int n) const
+QGaussLog<1>::get_quadrature_points(const unsigned int n)
{
std::vector<double> q_points(n);
template <>
std::vector<double>
-QGaussLog<1>::get_quadrature_weights(const unsigned int n) const
+QGaussLog<1>::get_quadrature_weights(const unsigned int n)
{
std::vector<double> quadrature_weights(n);