/usr/include/libmesh/face_inf_quad.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 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 | // $Id: face_inf_quad.h 4279 2011-03-21 17:01:31Z 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 __face_inf_quad_h__
#define __face_inf_quad_h__
// C++ includes
// Local includes
#include "libmesh_config.h"
#ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS
#include "elem.h"
namespace libMesh
{
// Forward declarations
/**
* The \p InfQuad is an abstract element type that lives in
* two dimensions. Here, an infinite face is @e always a quadrilateral,
* so this class is directly derived from \p Elem, without an intermediate
* \p InfFace class or so.
* It looks like this:
\verbatim
closer to infinity
| |
| |
side 2 | | side 1
| |
| |
----------- base side
side 0
\endverbatim
*/
// ------------------------------------------------------------
// InfQuad class definition
class InfQuad : public Elem
{
public:
/**
* Constructor. Derived classes implement 'true' elements.
*/
InfQuad (const unsigned int nn,
Elem* p) :
Elem(nn, InfQuad::n_sides(), p) {}
/**
* Constructor. Explicitly specifies the number of
* nodes and neighbors for which storage will be allocated.
*/
InfQuad (const unsigned int nn,
const unsigned int ns,
Elem* p) :
Elem(nn, ns, p) {}
/**
* @returns 2, the dimensionality of the object.
*/
unsigned int dim() const { return 2; }
// /**
// * @returns 2 for the base, 1 otherwise
// */
// unsigned int n_children_per_side(const unsigned int s) const;
/**
* @returns 3. Infinite faces have one side less
* than their conventional counterparts, since one
* side is supposed to be located at infinity.
*/
unsigned int n_sides() const { return 3; }
/**
* @returns 4. All infinite quads (in our setting) have 4 vertices.
*/
unsigned int n_vertices() const { return 4; }
/**
* @returns 3. All infinite quads have 1 edge in the
* base, and 2 perpendicular to the base.
*/
unsigned int n_edges() const { return 3; }
/**
* @returns 0. All 2D elements have no faces, just
* edges.
*/
unsigned int n_faces() const { return 0; }
/**
* @returns 2
*/
unsigned int n_children() const { return 2; }
/*
* @returns true iff the specified child is on the
* specified side
*/
virtual bool is_child_on_side(const unsigned int c,
const unsigned int s) const;
/**
* @returns an id associated with the \p s side of this element.
* The id is not necessariy unique, but should be close. This is
* particularly useful in the \p MeshBase::find_neighbors() routine.
*/
unsigned int key (const unsigned int s) const;
/**
* @returns a primitive (2-noded) edge or infedge for
* edge \p i.
*/
AutoPtr<Elem> side (const unsigned int i) const;
/**
* build_edge and build_side are identical in 2D
*/
AutoPtr<Elem> build_edge (const unsigned int i) const
{ return build_side(i); }
/**
* Based on the quality metric \p q specified by the user,
* returns a quantitative assessment of element quality.
*/
Real quality (const ElemQuality q) const;
/**
* Returns the suggested quality bounds for
* the hex based on quality measure q. These are
* the values suggested by the CUBIT User's Manual.
*/
std::pair<Real, Real> qual_bounds (const ElemQuality q) const;
/**
* @returns \p true. All classes derived from \p InfQuad
* are infinite elements.
*/
bool infinite () const { return true; }
/**
* @returns the origin of this infinite element.
*/
Point origin () const;
};
// ------------------------------------------------------------
// InfQuad class member functions
// inline
// unsigned int InfQuad::n_children_per_side(const unsigned int s) const
// {
// libmesh_assert (s < this->n_sides());
// switch (s)
// {
// case 0:
// // every infinite face has 2 children in the base edge
// return 2;
// default:
// // on infinite edges only 1 child
// return 1;
// }
// }
inline
Point InfQuad::origin () const
{
return ( this->point(0)*2 - this->point(this->n_vertices()/2) );
}
} // namespace libMesh
#endif // ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS
#endif
|