/usr/include/dune/grid/albertagrid/albertareader.hh is in libdune-grid-dev 2.2.1-2.
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 | #ifndef DUNE_ALBERTA_ALBERTAREADER_HH
#define DUNE_ALBERTA_ALBERTAREADER_HH
#include <dune/grid/common/grid.hh>
#include <dune/grid/common/gridfactory.hh>
#include <dune/grid/utility/grapedataioformattypes.hh>
#include <dune/grid/albertagrid/macrodata.hh>
#if HAVE_ALBERTA
namespace Dune
{
template< class Grid >
class AlbertaReader
{
typedef AlbertaReader< Grid > This;
public:
typedef Dune::GridFactory< Grid > GridFactory;
typedef typename Grid::ctype ctype;
static const int dimension = Grid::dimension;
static const int dimensionworld = Grid::dimensionworld;
private:
dune_static_assert( dimensionworld == Alberta::dimWorld,
"AlbertaReader: world dimension must match ALBERTA's world dimension." );
typedef Alberta::MacroData< dimension > MacroData;
MacroData macroData_;
AlbertaReader ( const This & );
This &operator= ( const This & );
public:
AlbertaReader ()
{}
template< GrapeIOFileFormatType type >
void readGrid ( const std::string &fileName, GridFactory &factory )
{
dune_static_assert( type != pgm, "AlbertaReader: reading pgm format is not supported." );
// read ALBERTA macro triangulation
macroData_.read( fileName, (type == xdr) );
// insert all vertices into the factory
const int numVertices = macroData_.vertexCount();
for( int i = 0; i < numVertices; ++i )
{
FieldVector< ctype, dimensionworld > v;
const Alberta::GlobalVector &coords = macroData_.vertex( i );
for( int j = 0; j < dimensionworld; ++j )
v[ j ] = coords[ j ];
factory.insertVertex( v );
}
// insert all elements into the factory
std::vector< unsigned int > vertices( dimension+1 );
const int numElements = macroData_.elementCount();
for( int i = 0; i < numElements; ++i )
{
const typename MacroData::ElementId &id = macroData_.element( i );
for( int j = 0; j <= dimension; ++j )
vertices[ j ] = id[ j ];
typedef typename GenericGeometry::SimplexTopology< dimension >::type Topology;
factory.insertElement( GeometryType( Topology() ), vertices );
}
// release ALBERTA macro data
macroData_.release();
}
void readGrid ( const std::string &filename, GridFactory &factory )
{
readGrid< ascii >( filename, factory );
}
};
}
#endif // #if HAVE_ALBERTA
#endif
|