normals_3.exe : normals_3.g.$(OBJEXT) $(libraries)
normals_4.exe : normals_4.g.$(OBJEXT) $(libraries)
mapping_cartesian_1.exe : mapping_cartesian_1.g.$(OBJEXT) $(libraries)
+mapping_q4_3d.exe : mapping_q4_3d.g.$(OBJEXT) $(libraries)
q_points.exe : q_points.g.$(OBJEXT) $(libraries)
q_point_sum_1.exe : q_point_sum_1.g.$(OBJEXT) $(libraries)
q_point_sum_2.exe : q_point_sum_2.g.$(OBJEXT) $(libraries)
q_point_sum_1 q_point_sum_2 q_point_sum_3 q_point_sum_4 \
volume_1 volume_2 volume_3 volume_4 \
mapping_cartesian_1 \
+ mapping_q4_3d \
q_points
############################################################
--- /dev/null
+//---------------------------- mapping_q4_3d.cc ---------------------------
+// $Id$
+// Version: $Name$
+//
+// Copyright (C) 2003 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//---------------------------- mapping_q4_3d.cc ---------------------------
+
+
+// this test used to fail somewhere in the 3d parts of MappingQ for
+// fourth order mappings
+
+#include <base/logstream.h>
+#include <base/quadrature_lib.h>
+#include <grid/tria.h>
+#include <grid/tria_accessor.h>
+#include <grid/tria_iterator.h>
+#include <grid/grid_generator.h>
+#include <dofs/dof_handler.h>
+#include <fe/fe_q.h>
+#include <fe/fe_values.h>
+#include <fe/mapping_q.h>
+
+#include <fstream>
+
+
+int main ()
+{
+ std::ofstream logfile("mapping_q4_3d.output");
+ deallog.attach(logfile);
+ deallog.depth_console(0);
+
+ Triangulation<3> coarse_grid;
+ GridGenerator::hyper_cube (coarse_grid, -1, 1);
+ coarse_grid.refine_global (1);
+
+ FE_Q<3> fe(1);
+ DoFHandler<3> dof_handler (coarse_grid);
+ dof_handler.distribute_dofs (fe);
+
+ MappingQ<3> mapping(4);
+ QGauss3<2> q_face;
+ FEFaceValues<3> fe_face_values (mapping, fe, q_face,
+ update_normal_vectors | update_JxW_values);
+ fe_face_values.reinit (dof_handler.begin_active(), 0);
+}
+
+