/usr/include/dolfin/common/IndexSet.h is in libdolfin1.0-dev 1.0.0-1.
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 | // Copyright (C) 2011 Anders Logg
//
// This file is part of DOLFIN.
//
// DOLFIN 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 3 of the License, or
// (at your option) any later version.
//
// DOLFIN 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 DOLFIN. If not, see <http://www.gnu.org/licenses/>.
//
// First added: 2011-02-07
// Last changed: 2011-08-28
#ifndef __INDEX_SET_H
#define __INDEX_SET_H
#include <vector>
#include "types.h"
namespace dolfin
{
/// This class provides an efficient data structure for index sets.
/// The cost of checking whether a given index is in the set is O(1)
/// and very very fast (optimal) at the cost of extra storage.
class IndexSet
{
public:
/// Create index set of given size
IndexSet(uint size) : _size(size), _has_index(size), _positions(size)
{
_indices.reserve(size);
clear();
}
/// Destructor
~IndexSet() {}
/// Return size of set
uint size() const
{ return _indices.size(); }
/// Check whether index is in set
bool has_index(uint index) const
{
dolfin_assert(index < _size);
return _has_index[index];
}
/// Return position (if any) for given index
uint find(uint index) const
{
dolfin_assert(index < _size);
if (!_has_index[index])
dolfin_error("IndexSet.h",
"locate position of index",
"Index %d is not in index set", index);
return _positions[index];
}
/// Return given index
uint& operator[] (uint i)
{
dolfin_assert(i < _indices.size());
return _indices[i];
}
/// Return given index (const version)
const uint& operator[] (uint i) const
{
dolfin_assert(i < _indices.size());
return _indices[i];
}
/// Insert index into set
void insert(uint index)
{
dolfin_assert(index < _size);
if (_has_index[index])
return;
_indices.push_back(index);
_has_index[index] = true;
_positions[index] = _indices.size() - 1;
}
/// Fill index set with indices 0, 1, 2, ..., size - 1
void fill()
{
_indices.clear();
for (uint i = 0; i < _size; i++)
_indices.push_back(i);
std::fill(_has_index.begin(), _has_index.end(), true);
}
/// Clear set
void clear()
{
_indices.clear();
std::fill(_has_index.begin(), _has_index.end(), false);
std::fill(_positions.begin(), _positions.end(), 0);
}
private:
// Size (maximum index + 1)
uint _size;
// Vector of indices
std::vector<uint> _indices;
// Indicators for which indices are in the set
std::vector<uint> _has_index;
// Mapping from indices to positions
std::vector<uint> _positions;
};
}
#endif
|