libMesh
quadrature_gauss.h
Go to the documentation of this file.
1 // The libMesh Finite Element Library.
2 // Copyright (C) 2002-2025 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_GAUSS_H
21 #define LIBMESH_QUADRATURE_GAUSS_H
22 
23 // Local includes
24 #include "libmesh/quadrature.h"
25 
26 namespace libMesh
27 {
28 
39 class QGauss : public QBase
40 {
41 public:
42 
52  QGauss (unsigned int dim,
53  Order order=INVALID_ORDER) :
54  QBase(dim, order)
55  {
56  if (dim == 1)
57  init(EDGE2);
58  }
59 
64  QGauss (const QGauss &) = default;
65  QGauss (QGauss &&) = default;
66  QGauss & operator= (const QGauss &) = default;
67  QGauss & operator= (QGauss &&) = default;
68  virtual ~QGauss() = default;
69 
73  virtual QuadratureType type() const override;
74 
75  virtual std::unique_ptr<QBase> clone() const override;
76 
77 private:
78 
79  virtual void init_1D () override;
80  virtual void init_2D () override;
81  virtual void init_3D () override;
82 
88  void dunavant_rule(const Real rule_data[][4],
89  const unsigned int n_pts);
90 
91  void dunavant_rule2(const Real * wts,
92  const Real * a,
93  const Real * b,
94  const unsigned int * permutation_ids,
95  const unsigned int n_wts);
96 
102  void keast_rule(const Real rule_data[][4],
103  const unsigned int n_pts);
104 };
105 
106 } // namespace libMesh
107 
108 #endif // LIBMESH_QUADRATURE_GAUSS_H
Order
defines an enum for polynomial orders.
Definition: enum_order.h:40
virtual QuadratureType type() const override
unsigned int dim
QuadratureType
Defines an enum for currently available quadrature rules.
The libMesh namespace provides an interface to certain functionality in the library.
void dunavant_rule2(const Real *wts, const Real *a, const Real *b, const unsigned int *permutation_ids, const unsigned int n_wts)
virtual std::unique_ptr< QBase > clone() const override
virtual ~QGauss()=default
void dunavant_rule(const Real rule_data[][4], const unsigned int n_pts)
The Dunavant rules are for triangles.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
void keast_rule(const Real rule_data[][4], const unsigned int n_pts)
The Keast rules are for tets.
This class implements specific orders of Gauss quadrature.
virtual void init(const Elem &e, unsigned int p_level=invalid_uint)
Initializes the data structures for a quadrature rule for the element e.
Definition: quadrature.C:65
virtual void init_3D() override
Initializes the 3D quadrature rule by filling the points and weights vectors with the appropriate val...
QGauss(unsigned int dim, Order order=INVALID_ORDER)
Constructor.
virtual void init_2D() override
Initializes the 2D quadrature rule by filling the points and weights vectors with the appropriate val...
virtual void init_1D() override
Initializes the 1D quadrature rule by filling the points and weights vectors with the appropriate val...
The QBase class provides the basic functionality from which various quadrature rules can be derived...
Definition: quadrature.h:61
QGauss & operator=(const QGauss &)=default