]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
More and more thorough testing.
authorwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 3 Sep 2002 17:26:13 +0000 (17:26 +0000)
committerwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 3 Sep 2002 17:26:13 +0000 (17:26 +0000)
git-svn-id: https://svn.dealii.org/trunk@6354 0785d39b-7218-0410-832d-ea1e28bc413d

tests/fe/nedelec.cc
tests/fe/nedelec.checked

index 88597efe93495fbcbdf872939d240d38985d3337..21ddad0ec241ba7f0f6c0c6f5ccade83e39ad928 100644 (file)
 
 #define PRECISION 2
 
-char fname[50];
+
+Point<2> stretch_coordinates (const Point<2> p)
+{
+  return Point<2>(2*p(0), p(1));
+};
+
+
+
+Point<2> tilt_coordinates (const Point<2> p)
+{
+  return Point<2>(p(0)+p(1), p(1));
+};
+
 
 void
 transform_grid (Triangulation<2> &tria,
@@ -28,20 +40,33 @@ transform_grid (Triangulation<2> &tria,
 {
   switch (transform)
     {
+                                       // first round: take
+                                       // original grid
       case 0:
-                                            // first round: take
-                                            // original grid
            break;
-      case 1:
+            
                                             // second round: rotate
                                             // triangulation
+      case 1:
            GridTools::rotate (3.14159265358/2, tria);
            break;
-      case 2:
+
                                             // third round: inflate
                                             // by a factor of 2
+      case 2:
            GridTools::scale (2, tria);
            break;
+
+                                             // third round: scale
+                                             // back, rotate back,
+                                             // stretch
+      case 3:
+            GridTools::scale (.5, tria);
+           GridTools::rotate (-3.14159265358/2, tria);
+            GridTools::transform (&stretch_coordinates, tria);
+            
+           break;
+            
       default:
            Assert (false, ExcNotImplemented());
     };
@@ -61,8 +86,10 @@ plot_shape_functions()
                                   // check the following with a
                                   // number of transformed
                                   // triangulations
-  for (unsigned int transform=0; transform<3; ++transform)
+  for (unsigned int transform=0; transform<4; ++transform)
     {
+      deallog << "GRID TRANSFORMATION " << transform << std::endl;
+      
       transform_grid (tr, transform);
 
       DoFHandler<dim> dof(tr);
@@ -72,26 +99,42 @@ plot_shape_functions()
       QTrapez<1> q_trapez;
       const unsigned int div=2;
       QIterated<dim> q(q_trapez, div);
-      FEValues<dim> fe(fe_ned, q, update_values);
+      FEValues<dim> fe(fe_ned, q, update_values|update_gradients|update_q_points);
       fe.reinit(c);
       
-      unsigned int k=0;
+      unsigned int q_point=0;
       for (unsigned int mz=0;mz<=((dim>2) ? div : 0) ;++mz)
        for (unsigned int my=0;my<=((dim>1) ? div : 0) ;++my)
          for (unsigned int mx=0;mx<=div;++mx)
            {
-             deallog << "q_point(" << k << ")=" << q.point(k)
+             deallog << "q_point(" << q_point << ")=" << fe.quadrature_point(q_point)
                      << std::endl;
              
              for (unsigned int i=0;i<fe_ned.dofs_per_cell;++i)
                {
-                 deallog << "  shape_function(" << i << ")=";
+                 deallog << "  phi(" << i << ") = [";
                  for (unsigned int c=0; c<fe.get_fe().n_components(); ++c)
-                   deallog << " " << fe.shape_value_component(i,k,c);
-                 deallog << std::endl;
+                   deallog << " " << fe.shape_value_component(i,q_point,c);
+                 deallog << " ]" << std::endl;
                };
-             k++;
+             for (unsigned int i=0;i<fe_ned.dofs_per_cell;++i)
+               {
+                 deallog << "  grad phi(" << i << ") = ";
+                 for (unsigned int c=0; c<dim; ++c)
+                    {
+                      deallog << "[";
+                      for (unsigned int d=0; d<dim; ++d)
+                        deallog << " " << fe.shape_grad_component(i,q_point,c)[d];
+                      deallog << " ]" << std::endl;
+                      if (c != dim-1)
+                        deallog << "                ";
+                    };
+               };
+
+              q_point++;
            }
+
+      deallog << std::endl;
     }
 };
 
index d81ee551e260baa7c6714cc2f8e6d2a4c305d74c..f9e34732b6a60444b854b68c48020a11c5cab062 100644 (file)
 
+DEAL::GRID TRANSFORMATION 0
 DEAL::q_point(0)=0.00 0.00
-DEAL::  shape_function(0)= 1.00 0.00
-DEAL::  shape_function(1)= 0.00 0.00
-DEAL::  shape_function(2)= 0.00 0.00
-DEAL::  shape_function(3)= 0.00 1.00
+DEAL::  phi(0) = [ 1.00 0.00 ]
+DEAL::  phi(1) = [ 0.00 0.00 ]
+DEAL::  phi(2) = [ 0.00 0.00 ]
+DEAL::  phi(3) = [ 0.00 1.00 ]
+DEAL::  grad phi(0) = [ 0.00 -1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 0.00 ]
+DEAL::                [ 1.00 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 0.00 ]
+DEAL::                [ -1.00 0.00 ]
 DEAL::q_point(1)=0.00 0.50
-DEAL::  shape_function(0)= 0.50 0.00
-DEAL::  shape_function(1)= 0.00 0.00
-DEAL::  shape_function(2)= 0.50 0.00
-DEAL::  shape_function(3)= 0.00 1.00
+DEAL::  phi(0) = [ 0.50 0.00 ]
+DEAL::  phi(1) = [ 0.00 0.00 ]
+DEAL::  phi(2) = [ 0.50 0.00 ]
+DEAL::  phi(3) = [ 0.00 1.00 ]
+DEAL::  grad phi(0) = [ 0.00 -1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 0.00 ]
+DEAL::                [ 1.00 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 0.00 ]
+DEAL::                [ -1.00 0.00 ]
 DEAL::q_point(2)=0.00 1.00
-DEAL::  shape_function(0)= 0.00 0.00
-DEAL::  shape_function(1)= 0.00 0.00
-DEAL::  shape_function(2)= 1.00 0.00
-DEAL::  shape_function(3)= 0.00 1.00
+DEAL::  phi(0) = [ 0.00 0.00 ]
+DEAL::  phi(1) = [ 0.00 0.00 ]
+DEAL::  phi(2) = [ 1.00 0.00 ]
+DEAL::  phi(3) = [ 0.00 1.00 ]
+DEAL::  grad phi(0) = [ 0.00 -1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 0.00 ]
+DEAL::                [ 1.00 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 0.00 ]
+DEAL::                [ -1.00 0.00 ]
 DEAL::q_point(3)=0.50 0.00
-DEAL::  shape_function(0)= 1.00 0.00
-DEAL::  shape_function(1)= 0.00 0.50
-DEAL::  shape_function(2)= 0.00 0.00
-DEAL::  shape_function(3)= 0.00 0.50
+DEAL::  phi(0) = [ 1.00 0.00 ]
+DEAL::  phi(1) = [ 0.00 0.50 ]
+DEAL::  phi(2) = [ 0.00 0.00 ]
+DEAL::  phi(3) = [ 0.00 0.50 ]
+DEAL::  grad phi(0) = [ 0.00 -1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 0.00 ]
+DEAL::                [ 1.00 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 0.00 ]
+DEAL::                [ -1.00 0.00 ]
 DEAL::q_point(4)=0.50 0.50
-DEAL::  shape_function(0)= 0.50 0.00
-DEAL::  shape_function(1)= 0.00 0.50
-DEAL::  shape_function(2)= 0.50 0.00
-DEAL::  shape_function(3)= 0.00 0.50
+DEAL::  phi(0) = [ 0.50 0.00 ]
+DEAL::  phi(1) = [ 0.00 0.50 ]
+DEAL::  phi(2) = [ 0.50 0.00 ]
+DEAL::  phi(3) = [ 0.00 0.50 ]
+DEAL::  grad phi(0) = [ 0.00 -1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 0.00 ]
+DEAL::                [ 1.00 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 0.00 ]
+DEAL::                [ -1.00 0.00 ]
 DEAL::q_point(5)=0.50 1.00
-DEAL::  shape_function(0)= 0.00 0.00
-DEAL::  shape_function(1)= 0.00 0.50
-DEAL::  shape_function(2)= 1.00 0.00
-DEAL::  shape_function(3)= 0.00 0.50
+DEAL::  phi(0) = [ 0.00 0.00 ]
+DEAL::  phi(1) = [ 0.00 0.50 ]
+DEAL::  phi(2) = [ 1.00 0.00 ]
+DEAL::  phi(3) = [ 0.00 0.50 ]
+DEAL::  grad phi(0) = [ 0.00 -1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 0.00 ]
+DEAL::                [ 1.00 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 0.00 ]
+DEAL::                [ -1.00 0.00 ]
 DEAL::q_point(6)=1.00 0.00
-DEAL::  shape_function(0)= 1.00 0.00
-DEAL::  shape_function(1)= 0.00 1.00
-DEAL::  shape_function(2)= 0.00 0.00
-DEAL::  shape_function(3)= 0.00 0.00
+DEAL::  phi(0) = [ 1.00 0.00 ]
+DEAL::  phi(1) = [ 0.00 1.00 ]
+DEAL::  phi(2) = [ 0.00 0.00 ]
+DEAL::  phi(3) = [ 0.00 0.00 ]
+DEAL::  grad phi(0) = [ 0.00 -1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 0.00 ]
+DEAL::                [ 1.00 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 0.00 ]
+DEAL::                [ -1.00 0.00 ]
 DEAL::q_point(7)=1.00 0.50
-DEAL::  shape_function(0)= 0.50 0.00
-DEAL::  shape_function(1)= 0.00 1.00
-DEAL::  shape_function(2)= 0.50 0.00
-DEAL::  shape_function(3)= 0.00 0.00
+DEAL::  phi(0) = [ 0.50 0.00 ]
+DEAL::  phi(1) = [ 0.00 1.00 ]
+DEAL::  phi(2) = [ 0.50 0.00 ]
+DEAL::  phi(3) = [ 0.00 0.00 ]
+DEAL::  grad phi(0) = [ 0.00 -1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 0.00 ]
+DEAL::                [ 1.00 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 0.00 ]
+DEAL::                [ -1.00 0.00 ]
 DEAL::q_point(8)=1.00 1.00
-DEAL::  shape_function(0)= 0.00 0.00
-DEAL::  shape_function(1)= 0.00 1.00
-DEAL::  shape_function(2)= 1.00 0.00
-DEAL::  shape_function(3)= 0.00 0.00
+DEAL::  phi(0) = [ 0.00 0.00 ]
+DEAL::  phi(1) = [ 0.00 1.00 ]
+DEAL::  phi(2) = [ 1.00 0.00 ]
+DEAL::  phi(3) = [ 0.00 0.00 ]
+DEAL::  grad phi(0) = [ 0.00 -1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 0.00 ]
+DEAL::                [ 1.00 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 0.00 ]
+DEAL::                [ -1.00 0.00 ]
+DEAL::
+DEAL::GRID TRANSFORMATION 1
 DEAL::q_point(0)=0.00 0.00
-DEAL::  shape_function(0)= 0.00 1.00
-DEAL::  shape_function(1)= 0.00 0.00
-DEAL::  shape_function(2)= 0.00 0.00
-DEAL::  shape_function(3)= -1.00 0.00
-DEAL::q_point(1)=0.00 0.50
-DEAL::  shape_function(0)= 0.00 0.50
-DEAL::  shape_function(1)= 0.00 0.00
-DEAL::  shape_function(2)= 0.00 0.50
-DEAL::  shape_function(3)= -1.00 0.00
-DEAL::q_point(2)=0.00 1.00
-DEAL::  shape_function(0)= 0.00 0.00
-DEAL::  shape_function(1)= 0.00 0.00
-DEAL::  shape_function(2)= 0.00 1.00
-DEAL::  shape_function(3)= -1.00 0.00
-DEAL::q_point(3)=0.50 0.00
-DEAL::  shape_function(0)= 0.00 1.00
-DEAL::  shape_function(1)= -0.50 0.00
-DEAL::  shape_function(2)= 0.00 0.00
-DEAL::  shape_function(3)= -0.50 0.00
-DEAL::q_point(4)=0.50 0.50
-DEAL::  shape_function(0)= 0.00 0.50
-DEAL::  shape_function(1)= -0.50 0.00
-DEAL::  shape_function(2)= 0.00 0.50
-DEAL::  shape_function(3)= -0.50 0.00
-DEAL::q_point(5)=0.50 1.00
-DEAL::  shape_function(0)= 0.00 0.00
-DEAL::  shape_function(1)= -0.50 0.00
-DEAL::  shape_function(2)= 0.00 1.00
-DEAL::  shape_function(3)= -0.50 0.00
-DEAL::q_point(6)=1.00 0.00
-DEAL::  shape_function(0)= 0.00 1.00
-DEAL::  shape_function(1)= -1.00 0.00
-DEAL::  shape_function(2)= 0.00 0.00
-DEAL::  shape_function(3)= 0.00 0.00
-DEAL::q_point(7)=1.00 0.50
-DEAL::  shape_function(0)= 0.00 0.50
-DEAL::  shape_function(1)= -1.00 0.00
-DEAL::  shape_function(2)= 0.00 0.50
-DEAL::  shape_function(3)= 0.00 0.00
-DEAL::q_point(8)=1.00 1.00
-DEAL::  shape_function(0)= 0.00 0.00
-DEAL::  shape_function(1)= -1.00 0.00
-DEAL::  shape_function(2)= 0.00 1.00
-DEAL::  shape_function(3)= 0.00 0.00
+DEAL::  phi(0) = [ 0.00 1.00 ]
+DEAL::  phi(1) = [ 0.00 0.00 ]
+DEAL::  phi(2) = [ 0.00 0.00 ]
+DEAL::  phi(3) = [ -1.00 0.00 ]
+DEAL::  grad phi(0) = [ 0.00 0.00 ]
+DEAL::                [ 1.00 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 -1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 0.00 ]
+DEAL::                [ -1.00 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::q_point(1)=-0.50 0.00
+DEAL::  phi(0) = [ 0.00 0.50 ]
+DEAL::  phi(1) = [ 0.00 0.00 ]
+DEAL::  phi(2) = [ 0.00 0.50 ]
+DEAL::  phi(3) = [ -1.00 0.00 ]
+DEAL::  grad phi(0) = [ 0.00 0.00 ]
+DEAL::                [ 1.00 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 -1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 0.00 ]
+DEAL::                [ -1.00 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::q_point(2)=-1.00 0.00
+DEAL::  phi(0) = [ 0.00 0.00 ]
+DEAL::  phi(1) = [ 0.00 0.00 ]
+DEAL::  phi(2) = [ 0.00 1.00 ]
+DEAL::  phi(3) = [ -1.00 0.00 ]
+DEAL::  grad phi(0) = [ 0.00 0.00 ]
+DEAL::                [ 1.00 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 -1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 0.00 ]
+DEAL::                [ -1.00 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::q_point(3)=0.00 0.50
+DEAL::  phi(0) = [ 0.00 1.00 ]
+DEAL::  phi(1) = [ -0.50 0.00 ]
+DEAL::  phi(2) = [ 0.00 0.00 ]
+DEAL::  phi(3) = [ -0.50 0.00 ]
+DEAL::  grad phi(0) = [ 0.00 0.00 ]
+DEAL::                [ 1.00 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 -1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 0.00 ]
+DEAL::                [ -1.00 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::q_point(4)=-0.50 0.50
+DEAL::  phi(0) = [ 0.00 0.50 ]
+DEAL::  phi(1) = [ -0.50 0.00 ]
+DEAL::  phi(2) = [ 0.00 0.50 ]
+DEAL::  phi(3) = [ -0.50 0.00 ]
+DEAL::  grad phi(0) = [ 0.00 0.00 ]
+DEAL::                [ 1.00 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 -1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 0.00 ]
+DEAL::                [ -1.00 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::q_point(5)=-1.00 0.50
+DEAL::  phi(0) = [ 0.00 0.00 ]
+DEAL::  phi(1) = [ -0.50 0.00 ]
+DEAL::  phi(2) = [ 0.00 1.00 ]
+DEAL::  phi(3) = [ -0.50 0.00 ]
+DEAL::  grad phi(0) = [ 0.00 0.00 ]
+DEAL::                [ 1.00 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 -1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 0.00 ]
+DEAL::                [ -1.00 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::q_point(6)=0.00 1.00
+DEAL::  phi(0) = [ 0.00 1.00 ]
+DEAL::  phi(1) = [ -1.00 0.00 ]
+DEAL::  phi(2) = [ 0.00 0.00 ]
+DEAL::  phi(3) = [ 0.00 0.00 ]
+DEAL::  grad phi(0) = [ 0.00 0.00 ]
+DEAL::                [ 1.00 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 -1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 0.00 ]
+DEAL::                [ -1.00 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::q_point(7)=-0.50 1.00
+DEAL::  phi(0) = [ 0.00 0.50 ]
+DEAL::  phi(1) = [ -1.00 0.00 ]
+DEAL::  phi(2) = [ 0.00 0.50 ]
+DEAL::  phi(3) = [ 0.00 0.00 ]
+DEAL::  grad phi(0) = [ 0.00 0.00 ]
+DEAL::                [ 1.00 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 -1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 0.00 ]
+DEAL::                [ -1.00 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::q_point(8)=-1.00 1.00
+DEAL::  phi(0) = [ 0.00 0.00 ]
+DEAL::  phi(1) = [ -1.00 0.00 ]
+DEAL::  phi(2) = [ 0.00 1.00 ]
+DEAL::  phi(3) = [ 0.00 0.00 ]
+DEAL::  grad phi(0) = [ 0.00 0.00 ]
+DEAL::                [ 1.00 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 -1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 0.00 ]
+DEAL::                [ -1.00 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 1.00 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::
+DEAL::GRID TRANSFORMATION 2
 DEAL::q_point(0)=0.00 0.00
-DEAL::  shape_function(0)= 0.00 2.00
-DEAL::  shape_function(1)= 0.00 0.00
-DEAL::  shape_function(2)= 0.00 0.00
-DEAL::  shape_function(3)= -2.00 0.00
+DEAL::  phi(0) = [ 0.00 0.50 ]
+DEAL::  phi(1) = [ 0.00 0.00 ]
+DEAL::  phi(2) = [ 0.00 0.00 ]
+DEAL::  phi(3) = [ -0.50 0.00 ]
+DEAL::  grad phi(0) = [ 0.00 0.00 ]
+DEAL::                [ 0.25 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 -0.25 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 0.00 ]
+DEAL::                [ -0.25 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 0.25 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::q_point(1)=-1.00 0.00
+DEAL::  phi(0) = [ 0.00 0.25 ]
+DEAL::  phi(1) = [ 0.00 0.00 ]
+DEAL::  phi(2) = [ 0.00 0.25 ]
+DEAL::  phi(3) = [ -0.50 0.00 ]
+DEAL::  grad phi(0) = [ 0.00 0.00 ]
+DEAL::                [ 0.25 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 -0.25 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 0.00 ]
+DEAL::                [ -0.25 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 0.25 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::q_point(2)=-2.00 0.00
+DEAL::  phi(0) = [ 0.00 0.00 ]
+DEAL::  phi(1) = [ 0.00 0.00 ]
+DEAL::  phi(2) = [ 0.00 0.50 ]
+DEAL::  phi(3) = [ -0.50 0.00 ]
+DEAL::  grad phi(0) = [ 0.00 0.00 ]
+DEAL::                [ 0.25 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 -0.25 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 0.00 ]
+DEAL::                [ -0.25 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 0.25 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::q_point(3)=0.00 1.00
+DEAL::  phi(0) = [ 0.00 0.50 ]
+DEAL::  phi(1) = [ -0.25 0.00 ]
+DEAL::  phi(2) = [ 0.00 0.00 ]
+DEAL::  phi(3) = [ -0.25 0.00 ]
+DEAL::  grad phi(0) = [ 0.00 0.00 ]
+DEAL::                [ 0.25 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 -0.25 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 0.00 ]
+DEAL::                [ -0.25 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 0.25 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::q_point(4)=-1.00 1.00
+DEAL::  phi(0) = [ 0.00 0.25 ]
+DEAL::  phi(1) = [ -0.25 0.00 ]
+DEAL::  phi(2) = [ 0.00 0.25 ]
+DEAL::  phi(3) = [ -0.25 0.00 ]
+DEAL::  grad phi(0) = [ 0.00 0.00 ]
+DEAL::                [ 0.25 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 -0.25 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 0.00 ]
+DEAL::                [ -0.25 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 0.25 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::q_point(5)=-2.00 1.00
+DEAL::  phi(0) = [ 0.00 0.00 ]
+DEAL::  phi(1) = [ -0.25 0.00 ]
+DEAL::  phi(2) = [ 0.00 0.50 ]
+DEAL::  phi(3) = [ -0.25 0.00 ]
+DEAL::  grad phi(0) = [ 0.00 0.00 ]
+DEAL::                [ 0.25 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 -0.25 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 0.00 ]
+DEAL::                [ -0.25 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 0.25 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::q_point(6)=0.00 2.00
+DEAL::  phi(0) = [ 0.00 0.50 ]
+DEAL::  phi(1) = [ -0.50 0.00 ]
+DEAL::  phi(2) = [ 0.00 0.00 ]
+DEAL::  phi(3) = [ 0.00 0.00 ]
+DEAL::  grad phi(0) = [ 0.00 0.00 ]
+DEAL::                [ 0.25 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 -0.25 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 0.00 ]
+DEAL::                [ -0.25 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 0.25 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::q_point(7)=-1.00 2.00
+DEAL::  phi(0) = [ 0.00 0.25 ]
+DEAL::  phi(1) = [ -0.50 0.00 ]
+DEAL::  phi(2) = [ 0.00 0.25 ]
+DEAL::  phi(3) = [ 0.00 0.00 ]
+DEAL::  grad phi(0) = [ 0.00 0.00 ]
+DEAL::                [ 0.25 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 -0.25 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 0.00 ]
+DEAL::                [ -0.25 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 0.25 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::q_point(8)=-2.00 2.00
+DEAL::  phi(0) = [ 0.00 0.00 ]
+DEAL::  phi(1) = [ -0.50 0.00 ]
+DEAL::  phi(2) = [ 0.00 0.50 ]
+DEAL::  phi(3) = [ 0.00 0.00 ]
+DEAL::  grad phi(0) = [ 0.00 0.00 ]
+DEAL::                [ 0.25 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 -0.25 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 0.00 ]
+DEAL::                [ -0.25 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 0.25 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::
+DEAL::GRID TRANSFORMATION 3
+DEAL::q_point(0)=0.00 0.00
+DEAL::  phi(0) = [ 0.50 0.00 ]
+DEAL::  phi(1) = [ 0.00 0.00 ]
+DEAL::  phi(2) = [ 0.00 0.00 ]
+DEAL::  phi(3) = [ 0.00 1.00 ]
+DEAL::  grad phi(0) = [ 0.00 -0.50 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 0.00 ]
+DEAL::                [ 0.50 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 0.50 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 0.00 ]
+DEAL::                [ -0.50 0.00 ]
 DEAL::q_point(1)=0.00 0.50
-DEAL::  shape_function(0)= 0.00 1.00
-DEAL::  shape_function(1)= 0.00 0.00
-DEAL::  shape_function(2)= 0.00 1.00
-DEAL::  shape_function(3)= -2.00 0.00
+DEAL::  phi(0) = [ 0.25 0.00 ]
+DEAL::  phi(1) = [ 0.00 0.00 ]
+DEAL::  phi(2) = [ 0.25 0.00 ]
+DEAL::  phi(3) = [ 0.00 1.00 ]
+DEAL::  grad phi(0) = [ 0.00 -0.50 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 0.00 ]
+DEAL::                [ 0.50 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 0.50 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 0.00 ]
+DEAL::                [ -0.50 0.00 ]
 DEAL::q_point(2)=0.00 1.00
-DEAL::  shape_function(0)= 0.00 0.00
-DEAL::  shape_function(1)= 0.00 0.00
-DEAL::  shape_function(2)= 0.00 2.00
-DEAL::  shape_function(3)= -2.00 0.00
-DEAL::q_point(3)=0.50 0.00
-DEAL::  shape_function(0)= 0.00 2.00
-DEAL::  shape_function(1)= -1.00 0.00
-DEAL::  shape_function(2)= 0.00 0.00
-DEAL::  shape_function(3)= -1.00 0.00
-DEAL::q_point(4)=0.50 0.50
-DEAL::  shape_function(0)= 0.00 1.00
-DEAL::  shape_function(1)= -1.00 0.00
-DEAL::  shape_function(2)= 0.00 1.00
-DEAL::  shape_function(3)= -1.00 0.00
-DEAL::q_point(5)=0.50 1.00
-DEAL::  shape_function(0)= 0.00 0.00
-DEAL::  shape_function(1)= -1.00 0.00
-DEAL::  shape_function(2)= 0.00 2.00
-DEAL::  shape_function(3)= -1.00 0.00
-DEAL::q_point(6)=1.00 0.00
-DEAL::  shape_function(0)= 0.00 2.00
-DEAL::  shape_function(1)= -2.00 0.00
-DEAL::  shape_function(2)= 0.00 0.00
-DEAL::  shape_function(3)= 0.00 0.00
-DEAL::q_point(7)=1.00 0.50
-DEAL::  shape_function(0)= 0.00 1.00
-DEAL::  shape_function(1)= -2.00 0.00
-DEAL::  shape_function(2)= 0.00 1.00
-DEAL::  shape_function(3)= 0.00 0.00
-DEAL::q_point(8)=1.00 1.00
-DEAL::  shape_function(0)= 0.00 0.00
-DEAL::  shape_function(1)= -2.00 0.00
-DEAL::  shape_function(2)= 0.00 2.00
-DEAL::  shape_function(3)= 0.00 0.00
+DEAL::  phi(0) = [ 0.00 0.00 ]
+DEAL::  phi(1) = [ 0.00 0.00 ]
+DEAL::  phi(2) = [ 0.50 0.00 ]
+DEAL::  phi(3) = [ 0.00 1.00 ]
+DEAL::  grad phi(0) = [ 0.00 -0.50 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 0.00 ]
+DEAL::                [ 0.50 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 0.50 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 0.00 ]
+DEAL::                [ -0.50 0.00 ]
+DEAL::q_point(3)=1.00 0.00
+DEAL::  phi(0) = [ 0.50 0.00 ]
+DEAL::  phi(1) = [ 0.00 0.50 ]
+DEAL::  phi(2) = [ 0.00 0.00 ]
+DEAL::  phi(3) = [ 0.00 0.50 ]
+DEAL::  grad phi(0) = [ 0.00 -0.50 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 0.00 ]
+DEAL::                [ 0.50 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 0.50 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 0.00 ]
+DEAL::                [ -0.50 0.00 ]
+DEAL::q_point(4)=1.00 0.50
+DEAL::  phi(0) = [ 0.25 0.00 ]
+DEAL::  phi(1) = [ 0.00 0.50 ]
+DEAL::  phi(2) = [ 0.25 0.00 ]
+DEAL::  phi(3) = [ 0.00 0.50 ]
+DEAL::  grad phi(0) = [ 0.00 -0.50 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 0.00 ]
+DEAL::                [ 0.50 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 0.50 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 0.00 ]
+DEAL::                [ -0.50 0.00 ]
+DEAL::q_point(5)=1.00 1.00
+DEAL::  phi(0) = [ 0.00 0.00 ]
+DEAL::  phi(1) = [ 0.00 0.50 ]
+DEAL::  phi(2) = [ 0.50 0.00 ]
+DEAL::  phi(3) = [ 0.00 0.50 ]
+DEAL::  grad phi(0) = [ 0.00 -0.50 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 0.00 ]
+DEAL::                [ 0.50 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 0.50 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 0.00 ]
+DEAL::                [ -0.50 0.00 ]
+DEAL::q_point(6)=2.00 0.00
+DEAL::  phi(0) = [ 0.50 0.00 ]
+DEAL::  phi(1) = [ 0.00 1.00 ]
+DEAL::  phi(2) = [ 0.00 0.00 ]
+DEAL::  phi(3) = [ 0.00 0.00 ]
+DEAL::  grad phi(0) = [ 0.00 -0.50 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 0.00 ]
+DEAL::                [ 0.50 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 0.50 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 0.00 ]
+DEAL::                [ -0.50 0.00 ]
+DEAL::q_point(7)=2.00 0.50
+DEAL::  phi(0) = [ 0.25 0.00 ]
+DEAL::  phi(1) = [ 0.00 1.00 ]
+DEAL::  phi(2) = [ 0.25 0.00 ]
+DEAL::  phi(3) = [ 0.00 0.00 ]
+DEAL::  grad phi(0) = [ 0.00 -0.50 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 0.00 ]
+DEAL::                [ 0.50 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 0.50 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 0.00 ]
+DEAL::                [ -0.50 0.00 ]
+DEAL::q_point(8)=2.00 1.00
+DEAL::  phi(0) = [ 0.00 0.00 ]
+DEAL::  phi(1) = [ 0.00 1.00 ]
+DEAL::  phi(2) = [ 0.50 0.00 ]
+DEAL::  phi(3) = [ 0.00 0.00 ]
+DEAL::  grad phi(0) = [ 0.00 -0.50 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(1) = [ 0.00 0.00 ]
+DEAL::                [ 0.50 0.00 ]
+DEAL::  grad phi(2) = [ 0.00 0.50 ]
+DEAL::                [ 0.00 0.00 ]
+DEAL::  grad phi(3) = [ 0.00 0.00 ]
+DEAL::                [ -0.50 0.00 ]
+DEAL::

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.