3#ifndef DUNE_GRID_IO_FILE_DGFPARSER_DGFUG_HH
4#define DUNE_GRID_IO_FILE_DGFPARSER_DGFUG_HH
13#include <dune/common/exceptions.hh>
14#include <dune/common/fvector.hh>
15#include <dune/common/parallel/mpihelper.hh>
79 struct DGFGridFactory< UGGrid< dim > >
82 typedef UGGrid< dim >
Grid;
93 dgf_( rank( comm ), size( comm ) )
103 dgf_( rank( comm ), size( comm ) )
105 std::ifstream input( filename.c_str() );
107 DUNE_THROW( DGFException,
"Error: Macrofile " << filename <<
" not found" );
118 template<
class GG,
class II >
121 return factory_.wasInserted( intersection );
125 template<
class GG,
class II >
132 template<
int codim >
136 return dgf_.nofelparams;
138 return dgf_.nofvtxparams;
144 template<
class Entity >
147 return numParameters< Entity::codimension >();
151 std::vector< double > &
parameter (
const typename Grid::template Codim< 0 >::Entity &element )
153 if( numParameters< 0 >() <= 0 )
155 DUNE_THROW( InvalidStateException,
156 "Calling DGFGridFactory::parameter is only allowed if there are parameters." );
158 return dgf_.elParams[ factory_.insertionIndex( element ) ];
162 std::vector< double > &
parameter (
const typename Grid::template Codim< dimension >::Entity &vertex )
164 if( numParameters< dimension >() <= 0 )
166 DUNE_THROW( InvalidStateException,
167 "Calling DGFGridFactory::parameter is only allowed if there are parameters." );
169 return dgf_.vtxParams[ factory_.insertionIndex( vertex ) ];
175 return dgf_.haveBndParameters;
179 template<
class GG,
class II >
186 const ReferenceElement< double, dimension > &refElem
187 = ReferenceElements< double, dimension >::general( entity.type() );
188 int corners = refElem.size( face, 1,
dimension );
189 std::vector< unsigned int > bound( corners );
190 for(
int i = 0; i < corners; ++i )
192 const int k = refElem.subEntity( face, 1, i,
dimension );
193 bound[ i ] = factory_.insertionIndex( entity.template subEntity< dimension >( k ) );
196 DuneGridFormatParser::facemap_t::key_type key( bound,
false );
197 const DuneGridFormatParser::facemap_t::const_iterator pos = dgf_.facemap.find( key );
198 if( pos != dgf_.facemap.end() )
199 return dgf_.facemap.find( key )->second.second;
206 void generate ( std::istream &input );
213 MPI_Comm_rank( MPICOMM, &rank );
223 MPI_Comm_size( MPICOMM, &size );
229 GridFactory< UGGrid< dim > > factory_;
230 DuneGridFormatParser dgf_;
Include standard header files.
Definition agrid.hh:60
int numParameters() const
Definition dgfgridfactory.hh:113
int boundaryId(const Intersection &intersection) const
Definition dgfgridfactory.hh:107
const DGFBoundaryParameter::type & boundaryParameter(const Intersection< GG, II > &intersection) const
Definition dgfgridfactory.hh:163
MPIHelper::MPICommunicator MPICommunicatorType
Definition dgfgridfactory.hh:39
bool wasInserted(const Intersection &intersection) const
Definition dgfgridfactory.hh:101
G Grid
Definition dgfgridfactory.hh:37
static const int dimension
Definition dgfgridfactory.hh:38
Grid * grid()
Definition dgfgridfactory.hh:95
std::vector< double > & parameter(const Element &element)
Definition dgfgridfactory.hh:129
bool haveBoundaryParameters() const
Definition dgfgridfactory.hh:156
DGFGridFactory(std::istream &input, MPICommunicatorType comm=MPIHelper::getCommunicator())
Definition dgfgridfactory.hh:47
Intersection of a mesh entity of codimension 0 ("element") with a "neighboring" element or with the d...
Definition common/intersection.hh:162
int indexInInside() const
Local index of codim 1 entity in the inside() entity where intersection is contained in.
Definition common/intersection.hh:374
size_t boundarySegmentIndex() const
index of the boundary segment within the macro grid
Definition common/intersection.hh:264
Entity inside() const
return Entity on the inside of this intersection. That is the Entity where we started this.
Definition common/intersection.hh:278
GridImp::template Codim< 0 >::Entity Entity
Type of entity that this Intersection belongs to.
Definition common/intersection.hh:186
Common Grid parameters.
Definition gridparameter.hh:33
Some simple static information for a given GridType.
Definition io/file/dgfparser/dgfparser.hh:55
static double refineWeight()
static int refineStepsForHalf()
number of globalRefine steps needed to refuce h by 0.5
size_t heapSize() const
returns heap size used on construction of the grid
Definition dgfug.hh:46
bool noCopy_
Definition dgfug.hh:50
size_t heapSize_
Definition dgfug.hh:51
bool noClosure_
Definition dgfug.hh:49
bool noCopy() const
returns true if no copies are made for UGGrid elements
Definition dgfug.hh:44
bool noClosure() const
returns true if no closure should be used for UGGrid
Definition dgfug.hh:42
static const type & defaultValue()
default constructor
Definition parser.hh:26
std::string type
type of additional boundary parameters
Definition parser.hh:23
Front-end for the grid manager of the finite element toolbox UG.
Definition uggrid.hh:231