4 #ifndef DUNE_GRID_IO_FILE_VTK_FUNCTION_HH
5 #define DUNE_GRID_IO_FILE_VTK_FUNCTION_HH
9 #include <dune/common/exceptions.hh>
10 #include <dune/common/fvector.hh>
12 #include <dune/geometry/type.hh>
13 #include <dune/geometry/referenceelements.hh>
14 #include <dune/geometry/multilineargeometry.hh>
38 template<
class Gr
idView >
44 typedef typename GridView::template Codim< 0 >::Entity
Entity;
59 const Dune::FieldVector<ctype,dim>& xi)
const = 0;
62 virtual std::string
name ()
const = 0;
92 template<
typename GV,
typename V>
128 const Dune::FieldVector<ctype,dim>&)
const override
130 return v[mapper.
index(e)*ncomps_+mycomp_];
134 std::string
name ()
const override
172 if (v.size()!=(
unsigned int)(mapper.
size()*ncomps_))
173 DUNE_THROW(IOError,
"P0VTKFunction: size mismatch");
200 template<
typename GV,
typename V>
236 const Dune::FieldVector<ctype,dim>& xi)
const override
239 const unsigned int nVertices = e.subEntities(
dim);
241 std::vector<FieldVector<ctype,1> > cornerValues(nVertices);
242 for (
unsigned i=0; i<nVertices; ++i)
243 cornerValues[i] = v[mapper.
subIndex(e,i,
dim)*ncomps_+mycomp_];
246 const MultiLinearGeometry<ctype,dim,1> interpolation(e.type(), cornerValues);
247 return interpolation.global(xi);
251 std::string
name ()
const override
289 if (v.size()!=(
unsigned int)(mapper.
size()*ncomps_))
290 DUNE_THROW(IOError,
"P1VTKFunction: size mismatch");
Mapper for multiple codim and multiple geometry types.
Common stuff for the VTKWriter.
Grid::ctype ctype
type used for coordinates in grid
Definition: common/gridview.hh:124
@ dimension
The dimension of the grid.
Definition: common/gridview.hh:127
MCMGLayout mcmgElementLayout()
layout for elements (codim-0 entities)
Definition: mcmgmapper.hh:150
MCMGLayout mcmgVertexLayout()
layout for vertices (dim-0 entities)
Definition: mcmgmapper.hh:160
Include standard header files.
Definition: agrid.hh:59
Precision
which precision to use when writing out data to vtk files
Definition: common.hh:319
@ mydimension
Dimensionality of the reference element of the entity.
Definition: common/entity.hh:113
Implementation class for a multiple codim and multiple geometry type mapper.
Definition: mcmgmapper.hh:200
Index subIndex(const typename GV::template Codim< 0 >::Entity &e, int i, unsigned int codim) const
Map subentity of codim 0 entity to starting index in array for dof block.
Definition: mcmgmapper.hh:268
size_type size() const
Return total number of entities in the entity set managed by the mapper.
Definition: mcmgmapper.hh:287
Index index(const EntityType &e) const
Map entity to starting index in array for dof block.
Definition: mcmgmapper.hh:254
A base class for grid functions with any return type and dimension.
Definition: function.hh:40
virtual double evaluate(int comp, const Entity &e, const Dune::FieldVector< ctype, dim > &xi) const =0
evaluate single component comp in the entity e at local coordinates xi
GridView::ctype ctype
Definition: function.hh:42
GridView::template Codim< 0 >::Entity Entity
Definition: function.hh:44
virtual std::string name() const =0
get name
virtual VTK::Precision precision() const
get output precision for the field
Definition: function.hh:65
@ dim
Definition: function.hh:43
virtual int ncomps() const =0
virtual ~VTKFunction()
virtual destructor
Definition: function.hh:69
Take a vector and interpret it as cell data for the VTKWriter.
Definition: function.hh:95
std::string name() const override
get name
Definition: function.hh:134
Base::ctype ctype
Definition: function.hh:117
double evaluate(int, const Entity &e, const Dune::FieldVector< ctype, dim > &) const override
evaluate
Definition: function.hh:127
P0VTKFunction(const GV &gv, const V &v_, const std::string &s_, int ncomps=1, int mycomp=0, VTK::Precision prec=VTK::Precision::float32)
construct from a vector and a name
Definition: function.hh:163
VTK::Precision precision() const override
get output precision for the field
Definition: function.hh:140
Base::Entity Entity
Definition: function.hh:116
virtual ~P0VTKFunction()
destructor
Definition: function.hh:177
int ncomps() const override
return number of components
Definition: function.hh:121
Take a vector and interpret it as point data for the VTKWriter.
Definition: function.hh:203
std::string name() const override
get name
Definition: function.hh:251
virtual ~P1VTKFunction()
destructor
Definition: function.hh:294
Base::Entity Entity
Definition: function.hh:224
Base::ctype ctype
Definition: function.hh:225
VTK::Precision precision() const override
get output precision for the field
Definition: function.hh:257
double evaluate(int comp, const Entity &e, const Dune::FieldVector< ctype, dim > &xi) const override
evaluate
Definition: function.hh:235
P1VTKFunction(const GV &gv, const V &v_, const std::string &s_, int ncomps=1, int mycomp=0, VTK::Precision prec=VTK::Precision::float32)
construct from a vector and a name
Definition: function.hh:280
int ncomps() const override
return number of components
Definition: function.hh:229