]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Implement some things with gradients in 3d for the Nedelec element.
authorwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 16 Sep 2002 22:00:11 +0000 (22:00 +0000)
committerwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 16 Sep 2002 22:00:11 +0000 (22:00 +0000)
git-svn-id: https://svn.dealii.org/trunk@6436 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/source/fe/fe_nedelec.cc
tests/fe/non_primitive_2.cc
tests/fe/non_primitive_2.checked [new file with mode: 0644]

index 4ac1f20bee91761e550c4917ba907c1fdbac7145..2f283743cf17eee068cf1405091ce4ecc39984f0 100644 (file)
@@ -296,7 +296,7 @@ FE_Nedelec<dim>::shape_value_component (const unsigned int i,
                  case 2: return (component == 0 ? p(1) : 0);
                                                         // (0, 1-x)
                  case 3: return (component == 0 ? 0 : 1-p(0));
-
+                        
                                                         // there are
                                                         // only four
                                                         // shape
@@ -402,7 +402,7 @@ FE_Nedelec<dim>::shape_value_component (const unsigned int i,
 template <int dim>
 Tensor<1,dim>
 FE_Nedelec<dim>::shape_grad_component (const unsigned int i,
-                                      const Point<dim> &/*p*/,
+                                      const Point<dim> &p,
                                       const unsigned int component) const
 {
   Assert (i<this->dofs_per_cell, ExcIndexRange(i,0,this->dofs_per_cell));
@@ -448,6 +448,79 @@ FE_Nedelec<dim>::shape_grad_component (const unsigned int i,
          };
       };
 
+      case 3:  // 3d
+      {
+       switch (degree)
+         {
+                                            // first order Nedelec
+                                            // elements
+           case 1:
+           {
+                                              // on the unit cell,
+                                              // the gradients of
+                                              // these shape
+                                              // functions are
+                                              // linear. we pack them
+                                              // into an array,
+                                              // knowing that it may
+                                              // be expensive to
+                                              // recompute the whole
+                                              // array each
+                                              // time. maybe some
+                                              // clever compiler can
+                                              // optimize this out,
+                                              // seeing that except
+                                              // for one element all
+                                              // the other ones are
+                                              // dead stores...
+                                              //
+                                              // the format is: first
+                                              // index=shape function
+                                              // number; second
+                                              // index=vector
+                                              // component, thrid
+                                              // index=component
+                                              // within gradient
+             const double x = p(0),
+                          y = p(1),
+                          z = p(2);
+             static const double unit_gradients[12][3][3]
+               = { { {0,-(1-z), -(1-y)}, {0,0,0}, {     0,      0, 0} },
+                   { {0,     0,      0}, {0,0,0}, { (1-y),     -x, 0} },
+                    { {0,    -z,  (1-y)}, {0,0,0}, {     0,      0, 0} },
+                   { {0,     0,      0}, {0,0,0}, {-(1-y), -(1-x), 0} },
+                    
+                    { {0, (1-z),     -y}, {0,0,0}, {     0,      0, 0} },
+                   { {0,     0,      0}, {0,0,0}, {     y,      x, 0} },
+                    { {0,     z,      y}, {0,0,0}, {     0,      0, 0} },
+                   { {0,     0,      0}, {0,0,0}, {    -y,  (1-x), 0} },
+                    
+                    { {0, 0, 0}, {-(1-z), 0, -(1-x)}, {0, 0, 0} },
+                    { {0, 0, 0}, { (1-z), 0,     -x}, {0, 0, 0} },
+                    { {0, 0, 0}, {     z, 0,      x}, {0, 0, 0} },
+                    { {0, 0, 0}, {    -z, 0,  (1-x)}, {0, 0, 0} } };
+                                               // note: simple check
+                                               // whether this can at
+                                               // all be: build the
+                                               // sum over all these
+                                               // tensors. since the
+                                               // sum of the shape
+                                               // functions is a
+                                               // constant, the
+                                               // gradient must
+                                               // necessarily be
+                                               // zero. this is in
+                                               // fact the case here,
+                                               // so test successfull
+             return Tensor<1,dim>(unit_gradients[i][component]);
+           };
+
+                                            // no other degrees
+                                            // implemented
+           default:
+                 Assert (false, ExcNotImplemented());
+         };
+      };
                                       // presently no other space
                                       // dimension implemented
       default:
index f8f0b884235ffe97df47dcea40e0112d1a203f1f..52533c176e080f5b6aad3996f3be49850d902f19 100644 (file)
@@ -226,10 +226,10 @@ test ()
                                   // for the stokes equation the
                                   // pressure does not couple to
                                   // itself
-  std::vector<std::vector<bool> > mask (dim+1, std::vector<bool> (dim+1, true));
-  mask[dim][dim] = false;
+//   std::vector<std::vector<bool> > mask (dim+1, std::vector<bool> (dim+1, true));
+//   mask[dim][dim] = false;
   
-  DoFTools::make_sparsity_pattern (dof_handler, mask, sparsity);
+  DoFTools::make_sparsity_pattern (dof_handler, /*mask,*/ sparsity);
   sparsity.compress ();
   
   SparseMatrix<double> A2 (sparsity);
@@ -257,7 +257,7 @@ test ()
 
 int main ()
 {
-  std::ofstream logfile("non_primitive_1.output");
+  std::ofstream logfile("non_primitive_2.output");
   deallog.attach(logfile);
   deallog.depth_console(0);
 
diff --git a/tests/fe/non_primitive_2.checked b/tests/fe/non_primitive_2.checked
new file mode 100644 (file)
index 0000000..9f97140
--- /dev/null
@@ -0,0 +1,1105 @@
+
+DEAL::dim=2, n_cells=100
+DEAL::0 0.00000
+DEAL::8 0.00000
+DEAL::16 0.00000
+DEAL::24 0.00000
+DEAL::32 0.00000
+DEAL::40 0.00000
+DEAL::48 0.00000
+DEAL::56 0.00000
+DEAL::64 0.00000
+DEAL::72 0.00000
+DEAL::80 -2.01062
+DEAL::88 0.00000
+DEAL::96 0.00000
+DEAL::104 0.00000
+DEAL::112 0.00000
+DEAL::120 0.00000
+DEAL::128 0.00000
+DEAL::136 0.00000
+DEAL::144 -2.01062
+DEAL::152 0.00000
+DEAL::160 0.00000
+DEAL::168 0.00000
+DEAL::176 0.00000
+DEAL::184 0.00000
+DEAL::192 0.00000
+DEAL::200 2.01062
+DEAL::208 4.02124
+DEAL::216 0.00000
+DEAL::224 0.00000
+DEAL::232 0.00000
+DEAL::240 0.00000
+DEAL::248 0.00000
+DEAL::256 0.00000
+DEAL::264 0.00000
+DEAL::272 2.01062
+DEAL::280 0.00000
+DEAL::288 0.00000
+DEAL::296 0.00000
+DEAL::304 0.00000
+DEAL::312 0.00000
+DEAL::320 0.00000
+DEAL::328 0.00000
+DEAL::336 0.00000
+DEAL::344 0.00000
+DEAL::352 2.01062
+DEAL::360 0.00000
+DEAL::368 0.00000
+DEAL::376 0.00000
+DEAL::384 -2.01062
+DEAL::392 -2.01062
+DEAL::400 0.00000
+DEAL::408 0.00000
+DEAL::416 0.00000
+DEAL::424 0.00000
+DEAL::432 0.00000
+DEAL::440 0.00000
+DEAL::448 0.00000
+DEAL::456 0.00000
+DEAL::464 0.00000
+DEAL::472 0.00000
+DEAL::480 0.00000
+DEAL::488 0.00000
+DEAL::496 0.00000
+DEAL::504 0.00000
+DEAL::512 0.00000
+DEAL::520 0.00000
+DEAL::528 0.00000
+DEAL::536 0.00000
+DEAL::544 0.00000
+DEAL::552 0.00000
+DEAL::560 0.00000
+DEAL::568 0.00000
+DEAL::576 0.00000
+DEAL::584 0.00000
+DEAL::592 -2.01062
+DEAL::600 0.00
+DEAL::608 0.00000
+DEAL::616 0.00000
+DEAL::624 0.00000
+DEAL::632 0.00000
+DEAL::640 0.00000
+DEAL::648 0.00000
+DEAL::656 -2.01062
+DEAL::664 0.00
+DEAL::672 0.00000
+DEAL::680 0.00
+DEAL::688 0.00000
+DEAL::696 0.00000
+DEAL::704 0.00000
+DEAL::712 0.00000
+DEAL::720 0.00000
+DEAL::728 0.00
+DEAL::736 -2.01062
+DEAL::744 0.00
+DEAL::752 0.00000
+DEAL::760 0.00
+DEAL::768 0.00000
+DEAL::776 0.00000
+DEAL::784 0.00
+DEAL::792 0.00000
+DEAL::800 0.00000
+DEAL::808 0.00000
+DEAL::816 0.00000
+DEAL::824 0.00000
+DEAL::832 2.01062
+DEAL::840 4.02124
+DEAL::848 0.00000
+DEAL::856 0.00
+DEAL::864 0.00
+DEAL::872 0.00000
+DEAL::880 0.00000
+DEAL::888 0.00000
+DEAL::896 0.00000
+DEAL::904 0.00000
+DEAL::912 0.00000
+DEAL::920 0.00000
+DEAL::928 0.00000
+DEAL::936 0.00000
+DEAL::944 0.00000
+DEAL::952 0.00000
+DEAL::960 0.00000
+DEAL::968 0.00000
+DEAL::976 0.00
+DEAL::984 4.02124
+DEAL::992 0.00000
+DEAL::1000 0.00000
+DEAL::1008 0.00000
+DEAL::1016 0.00000
+DEAL::1024 0.00000
+DEAL::1032 0.00000
+DEAL::1040 0.00000
+DEAL::1048 -2.01062
+DEAL::1056 0.00000
+DEAL::1064 0.00
+DEAL::1072 0.00000
+DEAL::1080 0.00000
+DEAL::1088 0.00000
+DEAL::1096 -2.01062
+DEAL::1104 0.00
+DEAL::1112 0.00000
+DEAL::1120 0.00000
+DEAL::1128 0.00000
+DEAL::1136 0.00000
+DEAL::1144 0.00000
+DEAL::1152 0.00
+DEAL::1160 0.00
+DEAL::1168 0.00000
+DEAL::1176 0.00
+DEAL::1184 0.00000
+DEAL::1192 0.00000
+DEAL::1200 0.00000
+DEAL::1208 0.00
+DEAL::1216 0.00000
+DEAL::1224 0.00000
+DEAL::1232 0.00000
+DEAL::1240 0.00000
+DEAL::1248 0.00000
+DEAL::1256 0.00000
+DEAL::1264 0.00
+DEAL::1272 -2.01062
+DEAL::1280 0.00000
+DEAL::1288 0.00000
+DEAL::1296 0.00000
+DEAL::1304 0.00
+DEAL::1312 0.00000
+DEAL::1320 0.00000
+DEAL::1328 0.00000
+DEAL::1336 0.00
+DEAL::1344 0.00000
+DEAL::1352 0.00
+DEAL::1360 0.00000
+DEAL::1368 2.01062
+DEAL::1376 0.00000
+DEAL::1384 0.00000
+DEAL::1392 0.00000
+DEAL::1400 0.00000
+DEAL::1408 0.00000
+DEAL::1416 0.00000
+DEAL::1424 0.00000
+DEAL::1432 0.00000
+DEAL::1440 0.00000
+DEAL::1448 0.00000
+DEAL::1456 -2.01062
+DEAL::1464 0.00
+DEAL::1472 0.00000
+DEAL::1480 0.00
+DEAL::1488 0.00000
+DEAL::1496 0.00000
+DEAL::1504 0.00000
+DEAL::1512 0.00000
+DEAL::1520 -2.01062
+DEAL::1528 0.00000
+DEAL::1536 0.00000
+DEAL::1544 0.00
+DEAL::1552 0.00000
+DEAL::1560 0.00000
+DEAL::1568 0.00000
+DEAL::1576 0.00000
+DEAL::1584 0.00000
+DEAL::1592 0.00
+DEAL::1600 -2.01062
+DEAL::1608 0.00
+DEAL::1616 0.00000
+DEAL::1624 0.00
+DEAL::1632 0.00000
+DEAL::1640 0.00000
+DEAL::1648 0.00
+DEAL::1656 0.00000
+DEAL::1664 0.00000
+DEAL::1672 0.00000
+DEAL::1680 0.00000
+DEAL::1688 0.00
+DEAL::1696 4.02124
+DEAL::1704 0.00000
+DEAL::1712 0.00000
+DEAL::1720 0.00000
+DEAL::1728 0.00000
+DEAL::1736 0.00000
+DEAL::1744 0.00000
+DEAL::1752 0.00000
+DEAL::1760 0.00000
+DEAL::1768 0.00000
+DEAL::1776 0.00000
+DEAL::1784 0.00000
+DEAL::1792 0.00000
+DEAL::1800 0.00000
+DEAL::1808 0.00000
+DEAL::1816 0.00000
+DEAL::1824 4.02124
+DEAL::1832 0.00000
+DEAL::1840 0.00
+DEAL::1848 -2.01062
+DEAL::1856 0.00000
+DEAL::1864 0.00000
+DEAL::1872 0.00000
+DEAL::1880 0.00000
+DEAL::1888 0.00000
+DEAL::1896 0.00
+DEAL::1904 0.00
+DEAL::1912 0.00
+DEAL::1920 0.00000
+DEAL::1928 0.00000
+DEAL::1936 0.00000
+DEAL::1944 0.00000
+DEAL::1952 4.02124
+DEAL::1960 0.00000
+DEAL::1968 -2.01062
+DEAL::1976 -2.01062
+DEAL::1984 0.00000
+DEAL::1992 0.00000
+DEAL::2000 0.00000
+DEAL::2008 0.00000
+DEAL::2016 0.00
+DEAL::2024 0.00
+DEAL::2032 0.00000
+DEAL::2040 0.00000
+DEAL::2048 0.00000
+DEAL::2056 0.00
+DEAL::2064 -2.01062
+DEAL::2072 0.00000
+DEAL::2080 0.00000
+DEAL::2088 0.00000
+DEAL::2096 0.00000
+DEAL::2104 0.00000
+DEAL::2112 0.00000
+DEAL::2120 0.00000
+DEAL::2128 0.00
+DEAL::2136 0.00000
+DEAL::2144 0.00000
+DEAL::2152 0.00
+DEAL::2160 0.00000
+DEAL::2168 0.00000
+DEAL::2176 0.00000
+DEAL::2184 0.00000
+DEAL::2192 0.00000
+DEAL::2200 0.00
+DEAL::2208 -2.01062
+DEAL::2216 0.00000
+DEAL::2224 0.00000
+DEAL::2232 4.02124
+DEAL::2240 0.00000
+DEAL::2248 -2.01062
+DEAL::2256 0.00
+DEAL::2264 0.00000
+DEAL::2272 0.00000
+DEAL::2280 0.00000
+DEAL::2288 0.00
+DEAL::2296 0.00000
+DEAL::2304 0.00
+DEAL::2312 0.00000
+DEAL::2320 0.00000
+DEAL::2328 0.00000
+DEAL::2336 0.00000
+DEAL::2344 0.00000
+DEAL::2352 0.00000
+DEAL::2360 4.02124
+DEAL::2368 0.00000
+DEAL::2376 -2.01062
+DEAL::2384 0.00
+DEAL::2392 0.00000
+DEAL::2400 0.00000
+DEAL::2408 0.00000
+DEAL::2416 0.00000
+DEAL::2424 0.00000
+DEAL::2432 0.00
+DEAL::2440 0.00
+DEAL::2448 0.00000
+DEAL::2456 0.00000
+DEAL::2464 -2.01062
+DEAL::2472 0.00
+DEAL::2480 0.00000
+DEAL::2488 0.00000
+DEAL::2496 0.00000
+DEAL::2504 0.00000
+DEAL::2512 0.00000
+DEAL::2520 0.00
+DEAL::2528 0.00
+DEAL::2536 0.00000
+DEAL::2544 0.00000
+DEAL::2552 0.00000
+DEAL::2560 0.00
+DEAL::2568 0.00000
+DEAL::2576 0.00000
+DEAL::2584 0.00000
+DEAL::2592 0.00000
+DEAL::2600 0.00000
+DEAL::2608 0.00
+DEAL::2616 4.02124
+DEAL::2624 0.00000
+DEAL::2632 0.00
+DEAL::2640 -2.01062
+DEAL::2648 0.00000
+DEAL::2656 0.00000
+DEAL::2664 0.00000
+DEAL::2672 0.00000
+DEAL::2680 0.00000
+DEAL::2688 0.00000
+DEAL::2696 0.00000
+DEAL::2704 -2.01062
+DEAL::2712 0.00
+DEAL::2720 0.00000
+DEAL::2728 0.00000
+DEAL::2736 0.00000
+DEAL::2744 0.00
+DEAL::2752 -2.01062
+DEAL::2760 0.00000
+DEAL::2768 0.00000
+DEAL::2776 0.00000
+DEAL::2784 -2.01062
+DEAL::2792 0.00
+DEAL::2800 0.00000
+DEAL::2808 0.00000
+DEAL::2816 0.00
+DEAL::2824 0.00000
+DEAL::2832 0.00000
+DEAL::2840 0.00
+DEAL::2848 0.00000
+DEAL::2856 0.00000
+DEAL::2864 0.00000
+DEAL::2872 0.00000
+DEAL::2880 0.00000
+DEAL::2888 0.00000
+DEAL::2896 0.00000
+DEAL::2904 0.00000
+DEAL::2912 0.00000
+DEAL::2920 0.00000
+DEAL::2928 0.00000
+DEAL::2936 -8.04248
+DEAL::2944 0.00000
+DEAL::2952 0.00000
+DEAL::2960 0.00000
+DEAL::2968 0.00000
+DEAL::2976 0.00000
+DEAL::2984 0.00000
+DEAL::2992 0.00000
+DEAL::3000 -8.04248
+DEAL::3008 0.00000
+DEAL::3016 0.00000
+DEAL::3024 0.00000
+DEAL::3032 0.00000
+DEAL::3040 0.00000
+DEAL::3048 0.00000
+DEAL::3056 0.00000
+DEAL::3064 0.00000
+DEAL::3072 0.00000
+DEAL::3080 0.00000
+DEAL::3088 -8.04248
+DEAL::3096 0.00000
+DEAL::3104 0.00000
+DEAL::3112 0.00000
+DEAL::3120 0.00000
+DEAL::3128 0.00000
+DEAL::3136 0.00000
+DEAL::3144 0.00000
+DEAL::3152 0.00000
+DEAL::3160 0.00000
+DEAL::3168 0.00000
+DEAL::3176 8.04248
+DEAL::3184 16.0850
+DEAL::3192 0.00000
+DEAL::3200 0.00000
+DEAL::3208 0.00000
+DEAL::3216 0.00000
+DEAL::3224 0.00000
+DEAL::3232 0.00000
+DEAL::3240 0.00000
+DEAL::3248 0.00000
+DEAL::3256 0.00000
+DEAL::3264 0.00000
+DEAL::3272 0.00000
+DEAL::3280 -8.04248
+DEAL::3288 0.00000
+DEAL::3296 0.00000
+DEAL::3304 -8.04248
+DEAL::3312 -8.04248
+DEAL::3320 0.00000
+DEAL::3328 0.00000
+DEAL::3336 0.00000
+DEAL::3344 0.00000
+DEAL::3352 0.00000
+DEAL::3360 0.00000
+DEAL::3368 0.00000
+DEAL::3376 0.00000
+DEAL::3384 0.00000
+DEAL::3392 8.04248
+DEAL::3400 0.00000
+DEAL::3408 0.00000
+DEAL::3416 0.00000
+DEAL::3424 0.00000
+DEAL::3432 -8.04248
+DEAL::3440 0.00000
+DEAL::3448 0.00000
+DEAL::3456 0.00000
+DEAL::3464 0.00000
+DEAL::3472 0.00
+DEAL::3480 0.00000
+DEAL::3488 -8.04248
+DEAL::3496 0.00000
+DEAL::3504 0.00000
+DEAL::3512 0.00000
+DEAL::3520 0.00000
+DEAL::3528 0.00000
+DEAL::3536 0.00000
+DEAL::3544 0.00000
+DEAL::3552 0.00000
+DEAL::3560 0.00000
+DEAL::3568 0.00000
+DEAL::3576 0.00000
+DEAL::3584 0.00000
+DEAL::3592 0.00000
+DEAL::3600 0.00000
+DEAL::3608 0.00000
+DEAL::3616 0.00
+DEAL::3624 0.00
+DEAL::3632 0.00000
+DEAL::3640 0.00000
+DEAL::3648 0.00
+DEAL::3656 0.00000
+DEAL::3664 0.00000
+DEAL::3672 8.04248
+DEAL::3680 16.0850
+DEAL::3688 0.00
+DEAL::3696 0.00000
+DEAL::3704 -8.04248
+DEAL::3712 0.00000
+DEAL::3720 0.00000
+DEAL::3728 0.00000
+DEAL::3736 0.00000
+DEAL::3744 0.00000
+DEAL::3752 0.00000
+DEAL::3760 0.00
+DEAL::3768 0.00
+DEAL::3776 0.00000
+DEAL::3784 0.00000
+DEAL::3792 0.00000
+DEAL::3800 16.0850
+DEAL::3808 0.00000
+DEAL::3816 0.00
+DEAL::3824 0.00000
+DEAL::3832 0.00000
+DEAL::3840 0.00000
+DEAL::3848 0.00000
+DEAL::3856 0.00000
+DEAL::3864 0.00
+DEAL::3872 0.00000
+DEAL::3880 0.00000
+DEAL::3888 0.00
+DEAL::3896 0.00000
+DEAL::3904 0.00000
+DEAL::3912 0.00000
+DEAL::3920 0.00000
+DEAL::3928 0.00000
+DEAL::3936 0.00000
+DEAL::3944 16.0850
+DEAL::3952 0.00
+DEAL::3960 0.00
+DEAL::3968 -8.04248
+DEAL::3976 0.00000
+DEAL::3984 0.00000
+DEAL::3992 0.00000
+DEAL::4000 0.00000
+DEAL::4008 0.00000
+DEAL::4016 0.00
+DEAL::4024 0.00
+DEAL::4032 0.00
+DEAL::4040 0.00000
+DEAL::4048 0.00000
+DEAL::4056 0.00000
+DEAL::4064 8.04248
+DEAL::4072 16.0850
+DEAL::4080 0.00
+DEAL::4088 0.00000
+DEAL::4096 0.00000
+DEAL::4104 0.00
+DEAL::4112 0.00
+DEAL::4120 0.00000
+DEAL::4128 0.00000
+DEAL::4136 0.00000
+DEAL::4144 0.00
+DEAL::4152 0.00000
+DEAL::4160 0.00000
+DEAL::4168 0.00000
+DEAL::4176 0.00000
+DEAL::4184 0.00000
+DEAL::4192 0.00
+DEAL::4200 0.00000
+DEAL::4208 8.04248
+DEAL::4216 16.0850
+DEAL::4224 0.00000
+DEAL::4232 0.00000
+DEAL::4240 0.00000
+DEAL::4248 0.00000
+DEAL::4256 0.00000
+DEAL::4264 0.00000
+DEAL::4272 0.00000
+DEAL::4280 0.00000
+DEAL::4288 0.00000
+DEAL::4296 0.00000
+DEAL::4304 0.00000
+DEAL::4312 0.00000
+DEAL::4320 0.00000
+DEAL::4328 0.00000
+DEAL::4336 0.00000
+DEAL::4344 0.00
+DEAL::4352 -32.1699
+DEAL::4360 32.1699
+DEAL::4368 0.00000
+DEAL::4376 0.00000
+DEAL::4384 0.00000
+DEAL::4392 0.00000
+DEAL::4400 0.00000
+DEAL::4408 -32.1699
+DEAL::4416 -32.1699
+DEAL::4424 0.00000
+DEAL::4432 0.00000
+DEAL::4440 0.00000
+DEAL::4448 0.00000
+DEAL::4456 0.00000
+DEAL::4464 0.00000
+DEAL::4472 -32.1699
+DEAL::4480 0.00000
+DEAL::4488 0.00000
+DEAL::4496 0.00
+DEAL::4504 0.00000
+DEAL::4512 0.00000
+DEAL::4520 -32.1699
+DEAL::4528 -32.1699
+DEAL::4536 0.00000
+DEAL::4544 0.00000
+DEAL::4552 0.00000
+DEAL::4560 0.00000
+DEAL::4568 0.00000
+DEAL::4576 0.00000
+DEAL::4584 -32.1699
+DEAL::4592 0.00
+DEAL::4600 0.00000
+DEAL::4608 0.00000
+DEAL::4616 0.00000
+DEAL::4624 0.00000
+DEAL::4632 0.00000
+DEAL::4640 0.00
+DEAL::4648 0.00000
+DEAL::4656 0.00000
+DEAL::4664 0.00000
+DEAL::4672 0.00000
+DEAL::4680 0.00000
+DEAL::4688 0.00
+DEAL::4696 0.00
+DEAL::4704 0.00000
+DEAL::4712 0.00000
+DEAL::4720 0.00000
+DEAL::4728 0.00000
+DEAL::4736 0.00000
+DEAL::4744 0.00000
+DEAL::4752 0.00000
+DEAL::4760 0.00000
+DEAL::4768 0.00000
+DEAL::4776 0.00000
+DEAL::4784 64.3398
+DEAL::4792 -32.1699
+DEAL::4800 0.00000
+DEAL::4808 0.00000
+DEAL::4816 0.00000
+DEAL::4824 0.00
+DEAL::4832 0.00000
+DEAL::4840 0.00000
+DEAL::4848 0.00
+DEAL::4856 -32.1699
+DEAL::4864 0.00000
+DEAL::4872 -32.1699
+DEAL::4880 0.00
+DEAL::4888 0.00
+DEAL::4896 0.00000
+DEAL::4904 0.00000
+DEAL::4912 0.00000
+DEAL::4920 0.00000
+DEAL::4928 64.3398
+DEAL::4936 0.00000
+DEAL::4944 0.00000
+DEAL::4952 -32.1699
+DEAL::4960 -32.1699
+DEAL::4968 0.00
+DEAL::4976 0.00000
+DEAL::4984 0.00000
+DEAL::4992 0.00000
+DEAL::5000 -32.1699
+DEAL::5008 -32.1699
+DEAL::5016 0.00
+DEAL::5024 0.00
+DEAL::dim=3, n_cells=43
+DEAL::0 0.00000
+DEAL::27 0.00000
+DEAL::54 0.00000
+DEAL::81 0.00000
+DEAL::108 0.00000
+DEAL::135 0.00000
+DEAL::162 0.00000
+DEAL::189 0.00000
+DEAL::216 0.00000
+DEAL::243 0.00000
+DEAL::270 0.00000
+DEAL::297 0.00000
+DEAL::324 0.00000
+DEAL::351 0.00000
+DEAL::378 0.00000
+DEAL::405 0.00000
+DEAL::432 -0.0431842
+DEAL::459 0.00
+DEAL::486 -0.0125664
+DEAL::513 0.00000
+DEAL::540 -0.0376991
+DEAL::567 0.00
+DEAL::594 0.00
+DEAL::621 0.00
+DEAL::648 0.00000
+DEAL::675 0.00
+DEAL::702 0.00000
+DEAL::729 -0.0431842
+DEAL::756 0.00000
+DEAL::783 0.00000
+DEAL::810 0.00000
+DEAL::837 0.00000
+DEAL::864 0.00000
+DEAL::891 0.00000
+DEAL::918 0.0502655
+DEAL::945 0.00
+DEAL::972 0.00
+DEAL::999 0.00000
+DEAL::1026 0.00000
+DEAL::1053 0.00
+DEAL::1080 0.00
+DEAL::1107 0.00000
+DEAL::1134 0.00000
+DEAL::1161 0.00000
+DEAL::1188 0.00000
+DEAL::1215 0.00000
+DEAL::1242 0.00000
+DEAL::1269 0.00000
+DEAL::1296 0.00
+DEAL::1323 0.00000
+DEAL::1350 0.00000
+DEAL::1377 -0.0431842
+DEAL::1404 -0.0376991
+DEAL::1431 0.00
+DEAL::1458 0.00000
+DEAL::1485 0.00000
+DEAL::1512 0.00000
+DEAL::1539 0.00
+DEAL::1566 0.00
+DEAL::1593 0.00000
+DEAL::1620 0.00000
+DEAL::1647 0.00000
+DEAL::1674 0.00000
+DEAL::1701 0.00000
+DEAL::1728 0.00000
+DEAL::1755 0.00000
+DEAL::1782 0.00000
+DEAL::1809 0.00000
+DEAL::1836 0.00000
+DEAL::1863 0.00000
+DEAL::1890 0.00000
+DEAL::1917 0.00000
+DEAL::1944 0.00000
+DEAL::1971 0.00000
+DEAL::1998 0.00
+DEAL::2025 -0.0863685
+DEAL::2052 0.00000
+DEAL::2079 -0.0863685
+DEAL::2106 0.00000
+DEAL::2133 -0.0753982
+DEAL::2160 0.00
+DEAL::2187 0.00000
+DEAL::2214 0.00000
+DEAL::2241 0.00
+DEAL::2268 0.00000
+DEAL::2295 0.402124
+DEAL::2322 0.00
+DEAL::2349 0.00000
+DEAL::2376 0.00
+DEAL::2403 0.00000
+DEAL::2430 0.00000
+DEAL::2457 0.00000
+DEAL::2484 0.00000
+DEAL::2511 0.00
+DEAL::2538 0.00000
+DEAL::2565 0.298401
+DEAL::2592 -0.0863685
+DEAL::2619 0.00000
+DEAL::2646 0.00
+DEAL::2673 0.00
+DEAL::2700 0.103723
+DEAL::2727 0.00
+DEAL::2754 0.00000
+DEAL::2781 0.00000
+DEAL::2808 0.00000
+DEAL::2835 0.00000
+DEAL::2862 0.00000
+DEAL::2889 0.00000
+DEAL::2916 0.00000
+DEAL::2943 0.00
+DEAL::2970 0.00000
+DEAL::2997 0.00000
+DEAL::3024 0.00
+DEAL::3051 0.00
+DEAL::3078 0.00000
+DEAL::3105 0.00
+DEAL::3132 0.00000
+DEAL::3159 0.00
+DEAL::3186 0.00
+DEAL::3213 0.00
+DEAL::3240 0.00000
+DEAL::3267 0.00000
+DEAL::3294 0.00000
+DEAL::3321 0.00000
+DEAL::3348 0.00
+DEAL::3375 0.00000
+DEAL::3402 0.00
+DEAL::3429 0.00000
+DEAL::3456 0.00000
+DEAL::3483 0.00000
+DEAL::3510 0.00000
+DEAL::3537 0.00000
+DEAL::3564 0.00000
+DEAL::3591 0.00000
+DEAL::3618 0.00000
+DEAL::3645 0.00000
+DEAL::3672 0.00000
+DEAL::3699 0.00000
+DEAL::3726 0.00
+DEAL::3753 0.00
+DEAL::3780 -0.0251327
+DEAL::3807 0.00
+DEAL::3834 0.00
+DEAL::3861 -0.0251327
+DEAL::3888 0.298401
+DEAL::3915 0.00
+DEAL::3942 0.00000
+DEAL::3969 0.00
+DEAL::3996 0.00000
+DEAL::4023 0.00000
+DEAL::4050 0.00000
+DEAL::4077 0.00000
+DEAL::4104 0.00
+DEAL::4131 0.00000
+DEAL::4158 0.00
+DEAL::4185 0.00
+DEAL::4212 0.00000
+DEAL::4239 0.00000
+DEAL::4266 0.00000
+DEAL::4293 0.00000
+DEAL::4320 0.00000
+DEAL::4347 0.00
+DEAL::4374 0.00000
+DEAL::4401 0.00
+DEAL::4428 0.00000
+DEAL::4455 0.00
+DEAL::4482 0.00
+DEAL::4509 0.00000
+DEAL::4536 0.00000
+DEAL::4563 0.00000
+DEAL::4590 0.00000
+DEAL::4617 0.00000
+DEAL::4644 0.00000
+DEAL::4671 0.00000
+DEAL::4698 0.00000
+DEAL::4725 0.00000
+DEAL::4752 0.00000
+DEAL::4779 0.00
+DEAL::4806 0.00
+DEAL::4833 0.00
+DEAL::4860 0.00000
+DEAL::4887 0.00000
+DEAL::4914 0.00000
+DEAL::4941 0.00
+DEAL::4968 0.00000
+DEAL::4995 0.00000
+DEAL::5022 0.00000
+DEAL::5049 0.00000
+DEAL::5076 0.00000
+DEAL::5103 0.00000
+DEAL::5130 0.00000
+DEAL::5157 0.00000
+DEAL::5184 0.402124
+DEAL::5211 0.00000
+DEAL::5238 0.00000
+DEAL::5265 0.0109702
+DEAL::5292 0.00000
+DEAL::5319 0.00000
+DEAL::5346 0.00000
+DEAL::5373 0.00000
+DEAL::5400 0.00
+DEAL::5427 0.00000
+DEAL::5454 0.103723
+DEAL::5481 0.0109702
+DEAL::5508 0.00
+DEAL::5535 0.00000
+DEAL::5562 -0.0753982
+DEAL::5589 0.103723
+DEAL::5616 0.00
+DEAL::5643 0.00000
+DEAL::5670 0.00000
+DEAL::5697 0.00000
+DEAL::5724 0.00
+DEAL::5751 0.00
+DEAL::5778 0.00000
+DEAL::5805 0.00000
+DEAL::5832 0.00000
+DEAL::5859 0.00000
+DEAL::5886 0.00000
+DEAL::5913 0.00000
+DEAL::5940 0.00000
+DEAL::5967 0.00000
+DEAL::5994 0.00000
+DEAL::6021 0.00000
+DEAL::6048 0.00000
+DEAL::6075 0.00
+DEAL::6102 0.00
+DEAL::6129 0.00
+DEAL::6156 0.00
+DEAL::6183 0.00000
+DEAL::6210 0.00000
+DEAL::6237 0.00000
+DEAL::6264 -0.0863685
+DEAL::6291 0.00000
+DEAL::6318 0.00000
+DEAL::6345 0.00000
+DEAL::6372 0.00000
+DEAL::6399 0.00
+DEAL::6426 0.00
+DEAL::6453 0.00000
+DEAL::6480 0.00000
+DEAL::6507 -0.0753982
+DEAL::6534 0.00000
+DEAL::6561 0.00000
+DEAL::6588 0.00
+DEAL::6615 0.00000
+DEAL::6642 0.00
+DEAL::6669 0.00000
+DEAL::6696 0.0109702
+DEAL::6723 0.00000
+DEAL::6750 0.00
+DEAL::6777 0.00000
+DEAL::6804 0.00000
+DEAL::6831 0.00000
+DEAL::6858 0.00000
+DEAL::6885 0.00000
+DEAL::6912 0.00000
+DEAL::6939 0.00000
+DEAL::6966 0.00000
+DEAL::6993 -0.0863685
+DEAL::7020 0.402124
+DEAL::7047 0.00000
+DEAL::7074 0.00000
+DEAL::7101 0.00
+DEAL::7128 -0.0251327
+DEAL::7155 0.00000
+DEAL::7182 0.00
+DEAL::7209 0.00000
+DEAL::7236 -0.0753982
+DEAL::7263 0.00
+DEAL::7290 0.00000
+DEAL::7317 0.00000
+DEAL::7344 0.00000
+DEAL::7371 0.00000
+DEAL::7398 0.00000
+DEAL::7425 0.00000
+DEAL::7452 0.00000
+DEAL::7479 0.00
+DEAL::7506 0.0109702
+DEAL::7533 0.00000
+DEAL::7560 0.00000
+DEAL::7587 0.00000
+DEAL::7614 0.00000
+DEAL::7641 0.00000
+DEAL::7668 0.00000
+DEAL::7695 0.00
+DEAL::7722 0.00
+DEAL::7749 0.00
+DEAL::7776 0.00000
+DEAL::7803 -0.0251327
+DEAL::7830 0.00000
+DEAL::7857 0.00
+DEAL::7884 0.00000
+DEAL::7911 0.00000
+DEAL::7938 0.00000
+DEAL::7965 0.00000
+DEAL::7992 0.00000
+DEAL::8019 0.00000
+DEAL::8046 0.00000
+DEAL::8073 -0.0863685
+DEAL::8100 0.00000
+DEAL::8127 0.00000
+DEAL::8154 0.00000
+DEAL::8181 0.00000
+DEAL::8208 0.00000
+DEAL::8235 0.00
+DEAL::8262 0.00000
+DEAL::8289 0.00000
+DEAL::8316 0.00000
+DEAL::8343 0.00000
+DEAL::8370 0.00000
+DEAL::8397 -0.0251327
+DEAL::8424 0.00
+DEAL::8451 0.00000
+DEAL::8478 0.00000
+DEAL::8505 0.00000
+DEAL::8532 0.00000
+DEAL::8559 0.00
+DEAL::8586 0.00
+DEAL::8613 0.00000
+DEAL::8640 0.00000
+DEAL::8667 0.00000
+DEAL::8694 0.00
+DEAL::8721 -0.0251327
+DEAL::8748 0.00
+DEAL::8775 0.00
+DEAL::8802 -0.0753982
+DEAL::8829 0.00000
+DEAL::8856 0.00000
+DEAL::8883 0.00
+DEAL::8910 0.00000
+DEAL::8937 -0.0863685
+DEAL::8964 0.00000
+DEAL::8991 0.00000
+DEAL::9018 0.00000
+DEAL::9045 0.00000
+DEAL::9072 0.00000
+DEAL::9099 0.00
+DEAL::9126 -0.0251327
+DEAL::9153 0.00
+DEAL::9180 -0.0863685
+DEAL::9207 0.00000
+DEAL::9234 0.00000
+DEAL::9261 0.00
+DEAL::9288 0.00000
+DEAL::9315 0.00000
+DEAL::9342 0.00
+DEAL::9369 0.00
+DEAL::9396 0.00
+DEAL::9423 0.00000
+DEAL::9450 0.00
+DEAL::9477 0.00000
+DEAL::9504 0.00000
+DEAL::9531 0.00000
+DEAL::9558 0.00
+DEAL::9585 0.00
+DEAL::9612 0.00000
+DEAL::9639 0.00
+DEAL::9666 0.00000
+DEAL::9693 0.00000
+DEAL::9720 0.00000
+DEAL::9747 0.00
+DEAL::9774 0.00
+DEAL::9801 0.00000
+DEAL::9828 0.00000
+DEAL::9855 0.00000
+DEAL::9882 0.00000
+DEAL::9909 0.00000
+DEAL::9936 0.00000
+DEAL::9963 0.0109702
+DEAL::9990 0.00000
+DEAL::10017 0.00000
+DEAL::10044 0.00
+DEAL::10071 0.00000
+DEAL::10098 0.00000
+DEAL::10125 0.00000
+DEAL::10152 0.00
+DEAL::10179 0.00000
+DEAL::10206 0.00000
+DEAL::10233 0.00
+DEAL::10260 0.00000
+DEAL::10287 0.0109702
+DEAL::10314 0.00000
+DEAL::10341 0.00000
+DEAL::10368 0.0109702
+DEAL::10395 0.00000
+DEAL::10422 0.00000
+DEAL::10449 0.00000
+DEAL::10476 0.00000
+DEAL::10503 0.00000
+DEAL::10530 0.00000
+DEAL::10557 0.00000
+DEAL::10584 0.00000
+DEAL::10611 0.00000
+DEAL::10638 0.00000
+DEAL::10665 0.00000
+DEAL::10692 0.00000
+DEAL::10719 0.00000
+DEAL::10746 0.00000
+DEAL::10773 0.00000
+DEAL::10800 0.00
+DEAL::10827 0.00
+DEAL::10854 0.00
+DEAL::10881 -0.0502655
+DEAL::10908 0.00
+DEAL::10935 0.00
+DEAL::10962 0.00000
+DEAL::10989 0.596801
+DEAL::11016 -0.172737
+DEAL::11043 0.00000
+DEAL::11070 0.00
+DEAL::11097 0.00000
+DEAL::11124 0.596801
+DEAL::11151 0.00
+DEAL::11178 0.00000
+DEAL::11205 0.00000
+DEAL::11232 0.00000
+DEAL::11259 0.00000
+DEAL::11286 -0.172737
+DEAL::11313 0.00000
+DEAL::11340 0.00
+DEAL::11367 0.00
+DEAL::11394 0.00
+DEAL::11421 0.00000
+DEAL::11448 0.00
+DEAL::11475 0.00
+DEAL::11502 0.00
+DEAL::11529 0.00000
+DEAL::11556 0.00000
+DEAL::11583 0.00000
+DEAL::11610 0.00000
+DEAL::11637 0.00000
+DEAL::11664 0.00
+DEAL::11691 0.00000
+DEAL::11718 0.00000
+DEAL::11745 0.00
+DEAL::11772 0.00
+DEAL::11799 0.00
+DEAL::11826 0.00
+DEAL::11853 0.00
+DEAL::11880 0.00000
+DEAL::11907 0.00000
+DEAL::11934 -0.172737
+DEAL::11961 0.00
+DEAL::11988 0.00
+DEAL::12015 0.00000
+DEAL::12042 0.00000
+DEAL::12069 0.00000
+DEAL::12096 0.00000
+DEAL::12123 0.00000
+DEAL::12150 0.00
+DEAL::12177 0.00
+DEAL::12204 0.00000
+DEAL::12231 0.00000
+DEAL::12258 -0.172737
+DEAL::12285 0.00000
+DEAL::12312 0.00
+DEAL::12339 0.00
+DEAL::12366 -0.0502655
+DEAL::12393 0.00000
+DEAL::12420 0.00
+DEAL::12447 0.00000
+DEAL::12474 0.0219405
+DEAL::12501 0.00000
+DEAL::12528 0.0219405
+DEAL::12555 0.00000
+DEAL::12582 0.00638455
+DEAL::12609 0.00
+DEAL::12636 0.00000
+DEAL::12663 0.00000
+DEAL::12690 -0.0502655
+DEAL::12717 0.00000
+DEAL::12744 0.00

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.