]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Since output formats cant write more than Q1, make that explicit in the creation...
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Wed, 11 Apr 2001 13:53:08 +0000 (13:53 +0000)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Wed, 11 Apr 2001 13:53:08 +0000 (13:53 +0000)
git-svn-id: https://svn.dealii.org/trunk@4439 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/source/numerics/data_out.cc
deal.II/deal.II/source/numerics/data_out_faces.cc
deal.II/deal.II/source/numerics/data_out_rotation.cc
deal.II/deal.II/source/numerics/data_out_stack.cc

index 396b35d4cd7b3eb254aacbef2a42564199f31b09..6593b8c135e0e94e1bd0a98e18d0be8474d0941a 100644 (file)
@@ -22,6 +22,7 @@
 #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>
@@ -302,11 +303,14 @@ void DataOut<dim>::build_some_patches (Data data)
   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;
   
index ae028c8140384705ed4deea9fd113b4fa352a983..d50d93209555ecc03962f1bb892ca6b5f4409c27 100644 (file)
@@ -21,6 +21,7 @@
 #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>
@@ -36,11 +37,14 @@ void DataOutFaces<dim>::build_some_patches (Data data)
   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;
   
index 17f6b43f33d8580b62b6f113722de19eb285987c..4aa135e855406ce91f68fa889c2e25fd71fec0c0 100644 (file)
@@ -21,6 +21,7 @@
 #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>
@@ -41,11 +42,14 @@ void DataOutRotation<dim>::build_some_patches (Data data)
   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;
 
index 4a2712bfe89d476987dc5bdd0f238afa75c8baaa..62de438914d53e7314c9e096e6d27e5e1ef6a9de 100644 (file)
@@ -21,6 +21,8 @@
 #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
@@ -219,11 +221,14 @@ void DataOutStack<dim>::build_patches (const unsigned int n_subdivisions)
   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,

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.