/usr/include/rheolef/geo_element_indirect.h is in librheolef-dev 6.5-1+b1.
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 | #ifndef _RHEOLEF_GEO_ELEMENT_INDIRECT_H
#define _RHEOLEF_GEO_ELEMENT_INDIRECT_H
///
/// This file is part of Rheolef.
///
/// Copyright (C) 2000-2009 Pierre Saramito <Pierre.Saramito@imag.fr>
///
/// Rheolef is free software; you can redistribute it and/or modify
/// it under the terms of the GNU General Public License as published by
/// the Free Software Foundation; either version 2 of the License, or
/// (at your option) any later version.
///
/// Rheolef 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 General Public License for more details.
///
/// You should have received a copy of the GNU General Public License
/// along with Rheolef; if not, write to the Free Software
/// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
///
/// =========================================================================
#include "rheolef/compiler.h"
namespace rheolef {
// =====================================================================
// geo_element_indirect = index + orient + face shift
// =====================================================================
// TODO: compact the sign bit and size_t in a long int
class geo_element_indirect {
public:
// typedefs:
typedef size_t size_type;
typedef short int orientation_type; // for sign (+1,-1)
typedef short int shift_type; // for face [0:4[ shift
// allocators:
geo_element_indirect ()
: _all(std::numeric_limits<size_type>::max()) {}
geo_element_indirect (orientation_type orient, size_type ige, size_type shift = 0)
: _all(std::numeric_limits<size_type>::max())
{ set (orient, ige, shift); }
// accessors:
size_type index() const { return _all & index_mask; }
orientation_type orientation() const { return (_all & orient_mask) ? -1 : 1; }
shift_type shift() const { return (_all & shift_mask) >> shift_position; }
// modifiers:
void set_orientation (orientation_type orient) {
_all = (_all & (~orient_mask));
if (orient < 0) _all = _all | orient_mask;
}
void set_shift (size_type shift) {
_all = (_all & (~shift_mask));
_all = _all | (shift_mask & (shift << shift_position));
}
void set_index (size_type ige) {
_all = (_all & (~index_mask));
_all = _all | (index_mask & ige);
}
void set (orientation_type orient, size_type ige, size_type shift = 0) {
set_orientation (orient);
set_index (ige);
set_shift (shift);
}
// i/o:
void dump(std::ostream& out = std::clog) const {
for (int i = 8*sizeof(size_type) - 1; i >= 0; i--) {
size_type bit = (_all & (size_type(1) << i)) >> i;
out << bit;
if (i % 8 == 0 && i != 0) out << " ";
}
}
template<class Archive>
void serialize (Archive& ar, const unsigned int version) { ar & _all; }
protected:
// implementation:
// left bit(0) = first left bit: orientation
static const size_type orient_position = 8*sizeof(size_type)-1;
static const size_type orient_mask = (size_type(1) << orient_position);
// left bits(1:3) = face shift (for 3d faces)
static const size_type shift_position = 8*sizeof(size_type)-4;
static const size_type shift_mask = (size_type(1) << shift_position) |
(size_type(1) << (shift_position+1)) |
(size_type(1) << (shift_position+2));
// others right bits = index of the subgeo (face or edge)
static const size_type index_mask = ~(orient_mask | shift_mask);
// data:
size_type _all;
};
// i/o:
std::istream& operator>> (std::istream& is, geo_element_indirect& x);
std::ostream& operator<< (std::ostream& os, const geo_element_indirect& x);
} // namespace rheolef
#ifdef _RHEOLEF_HAVE_MPI
#include "rheolef/distributed.h"
// =====================================================================
// Some serializable types, have a fixed amount of data stored at fixed field positions.
// When this is the case, boost::mpi can optimize their serialization and transmission to avoid extraneous
// copy operations.
// To enable this optimization, we specialize the type trait is_mpi_datatype, e.g.:
namespace boost {
namespace mpi {
template <> struct is_mpi_datatype<rheolef::geo_element_indirect> : mpl::true_ { };
} // namespace mpi
} // namespace boost
#endif // _RHEOLEF_HAVE_MPI
#endif // _RHEOLEF_GEO_ELEMENT_INDIRECT_H
|