5#ifndef DUNE_GRID_IO_FILE_GMSHWRITER_HH
6#define DUNE_GRID_IO_FILE_GMSHWRITER_HH
14#include <dune/common/exceptions.hh>
15#include <dune/geometry/type.hh>
16#include <dune/geometry/referenceelements.hh>
35 template <
class Gr
idView>
44 static_assert( (dimWorld <= 3),
"GmshWriter requires dimWorld <= 3." );
47 template<
typename Entity>
48 std::size_t nodeIndexFromEntity(
const Entity& entity,
int i)
const {
49 return gv.
indexSet().subIndex(entity, i, dim)+1;
55 static std::size_t translateDuneToGmshType(
const GeometryType& type) {
56 std::size_t element_type;
60 else if (type.isTriangle())
62 else if (type.isQuadrilateral())
64 else if (type.isTetrahedron())
66 else if (type.isHexahedron())
68 else if (type.isPrism())
70 else if (type.isPyramid())
72 else if (type.isVertex())
75 DUNE_THROW(Dune::IOError,
"GeometryType " << type <<
" is not supported by gmsh.");
94 void outputElements(std::ofstream& file,
const std::vector<int>& physicalEntities,
const std::vector<int>& physicalBoundaries)
const {
96 std::size_t counter(1);
97 for (
const auto& entity : elements(gv)) {
100 std::size_t element_type = translateDuneToGmshType(entity.
type());
102 file << counter <<
" " << element_type;
104 if (!physicalEntities.empty())
105 file <<
" " << 1 <<
" " << physicalEntities[elementMapper.
index(entity)];
111 if (3 == element_type)
113 << nodeIndexFromEntity(entity, 0) <<
" " << nodeIndexFromEntity(entity, 1) <<
" "
114 << nodeIndexFromEntity(entity, 3) <<
" " << nodeIndexFromEntity(entity, 2);
115 else if (5 == element_type)
117 << nodeIndexFromEntity(entity, 0) <<
" " << nodeIndexFromEntity(entity, 1) <<
" "
118 << nodeIndexFromEntity(entity, 3) <<
" " << nodeIndexFromEntity(entity, 2) <<
" "
119 << nodeIndexFromEntity(entity, 4) <<
" " << nodeIndexFromEntity(entity, 5) <<
" "
120 << nodeIndexFromEntity(entity, 7) <<
" " << nodeIndexFromEntity(entity, 6);
121 else if (7 == element_type)
123 << nodeIndexFromEntity(entity, 0) <<
" " << nodeIndexFromEntity(entity, 1) <<
" "
124 << nodeIndexFromEntity(entity, 3) <<
" " << nodeIndexFromEntity(entity, 2) <<
" "
125 << nodeIndexFromEntity(entity, 4);
127 for (
int k = 0; k < entity.
geometry().corners(); ++k)
128 file <<
" " << nodeIndexFromEntity(entity, k);
135 if (!physicalBoundaries.empty()) {
136 auto refElement = referenceElement<typename GridView::ctype,dim>(entity.
type());
137 for(
const auto& intersection : intersections(gv, entity)) {
138 if(intersection.boundary()) {
139 const auto faceLocalIndex(intersection.indexInInside());
140 file << counter <<
" " << translateDuneToGmshType(intersection.type())
141 <<
" " << 1 <<
" " << physicalBoundaries[intersection.boundarySegmentIndex()];
142 for (
int k = 0; k < intersection.geometry().corners(); ++k)
144 const auto vtxLocalIndex(refElement.subEntity(faceLocalIndex, 1, k, dim));
145 file <<
" " << nodeIndexFromEntity(entity, vtxLocalIndex);
153 }
catch(Exception& e) {
167 void outputNodes(std::ofstream& file)
const {
168 for (
const auto&
vertex : vertices(gv)) {
169 const auto globalCoord =
vertex.geometry().center();
173 file << nodeIndex <<
" " << globalCoord[0] <<
" " << 0 <<
" " << 0 << std::endl;
174 else if (2 == dimWorld)
175 file << nodeIndex <<
" " << globalCoord[0] <<
" " << globalCoord[1] <<
" " << 0 << std::endl;
177 file << nodeIndex <<
" " << globalCoord[0] <<
" " << globalCoord[1] <<
" " << globalCoord[2] << std::endl;
194 precision = numDigits;
218 void write(
const std::string& fileName,
219 const std::vector<int>& physicalEntities=std::vector<int>(),
220 const std::vector<int>& physicalBoundaries=std::vector<int>())
const {
222 std::ofstream file(fileName.c_str());
224 DUNE_THROW(Dune::IOError,
"Could not open " << fileName <<
" with write access.");
227 file << std::setprecision( precision );
230 file <<
"$MeshFormat" << std::endl
231 <<
"2.0 0 " <<
sizeof(double) << std::endl
232 <<
"$EndMeshFormat" << std::endl;
235 file <<
"$Nodes" << std::endl
236 << gv.
size(dim) << std::endl;
240 file <<
"$EndNodes" << std::endl;
243 int boundariesSize(0);
244 if(!physicalBoundaries.empty())
245 for(
const auto& entity : elements(gv))
246 for(
const auto& intersection : intersections(gv, entity))
247 if(intersection.boundary())
250 file <<
"$Elements" << std::endl
251 << gv.
size(0) + boundariesSize<< std::endl;
253 outputElements(file, physicalEntities, physicalBoundaries);
255 file <<
"$EndElements" << std::endl;
Mapper for multiple codim and multiple geometry types.
const IndexSet & indexSet() const
obtain the index set
Definition: common/gridview.hh:191
int size(int codim) const
obtain number of entities in a given codimension
Definition: common/gridview.hh:197
static constexpr int dimension
The dimension of the grid.
Definition: common/gridview.hh:148
static constexpr int dimensionworld
The dimension of the world the grid lives in.
Definition: common/gridview.hh:151
MCMGLayout mcmgElementLayout()
layout for elements (codim-0 entities)
Definition: mcmgmapper.hh:97
Include standard header files.
Definition: agrid.hh:60
GeometryType
Type representing VTK's entity geometry types.
Definition: common.hh:132
@ vertex
Definition: common.hh:133
Wrapper class for entities.
Definition: common/entity.hh:66
Geometry geometry() const
obtain geometric realization of the entity
Definition: common/entity.hh:141
GeometryType type() const
Return the name of the reference element. The type can be used to access the Dune::ReferenceElement.
Definition: common/entity.hh:146
Grid view abstract base class.
Definition: common/gridview.hh:66
Implementation class for a multiple codim and multiple geometry type mapper.
Definition: mcmgmapper.hh:129
Index index(const EntityType &e) const
Map entity to starting index in array for dof block.
Definition: mcmgmapper.hh:171
Write Gmsh mesh file.
Definition: gmshwriter.hh:37
GmshWriter(const GridView &gridView, int numDigits=6)
Constructor expecting GridView of Grid to be written.
Definition: gmshwriter.hh:187
void setPrecision(int numDigits)
Set the number of digits to be used when writing the vertices. By default is 6.
Definition: gmshwriter.hh:193
void write(const std::string &fileName, const std::vector< int > &physicalEntities=std::vector< int >(), const std::vector< int > &physicalBoundaries=std::vector< int >()) const
Write given grid in Gmsh 2.0 compatible ASCII file.
Definition: gmshwriter.hh:218
Different resources needed by all grid implementations.