1 #ifndef DUNE_ALBERTA_ALBERTAREADER_HH
2 #define DUNE_ALBERTA_ALBERTAREADER_HH
16 template<
class Gr
id >
31 "AlbertaReader: world dimension must match ALBERTA's world dimension." );
44 template< GrapeIOFileFormatType type >
47 dune_static_assert( type !=
pgm,
"AlbertaReader: reading pgm format is not supported." );
50 macroData_.
read( fileName, (type ==
xdr) );
54 for(
int i = 0; i < numVertices; ++i )
56 FieldVector< ctype, dimensionworld > v;
64 std::vector< unsigned int > vertices(
dimension+1 );
66 for(
int i = 0; i < numElements; ++i )
70 vertices[ j ] =
id[ j ];
71 typedef typename GenericGeometry::SimplexTopology< dimension >::type Topology;
81 readGrid< ascii >( filename, factory );
87 #endif // #if HAVE_ALBERTA