/usr/include/libmesh/discontinuity_measure.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 | // $Id: discontinuity_measure.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 __discontinuity_measure_h___
#define __discontinuity_measure_h__
// C++ includes
#include <vector>
#include <string>
// Local Includes
#include "jump_error_estimator.h"
namespace libMesh
{
// Forward Declarations
class Point;
/**
* This class measures discontinuities between elements
* for debugging purposes. It derives from ErrorEstimator
* just in case someone finds it useful in a DG framework.
*
* @author Roy H. Stogner, 2006.
*/
class DiscontinuityMeasure : public JumpErrorEstimator
{
public:
/**
* Constructor. Responsible for initializing the _bc_function function
* pointer to NULL. Defaults to L2 norm; changes to system norm are
* ignored.
*/
DiscontinuityMeasure() : _bc_function(NULL) { error_norm = L2; }
/**
* Destructor.
*/
~DiscontinuityMeasure() {}
/**
* Register a user function to use in computing the essential BCs.
* The return value is std::pair<bool, Real>
*/
void attach_essential_bc_function (std::pair<bool,Real> fptr(const System& system,
const Point& p,
const std::string& var_name));
protected:
/**
* An initialization function, for requesting specific data from the FE
* objects
*/
virtual void initialize(const System& system,
ErrorVector& error_per_cell,
bool estimate_parent_error);
/**
* The function which calculates a normal derivative jump based error
* term on an internal side
*/
virtual void internal_side_integration();
/**
* The function which calculates a normal derivative jump based error
* term on a boundary side.
* Returns true if the flux bc function is in fact defined on the current side.
*/
virtual bool boundary_side_integration();
/**
* A pointer to the current System
*/
const System *my_system;
/**
* Pointer to function that returns BC information.
*/
std::pair<bool,Real> (* _bc_function) (const System& system,
const Point& p,
const std::string& var_name);
};
} // namespace libMesh
#endif
|