grid/grid_test: grid/grid_test.o $(LIBFILES)
@echo ================= Linking $@
- @$(CXX) $(CXXFLAGS) -o $@ $< $(LIBS)
+ @$(CXX) $(CXXFLAGS) -o $@ $< $(LIBS.g)
dof/dof_test: dof/dof_test.o $(LIBFILES.g)
@echo ================= Linking $@
+
DoFHandler<1>::raw_cell_iterator
DoFHandler<1>::begin_raw (const unsigned int level) const {
return begin_raw_line (level);
//------------------------------------------------------------------
-DoFHandler<1>::raw_cell_iterator
-DoFHandler<1>::begin_raw_line (const unsigned int level) const {
- return raw_line_iterator (tria,
- tria->begin_raw_line(level)->level(),
- tria->begin_raw_line(level)->index(),
- this);
-};
-
-
-DoFHandler<2>::raw_line_iterator
-DoFHandler<2>::begin_raw_line (const unsigned int level) const {
+template <int dim>
+typename DoFHandler<dim>::raw_line_iterator
+DoFHandler<dim>::begin_raw_line (const unsigned int level) const {
return raw_line_iterator (tria,
tria->begin_raw_line(level)->level(),
tria->begin_raw_line(level)->index(),
-DoFHandler<1>::line_iterator
-DoFHandler<1>::begin_line (const unsigned int level) const {
- return line_iterator (tria,
- tria->begin_line(level)->level(),
- tria->begin_line(level)->index(),
- this);
-};
-
-
-
-DoFHandler<2>::line_iterator
-DoFHandler<2>::begin_line (const unsigned int level) const {
+template <int dim>
+typename DoFHandler<dim>::line_iterator
+DoFHandler<dim>::begin_line (const unsigned int level) const {
return line_iterator (tria,
tria->begin_line(level)->level(),
tria->begin_line(level)->index(),
-DoFHandler<1>::active_line_iterator
-DoFHandler<1>::begin_active_line (const unsigned int level) const {
+template <int dim>
+typename DoFHandler<dim>::active_line_iterator
+DoFHandler<dim>::begin_active_line (const unsigned int level) const {
return active_line_iterator (tria,
tria->begin_active_line(level)->level(),
tria->begin_active_line(level)->index(),
-DoFHandler<2>::active_line_iterator
-DoFHandler<2>::begin_active_line (const unsigned int level) const {
- return active_line_iterator (tria,
- tria->begin_active_line(level)->level(),
- tria->begin_active_line(level)->index(),
- this);
+DoFHandler<1>::raw_quad_iterator
+DoFHandler<1>::begin_raw_quad (const unsigned int) const {
+ Assert (false, ExcNotImplemented());
+ return 0;
};
-DoFHandler<2>::raw_quad_iterator
-DoFHandler<2>::begin_raw_quad (const unsigned int level) const {
+template <int dim>
+typename DoFHandler<dim>::raw_quad_iterator
+DoFHandler<dim>::begin_raw_quad (const unsigned int level) const {
return raw_quad_iterator (tria,
tria->begin_raw_quad(level)->level(),
tria->begin_raw_quad(level)->index(),
-DoFHandler<2>::quad_iterator
-DoFHandler<2>::begin_quad (const unsigned int level) const {
+DoFHandler<1>::quad_iterator
+DoFHandler<1>::begin_quad (const unsigned int) const {
+ Assert (false, ExcNotImplemented());
+ return 0;
+};
+
+
+
+template <int dim>
+typename DoFHandler<dim>::quad_iterator
+DoFHandler<dim>::begin_quad (const unsigned int level) const {
return quad_iterator (tria,
tria->begin_quad(level)->level(),
tria->begin_quad(level)->index(),
-DoFHandler<2>::active_quad_iterator
-DoFHandler<2>::begin_active_quad (const unsigned int level) const {
+DoFHandler<1>::active_quad_iterator
+DoFHandler<1>::begin_active_quad (const unsigned int) const {
+ Assert (false, ExcNotImplemented());
+ return 0;
+};
+
+
+
+template <int dim>
+typename DoFHandler<dim>::active_quad_iterator
+DoFHandler<dim>::begin_active_quad (const unsigned int level) const {
return active_quad_iterator (tria,
tria->begin_active_quad(level)->level(),
tria->begin_active_quad(level)->index(),
-DoFHandler<1>::raw_line_iterator
-DoFHandler<1>::end_line () const {
+template <int dim>
+typename DoFHandler<dim>::raw_line_iterator
+DoFHandler<dim>::end_line () const {
return raw_line_iterator (tria, -1, -1, this);
};
-DoFHandler<2>::raw_line_iterator
-DoFHandler<2>::end_line () const {
- return raw_line_iterator (tria, -1, -1, this);
+DoFHandler<1>::raw_quad_iterator
+DoFHandler<1>::end_quad () const {
+ Assert (false, ExcNotImplemented());
+ return 0;
};
-DoFHandler<2>::raw_quad_iterator
-DoFHandler<2>::end_quad () const {
+template <int dim>
+typename DoFHandler<dim>::raw_quad_iterator
+DoFHandler<dim>::end_quad () const {
return raw_quad_iterator (tria, -1, -1, this);
};
-DoFHandler<1>::raw_line_iterator
-DoFHandler<1>::last_raw_line (const unsigned int level) const {
- return raw_line_iterator (tria,
- tria->last_raw_line(level)->level(),
- tria->last_raw_line(level)->index(),
- this);
-};
-
-
-
-DoFHandler<2>::raw_line_iterator
-DoFHandler<2>::last_raw_line (const unsigned int level) const {
+template <int dim>
+typename DoFHandler<dim>::raw_line_iterator
+DoFHandler<dim>::last_raw_line (const unsigned int level) const {
return raw_line_iterator (tria,
tria->last_raw_line(level)->level(),
tria->last_raw_line(level)->index(),
-DoFHandler<1>::line_iterator
-DoFHandler<1>::last_line (const unsigned int level) const {
- return line_iterator (tria,
- tria->last_line(level)->level(),
- tria->last_line(level)->index(),
- this);
-};
-
-
-
-DoFHandler<2>::line_iterator
-DoFHandler<2>::last_line (const unsigned int level) const {
+template <int dim>
+typename DoFHandler<dim>::line_iterator
+DoFHandler<dim>::last_line (const unsigned int level) const {
return line_iterator (tria,
tria->last_line(level)->level(),
tria->last_line(level)->index(),
-DoFHandler<1>::active_line_iterator
-DoFHandler<1>::last_active_line (const unsigned int level) const {
+template <int dim>
+typename DoFHandler<dim>::active_line_iterator
+DoFHandler<dim>::last_active_line (const unsigned int level) const {
return active_line_iterator (tria,
tria->last_active_line(level)->level(),
tria->last_active_line(level)->index(),
-DoFHandler<2>::active_line_iterator
-DoFHandler<2>::last_active_line (const unsigned int level) const {
- return active_line_iterator (tria,
- tria->last_active_line(level)->level(),
- tria->last_active_line(level)->index(),
- this);
+DoFHandler<1>::raw_quad_iterator
+DoFHandler<1>::last_raw_quad (const unsigned int) const {
+ Assert (false, ExcNotImplemented());
+ return 0;
};
-
-DoFHandler<2>::raw_quad_iterator
-DoFHandler<2>::last_raw_quad (const unsigned int level) const {
+template <int dim>
+typename DoFHandler<dim>::raw_quad_iterator
+DoFHandler<dim>::last_raw_quad (const unsigned int level) const {
return raw_quad_iterator (tria,
tria->last_raw_quad(level)->level(),
tria->last_raw_quad(level)->index(),
-DoFHandler<2>::quad_iterator
-DoFHandler<2>::last_quad (const unsigned int level) const {
+DoFHandler<1>::quad_iterator
+DoFHandler<1>::last_quad (const unsigned int) const {
+ Assert (false, ExcNotImplemented());
+ return 0;
+};
+
+
+
+template <int dim>
+typename DoFHandler<dim>::quad_iterator
+DoFHandler<dim>::last_quad (const unsigned int level) const {
return quad_iterator (tria,
tria->last_quad(level)->level(),
tria->last_quad(level)->index(),
-DoFHandler<2>::active_quad_iterator
-DoFHandler<2>::last_active_quad (const unsigned int level) const {
+DoFHandler<1>::active_quad_iterator
+DoFHandler<1>::last_active_quad (const unsigned int) const {
+ Assert (false, ExcNotImplemented());
+ return 0;
+};
+
+
+
+template <int dim>
+typename DoFHandler<dim>::active_quad_iterator
+DoFHandler<dim>::last_active_quad (const unsigned int level) const {
return active_quad_iterator (tria,
tria->last_active_quad(level)->level(),
tria->last_active_quad(level)->index(),
-DoFHandler<1>::raw_line_iterator
-DoFHandler<1>::last_raw_line () const {
+template <int dim>
+typename DoFHandler<dim>::raw_line_iterator
+DoFHandler<dim>::last_raw_line () const {
return last_raw_line (levels.size()-1);
};
-DoFHandler<2>::raw_line_iterator
-DoFHandler<2>::last_raw_line () const {
- return last_raw_line (levels.size()-1);
+DoFHandler<1>::raw_quad_iterator
+DoFHandler<1>::last_raw_quad () const {
+ Assert (false, ExcNotImplemented());
+ return 0;
};
-DoFHandler<2>::raw_quad_iterator
-DoFHandler<2>::last_raw_quad () const {
+template <int dim>
+typename DoFHandler<dim>::raw_quad_iterator
+DoFHandler<dim>::last_raw_quad () const {
return last_raw_quad (levels.size()-1);
};
-DoFHandler<1>::line_iterator
-DoFHandler<1>::last_line () const {
+template <int dim>
+typename DoFHandler<dim>::line_iterator
+DoFHandler<dim>::last_line () const {
return last_line (levels.size()-1);
};
-DoFHandler<2>::line_iterator
-DoFHandler<2>::last_line () const {
- return last_line (levels.size()-1);
+DoFHandler<1>::quad_iterator
+DoFHandler<1>::last_quad () const {
+ Assert (false, ExcNotImplemented());
+ return 0;
};
-DoFHandler<2>::quad_iterator
-DoFHandler<2>::last_quad () const {
+template <int dim>
+typename DoFHandler<dim>::quad_iterator
+DoFHandler<dim>::last_quad () const {
return last_quad (levels.size()-1);
};
-DoFHandler<1>::active_line_iterator
-DoFHandler<1>::last_active_line () const {
+template <int dim>
+typename DoFHandler<dim>::active_line_iterator
+DoFHandler<dim>::last_active_line () const {
return last_active_line (levels.size()-1);
};
-DoFHandler<2>::active_line_iterator
-DoFHandler<2>::last_active_line () const {
- return last_active_line (levels.size()-1);
+
+DoFHandler<1>::active_quad_iterator
+DoFHandler<1>::last_active_quad () const {
+ Assert (false, ExcNotImplemented());
+ return 0;
};
-DoFHandler<2>::active_quad_iterator
-DoFHandler<2>::last_active_quad () const {
+template <int dim>
+typename DoFHandler<dim>::active_quad_iterator
+DoFHandler<dim>::last_active_quad () const {
return last_active_quad (levels.size()-1);
};
-void FEValues<1>::reinit (const Triangulation<1>::cell_iterator &cell,
- const FiniteElement<1> &fe) {
- const unsigned int dim=1;
+template <int dim>
+void FEValues<dim>::reinit (const typename Triangulation<dim>::cell_iterator &cell,
+ const FiniteElement<dim> &fe) {
// fill jacobi matrices and real
//quadrature points
fe.fill_fe_values (cell,
-void FEValues<2>::reinit (const Triangulation<2>::cell_iterator &cell,
- const FiniteElement<2> &fe) {
- const unsigned int dim=2;
- // fill jacobi matrices and real
- //quadrature points
- fe.fill_fe_values (cell,
- unit_quadrature_points,
- jacobi_matrices,
- quadrature_points);
- // compute gradients on real element
- for (unsigned int i=0; i<fe.total_dofs; ++i)
- for (unsigned int j=0; j<n_quadrature_points; ++j)
- for (unsigned int s=0; s<dim; ++s)
- {
- shape_gradients[i][j](s) = 0;
-
- // (grad psi)_s =
- // (grad_{\xi\eta})_b J_{bs}
- // with J_{bs}=(d\xi_b)/(dx_s)
- for (unsigned int b=0; b<dim; ++b)
- shape_gradients[i][j](s)
- +=
- unit_shape_gradients[i][j](b) * jacobi_matrices[j](b,s);
- };
-
- // refer to the general doc for
- // why we take the inverse of the
- // determinant
- for (unsigned int i=0; i<n_quadrature_points; ++i)
- JxW_values[i] = weights[i] / jacobi_matrices[i].determinant();
-};
-void FiniteElementBase<1>::fill_fe_values (const Triangulation<1>::cell_iterator &,
- const vector<Point<1> > &,
- vector<dFMatrix> &,
- vector<Point<1> > &) const {
- Assert (false, ExcPureFunctionCalled());
-};
-
-
-void FiniteElementBase<2>::fill_fe_values (const Triangulation<2>::cell_iterator &,
- const vector<Point<2> > &,
- vector<dFMatrix> &,
- vector<Point<2> > &) const {
+template <int dim>
+void FiniteElementBase<dim>::fill_fe_values (const typename Triangulation<dim>::cell_iterator &,
+ const vector<Point<dim> > &,
+ vector<dFMatrix> &,
+ vector<Point<dim> > &) const {
Assert (false, ExcPureFunctionCalled());
};
vector<dFMatrix> &jacobians,
vector<Point<1> > &points) const {
// local mesh width
- double h=(cell->vertex(1)(0) - cell->vertex(0)(0));
+ const double h=(cell->vertex(1)(0) - cell->vertex(0)(0));
for (unsigned int i=0; i<points.size(); ++i)
{
for (unsigned int j=0; j<4; ++j)
cells[0].push_back (cell_vertices[0][j]);
- create_triangulation (vector<Point<2> >(&vertices[0], &vertices[8]),
+ create_triangulation (vector<Point<2> >(&vertices[0], &vertices[4]),
cells);
};
-TriaDimensionInfo<1>::raw_line_iterator
-Triangulation<1>::begin_raw_line (unsigned int level) const {
+template <int dim>
+typename TriaDimensionInfo<dim>::raw_line_iterator
+Triangulation<dim>::begin_raw_line (unsigned int level) const {
Assert (level<levels.size(),
ExcInvalidLevel(level));
// does this level have cells?
if (levels[level]->lines.lines.size() == 0)
- return end();
+ return end_line();
- return raw_cell_iterator ((Triangulation<1>*)this,
+ return raw_line_iterator ((Triangulation<dim>*)this,
level,
0);
};
-TriaDimensionInfo<2>::raw_line_iterator
-Triangulation<2>::begin_raw_line (unsigned int level) const {
- Assert (level<levels.size(),
- ExcInvalidLevel(level));
- // does this level have cells?
- if (levels[level]->lines.lines.size() == 0)
- return end_line();
-
- return raw_line_iterator ((Triangulation<2>*)this,
- level,
- 0);
+TriaDimensionInfo<1>::raw_quad_iterator
+Triangulation<1>::begin_raw_quad (unsigned int) const {
+ Assert (false, ExcFunctionNotUseful());
+ return 0;
};
-TriaDimensionInfo<2>::raw_quad_iterator
-Triangulation<2>::begin_raw_quad (unsigned int level) const {
+template <int dim>
+typename TriaDimensionInfo<dim>::raw_quad_iterator
+Triangulation<dim>::begin_raw_quad (unsigned int level) const {
Assert (level<levels.size(),
ExcInvalidLevel(level));
if (levels[level]->quads.quads.size() == 0)
return end_quad();
- return raw_quad_iterator ((Triangulation<2>*)this,
+ return raw_quad_iterator ((Triangulation<dim>*)this,
level,
0);
};
-TriaDimensionInfo<1>::line_iterator
-Triangulation<1>::begin_line (unsigned int level) const {
+
+template <int dim>
+typename TriaDimensionInfo<dim>::line_iterator
+Triangulation<dim>::begin_line (unsigned int level) const {
// level is checked in begin_raw
raw_line_iterator ri = begin_raw_line (level);
if (ri.state() != valid)
-TriaDimensionInfo<2>::line_iterator
-Triangulation<2>::begin_line (unsigned int level) const {
- // level is checked in begin_raw
- raw_line_iterator ri = begin_raw_line (level);
- if (ri.state() != valid)
- return ri;
- while (ri->used() == false)
- if ((++ri).state() != valid)
- return ri;
- return ri;
+
+TriaDimensionInfo<1>::quad_iterator
+Triangulation<1>::begin_quad (unsigned int) const {
+ Assert (false, ExcFunctionNotUseful());
+ return 0;
};
-TriaDimensionInfo<2>::quad_iterator
-Triangulation<2>::begin_quad (unsigned int level) const {
+template <int dim>
+typename TriaDimensionInfo<dim>::quad_iterator
+Triangulation<dim>::begin_quad (unsigned int level) const {
// level is checked in begin_raw
raw_quad_iterator ri = begin_raw_quad (level);
if (ri.state() != valid)
-TriaDimensionInfo<1>::active_line_iterator
-Triangulation<1>::begin_active_line (unsigned int level) const {
+
+template <int dim>
+typename TriaDimensionInfo<dim>::active_line_iterator
+Triangulation<dim>::begin_active_line (unsigned int level) const {
// level is checked in begin_raw
line_iterator i = begin_line (level);
if (i.state() != valid)
-TriaDimensionInfo<2>::active_line_iterator
-Triangulation<2>::begin_active_line (unsigned int level) const {
- // level is checked in begin_raw
- line_iterator i = begin_line (level);
- if (i.state() != valid)
- return i;
- while (i->has_children())
- if ((++i).state() != valid)
- return i;
- return i;
+TriaDimensionInfo<1>::active_quad_iterator
+Triangulation<1>::begin_active_quad (unsigned int) const {
+ Assert (false, ExcFunctionNotUseful());
+ return 0;
};
-TriaDimensionInfo<2>::active_quad_iterator
-Triangulation<2>::begin_active_quad (unsigned int level) const {
+template <int dim>
+typename TriaDimensionInfo<dim>::active_quad_iterator
+Triangulation<dim>::begin_active_quad (unsigned int level) const {
// level is checked in begin_raw
quad_iterator i = begin(level);
if (i.state() != valid)
-TriaDimensionInfo<1>::raw_line_iterator
-Triangulation<1>::end_line () const {
- return raw_line_iterator ((Triangulation<1>*)this,
+template <int dim>
+typename TriaDimensionInfo<dim>::raw_line_iterator
+Triangulation<dim>::end_line () const {
+ return raw_line_iterator ((Triangulation<dim>*)this,
-1,
-1);
};
-TriaDimensionInfo<2>::raw_line_iterator
-Triangulation<2>::end_line () const {
- return raw_line_iterator ((Triangulation<2>*)this,
- -1,
- -1);
+
+TriaDimensionInfo<1>::raw_quad_iterator
+Triangulation<1>::end_quad () const {
+ Assert (false, ExcFunctionNotUseful());
+ return 0;
};
-TriaDimensionInfo<2>::raw_quad_iterator
-Triangulation<2>::end_quad () const {
- return raw_quad_iterator ((Triangulation<2>*)this,
+template <int dim>
+typename TriaDimensionInfo<dim>::raw_quad_iterator
+Triangulation<dim>::end_quad () const {
+ return raw_quad_iterator ((Triangulation<dim>*)this,
-1,
-1);
};
-TriaDimensionInfo<1>::raw_line_iterator
-Triangulation<1>::last_raw_line (const unsigned int level) const {
+template <int dim>
+typename TriaDimensionInfo<dim>::raw_line_iterator
+Triangulation<dim>::last_raw_line (const unsigned int level) const {
Assert (level<levels.size(),
ExcInvalidLevel(level));
- return raw_line_iterator ((Triangulation<1>*)this,
+ return raw_line_iterator ((Triangulation<dim>*)this,
level,
levels[level]->lines.lines.size()-1);
};
-TriaDimensionInfo<2>::raw_line_iterator
-Triangulation<2>::last_raw_line (const unsigned int level) const {
- Assert (level<levels.size(),
- ExcInvalidLevel(level));
-
- return raw_line_iterator ((Triangulation<2>*)this,
- level,
- levels[level]->lines.lines.size()-1);
+
+
+TriaDimensionInfo<1>::raw_quad_iterator
+Triangulation<1>::last_raw_quad (const unsigned int) const {
+ Assert (false, ExcFunctionNotUseful());
+ return 0;
};
-TriaDimensionInfo<2>::raw_quad_iterator
-Triangulation<2>::last_raw_quad (const unsigned int level) const {
+template <int dim>
+typename TriaDimensionInfo<dim>::raw_quad_iterator
+Triangulation<dim>::last_raw_quad (const unsigned int level) const {
Assert (level<levels.size(),
ExcInvalidLevel(level));
- return raw_quad_iterator ((Triangulation<2>*)this,
+ return raw_quad_iterator ((Triangulation<dim>*)this,
level,
levels[level]->quads.quads.size()-1);
};
-TriaDimensionInfo<1>::raw_line_iterator
-Triangulation<1>::last_raw_line () const {
+template <int dim>
+typename TriaDimensionInfo<dim>::raw_line_iterator
+Triangulation<dim>::last_raw_line () const {
return last_raw_line (levels.size()-1);
};
-TriaDimensionInfo<2>::raw_line_iterator
-Triangulation<2>::last_raw_line () const {
- return last_raw_line (levels.size()-1);
+
+TriaDimensionInfo<1>::raw_quad_iterator
+Triangulation<1>::last_raw_quad () const {
+ Assert (false, ExcFunctionNotUseful());
+ return 0;
};
-TriaDimensionInfo<2>::raw_quad_iterator
-Triangulation<2>::last_raw_quad () const {
+template <int dim>
+typename TriaDimensionInfo<dim>::raw_quad_iterator
+Triangulation<dim>::last_raw_quad () const {
return last_raw_quad (levels.size()-1);
};
-TriaDimensionInfo<1>::line_iterator
-Triangulation<1>::last_line (const unsigned int level) const {
+template <int dim>
+typename TriaDimensionInfo<dim>::line_iterator
+Triangulation<dim>::last_line (const unsigned int level) const {
// level is checked in begin_raw
raw_line_iterator ri = last_raw_line(level);
if (ri->used()==true)
-TriaDimensionInfo<2>::line_iterator
-Triangulation<2>::last_line (const unsigned int level) const {
- // level is checked in begin_raw
- raw_line_iterator ri = last_raw_line(level);
- if (ri->used()==true)
- return ri;
- while ((--ri).state() == valid)
- if (ri->used()==true)
- return ri;
- return ri;
+
+TriaDimensionInfo<1>::quad_iterator
+Triangulation<1>::last_quad (const unsigned int) const {
+ Assert (false, ExcFunctionNotUseful());
+ return 0;
};
-TriaDimensionInfo<2>::quad_iterator
-Triangulation<2>::last_quad (const unsigned int level) const {
+template <int dim>
+typename TriaDimensionInfo<dim>::quad_iterator
+Triangulation<dim>::last_quad (const unsigned int level) const {
// level is checked in begin_raw
raw_quad_iterator ri = last_raw_quad(level);
if (ri->used()==true)
-TriaDimensionInfo<1>::line_iterator
-Triangulation<1>::last_line () const {
+template <int dim>
+typename TriaDimensionInfo<dim>::line_iterator
+Triangulation<dim>::last_line () const {
return last_line (levels.size()-1);
};
-TriaDimensionInfo<2>::line_iterator
-Triangulation<2>::last_line () const {
- return last_line (levels.size()-1);
+TriaDimensionInfo<1>::quad_iterator
+Triangulation<1>::last_quad () const {
+ Assert (false, ExcFunctionNotUseful());
+ return 0;
};
-TriaDimensionInfo<2>::quad_iterator
-Triangulation<2>::last_quad () const {
+
+template <int dim>
+typename TriaDimensionInfo<dim>::quad_iterator
+Triangulation<dim>::last_quad () const {
return last_quad (levels.size()-1);
};
-TriaDimensionInfo<1>::active_line_iterator
-Triangulation<1>::last_active_line (const unsigned int level) const {
+
+template <int dim>
+typename TriaDimensionInfo<dim>::active_line_iterator
+Triangulation<dim>::last_active_line (const unsigned int level) const {
// level is checked in begin_raw
line_iterator i = last_line(level);
if (i->has_children()==false)
-TriaDimensionInfo<2>::active_line_iterator
-Triangulation<2>::last_active_line (const unsigned int level) const {
- // level is checked in begin_raw
- line_iterator i = last_line(level);
- if (i->has_children()==false)
- return i;
- while ((--i).state() == valid)
- if (i->has_children()==false)
- return i;
- return i;
+TriaDimensionInfo<1>::active_quad_iterator
+Triangulation<1>::last_active_quad (const unsigned int) const {
+ Assert (false, ExcFunctionNotUseful());
+ return 0;
};
-TriaDimensionInfo<2>::active_quad_iterator
-Triangulation<2>::last_active_quad (const unsigned int level) const {
+
+template <int dim>
+typename TriaDimensionInfo<dim>::active_quad_iterator
+Triangulation<dim>::last_active_quad (const unsigned int level) const {
// level is checked in begin_raw
quad_iterator i = last_quad(level);
if (i->has_children()==false)
-TriaDimensionInfo<1>::active_line_iterator
-Triangulation<1>::last_active_line () const {
+template <int dim>
+typename TriaDimensionInfo<dim>::active_line_iterator
+Triangulation<dim>::last_active_line () const {
return last_active_line (levels.size()-1);
};
-TriaDimensionInfo<2>::active_line_iterator
-Triangulation<2>::last_active_line () const {
- return last_active_line (levels.size()-1);
+TriaDimensionInfo<1>::active_quad_iterator
+Triangulation<1>::last_active_quad () const {
+ Assert (false, ExcFunctionNotUseful());
+ return 0;
};
-TriaDimensionInfo<2>::active_quad_iterator
-Triangulation<2>::last_active_quad () const {
+
+template <int dim>
+typename TriaDimensionInfo<dim>::active_quad_iterator
+Triangulation<dim>::last_active_quad () const {
return last_active_quad (levels.size()-1);
};
-void Equation<1>::assemble (dFMatrix &,
- dVector &,
- const FEValues<1> &,
- const Triangulation<1>::cell_iterator &) const {
- Assert (false, ExcPureVirtualFunctionCalled());
-};
-
-
-
-void Equation<2>::assemble (dFMatrix &,
- dVector &,
- const FEValues<2> &,
- const Triangulation<2>::cell_iterator &) const {
- Assert (false, ExcPureVirtualFunctionCalled());
-};
-
-
-
-void Equation<1>::assemble (dFMatrix &,
- const FEValues<1> &,
- const Triangulation<1>::cell_iterator &) const {
- Assert (false, ExcPureVirtualFunctionCalled());
-};
-
-
-
-void Equation<2>::assemble (dFMatrix &,
- const FEValues<2> &,
- const Triangulation<2>::cell_iterator &) const {
+template <int dim>
+void Equation<dim>::assemble (dFMatrix &,
+ dVector &,
+ const FEValues<dim> &,
+ const typename Triangulation<dim>::cell_iterator &) const {
Assert (false, ExcPureVirtualFunctionCalled());
};
-void Equation<1>::assemble (dVector &,
- const FEValues<1> &,
- const Triangulation<1>::cell_iterator &) const {
+template <int dim>
+void Equation<dim>::assemble (dFMatrix &,
+ const FEValues<dim> &,
+ const typename Triangulation<dim>::cell_iterator &) const {
Assert (false, ExcPureVirtualFunctionCalled());
};
-void Equation<2>::assemble (dVector &,
- const FEValues<2> &,
- const Triangulation<2>::cell_iterator &) const {
+template <int dim>
+void Equation<dim>::assemble (dVector &,
+ const FEValues<dim> &,
+ const typename Triangulation<dim>::cell_iterator &) const {
Assert (false, ExcPureVirtualFunctionCalled());
};
grid/grid_test: grid/grid_test.o $(LIBFILES)
@echo ================= Linking $@
- @$(CXX) $(CXXFLAGS) -o $@ $< $(LIBS)
+ @$(CXX) $(CXXFLAGS) -o $@ $< $(LIBS.g)
dof/dof_test: dof/dof_test.o $(LIBFILES.g)
@echo ================= Linking $@