#include <grid/tria_iterator.h>
#include <fe/fe.h>
#include <fe/fe_values.h>
-#include <fe/mapping_q1.h>
-
-//TODO:[RH,GK] Replace global by local object; better: have two functions, or by default arg
-static const MappingQ1<deal_II_dimension> mapping;
#ifdef DEAL_II_USE_MT
#include <base/thread_management.h>
{
QTrapez<1> q_trapez;
QIterated<dim> patch_points (q_trapez, data.n_subdivisions);
-
- FEValues<dim> fe_patch_values(mapping, dofs->get_fe(),
+
+ // this constructor implicitely
+ // uses a MappingQ1 mapping
+ FEValues<dim> fe_patch_values(dofs->get_fe(),
patch_points,
update_values);
#include <fe/fe.h>
#include <fe/fe_values.h>
-#include <fe/mapping_q1.h>
-
-//TODO:[RH,GK] Replace global by local object; better: have two functions, or by default arg
-static const MappingQ1<deal_II_dimension> mapping;
-
#ifdef DEAL_II_USE_MT
#include <base/thread_management.h>
#endif
QTrapez<1> q_trapez;
QIterated<dim-1> patch_points (q_trapez, data.n_subdivisions);
- FEFaceValues<dim> fe_patch_values(mapping,
- dofs->get_fe(),
+ // this constructor implicitely
+ // uses a MappingQ1 mapping
+ FEFaceValues<dim> fe_patch_values(dofs->get_fe(),
patch_points,
update_values);
#include <dofs/dof_handler.h>
#include <dofs/dof_accessor.h>
#include <grid/tria_iterator.h>
-#include <fe/mapping_q1.h>
#include <fe/fe.h>
#include <fe/fe_values.h>
-//TODO:[RH,GK] Replace global by local object; better: have two functions, or by default arg
-static const MappingQ1<deal_II_dimension> mapping;
-
#ifdef DEAL_II_USE_MT
#include <base/thread_management.h>
#endif
{
QTrapez<1> q_trapez;
QIterated<dim> patch_points (q_trapez, data.n_subdivisions);
-
- FEValues<dim> fe_patch_values(mapping,
- dofs->get_fe(),
+
+ // this constructor implicitely
+ // uses a MappingQ1 mapping
+ FEValues<dim> fe_patch_values(dofs->get_fe(),
patch_points,
update_values);
#include <grid/tria_iterator.h>
#include <fe/fe.h>
#include <fe/fe_values.h>
-#include <fe/mapping_q1.h>
-
// if necessary try to work around a bug in the IBM xlC compiler
#endif
-//TODO:[RH,GK] Replace global by local object; better: have two functions, or by default arg
-static const MappingQ1<deal_II_dimension> mapping;
-
template <int dim>
unsigned int
// cell to these points
QTrapez<1> q_trapez;
QIterated<dim> patch_points (q_trapez, n_subdivisions);
- FEValues<dim> fe_patch_values (mapping, dof_handler->get_fe(),
+
+ // this constructor implicitely
+ // uses a MappingQ1 mapping
+ FEValues<dim> fe_patch_values (dof_handler->get_fe(),
patch_points,
update_values);
const unsigned int n_q_points = patch_points.n_quadrature_points;