quadrature_jacobi.h
Go to the documentation of this file.
1 // The libMesh Finite Element Library.
2 // Copyright (C) 2002-2018 Benjamin S. Kirk, John W. Peterson, Roy H. Stogner
3 
4 // This library is free software; you can redistribute it and/or
5 // modify it under the terms of the GNU Lesser General Public
6 // License as published by the Free Software Foundation; either
7 // version 2.1 of the License, or (at your option) any later version.
8 
9 // This library is distributed in the hope that it will be useful,
10 // but WITHOUT ANY WARRANTY; without even the implied warranty of
11 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
12 // Lesser General Public License for more details.
13 
14 // You should have received a copy of the GNU Lesser General Public
15 // License along with this library; if not, write to the Free Software
16 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
17 
18 
19 
20 #ifndef LIBMESH_QUADRATURE_JACOBI_H
21 #define LIBMESH_QUADRATURE_JACOBI_H
22 
23 // Local includes
24 #include "libmesh/quadrature.h"
25 
26 namespace libMesh
27 {
28 
49 class QJacobi final : public QBase
50 {
51 public:
52 
56  QJacobi (unsigned int dim,
57  Order order=INVALID_ORDER,
58  unsigned int a=1,
59  unsigned int b=0) :
60  QBase(dim, order),
61  _alpha(a),
62  _beta(b)
63  {
64  if (_dim == 1)
65  init(EDGE2);
66  }
67 
72  QJacobi (const QJacobi &) = default;
73  QJacobi (QJacobi &&) = default;
74  QJacobi & operator= (const QJacobi &) = default;
75  QJacobi & operator= (QJacobi &&) = default;
76  virtual ~QJacobi() = default;
77 
82  virtual QuadratureType type() const override;
83 
84 private:
85  const unsigned int _alpha;
86  const unsigned int _beta;
87 
88  virtual void init_1D (const ElemType _type=INVALID_ELEM,
89  unsigned int p_level=0) override;
90 };
91 
92 } // namespace libMesh
93 
94 #endif // LIBMESH_QUADRATURE_JACOBI_H
virtual void init_1D(const ElemType _type=INVALID_ELEM, unsigned int p_level=0) override
virtual void init(const ElemType type=INVALID_ELEM, unsigned int p_level=0)
Definition: quadrature.C:28
QJacobi(unsigned int dim, Order order=INVALID_ORDER, unsigned int a=1, unsigned int b=0)
const unsigned int _alpha
ElemType _type
Definition: quadrature.h:337
unsigned int _dim
Definition: quadrature.h:325
const unsigned int _beta
QJacobi & operator=(const QJacobi &)=default
Implements 1D Gauss-Jacobi quadrature rules of various orders.
virtual QuadratureType type() const override
virtual ~QJacobi()=default
Base class for all quadrature families and orders.
Definition: quadrature.h:62