/usr/include/dolfin/io/XMLMeshValueCollection.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 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 | // 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-06-30
// Last changed: 2011-11-15
#ifndef __XML_MESH_VALUE_COLLECTION_H
#define __XML_MESH_VALUE_COLLECTION_H
#include <dolfin/mesh/MeshValueCollection.h>
#include "pugixml.hpp"
#include "xmlutils.h"
#include "XMLMesh.h"
namespace dolfin
{
class XMLMeshValueCollection
{
public:
// Read mesh value collection from XML file
template <typename T>
static void read(MeshValueCollection<T>& mesh_value_collection,
const std::string type,
const pugi::xml_node xml_node);
/// Write mesh value collection to XML file
template<typename T>
static void write(const MeshValueCollection<T>& mesh_value_collection,
const std::string type,
pugi::xml_node xml_node);
};
//---------------------------------------------------------------------------
template <typename T>
void XMLMeshValueCollection::read(MeshValueCollection<T>& mesh_value_collection,
const std::string type,
const pugi::xml_node xml_node)
{
// Get node
boost::shared_ptr<const pugi::xml_node>
mvc_node = get_node(xml_node, "mesh_value_collection");
dolfin_assert(mvc_node);
// Get attributes
const std::string type_file = mvc_node->attribute("type").value();
const uint dim = mvc_node->attribute("dim").as_uint();
// Set dim
mesh_value_collection.set_dim(dim);
// Check that types match
if (type != type_file)
{
dolfin_error("XMLMeshValueCollection.h",
"read mesh value collection from XML file",
"Type mismatch, found \"%s\" but expecting \"%s\"",
type_file.c_str(), type.c_str());
}
// Check that dimension matches
if (mesh_value_collection.dim() != dim)
{
dolfin_error("XMLMeshValueCollection.h",
"read mesh value collection from XML file",
"Dimension mismatch, found %d but expecting %d",
dim, mesh_value_collection.dim());
}
// Clear old values
mesh_value_collection.clear();
// Choose data type
if (type == "uint")
{
for (pugi::xml_node_iterator it = mvc_node->begin();
it != mvc_node->end(); ++it)
{
const uint cell_index = it->attribute("cell_index").as_uint();
const uint local_entity = it->attribute("local_entity").as_uint();
const uint value = it->attribute("value").as_uint();
mesh_value_collection.set_value(cell_index, local_entity, value);
}
}
else if (type == "int")
{
for (pugi::xml_node_iterator it = mvc_node->begin();
it != mvc_node->end(); ++it)
{
const uint cell_index = it->attribute("cell_index").as_uint();
const uint local_entity = it->attribute("local_entity").as_uint();
const int value = it->attribute("value").as_int();
mesh_value_collection.set_value(cell_index, local_entity, value);
}
}
else if (type == "double")
{
for (pugi::xml_node_iterator it = mvc_node->begin();
it != mvc_node->end(); ++it)
{
const uint cell_index = it->attribute("cell_index").as_uint();
const uint local_entity = it->attribute("local_entity").as_uint();
const double value = it->attribute("value").as_double();
mesh_value_collection.set_value(cell_index, local_entity, value);
}
}
else if (type == "bool")
{
for (pugi::xml_node_iterator it = mvc_node->begin();
it != mvc_node->end(); ++it)
{
const uint cell_index = it->attribute("cell_index").as_uint();
const uint local_entity = it->attribute("local_entity").as_uint();
const bool value = it->attribute("value").as_bool();
mesh_value_collection.set_value(cell_index, local_entity, value);
}
}
else
{
dolfin_error("XMLValueCollection.h",
"read mesh value collection from XML file",
"Unhandled value type \"%s\"", type.c_str());
}
}
//---------------------------------------------------------------------------
template<typename T>
void XMLMeshValueCollection::write(const MeshValueCollection<T>& mesh_value_collection,
const std::string type,
pugi::xml_node xml_node)
{
not_working_in_parallel("Writing XML MeshValueCollection");
// Add mesh function node and attributes
pugi::xml_node mf_node = xml_node.append_child("mesh_value_collection");
mf_node.append_attribute("type") = type.c_str();
mf_node.append_attribute("dim") = mesh_value_collection.dim();
mf_node.append_attribute("size") = mesh_value_collection.size();
// Add data
const std::map<std::pair<uint, uint>, T>&
values = mesh_value_collection.values();
typename std::map<std::pair<uint, uint>, T>::const_iterator it;
for (it = values.begin(); it != values.end(); ++it)
{
pugi::xml_node entity_node = mf_node.append_child("value");
entity_node.append_attribute("cell_index") = it->first.first;
entity_node.append_attribute("local_entity") = it->first.second;
entity_node.append_attribute("value") = it->second;
}
}
//---------------------------------------------------------------------------
}
#endif
|