]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Make a few parameters and variables constant.
authorbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 6 Nov 2006 05:40:12 +0000 (05:40 +0000)
committerbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 6 Nov 2006 05:40:12 +0000 (05:40 +0000)
git-svn-id: https://svn.dealii.org/trunk@14164 0785d39b-7218-0410-832d-ea1e28bc413d

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

index ac77bcadedd85e650d908e6e246e85c4934ecd5a..a1362d1078d158f33d17b3d3a2b6beb449a89d88 100644 (file)
@@ -89,9 +89,10 @@ class QGaussLobatto : public Quadrature<dim>
                                      * @param q  number of points.
                                      * @return vector containing nodes.
                                      */
-    std::vector<double> compute_quadrature_points(const unsigned int q,
-                                                  const int alpha,
-                                                  const int beta) const;
+    std::vector<double>
+    compute_quadrature_points (const unsigned int q,
+                              const int alpha,
+                              const int beta) const;
 
                                     /**
                                      * Compute Legendre-Gauss-Lobatto quadrature
@@ -99,9 +100,10 @@ class QGaussLobatto : public Quadrature<dim>
                                      * @param x  quadrature points.
                                      * @return vector containing weights.
                                      */
-    std::vector<double> compute_quadrature_weights(std::vector<double>& x,
-                                                   const int alpha,
-                                                   const int beta) const;
+    std::vector<double>
+    compute_quadrature_weights (const std::vector<double> &x,
+                               const int alpha,
+                               const int beta) const;
     
                                     /**
                                      * Evaluate a Jacobi polynomial
@@ -304,7 +306,7 @@ std::vector<double> QGaussLobatto<1>::
 compute_quadrature_points(const unsigned int, const int, const int) const;
 template <>
 std::vector<double> QGaussLobatto<1>::
-compute_quadrature_weights(std::vector<double>&, const int, const int) const;
+compute_quadrature_weights(const std::vector<double>&, const int, const int) const;
 template <>
 double QGaussLobatto<1>::
 JacobiP(const double, const int, const int, const unsigned int) const;
index 49a41311dccffad18b0abdcb2f62b4031c89572c..b556166fc4b9d45a7c061997b9f84fff40163671 100644 (file)
@@ -223,7 +223,7 @@ compute_quadrature_points(const unsigned int q,
 
 template <>
 std::vector<double> QGaussLobatto<1>::
-compute_quadrature_weights(std::vector<double>& x,
+compute_quadrature_weights(const std::vector<double> &x,
                           const int alpha,
                           const int beta) const
 {
@@ -231,8 +231,10 @@ compute_quadrature_weights(std::vector<double>& x,
   std::vector<double> w(q);
   double s = 0;
   
-  double factor = std::pow(2., alpha+beta+1)*gamma(alpha+q)*gamma(beta+q)/
-                 ((q-1)*gamma(q)*gamma(alpha+beta+q+1));
+  const double factor = std::pow(2., alpha+beta+1) *
+                       gamma(alpha+q) *
+                       gamma(beta+q) /
+                       ((q-1)*gamma(q)*gamma(alpha+beta+q+1));
   for (unsigned int i=0; i<q; ++i)
     {
       s = JacobiP(x[i], alpha, beta, q-1);

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.