Page MenuHomec4science
Files F63403350

dumper_iohelper_tmpl_quadrature_points_field.hh
No OneTemporary
Actions

File Metadata

Created
Sun, May 19, 21:30

dumper_iohelper_tmpl_quadrature_points_field.hh
View Options

/**
* @file dumper_iohelper_tmpl_quadrature_points_field.hh
* @author Nicolas Richart <nicolas.richart@epfl.ch>
* @date Thu Oct 25 14:40:14 2012
*
* @brief Description of quadrature points fields
*
* @section LICENSE
*
* Copyright (©) 2010-2011 EPFL (Ecole Polytechnique Fédérale de Lausanne)
* Laboratory (LSMS - Laboratoire de Simulation en Mécanique des Solides)
*
* Akantu is free software: you can redistribute it and/or modify it under the
* terms of the GNU Lesser General Public License as published by the Free
* Software Foundation, either version 3 of the License, or (at your option) any
* later version.
*
* Akantu is distributed in the hope that it will be useful, but WITHOUT ANY
* WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
* A PARTICULAR PURPOSE. See the GNU Lesser General Public License for more
* details.
*
* You should have received a copy of the GNU Lesser General Public License
* along with Akantu. If not, see <http://www.gnu.org/licenses/>.
*
*/
/* -------------------------------------------------------------------------- */
__END_AKANTU__
#include "fem.hh"
__BEGIN_AKANTU__
/* -------------------------------------------------------------------------- */
template<typename i_type, typename d_type,
template<typename> class ret_type, class daughter>
class DumperIOHelper::quadrature_point_iterator : public element_iterator<i_type, d_type,
ret_type, daughter> {
public:
typedef element_iterator<i_type, d_type, ret_type, daughter> parent;
typedef typename parent::it_type it_type;
typedef typename parent::data_type data_type;
typedef typename parent::return_type return_type;
typedef typename parent::field_type field_type;
typedef typename parent::internal_iterator internal_iterator;
public:
quadrature_point_iterator(const field_type & field,
UInt n,
const typename field_type::type_iterator & t_it,
const typename field_type::type_iterator & t_it_end,
const internal_iterator & it,
ElementType element_type,
const GhostType ghost_type = _not_ghost) :
parent(field, n, t_it, t_it_end, it, element_type, ghost_type) { }
void setFEM(const FEM & fem) { this->fem = &fem; }
protected:
virtual UInt getNbDataPerElem(const ElementType & type) { return fem->getNbQuadraturePoints(type); }
protected:
const FEM * fem;
};
/* -------------------------------------------------------------------------- */
/* Fields type description */
/* -------------------------------------------------------------------------- */
template<typename T, class iterator_type, template<class> class ret_type>
class DumperIOHelper::GenericQuadraturePointsField : public GenericElementalField<T, iterator_type, ret_type> {
typedef GenericElementalField<T, iterator_type, ret_type> parent;
public:
GenericQuadraturePointsField(const FEM & fem,
const ByElementTypeVector<T> & field,
UInt spatial_dimension = 0,
GhostType ghost_type = _not_ghost,
ElementKind element_kind = _ek_not_defined) :
parent(field, spatial_dimension,
ghost_type, element_kind), fem(fem) { }
GenericQuadraturePointsField(const FEM & fem,
const ByElementTypeVector<T> & field,
UInt n,
UInt spatial_dimension = 0,
GhostType ghost_type = _not_ghost,
ElementKind element_kind = _ek_not_defined) :
parent(field, n, spatial_dimension,
ghost_type, element_kind), fem(fem) { }
typedef iterator_type iterator;
/* ------------------------------------------------------------------------ */
virtual iterator begin() {
iterator it = parent::begin();
it.setFEM(fem);
return it;
}
virtual iterator end () {
iterator it = parent::end();
it.setFEM(fem);
return it;
}
const FEM & getFEM() { return fem; }
protected:
virtual UInt getNbDataPerElem(const ElementType & type) { return fem.getNbQuadraturePoints(type); }
protected:
const FEM & fem;
};

Event Timeline

c4science · Help