/usr/include/libmesh/quadrature_trap.h is in libmesh-dev 0.7.1-2ubuntu1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 | // $Id: quadrature_trap.h 3874 2010-07-02 21:57:26Z roystgnr $
// The libMesh Finite Element Library.
// Copyright (C) 2002-2008 Benjamin S. Kirk, John W. Peterson, Roy H. Stogner
// This library is free software; you can redistribute it and/or
// modify it under the terms of the GNU Lesser General Public
// License as published by the Free Software Foundation; either
// version 2.1 of the License, or (at your option) any later version.
// This library is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
// Lesser General Public License for more details.
// You should have received a copy of the GNU Lesser General Public
// License along with this library; if not, write to the Free Software
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
#ifndef __quadrature_trap_h__
#define __quadrature_trap_h__
// C++ includes
// Local includes
#include "quadrature.h"
namespace libMesh
{
// ------------------------------------------------------------
// QTrap class definition
/**
* This class implemenets trapezoidal quadratue. Sometimes
* also known as Newton-Cotes quadrature with two points. These rules
* sample at the corners and will integrate linears exactly.
*/
class QTrap : public QBase
{
public:
/**
* Constructor. Declares the dimension of the quadrature rule.
*/
QTrap (const unsigned int _dim);
/**
* Destructor. Empty.
*/
~QTrap() {}
/**
* @returns \p QTRAP
*/
QuadratureType type() const { return QTRAP; }
private:
void init_1D (const ElemType _type=INVALID_ELEM,
unsigned int p_level=0);
void init_2D (const ElemType _type=INVALID_ELEM,
unsigned int p_level=0);
void init_3D (const ElemType _type=INVALID_ELEM,
unsigned int p_level=0);
};
// ------------------------------------------------------------
// QTrap class members
inline
QTrap::QTrap(const unsigned int d) : QBase(d,FIRST)
{
// explicitly call the init function in 1D since the
// other tensor-product rules require this one.
// note that EDGE will not be used internally, however
// if we called the function with INVALID_ELEM it would try to
// be smart and return, thinking it had already done the work.
if (_dim == 1)
init(EDGE2);
}
} // namespace libMesh
#endif
|