#include <grid/tria_iterator.h>
#include <fe/fe.h>
#include <fe/fe_values.h>
+#include <fe/mapping_q1.h>
#ifdef DEAL_II_USE_MT
#include <base/thread_management.h>
QTrapez<1> q_trapez;
QIterated<dim> patch_points (q_trapez, data.n_subdivisions);
- // this constructor implicitely
- // uses a MappingQ1 mapping
- FEValues<dim> fe_patch_values(dofs->get_fe(),
- patch_points,
- update_values);
+ // since most output formats can't
+ // handle cells that are not
+ // transformed using a Q1 mapping,
+ // we don't support anything else
+ // as well
+ static const MappingQ1<dim> mapping;
+ FEValues<dim> fe_patch_values (mapping, dofs->get_fe(),
+ patch_points, update_values);
const unsigned int n_q_points = patch_points.n_quadrature_points;
#include <grid/tria_iterator.h>
#include <fe/fe.h>
#include <fe/fe_values.h>
+#include <fe/mapping_q1.h>
#ifdef DEAL_II_USE_MT
#include <base/thread_management.h>
QTrapez<1> q_trapez;
QIterated<dim-1> patch_points (q_trapez, data.n_subdivisions);
- // this constructor implicitely
- // uses a MappingQ1 mapping
- FEFaceValues<dim> fe_patch_values(dofs->get_fe(),
- patch_points,
- update_values);
+ // since most output formats can't
+ // handle cells that are not
+ // transformed using a Q1 mapping,
+ // we don't support anything else
+ // as well
+ static const MappingQ1<dim> mapping;
+ FEFaceValues<dim> fe_patch_values (mapping, dofs->get_fe(),
+ patch_points, update_values);
const unsigned int n_q_points = patch_points.n_quadrature_points;
#include <grid/tria_iterator.h>
#include <fe/fe.h>
#include <fe/fe_values.h>
+#include <fe/mapping_q1.h>
#ifdef DEAL_II_USE_MT
#include <base/thread_management.h>
QTrapez<1> q_trapez;
QIterated<dim> patch_points (q_trapez, data.n_subdivisions);
- // this constructor implicitely
- // uses a MappingQ1 mapping
- FEValues<dim> fe_patch_values(dofs->get_fe(),
- patch_points,
- update_values);
+ // since most output formats can't
+ // handle cells that are not
+ // transformed using a Q1 mapping,
+ // we don't support anything else
+ // as well
+ static const MappingQ1<dim> mapping;
+ FEValues<dim> fe_patch_values (mapping, dofs->get_fe(),
+ patch_points, update_values);
const unsigned int n_patches_per_circle = data.n_patches_per_circle;
#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
QTrapez<1> q_trapez;
QIterated<dim> patch_points (q_trapez, n_subdivisions);
- // this constructor implicitely
- // uses a MappingQ1 mapping
- FEValues<dim> fe_patch_values (dof_handler->get_fe(),
- patch_points,
- update_values);
+ // since most output formats can't
+ // handle cells that are not
+ // transformed using a Q1 mapping,
+ // we don't support anything else
+ // as well
+ static const MappingQ1<dim> mapping;
+ FEValues<dim> fe_patch_values (mapping, dof_handler->get_fe(),
+ patch_points, update_values);
const unsigned int n_q_points = patch_points.n_quadrature_points;
std::vector<double> patch_values (n_q_points);
std::vector<Vector<double> > patch_values_system (n_q_points,