Dune Core Modules (2.4.1)

albertareader.hh
1// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
2// vi: set et ts=4 sw=2 sts=2:
3#ifndef DUNE_ALBERTA_ALBERTAREADER_HH
4#define DUNE_ALBERTA_ALBERTAREADER_HH
5
8
9#include <dune/grid/utility/grapedataioformattypes.hh>
10
12
13#if HAVE_ALBERTA
14
15namespace Dune
16{
17
18 template< class Grid >
19 class AlbertaReader
20 {
21 typedef AlbertaReader< Grid > This;
22
23 public:
24 typedef Dune::GridFactory< Grid > GridFactory;
25
26 typedef typename Grid::ctype ctype;
27
28 static const int dimension = Grid::dimension;
29 static const int dimensionworld = Grid::dimensionworld;
30
31 private:
32 static_assert(dimensionworld == Alberta::dimWorld,
33 "AlbertaReader: world dimension must match ALBERTA's world dimension.");
34
35 typedef Alberta::MacroData< dimension > MacroData;
36
37 MacroData macroData_;
38
39 AlbertaReader ( const This & );
40 This &operator= ( const This & );
41
42 public:
43 AlbertaReader ()
44 {}
45
46 template< GrapeIOFileFormatType type >
47 void readGrid ( const std::string &fileName, GridFactory &factory )
48 {
49 static_assert(type != pgm, "AlbertaReader: reading pgm format is not supported.");
50
51 // read ALBERTA macro triangulation
52 macroData_.read( fileName, (type == xdr) );
53
54 // insert all vertices into the factory
55 const int numVertices = macroData_.vertexCount();
56 for( int i = 0; i < numVertices; ++i )
57 {
58 FieldVector< ctype, dimensionworld > v;
59 const Alberta::GlobalVector &coords = macroData_.vertex( i );
60 for( int j = 0; j < dimensionworld; ++j )
61 v[ j ] = coords[ j ];
62 factory.insertVertex( v );
63 }
64
65 // insert all elements into the factory
66 std::vector< unsigned int > vertices( dimension+1 );
67 const int numElements = macroData_.elementCount();
68 for( int i = 0; i < numElements; ++i )
69 {
70 const typename MacroData::ElementId &id = macroData_.element( i );
71 for( int j = 0; j <= dimension; ++j )
72 vertices[ j ] = id[ j ];
73 typedef typename GenericGeometry::SimplexTopology< dimension >::type Topology;
74 factory.insertElement( GeometryType( Topology() ), vertices );
75 }
76
77 // release ALBERTA macro data
78 macroData_.release();
79 }
80
81 void readGrid ( const std::string &filename, GridFactory &factory )
82 {
83 readGrid< ascii >( filename, factory );
84 }
85 };
86
87}
88
89#endif // #if HAVE_ALBERTA
90
91#endif
@ dimensionworld
The dimension of the world the grid lives in.
Definition: grid.hh:408
@ dimension
The dimension of the grid.
Definition: grid.hh:402
ct ctype
Define type used for coordinates in grid module.
Definition: grid.hh:548
Different resources needed by all grid implementations.
Provide a generic factory class for unstructured grids.
GeometryType
Type representing VTK's entity geometry types.
Definition: common.hh:178
provides a wrapper for ALBERTA's macro_data structure
Dune namespace.
Definition: alignment.hh:10
@ xdr
Definition: grapedataioformattypes.hh:16
Creative Commons License   |  Legal Statements / Impressum  |  Hosted by TU Dresden  |  generated with Hugo v0.111.3 (Nov 12, 23:30, 2024)