From: bangerth Date: Mon, 6 Nov 2006 05:40:12 +0000 (+0000) Subject: Make a few parameters and variables constant. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=5742b29faeb6c4d26ca14681e18f4c5da5179920;p=dealii-svn.git Make a few parameters and variables constant. git-svn-id: https://svn.dealii.org/trunk@14164 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/base/include/base/quadrature_lib.h b/deal.II/base/include/base/quadrature_lib.h index ac77bcaded..a1362d1078 100644 --- a/deal.II/base/include/base/quadrature_lib.h +++ b/deal.II/base/include/base/quadrature_lib.h @@ -89,9 +89,10 @@ class QGaussLobatto : public Quadrature * @param q number of points. * @return vector containing nodes. */ - std::vector compute_quadrature_points(const unsigned int q, - const int alpha, - const int beta) const; + std::vector + 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 * @param x quadrature points. * @return vector containing weights. */ - std::vector compute_quadrature_weights(std::vector& x, - const int alpha, - const int beta) const; + std::vector + compute_quadrature_weights (const std::vector &x, + const int alpha, + const int beta) const; /** * Evaluate a Jacobi polynomial @@ -304,7 +306,7 @@ std::vector QGaussLobatto<1>:: compute_quadrature_points(const unsigned int, const int, const int) const; template <> std::vector QGaussLobatto<1>:: -compute_quadrature_weights(std::vector&, const int, const int) const; +compute_quadrature_weights(const std::vector&, const int, const int) const; template <> double QGaussLobatto<1>:: JacobiP(const double, const int, const int, const unsigned int) const; diff --git a/deal.II/base/source/quadrature_lib.cc b/deal.II/base/source/quadrature_lib.cc index 49a41311dc..b556166fc4 100644 --- a/deal.II/base/source/quadrature_lib.cc +++ b/deal.II/base/source/quadrature_lib.cc @@ -223,7 +223,7 @@ compute_quadrature_points(const unsigned int q, template <> std::vector QGaussLobatto<1>:: -compute_quadrature_weights(std::vector& x, +compute_quadrature_weights(const std::vector &x, const int alpha, const int beta) const { @@ -231,8 +231,10 @@ compute_quadrature_weights(std::vector& x, std::vector 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