#include <numerics/error_estimator.h>
// Now for the C++ standard headers:
+#include <iostream>
#include <fstream>
#include <list>
void assemble_linear_system (LinearSystem &linear_system);
- void assemble_matrix (LinearSystem &linear_system,
- const DoFHandler<dim>::active_cell_iterator &begin_cell,
- const DoFHandler<dim>::active_cell_iterator &end_cell,
- Threads::ThreadMutex &mutex) const ;
+ void assemble_matrix (LinearSystem &linear_system,
+ const typename DoFHandler<dim>::active_cell_iterator &begin_cell,
+ const typename DoFHandler<dim>::active_cell_iterator &end_cell,
+ Threads::ThreadMutex &mutex) const ;
};
template <int dim>
void
- Solver<dim>::assemble_matrix (LinearSystem &linear_system,
- const DoFHandler<dim>::active_cell_iterator &begin_cell,
- const DoFHandler<dim>::active_cell_iterator &end_cell,
- Threads::ThreadMutex &mutex) const
+ Solver<dim>::assemble_matrix (LinearSystem &linear_system,
+ const typename DoFHandler<dim>::active_cell_iterator &begin_cell,
+ const typename DoFHandler<dim>::active_cell_iterator &end_cell,
+ Threads::ThreadMutex &mutex) const
{
//TODO: adaptive
QGauss4<dim> quadrature_formula;