/usr/include/dolfin/mesh/LocalMeshData.h is in libdolfin-dev 2016.2.0-2.
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 | // Copyright (C) 2008 Ola Skavhaug
//
// 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/>.
//
// Modified by Anders Logg, 2008-2009.
//
// First added: 2008-11-28
// Last changed: 2013-02-18
//
// Modified by Anders Logg, 2008-2009.
// Modified by Kent-Andre Mardal, 2011.
#ifndef __LOCAL_MESH_DATA_H
#define __LOCAL_MESH_DATA_H
#include <cstdint>
#include <map>
#include <vector>
#include <boost/multi_array.hpp>
#include <dolfin/common/MPI.h>
#include <dolfin/common/Variable.h>
#include "CellType.h"
namespace dolfin
{
class Mesh;
/// This class stores mesh data on a local processor corresponding
/// to a portion of a (larger) global mesh.
///
/// Note that the data stored in this class does typically not
/// correspond to a topologically connected mesh; it merely stores a
/// list of vertex coordinates, a list of cell-vertex mappings and a
/// list of global vertex numbers for the locally stored vertices.
///
/// It is typically used for parsing meshes in parallel from mesh
/// XML files. After local mesh data has been parsed on each
/// processor, a subsequent repartitioning takes place: first a
/// geometric partitioning of the vertices followed by a
/// redistribution of vertex and cell data, and then a topological
/// partitioning again followed by redistribution of vertex and cell
/// data, at that point corresponding to topologically connected
/// meshes instead of local mesh data.
class LocalMeshData : public Variable
{
public:
/// Create empty local mesh data
explicit LocalMeshData(const MPI_Comm mpi_comm);
/// Create local mesh data for given mesh
explicit LocalMeshData(const Mesh& mesh);
/// Destructor
~LocalMeshData();
/// Check that all essential data has been intialized, and throw error
/// if there is a problem
void check() const;
/// Return informal string representation (pretty-print)
std::string str(bool verbose) const;
// Clear all data
void clear();
// Copy data from mesh
void extract_mesh_data(const Mesh& mesh);
// Broadcast mesh data from main process (used when Mesh is created
// on one process)
void broadcast_mesh_data(const MPI_Comm mpi_comm);
// Receive mesh data from main process
void receive_mesh_data(const MPI_Comm mpi_comm);
// Reorder cell data
void reorder();
// Holder for geometry data
struct Geometry
{
Geometry() : dim(-1), num_global_vertices(-1) {}
// Geometric dimension
int dim;
// Global number of vertices
std::int64_t num_global_vertices;
// Coordinates for all vertices stored on local processor
boost::multi_array<double, 2> vertex_coordinates;
// Global vertex indices for all vertices stored on local processor
std::vector<std::int64_t> vertex_indices;
void clear()
{
dim = -1;
num_global_vertices = -1;
vertex_coordinates.resize(boost::extents[0][0]);
vertex_indices.clear();
}
// Unpack received vertex coordinates
void unpack_vertex_coordinates(const std::vector<double>& values);
};
Geometry geometry;
// Holder for topology data
struct Topology
{
Topology() : dim(-1), num_global_cells(-1) {}
// Topological dimension
int dim;
// Global number of cells
std::int64_t num_global_cells;
// Number of vertices per cell
int num_vertices_per_cell;
// Global vertex indices for all cells stored on local processor
boost::multi_array<std::int64_t, 2> cell_vertices;
// Global cell numbers for all cells stored on local processor
std::vector<std::int64_t> global_cell_indices;
// Optional process owner for each cell in global_cell_indices
std::vector<int> cell_partition;
// Optional weight for each cell for partitioning
std::vector<std::size_t> cell_weight;
// Cell type
// FIXME: this should replace the need for num_vertices_per_cell and tdim
CellType::Type cell_type;
void clear()
{
dim = -1;
num_global_cells = -1;
num_vertices_per_cell = -1;
cell_vertices.resize(boost::extents[0][0]);
global_cell_indices.clear();
cell_partition.clear();
cell_weight.clear();
}
// Unpack received cell vertices
void unpack_cell_vertices(const std::vector<std::int64_t>& values);
};
Topology topology;
// Mesh domain data [dim](line, (cell_index, local_index, value))
std::map<std::size_t, std::vector<std::pair<std::pair<std::size_t,
std::size_t>, std::size_t>>> domain_data;
// Return MPI communicator
MPI_Comm mpi_comm() const
{ return _mpi_comm; }
private:
// MPI communicator
MPI_Comm _mpi_comm;
};
}
#endif
|