]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Let the gamma function return a long double instead of an unsigned int, otherwise...
authorRalf Hartmann <Ralf.Hartmann@dlr.de>
Fri, 5 Jun 2009 03:56:02 +0000 (03:56 +0000)
committerRalf Hartmann <Ralf.Hartmann@dlr.de>
Fri, 5 Jun 2009 03:56:02 +0000 (03:56 +0000)
git-svn-id: https://svn.dealii.org/trunk@18907 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/base/include/base/quadrature_lib.h
deal.II/base/source/quadrature_lib.cc

index 300c7c03f7d2f5e51c557e30eb7c6e6af88171e2..ffc6ee4ee118a574d8c650cf50eb7d343c71650b 100644 (file)
@@ -134,7 +134,7 @@ class QGaussLobatto : public Quadrature<dim>
                                      * $ \Gamma(n) = (n-1)! $. 
                                      * @param n  point of evaluation (integer).
                                      */
-    unsigned int gamma(const unsigned int n) const;
+    long double gamma(const unsigned int n) const;
 };
 
   
@@ -513,7 +513,7 @@ template <>
 long double QGaussLobatto<1>::
 JacobiP(const long double, const int, const int, const unsigned int) const;
 template <>
-unsigned int 
+long double 
 QGaussLobatto<1>::gamma(const unsigned int n) const;
 
 template <> std::vector<double> QGaussLog<1>::set_quadrature_points(const unsigned int) const;
index b6a9ba31a072fc01bcf04b4c28cd5c019d006b25..954b9f010a1a3ebb377577d10cf09e070f64d537 100644 (file)
@@ -328,9 +328,9 @@ long double QGaussLobatto<1>::JacobiP(const long double x,
 
 
 template <>
-unsigned int QGaussLobatto<1>::gamma(const unsigned int n) const
+long double QGaussLobatto<1>::gamma(const unsigned int n) const
 {
-  unsigned int result = n - 1;
+  long double result = n - 1;
   for (int i=n-2; i>1; --i)
     result *= i;
   return result;

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.