quadrature_grid_1D.C
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 // C++ includes
21 
22 // Local includes
24 
25 namespace libMesh
26 {
27 
28 
29 
31  unsigned int)
32 {
33  //----------------------------------------------------------------------
34  // 1D quadrature rules
35 
36  // We ignore p - the grid rule is just for experimentation
37 
38  _points.resize(_order + 1);
39  _weights.resize(_order + 1);
40  const Real dx = Real(2)/(_order+1);
41  for (int i = 0; i != _order + 1; ++i)
42  {
43  _points[i](0) = (i+0.5)*dx-1;
44  _weights[i] = dx;
45  }
46  return;
47 }
48 
49 } // namespace libMesh
virtual void init_1D(const ElemType _type=INVALID_ELEM, unsigned int p_level=0) override
std::vector< Point > _points
Definition: quadrature.h:349
std::vector< Real > _weights
Definition: quadrature.h:355
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real