3 #ifndef DUNE_ALBERTA_ALBERTAREADER_HH
4 #define DUNE_ALBERTA_ALBERTAREADER_HH
18 template<
class Gr
id >
33 "AlbertaReader: world dimension must match ALBERTA's world dimension." );
46 template< GrapeIOFileFormatType type >
49 dune_static_assert( type !=
pgm,
"AlbertaReader: reading pgm format is not supported." );
52 macroData_.
read( fileName, (type ==
xdr) );
56 for(
int i = 0; i < numVertices; ++i )
58 FieldVector< ctype, dimensionworld > v;
66 std::vector< unsigned int > vertices(
dimension+1 );
68 for(
int i = 0; i < numElements; ++i )
72 vertices[ j ] =
id[ j ];
73 typedef typename GenericGeometry::SimplexTopology< dimension >::type Topology;
83 readGrid< ascii >( filename, factory );
89 #endif // #if HAVE_ALBERTA