*/
virtual void
fill_fe_values (const Mapping<dim> &mapping,
- const DoFHandler<dim>::cell_iterator &cell,
+ const typename DoFHandler<dim>::cell_iterator &cell,
const Quadrature<dim> &quadrature,
- Mapping<dim>::InternalDataBase &mapping_internal,
- Mapping<dim>::InternalDataBase &fe_internal,
+ typename Mapping<dim>::InternalDataBase &mapping_internal,
+ typename Mapping<dim>::InternalDataBase &fe_internal,
FEValuesData<dim> &data) const = 0;
/**
*/
virtual void
fill_fe_face_values (const Mapping<dim> &mapping,
- const DoFHandler<dim>::cell_iterator &cell,
+ const typename DoFHandler<dim>::cell_iterator &cell,
const unsigned int face_no,
const Quadrature<dim-1> &quadrature,
- Mapping<dim>::InternalDataBase &mapping_internal,
- Mapping<dim>::InternalDataBase &fe_internal,
+ typename Mapping<dim>::InternalDataBase &mapping_internal,
+ typename Mapping<dim>::InternalDataBase &fe_internal,
FEValuesData<dim> &data) const = 0;
/**
*/
virtual void
fill_fe_subface_values (const Mapping<dim> &mapping,
- const DoFHandler<dim>::cell_iterator &cell,
+ const typename DoFHandler<dim>::cell_iterator &cell,
const unsigned int face_no,
const unsigned int sub_no,
const Quadrature<dim-1> &quadrature,
* gradients.
*/
void compute_2nd (const Mapping<dim> &mapping,
- const DoFHandler<dim>::cell_iterator &cell,
+ const typename DoFHandler<dim>::cell_iterator &cell,
const unsigned int offset,
typename Mapping<dim>::InternalDataBase &mapping_internal,
InternalDataBase &fe_internal,
*/
virtual void
fill_fe_values (const Mapping<dim> &mapping,
- const DoFHandler<dim>::cell_iterator &cell,
+ const typename DoFHandler<dim>::cell_iterator &cell,
const Quadrature<dim> &quadrature,
- Mapping<dim>::InternalDataBase &mapping_internal,
- Mapping<dim>::InternalDataBase &fe_internal,
+ typename Mapping<dim>::InternalDataBase &mapping_internal,
+ typename Mapping<dim>::InternalDataBase &fe_internal,
FEValuesData<dim>& data) const;
/**
*/
virtual void
fill_fe_face_values (const Mapping<dim> &mapping,
- const DoFHandler<dim>::cell_iterator &cell,
+ const typename DoFHandler<dim>::cell_iterator &cell,
const unsigned int face_no,
const Quadrature<dim-1> &quadrature,
- Mapping<dim>::InternalDataBase &mapping_internal,
- Mapping<dim>::InternalDataBase &fe_internal,
+ typename Mapping<dim>::InternalDataBase &mapping_internal,
+ typename Mapping<dim>::InternalDataBase &fe_internal,
FEValuesData<dim>& data) const ;
/**
*/
virtual void
fill_fe_subface_values (const Mapping<dim> &mapping,
- const DoFHandler<dim>::cell_iterator &cell,
+ const typename DoFHandler<dim>::cell_iterator &cell,
const unsigned int face_no,
const unsigned int sub_no,
const Quadrature<dim-1> &quadrature,
- Mapping<dim>::InternalDataBase &mapping_internal,
- Mapping<dim>::InternalDataBase &fe_internal,
+ typename Mapping<dim>::InternalDataBase &mapping_internal,
+ typename Mapping<dim>::InternalDataBase &fe_internal,
FEValuesData<dim>& data) const ;
private:
*/
virtual void
fill_fe_values (const Mapping<dim> &mapping,
- const DoFHandler<dim>::cell_iterator &cell,
+ const typename DoFHandler<dim>::cell_iterator &cell,
const Quadrature<dim> &quadrature,
- Mapping<dim>::InternalDataBase &mapping_internal,
- Mapping<dim>::InternalDataBase &fe_internal,
+ typename Mapping<dim>::InternalDataBase &mapping_internal,
+ typename Mapping<dim>::InternalDataBase &fe_internal,
FEValuesData<dim>& data) const;
/**
*/
virtual void
fill_fe_face_values (const Mapping<dim> &mapping,
- const DoFHandler<dim>::cell_iterator &cell,
+ const typename DoFHandler<dim>::cell_iterator &cell,
const unsigned int face_no,
const Quadrature<dim-1> &quadrature,
- Mapping<dim>::InternalDataBase &mapping_internal,
- Mapping<dim>::InternalDataBase &fe_internal,
+ typename Mapping<dim>::InternalDataBase &mapping_internal,
+ typename Mapping<dim>::InternalDataBase &fe_internal,
FEValuesData<dim>& data) const ;
/**
*/
virtual void
fill_fe_subface_values (const Mapping<dim> &mapping,
- const DoFHandler<dim>::cell_iterator &cell,
+ const typename DoFHandler<dim>::cell_iterator &cell,
const unsigned int face_no,
const unsigned int sub_no,
const Quadrature<dim-1> &quadrature,
- Mapping<dim>::InternalDataBase &mapping_internal,
- Mapping<dim>::InternalDataBase &fe_internal,
+ typename Mapping<dim>::InternalDataBase &mapping_internal,
+ typename Mapping<dim>::InternalDataBase &fe_internal,
FEValuesData<dim>& data) const ;
private:
* to the @p{get_function_*} functions
* safer.
*/
- DoFHandler<dim>::cell_iterator present_cell;
+ typename DoFHandler<dim>::cell_iterator present_cell;
/**
* Storage for the mapping object.
* @ref{Mapping}.
*/
virtual void
- fill_fe_values (const DoFHandler<dim>::cell_iterator &cell,
+ fill_fe_values (const typename DoFHandler<dim>::cell_iterator &cell,
const Quadrature<dim>& quadrature,
- Mapping<dim>::InternalDataBase &mapping_data,
+ typename Mapping<dim>::InternalDataBase &mapping_data,
std::vector<Point<dim> > &quadrature_points,
std::vector<double> &JxW_values) const ;
* @ref{Mapping}.
*/
virtual void
- fill_fe_values (const DoFHandler<dim>::cell_iterator &cell,
+ fill_fe_values (const typename DoFHandler<dim>::cell_iterator &cell,
const Quadrature<dim> &quadrature,
- Mapping<dim>::InternalDataBase &mapping_data,
+ typename Mapping<dim>::InternalDataBase &mapping_data,
std::vector<Point<dim> > &quadrature_points,
std::vector<double> &JxW_values) const ;
* @ref{Mapping}.
*/
virtual void
- fill_fe_face_values (const DoFHandler<dim>::cell_iterator &cell,
+ fill_fe_face_values (const typename DoFHandler<dim>::cell_iterator &cell,
const unsigned int face_no,
const Quadrature<dim-1>& quadrature,
- Mapping<dim>::InternalDataBase &mapping_data,
+ typename Mapping<dim>::InternalDataBase &mapping_data,
std::vector<Point<dim> > &quadrature_points,
std::vector<double> &JxW_values,
std::vector<Tensor<1,dim> > &exterior_form,
* @p{compute_support_points_simple}
*/
//TODO:[RH] (later) remove this function altogether?
- void fill_quad_support_points_simple (const Triangulation<dim>::cell_iterator &cell,
+ void fill_quad_support_points_simple (const typename Triangulation<dim>::cell_iterator &cell,
std::vector<Point<dim> > &a) const;
/**
* @ref{Mapping}.
*/
virtual void
- fill_fe_values (const DoFHandler<dim>::cell_iterator &cell,
+ fill_fe_values (const typename DoFHandler<dim>::cell_iterator &cell,
const Quadrature<dim>& quadrature,
- Mapping<dim>::InternalDataBase &mapping_data,
+ typename Mapping<dim>::InternalDataBase &mapping_data,
std::vector<Point<dim> > &quadrature_points,
std::vector<double> &JxW_values) const ;
// already have built.
for (int level=maxlevel; level>=static_cast<int>(minlevel); --level)
{
- MGDoFHandler<dim>::active_cell_iterator
+ typename MGDoFHandler<dim>::active_cell_iterator
level_cell = mg_dof_handler->begin_active(level);
- const MGDoFHandler<dim>::active_cell_iterator
+ const typename MGDoFHandler<dim>::active_cell_iterator
level_end = mg_dof_handler->end_active(level);
//TODO:[?] Treat hanging nodes properly
for (unsigned int face_n=0; face_n<GeometryInfo<dim>::faces_per_cell; ++face_n)
{
- const MGDoFHandler<dim>::face_iterator face = level_cell->face(face_n);
+ const typename MGDoFHandler<dim>::face_iterator
+ face = level_cell->face(face_n);
if (face->has_children())
{
face->get_mg_dof_indices(level_face_indices);
std::vector<unsigned int> global_dof_indices (dofs_per_cell);
std::vector<unsigned int> level_dof_indices (dofs_per_cell);
- MGDoFHandler<dim>::active_cell_iterator
+ typename MGDoFHandler<dim>::active_cell_iterator
level_cell = mg_dof_handler->begin_active();
- const MGDoFHandler<dim>::active_cell_iterator
+ const typename MGDoFHandler<dim>::active_cell_iterator
endc = mg_dof_handler->end();
// traverse all cells and copy the
std::vector<unsigned int> global_dof_indices (dofs_per_cell);
std::vector<unsigned int> level_dof_indices (dofs_per_cell);
- DoFHandler<dim>::active_cell_iterator
+ typename DoFHandler<dim>::active_cell_iterator
global_cell = dof->begin_active(level);
- MGDoFHandler<dim>::active_cell_iterator
+ typename MGDoFHandler<dim>::active_cell_iterator
level_cell = mg_dof_handler->begin_active(level);
- const MGDoFHandler<dim>::cell_iterator
+ const typename MGDoFHandler<dim>::cell_iterator
endc = mg_dof_handler->end(level);
// traverse all cells and copy the
const unsigned int dofs_per_cell = dof.get_fe().dofs_per_cell;
std::vector<unsigned int> dofs_on_this_cell(dofs_per_cell);
- DoFHandler<dim>::active_cell_iterator cell = dof.begin_active(),
- endc = dof.end();
+ typename DoFHandler<dim>::active_cell_iterator cell = dof.begin_active(),
+ endc = dof.end();
for (; cell!=endc; ++cell)
{
cell->get_dof_indices (dofs_on_this_cell);
std::vector<unsigned int> dofs_on_this_cell(dofs_per_cell);
- DoFHandler<dim>::active_cell_iterator cell = dof.begin_active(),
- endc = dof.end();
+ typename DoFHandler<dim>::active_cell_iterator cell = dof.begin_active(),
+ endc = dof.end();
for (; cell!=endc; ++cell)
{
cell->get_dof_indices (dofs_on_this_cell);
const unsigned int total_dofs = dof.get_fe().dofs_per_cell;
std::vector<unsigned int> dofs_on_this_cell(total_dofs);
std::vector<unsigned int> dofs_on_other_cell(total_dofs);
- DoFHandler<dim>::active_cell_iterator cell = dof.begin_active(),
- endc = dof.end();
+ typename DoFHandler<dim>::active_cell_iterator cell = dof.begin_active(),
+ endc = dof.end();
// clear user flags for further use
const_cast<Triangulation<dim>&>(dof.get_tria()).clear_user_flags();
face < GeometryInfo<dim>::faces_per_cell;
++face)
{
- DoFHandler<dim>::face_iterator cell_face = cell->face(face);
+ typename DoFHandler<dim>::face_iterator cell_face = cell->face(face);
if (cell_face->user_flag_set ())
continue;
if (! cell->at_boundary (face) )
{
- DoFHandler<dim>::cell_iterator neighbor = cell->neighbor(face);
+ typename DoFHandler<dim>::cell_iterator
+ neighbor = cell->neighbor(face);
// Refinement edges are taken care of
// by coarser cells
if (neighbor->level() < cell->level())
sub_nr != GeometryInfo<dim>::subfaces_per_face;
++sub_nr)
{
- DoFHandler<dim>::cell_iterator sub_neighbor
+ typename DoFHandler<dim>::cell_iterator sub_neighbor
= neighbor->
child(GeometryInfo<dim>::child_cell_on_face(neighbor_face, sub_nr));
const unsigned int total_dofs = dof.get_fe().dofs_per_cell;
std::vector<unsigned int> dofs_on_this_cell(total_dofs);
std::vector<unsigned int> dofs_on_other_cell(total_dofs);
- DoFHandler<dim>::active_cell_iterator cell = dof.begin_active(),
- endc = dof.end();
+ typename DoFHandler<dim>::active_cell_iterator cell = dof.begin_active(),
+ endc = dof.end();
std::vector<std::vector<bool> > int_dof_mask(total_dofs,
face < GeometryInfo<dim>::faces_per_cell;
++face)
{
- DoFHandler<dim>::face_iterator cell_face = cell->face(face);
+ typename DoFHandler<dim>::face_iterator cell_face = cell->face(face);
if (cell_face->user_flag_set ())
continue;
if (! cell->at_boundary (face) )
{
- DoFHandler<dim>::cell_iterator neighbor = cell->neighbor(face);
+ typename DoFHandler<dim>::cell_iterator
+ neighbor = cell->neighbor(face);
// Refinement edges are taken care of
// by coarser cells
if (neighbor->level() < cell->level())
sub_nr != GeometryInfo<dim>::subfaces_per_face;
++sub_nr)
{
- DoFHandler<dim>::cell_iterator sub_neighbor
+ typename DoFHandler<dim>::cell_iterator sub_neighbor
= neighbor->
child(GeometryInfo<dim>::child_cell_on_face(neighbor_face, sub_nr));
// in the sum for each dof
std::vector<unsigned char> touch_count (dof_handler.n_dofs(), 0);
- DoFHandler<dim>::active_cell_iterator cell = dof_handler.begin_active(),
- endc = dof_handler.end();
+ typename DoFHandler<dim>::active_cell_iterator
+ cell = dof_handler.begin_active(),
+ endc = dof_handler.end();
unsigned int present_cell = 0;
const unsigned int dofs_per_cell = fe.dofs_per_cell;
std::vector<unsigned int> dof_indices (dofs_per_cell);
std::vector<unsigned int> indices(fe.dofs_per_cell);
- DoFHandler<dim>::active_cell_iterator c;
+ typename DoFHandler<dim>::active_cell_iterator c;
for (c = dof.begin_active() ; c != dof.end() ; ++ c)
{
c->get_dof_indices(indices);
std::vector<unsigned int> indices(fe.dofs_per_cell);
- MGDoFHandler<dim>::cell_iterator c;
+ typename MGDoFHandler<dim>::cell_iterator c;
for (c = dof.begin(level) ; c != dof.end(level) ; ++ c)
{
c->get_mg_dof_indices(indices);
// line is also part of a boundary
// face which we will be visiting
// sooner or later
- for (DoFHandler<dim>::active_cell_iterator cell=dof_handler.begin_active();
+ for (typename DoFHandler<dim>::active_cell_iterator
+ cell=dof_handler.begin_active();
cell!=dof_handler.end(); ++cell)
for (unsigned int face=0; face<GeometryInfo<dim>::faces_per_cell; ++face)
if (cell->at_boundary(face))
std::vector<bool> dof_is_interesting (fine_grid.n_dofs(), false);
std::vector<unsigned int> local_dof_indices (fine_fe.dofs_per_cell);
- for (DoFHandler<dim>::active_cell_iterator cell=fine_grid.begin_active();
+ for (typename DoFHandler<dim>::active_cell_iterator
+ cell=fine_grid.begin_active();
cell!=fine_grid.end(); ++cell)
{
cell->get_dof_indices (local_dof_indices);
void
FiniteElementBase<dim>::
compute_2nd (const Mapping<dim> &mapping,
- const DoFHandler<dim>::cell_iterator &cell,
+ const typename DoFHandler<dim>::cell_iterator &cell,
const unsigned int offset,
- Mapping<dim>::InternalDataBase &mapping_internal,
+ typename Mapping<dim>::InternalDataBase &mapping_internal,
InternalDataBase &fe_internal,
FEValuesData<dim> &data) const
{
template <int dim>
void
FE_DGQ<dim>::fill_fe_values (const Mapping<dim> &mapping,
- const DoFHandler<dim>::cell_iterator &cell,
+ const typename DoFHandler<dim>::cell_iterator &cell,
const Quadrature<dim> &quadrature,
- Mapping<dim>::InternalDataBase &mapping_data,
- Mapping<dim>::InternalDataBase &fedata,
+ typename Mapping<dim>::InternalDataBase &mapping_data,
+ typename Mapping<dim>::InternalDataBase &fedata,
FEValuesData<dim> &data) const
{
// convert data object to internal
template <int dim>
void
FE_DGQ<dim>::fill_fe_face_values (const Mapping<dim> &mapping,
- const DoFHandler<dim>::cell_iterator &cell,
+ const typename DoFHandler<dim>::cell_iterator &cell,
const unsigned int face,
const Quadrature<dim-1> &quadrature,
- Mapping<dim>::InternalDataBase &mapping_data,
- Mapping<dim>::InternalDataBase &fedata,
+ typename Mapping<dim>::InternalDataBase &mapping_data,
+ typename Mapping<dim>::InternalDataBase &fedata,
FEValuesData<dim> &data) const
{
// convert data object to internal
template <int dim>
void
FE_DGQ<dim>::fill_fe_subface_values (const Mapping<dim> &mapping,
- const DoFHandler<dim>::cell_iterator &cell,
+ const typename DoFHandler<dim>::cell_iterator &cell,
const unsigned int face,
const unsigned int subface,
const Quadrature<dim-1> &quadrature,
- Mapping<dim>::InternalDataBase &mapping_data,
- Mapping<dim>::InternalDataBase &fedata,
+ typename Mapping<dim>::InternalDataBase &mapping_data,
+ typename Mapping<dim>::InternalDataBase &fedata,
FEValuesData<dim> &data) const
{
// convert data object to internal
template <int dim>
void
FE_Q<dim>::fill_fe_values (const Mapping<dim> &mapping,
- const DoFHandler<dim>::cell_iterator &cell,
+ const typename DoFHandler<dim>::cell_iterator &cell,
const Quadrature<dim> &quadrature,
- Mapping<dim>::InternalDataBase &mapping_data,
- Mapping<dim>::InternalDataBase &fedata,
+ typename Mapping<dim>::InternalDataBase &mapping_data,
+ typename Mapping<dim>::InternalDataBase &fedata,
FEValuesData<dim> &data) const
{
// convert data object to internal
template <int dim>
void
FE_Q<dim>::fill_fe_face_values (const Mapping<dim> &mapping,
- const DoFHandler<dim>::cell_iterator &cell,
+ const typename DoFHandler<dim>::cell_iterator &cell,
const unsigned int face,
const Quadrature<dim-1> &quadrature,
- Mapping<dim>::InternalDataBase &mapping_data,
- Mapping<dim>::InternalDataBase &fedata,
+ typename Mapping<dim>::InternalDataBase &mapping_data,
+ typename Mapping<dim>::InternalDataBase &fedata,
FEValuesData<dim> &data) const
{
// convert data object to internal
template <int dim>
void
FE_Q<dim>::fill_fe_subface_values (const Mapping<dim> &mapping,
- const DoFHandler<dim>::cell_iterator &cell,
+ const typename DoFHandler<dim>::cell_iterator &cell,
const unsigned int face,
const unsigned int subface,
const Quadrature<dim-1> &quadrature,
- Mapping<dim>::InternalDataBase &mapping_data,
- Mapping<dim>::InternalDataBase &fedata,
+ typename Mapping<dim>::InternalDataBase &mapping_data,
+ typename Mapping<dim>::InternalDataBase &fedata,
FEValuesData<dim> &data) const
{
// convert data object to internal
template <int dim>
void
FESystem<dim>::fill_fe_values (const Mapping<dim> &mapping,
- const DoFHandler<dim>::cell_iterator &cell,
+ const typename DoFHandler<dim>::cell_iterator &cell,
const Quadrature<dim> &quadrature,
- Mapping<dim>::InternalDataBase &mapping_data,
- Mapping<dim>::InternalDataBase &fe_data,
+ typename Mapping<dim>::InternalDataBase &mapping_data,
+ typename Mapping<dim>::InternalDataBase &fe_data,
FEValuesData<dim> &data) const
{
compute_fill(mapping, cell, invalid_face_number, invalid_face_number,
template <int dim>
void
FESystem<dim>::fill_fe_face_values (const Mapping<dim> &mapping,
- const DoFHandler<dim>::cell_iterator &cell,
+ const typename DoFHandler<dim>::cell_iterator &cell,
const unsigned int face_no,
const Quadrature<dim-1> &quadrature,
- Mapping<dim>::InternalDataBase &mapping_data,
- Mapping<dim>::InternalDataBase &fe_data,
+ typename Mapping<dim>::InternalDataBase &mapping_data,
+ typename Mapping<dim>::InternalDataBase &fe_data,
FEValuesData<dim> &data) const
{
compute_fill(mapping, cell, face_no, invalid_face_number,
template <int dim>
void
FESystem<dim>::fill_fe_subface_values (const Mapping<dim> &mapping,
- const DoFHandler<dim>::cell_iterator &cell,
+ const typename DoFHandler<dim>::cell_iterator &cell,
const unsigned int face_no,
const unsigned int sub_no,
const Quadrature<dim-1> &quadrature,
- Mapping<dim>::InternalDataBase &mapping_data,
- Mapping<dim>::InternalDataBase &fe_data,
+ typename Mapping<dim>::InternalDataBase &mapping_data,
+ typename Mapping<dim>::InternalDataBase &fe_data,
FEValuesData<dim> &data) const
{
compute_fill(mapping, cell, face_no, sub_no,
template <int dim_1>
void
FESystem<dim>::compute_fill (const Mapping<dim> &mapping,
- const DoFHandler<dim>::cell_iterator &cell,
+ const typename DoFHandler<dim>::cell_iterator &cell,
const unsigned int face_no,
const unsigned int sub_no,
const Quadrature<dim_1> &quadrature,
- Mapping<dim>::InternalDataBase &mapping_data,
- Mapping<dim>::InternalDataBase &fedata,
+ typename Mapping<dim>::InternalDataBase &mapping_data,
+ typename Mapping<dim>::InternalDataBase &fedata,
FEValuesData<dim> &data) const
{
// convert data object to internal
FETools::get_interpolation_matrix(dof1.get_fe(), dof2.get_fe(),
interpolation_matrix);
- DoFHandler<dim>::active_cell_iterator cell1 = dof1.begin_active(),
- endc1 = dof1.end(),
- cell2 = dof2.begin_active(),
- endc2 = dof2.end();
+ typename DoFHandler<dim>::active_cell_iterator cell1 = dof1.begin_active(),
+ endc1 = dof1.end(),
+ cell2 = dof2.begin_active(),
+ endc2 = dof2.end();
std::vector<unsigned int> index_multiplicity(dof2.n_dofs(),0);
std::vector<unsigned int> dofs (dofs_per_cell2);
FETools::get_back_interpolation_matrix(dof1.get_fe(), fe2,
interpolation_matrix);
- DoFHandler<dim>::active_cell_iterator cell1 = dof1.begin_active(),
- endc1 = dof1.end();
+ typename DoFHandler<dim>::active_cell_iterator cell1 = dof1.begin_active(),
+ endc1 = dof1.end();
for (; cell1!=endc1; ++cell1)
{
FETools::get_interpolation_difference_matrix(dof1.get_fe(), fe2,
difference_matrix);
- DoFHandler<dim>::active_cell_iterator cell1 = dof1.begin_active(),
- endc1 = dof1.end();
+ typename DoFHandler<dim>::active_cell_iterator cell1 = dof1.begin_active(),
+ endc1 = dof1.end();
for (; cell1!=endc1; ++cell1)
{
template <int dim>
void
-MappingCartesian<dim>::fill_fe_values (const DoFHandler<dim>::cell_iterator& cell,
- const Quadrature<dim>& q,
- Mapping<dim>::InternalDataBase& mapping_data,
- std::vector<Point<dim> >& quadrature_points,
- std::vector<double>& JxW_values) const
+MappingCartesian<dim>::fill_fe_values (const typename DoFHandler<dim>::cell_iterator& cell,
+ const Quadrature<dim>& q,
+ typename Mapping<dim>::InternalDataBase& mapping_data,
+ std::vector<Point<dim> >& quadrature_points,
+ std::vector<double>& JxW_values) const
{
// convert data object to internal
// data for this class. fails with
template <int dim>
void
-MappingQ<dim>::fill_fe_values (const DoFHandler<dim>::cell_iterator &cell,
+MappingQ<dim>::fill_fe_values (const typename DoFHandler<dim>::cell_iterator &cell,
const Quadrature<dim> &q,
- Mapping<dim>::InternalDataBase &mapping_data,
+ typename Mapping<dim>::InternalDataBase &mapping_data,
std::vector<Point<dim> > &quadrature_points,
std::vector<double> &JxW_values) const
{
template <int dim>
void
-MappingQ<dim>::add_line_support_points (const Triangulation<dim>::cell_iterator &cell,
+MappingQ<dim>::add_line_support_points (const typename Triangulation<dim>::cell_iterator &cell,
std::vector<Point<dim> > &a) const
{
static const StraightBoundary<dim> straight_boundary;
template <int dim>
void
-MappingQ<dim>::fill_quad_support_points_simple (const Triangulation<dim>::cell_iterator &,
+MappingQ<dim>::fill_quad_support_points_simple (const typename Triangulation<dim>::cell_iterator &,
std::vector<Point<dim> > &) const
{}
template <int dim>
void
-MappingQ1<dim>::fill_fe_values (const DoFHandler<dim>::cell_iterator &cell,
+MappingQ1<dim>::fill_fe_values (const typename DoFHandler<dim>::cell_iterator &cell,
const Quadrature<dim> &q,
- Mapping<dim>::InternalDataBase &mapping_data,
+ typename Mapping<dim>::InternalDataBase &mapping_data,
std::vector<Point<dim> > &quadrature_points,
std::vector<double> &JxW_values) const
{
case 2:
{
- Triangulation<dim>::active_line_iterator line =tria.begin_active_line ();
- Triangulation<dim>::active_line_iterator endline=tria.end_line ();
+ typename Triangulation<dim>::active_line_iterator
+ line =tria.begin_active_line (),
+ endline=tria.end_line ();
// first treat all interior
// lines and make up a list
//TODO:[RH] curved boundaries in eps for 3d
Assert (mapping == 0, ExcNotImplemented());
- Triangulation<dim>::active_line_iterator line =tria.begin_active_line ();
- Triangulation<dim>::active_line_iterator endline=tria.end_line ();
+ typename Triangulation<dim>::active_line_iterator
+ line =tria.begin_active_line (),
+ endline=tria.end_line ();
// loop over all lines and compute their
// projection on the plane perpendicular
if (threshold==0)
return;
- Triangulation<dim>::active_cell_iterator cell = tria.begin_active();
+ typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active();
const unsigned int n_cells = criteria.size();
for (unsigned int index=0; index<n_cells; ++cell, ++index)
Assert (*std::min_element(criteria.begin(), criteria.end()) >= 0,
ExcInvalidParameterValue());
- Triangulation<dim>::active_cell_iterator cell = tria.begin_active();
+ typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active();
const unsigned int n_cells = criteria.size();
for (unsigned int index=0; index<n_cells; ++cell, ++index)
template <int dim>
void
HyperBallBoundary<dim>::get_intermediate_points_on_quad (
- const Triangulation<dim>::quad_iterator &,
+ const typename Triangulation<dim>::quad_iterator &,
typename std::vector<Point<dim> > &) const
{
Assert(false, Boundary<dim>::ExcFunctionNotUseful(dim));
template <int dim>
void
HyperBallBoundary<dim>::
-get_normals_at_vertices (const Triangulation<dim>::face_iterator &face,
+get_normals_at_vertices (const typename Triangulation<dim>::face_iterator &face,
typename Boundary<dim>::FaceVertexNormals &face_vertex_normals) const
{
for (unsigned int vertex=0; vertex<GeometryInfo<dim>::vertices_per_face; ++vertex)
const unsigned int dofs_per_cell = dof.get_fe().dofs_per_cell;
std::vector<unsigned int> dofs_on_this_cell(dofs_per_cell);
- MGDoFHandler<dim>::cell_iterator cell = dof.begin(level),
- endc = dof.end(level);
+ typename MGDoFHandler<dim>::cell_iterator cell = dof.begin(level),
+ endc = dof.end(level);
for (; cell!=endc; ++cell)
{
cell->get_mg_dof_indices (dofs_on_this_cell);
const unsigned int dofs_per_cell = dof.get_fe().dofs_per_cell;
std::vector<unsigned int> dofs_on_this_cell(dofs_per_cell);
std::vector<unsigned int> dofs_on_other_cell(dofs_per_cell);
- MGDoFHandler<dim>::cell_iterator cell = dof.begin(level),
- endc = dof.end(level);
+ typename MGDoFHandler<dim>::cell_iterator cell = dof.begin(level),
+ endc = dof.end(level);
for (; cell!=endc; ++cell)
{
cell->get_mg_dof_indices (dofs_on_this_cell);
if ( (! cell->at_boundary(face)) &&
(static_cast<unsigned int>(cell->neighbor_level(face)) == level) )
{
- MGDoFHandler<dim>::cell_iterator neighbor = cell->neighbor(face);
+ typename MGDoFHandler<dim>::cell_iterator
+ neighbor = cell->neighbor(face);
neighbor->get_mg_dof_indices (dofs_on_other_cell);
// only add one direction
// The other is taken care of
const unsigned int dofs_per_cell = dof.get_fe().dofs_per_cell;
std::vector<unsigned int> dofs_on_this_cell(dofs_per_cell);
std::vector<unsigned int> dofs_on_other_cell(dofs_per_cell);
- MGDoFHandler<dim>::cell_iterator cell = dof.begin(level),
- endc = dof.end(level);
+ typename MGDoFHandler<dim>::cell_iterator cell = dof.begin(level),
+ endc = dof.end(level);
for (; cell!=endc; ++cell)
{
cell->get_mg_dof_indices (dofs_on_this_cell);
if ( (! cell->at_boundary(face)) &&
(static_cast<unsigned int>(cell->neighbor_level(face)) != level) )
{
- MGDoFHandler<dim>::cell_iterator neighbor = cell->neighbor(face);
+ typename MGDoFHandler<dim>::cell_iterator
+ neighbor = cell->neighbor(face);
neighbor->get_mg_dof_indices (dofs_on_other_cell);
for (unsigned int i=0; i<dofs_per_cell; ++i)
const unsigned int total_dofs = dof.get_fe().dofs_per_cell;
std::vector<unsigned int> dofs_on_this_cell(total_dofs);
std::vector<unsigned int> dofs_on_other_cell(total_dofs);
- MGDoFHandler<dim>::cell_iterator cell = dof.begin(level),
- endc = dof.end(level);
+ typename MGDoFHandler<dim>::cell_iterator cell = dof.begin(level),
+ endc = dof.end(level);
std::vector<std::vector<bool> > int_dof_mask(total_dofs,
face < GeometryInfo<dim>::faces_per_cell;
++face)
{
- MGDoFHandler<dim>::face_iterator cell_face = cell->face(face);
+ typename MGDoFHandler<dim>::face_iterator cell_face = cell->face(face);
if (cell_face->user_flag_set ())
continue;
if (! cell->at_boundary (face) )
{
- MGDoFHandler<dim>::cell_iterator neighbor = cell->neighbor(face);
+ typename MGDoFHandler<dim>::cell_iterator
+ neighbor = cell->neighbor(face);
if (neighbor->level() < cell->level())
continue;
// at the boundary of this level's
// cells and if so add the dofs
// to the interior boundary dofs
- for (MGDoFHandler<dim>::cell_iterator cell=mg_dof.begin(level);
+ for (typename MGDoFHandler<dim>::cell_iterator cell=mg_dof.begin(level);
cell != mg_dof.end(level); ++cell)
for (unsigned int face=0; face<GeometryInfo<dim>::faces_per_cell; ++face)
if ((cell->neighbor(face).state() == valid) &&
prolongation_matrices.push_back (SparseMatrix<double>());
prolongation_matrices[level].reinit (prolongation_sparsities[level]);
-// now actually build the matrices
- for (MGDoFHandler<dim>::cell_iterator cell=mg_dof.begin(level);
+ // now actually build the matrices
+ for (typename MGDoFHandler<dim>::cell_iterator cell=mg_dof.begin(level);
cell != mg_dof.end(level); ++cell)
if (cell->has_children())
{
unsigned int cell_number = 0;
typename std::vector<DataOutBase::Patch<dim> >::iterator patch = patches.begin();
- DoFHandler<dim>::cell_iterator cell=first_cell();
+ typename DoFHandler<dim>::cell_iterator cell=first_cell();
// get first cell in this thread
for (unsigned int i=0; (i<data.this_thread)&&(cell != dofs->end()); ++i)
// create patches of and make sure
// there is enough memory for that
unsigned int n_patches = 0;
- for (DoFHandler<dim>::cell_iterator cell=first_cell();
+ for (typename DoFHandler<dim>::cell_iterator cell=first_cell();
cell != dofs->end();
cell = next_cell(cell))
++n_patches;
// create patches of and make sure
// there is enough memory for that
unsigned int n_patches = 0;
- for (DoFHandler<dim>::cell_iterator cell=first_cell();
+ for (typename DoFHandler<dim>::cell_iterator cell=first_cell();
cell != dofs->end();
cell = next_cell(cell))
++n_patches;
// create patches of and make sure
// there is enough memory for that
unsigned int n_patches = 0;
- for (DoFHandler<dim>::active_cell_iterator cell=dof_handler->begin_active();
+ for (typename DoFHandler<dim>::active_cell_iterator
+ cell=dof_handler->begin_active();
cell != dof_handler->end(); ++cell)
++n_patches;
indices_on_cell=std::vector<std::vector<unsigned int> > (n_active_cells,
std::vector<unsigned int> (dofs_per_cell));
- DoFHandler<dim>::cell_iterator cell = dof_handler->begin(),
- endc = dof_handler->end();
+ typename DoFHandler<dim>::cell_iterator cell = dof_handler->begin(),
+ endc = dof_handler->end();
for (unsigned int i=0; cell!=endc; ++cell)
{
unsigned int dofs_per_cell=dof_handler->get_fe().dofs_per_cell;
Vector<number> local_values(dofs_per_cell);
- DoFHandler<dim>::cell_iterator cell = dof_handler->begin(),
- endc = dof_handler->end();
+ typename DoFHandler<dim>::cell_iterator cell = dof_handler->begin(),
+ endc = dof_handler->end();
std::vector<unsigned int> *indexptr;
// and that'll stay or be refined
unsigned int n_cells_to_coarsen=0;
unsigned int n_cells_to_stay_or_refine=0;
- DoFHandler<dim>::active_cell_iterator act_cell = dof_handler->begin_active(),
- endc = dof_handler->end();
+ typename DoFHandler<dim>::active_cell_iterator
+ act_cell = dof_handler->begin_active(),
+ endc = dof_handler->end();
for (; act_cell!=endc; ++act_cell)
{
if (act_cell->coarsen_flag_set())
// and all the cells where a
// @p{Pointerstruct} is needed 'n'
unsigned int n_sr=0, n_cf=0, n=0;
- DoFHandler<dim>::cell_iterator cell = dof_handler->begin();
+ typename DoFHandler<dim>::cell_iterator cell = dof_handler->begin();
for (; cell!=endc; ++cell)
{
if (cell->active() && !cell->coarsen_flag_set())
std::vector<Vector<number> > *valuesptr;
std::vector<unsigned int> dofs(dofs_per_cell);
- DoFHandler<dim>::cell_iterator cell = dof_handler->begin(),
- endc = dof_handler->end();
+ typename DoFHandler<dim>::cell_iterator cell = dof_handler->begin(),
+ endc = dof_handler->end();
for (; cell!=endc; ++cell)
{
// desired so
// if (flags.max_refinement_level != 0)
// {
-// Triangulation<dim>::active_cell_iterator cell, endc;
+// typename Triangulation<dim>::active_cell_iterator cell, endc;
// for (cell = tria->begin_active(),
// endc = tria->end();
// cell!=endc; ++cell)
{
bool grids_changed = false;
- Triangulation<dim>::cell_iterator cell1 = tria1.begin(),
- cell2 = tria2.begin();
- Triangulation<dim>::cell_iterator endc;
+ typename Triangulation<dim>::cell_iterator cell1 = tria1.begin(),
+ cell2 = tria2.begin();
+ typename Triangulation<dim>::cell_iterator endc;
endc = (tria1.n_levels() == 1 ?
- Triangulation<dim>::cell_iterator(tria1.end()) :
+ typename Triangulation<dim>::cell_iterator(tria1.end()) :
tria1.begin(1));
for (; cell1!=endc; ++cell1, ++cell2)
grids_changed |= adapt_grids<dim> (cell1, cell2);
// cells, we have to subtract 3/4 of
// a cell for each flagged cell
double previous_cells = previous_tria->n_active_cells();
- Triangulation<dim>::active_cell_iterator cell, endc;
+ typename Triangulation<dim>::active_cell_iterator cell, endc;
cell = previous_tria->begin_active();
endc = previous_tria->end();
for (; cell!=endc; ++cell)
{
adapt_grids (*previous_tria, *tria);
- Triangulation<dim>::cell_iterator old_cell, new_cell, endc;
+ typename Triangulation<dim>::cell_iterator old_cell, new_cell, endc;
old_cell = previous_tria->begin(0);
new_cell = tria->begin(0);
endc = tria->end(0);
const unsigned int n_components = fe.n_components();
const bool fe_is_system = (n_components != 1);
- DoFHandler<dim>::active_cell_iterator cell = dof.begin_active(),
- endc = dof.end();
+ typename DoFHandler<dim>::active_cell_iterator cell = dof.begin_active(),
+ endc = dof.end();
// For FESystems many of the
// unit_support_points will
std::vector<short unsigned int> touch_count (dof_2.n_dofs(), 0);
std::vector<unsigned int> local_dof_indices (dof_2.get_fe().dofs_per_cell);
- DoFHandler<dim>::active_cell_iterator h = dof_1.begin_active();
- DoFHandler<dim>::active_cell_iterator l = dof_2.begin_active();
- const DoFHandler<dim>::cell_iterator endh = dof_1.end();
+ typename DoFHandler<dim>::active_cell_iterator h = dof_1.begin_active();
+ typename DoFHandler<dim>::active_cell_iterator l = dof_2.begin_active();
+ const typename DoFHandler<dim>::cell_iterator endh = dof_1.end();
for(; h != endh; ++h, ++l)
{
std::vector<Tensor<1,dim> > tmp_gradients (fe_values.n_quadrature_points);
// loop over all cells
- DoFHandler<dim>::active_cell_iterator cell = dof.begin_active(),
- endc = dof.end();
+ typename DoFHandler<dim>::active_cell_iterator cell = dof.begin_active(),
+ endc = dof.end();
for (unsigned int index=0; cell != endc; ++cell, ++index)
{
double diff=0;