4 #ifndef DUNE_GRID_UTILITY_GRIDINFO_HH 5 #define DUNE_GRID_UTILITY_GRIDINFO_HH 16 #include <dune/common/classname.hh> 17 #include <dune/common/exceptions.hh> 18 #include <dune/common/fvector.hh> 19 #include <dune/common/hybridutilities.hh> 20 #include <dune/common/std/utility.hh> 22 #include <dune/geometry/multilineargeometry.hh> 23 #include <dune/geometry/referenceelements.hh> 24 #include <dune/geometry/type.hh> 64 count(0), volumeMin(
std::numeric_limits<ctype>::infinity()),
65 volumeMax(-
std::numeric_limits<ctype>::infinity()), volumeSum(0)
77 public std::binary_function<GeometryType, GeometryType, bool>
82 return a.dim() < b.dim() ||
83 (a.dim() == b.dim() && (a.isNone() < b.isNone() ||
84 (a.isNone() == b.isNone() && (a.id() >> 1) < (b.id() >> 1))));
98 public std::map<GeometryType, EntityInfo<ctype>, GridViewInfoGTCompare>
125 void print(std::ostream &stream, std::string prefix)
const {
126 if(!gridName.empty()) {
127 stream << prefix << gridName <<
":\n";
130 if(!gridViewName.empty()) {
131 stream << prefix << gridViewName <<
":\n";
134 if(!partitionName.empty()) {
135 stream << prefix << partitionName <<
":\n";
139 typedef typename GridViewInfo::const_iterator Iterator;
140 std::size_t dim = ~0;
141 const Iterator &end = this->end();
142 for(Iterator it = this->begin(); it != end; ++it) {
143 if(it->first.dim() != dim) {
144 dim = it->first.dim();
145 stream << prefix <<
"Dim = " << dim <<
":\n";
147 stream << prefix <<
" " << it->first <<
": Count = " 148 << it->second.count <<
", Volume range = " 149 <<
"(" << it->second.volumeMin <<
".." 150 << it->second.volumeMax <<
"), Total volume = " 151 << it->second.volumeSum <<
"\n";
162 template<
class ctype>
166 info.
print(stream,
"");
173 struct FillGridInfoOperation {
174 template<
class Entity,
class Mapper,
class Visited,
class RefElem>
175 static void apply(
const Entity &e,
const Mapper &mapper, Visited &visited,
180 typedef typename Entity::Geometry::ctype ctype;
181 static const std::size_t dimw = Entity::Geometry::coorddimension;
183 std::vector<FieldVector<ctype, dimw> > coords;
184 for(
int i = 0; i < refelem.size(codim); ++i) {
185 int index = mapper.map(e, i, codim);
188 visited[index] =
true;
192 coords.resize( refelem.size(i, codim, dim) );
193 for(std::size_t corner = 0; corner < coords.size(); ++corner)
194 coords[ corner ] = geo.corner( refelem.subEntity( i, codim, corner, dim ) );
195 MultiLinearGeometry<ctype, dim-codim, dimw> mygeo(gt, coords);
197 ctype volume = mygeo.volume();
206 template<
int dimgr
id>
207 struct MCMGNonElementLayout {
208 bool contains(
GeometryType gt)
const {
return gt.dim() < dimgrid; }
221 typedef typename GV::ctype ctype;
222 static const std::size_t dim = GV::dimension;
223 typedef typename GV::template Codim<0>::Iterator EIterator;
224 typedef typename GV::template Codim<0>::Geometry EGeometry;
225 typedef typename GV::IndexSet
IndexSet;
229 typedef ReferenceElements<ctype, dim> RefElems;
232 std::vector<bool> visited(mapper.
size(),
false);
234 gridViewInfo.
gridName = className<typename GV::Grid>();
237 gridViewInfo.clear();
239 const EIterator &eend = gv.template end<0>();
240 for(EIterator eit = gv.template begin<0>(); eit != eend; ++eit) {
241 ctype volume = eit->geometry().volume();
247 if(!eit->type().isNone()) {
248 const EGeometry &geo = eit->geometry();
249 Hybrid::forEach(Std::make_index_sequence< dim >{},
250 [ & ](
auto i){ FillGridInfoOperation< i+1 >::apply(*eit, mapper, visited, geo, RefElems::general(eit->type()), gridViewInfo); } );
256 if(gridViewInfo.count(gt) > 0) {
257 for(std::size_t codim = 0; codim < dim; ++codim) {
258 gt.makeNone(dim-codim);
261 std::numeric_limits<ctype>::quiet_NaN();
268 const InfoIterator &end = gridViewInfo.end();
269 const IndexSet &is = gv.indexSet();
270 for(InfoIterator it = gridViewInfo.begin(); it != end; ++it) {
271 it->second.count = is.
size(it->first);
272 if(it->second.count == 0)
273 DUNE_THROW(Exception,
"Found Entities of geomentry type " <<
274 it->first <<
" while iterating through the grid, but " 275 "indexSet.size() == 0 for that geometry type");
283 #endif // DUNE_GRID_UTILITY_GRIDINFO_HH std::string partitionName
name of the partition this information was extracted from
Definition: utility/gridinfo.hh:109
std::string gridViewName
name of the class of the GridView this information was extracted from
Definition: utility/gridinfo.hh:103
ctype volumeMin
minimum volume of all entities in the set.
Definition: utility/gridinfo.hh:41
Index Set Interface base class.
Definition: common/grid.hh:346
std::string gridName
name of the grid class this information was extracted from
Definition: utility/gridinfo.hh:101
structure to hold information about a certain GridView.
Definition: utility/gridinfo.hh:97
GridImp::template Codim< cd >::Geometry Geometry
The corresponding geometry type.
Definition: common/entity.hh:98
Implementation class for a multiple codim and multiple geometry type mapper.
Definition: mcmgmapper.hh:198
IndexType size(GeometryType type) const
Return total number of entities of given geometry type in entity set .
Definition: indexidset.hh:220
Wrapper class for entities.
Definition: common/entity.hh:63
std::ostream & operator<<(std::ostream &stream, const GridViewInfo< ctype > &info)
write a GridViewInfo object
Definition: utility/gridinfo.hh:163
void fillGridViewInfoSerial(const GV &gv, GridViewInfo< typename GV::ctype > &gridViewInfo)
fill a GridViewInfo structure from a serial grid
Definition: utility/gridinfo.hh:218
bool operator()(const GeometryType &a, const GeometryType &b) const
compare two GeometryTypes
Definition: utility/gridinfo.hh:80
Comparison object to sort GeometryType by majorly dimension.
Definition: utility/gridinfo.hh:76
int max(const DofVectorPointer< int > &dofVector)
Definition: dofvector.hh:335
ctype volumeSum
sum of volumes of all entities in the set.
Definition: utility/gridinfo.hh:55
Know the grid dimension.
Definition: common/entity.hh:109
GeometryType
Type representing VTK's entity geometry types.
Definition: common.hh:178
Mapper for multiple codim and multiple geometry types.
int min(const DofVectorPointer< int > &dofVector)
Definition: dofvector.hh:346
Structure to hold statistical information about one type of entity.
Definition: utility/gridinfo.hh:32
ctype volumeMax
maximum volume of all entities in the set.
Definition: utility/gridinfo.hh:48
std::size_t count
number of entities in the set
Definition: utility/gridinfo.hh:34
void print(std::ostream &stream, std::string prefix) const
print the information contained in this object
Definition: utility/gridinfo.hh:125
Mapper interface.
Definition: mapper.hh:107
size_type size() const
Return total number of entities in the entity set managed by the mapper.
Definition: mcmgmapper.hh:287
Include standard header files.
Definition: agrid.hh:58
EntityInfo()
initialize the structure
Definition: utility/gridinfo.hh:63