]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Eliminate stores of zeroes.
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Wed, 24 May 2000 12:47:41 +0000 (12:47 +0000)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Wed, 24 May 2000 12:47:41 +0000 (12:47 +0000)
git-svn-id: https://svn.dealii.org/trunk@2937 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/source/fe/fe_lib.dg.cc

index 8849e667487380f3d075aa691630e39428105754..0a3386f340475d8dde26f5e27b474b1c1a322bdc 100644 (file)
@@ -34,8 +34,6 @@ FEDG_Q1<1>::FEDG_Q1():
   restriction[0](0,0) = 3.0/4.0;
   restriction[0](0,1) = 1.0/2.0;
   restriction[0](1,0) = -1.0/4.0;
-  restriction[0](1,1) = 0.0;
-  restriction[1](0,0) = 0.0;
   restriction[1](0,1) = -1.0/4.0;
   restriction[1](1,0) = 1.0/2.0;
   restriction[1](1,1) = 3.0/4.0;
@@ -218,60 +216,32 @@ FEDG_Q1<2>::FEDG_Q1():
   restriction[0](0,2) = 1.0/4.0;
   restriction[0](0,3) = 3.0/8.0;
   restriction[0](1,0) = -3.0/16.0;
-  restriction[0](1,1) = 0.0;
-  restriction[0](1,2) = 0.0;
   restriction[0](1,3) = -1.0/8.0;
   restriction[0](2,0) = 1.0/16.0;
-  restriction[0](2,1) = 0.0;
-  restriction[0](2,2) = 0.0;
-  restriction[0](2,3) = 0.0;
   restriction[0](3,0) = -3.0/16.0;
   restriction[0](3,1) = -1.0/8.0;
-  restriction[0](3,2) = 0.0;
-  restriction[0](3,3) = 0.0;
-  restriction[1](0,0) = 0.0;
   restriction[1](0,1) = -3.0/16.0;
   restriction[1](0,2) = -1.0/8.0;
-  restriction[1](0,3) = 0.0;
   restriction[1](1,0) = 3.0/8.0;
   restriction[1](1,1) = 9.0/16.0;
   restriction[1](1,2) = 3.0/8.0;
   restriction[1](1,3) = 1.0/4.0;
   restriction[1](2,0) = -1.0/8.0;
   restriction[1](2,1) = -3.0/16.0;
-  restriction[1](2,2) = 0.0;
-  restriction[1](2,3) = 0.0;
-  restriction[1](3,0) = 0.0;
   restriction[1](3,1) = 1.0/16.0;
-  restriction[1](3,2) = 0.0;
-  restriction[1](3,3) = 0.0;
-  restriction[2](0,0) = 0.0;
-  restriction[2](0,1) = 0.0;
   restriction[2](0,2) = 1.0/16.0;
-  restriction[2](0,3) = 0.0;
-  restriction[2](1,0) = 0.0;
-  restriction[2](1,1) = 0.0;
   restriction[2](1,2) = -3.0/16.0;
   restriction[2](1,3) = -1.0/8.0;
   restriction[2](2,0) = 1.0/4.0;
   restriction[2](2,1) = 3.0/8.0;
   restriction[2](2,2) = 9.0/16.0;
   restriction[2](2,3) = 3.0/8.0;
-  restriction[2](3,0) = 0.0;
   restriction[2](3,1) = -1.0/8.0;
   restriction[2](3,2) = -3.0/16.0;
-  restriction[2](3,3) = 0.0;
-  restriction[3](0,0) = 0.0;
-  restriction[3](0,1) = 0.0;
   restriction[3](0,2) = -1.0/8.0;
   restriction[3](0,3) = -3.0/16.0;
-  restriction[3](1,0) = 0.0;
-  restriction[3](1,1) = 0.0;
-  restriction[3](1,2) = 0.0;
   restriction[3](1,3) = 1.0/16.0;
   restriction[3](2,0) = -1.0/8.0;
-  restriction[3](2,1) = 0.0;
-  restriction[3](2,2) = 0.0;
   restriction[3](2,3) = -3.0/16.0;
   restriction[3](3,0) = 3.0/8.0;
   restriction[3](3,1) = 1.0/4.0;
@@ -1675,7 +1645,15 @@ FEDG_Q4<2>::FEDG_Q4():
   for (unsigned int child=0;
        child<GeometryInfo<2>::children_per_cell; ++ child)
     restriction[child].clear();
-  
+
+#if 1
+                                  // compiling the following function
+                                  // takes way too long. if you
+                                  // really need it, then change the
+                                  // preprocessor constant to let the
+                                  // following be compiled.
+  Assert (false, ExcNotImplemented());
+#else
   restriction[0](0,0) = 43681.0/82944.0;
   restriction[0](0,1) = 1045.0/82944.0;
   restriction[0](0,2) = 25.0/82944.0;
@@ -4176,6 +4154,7 @@ FEDG_Q4<2>::FEDG_Q4():
   restriction[3](24,22) = 25.0/6144.0;
   restriction[3](24,23) = 25.0/6144.0;
   restriction[3](24,24) = 25.0/4096.0;
+#endif
 };
 
 
@@ -4207,69 +4186,28 @@ FEDG_Q1<3>::FEDG_Q1():
   restriction[0](0,6) = 1.0/8.0;
   restriction[0](0,7) = 3.0/16.0;
   restriction[0](1,0) = -9.0/64.0;
-  restriction[0](1,1) = 0.0;
-  restriction[0](1,2) = 0.0;
   restriction[0](1,3) = -3.0/32.0;
   restriction[0](1,4) = -3.0/32.0;
-  restriction[0](1,5) = 0.0;
-  restriction[0](1,6) = 0.0;
   restriction[0](1,7) = -1.0/16.0;
   restriction[0](2,0) = 3.0/64.0;
-  restriction[0](2,1) = 0.0;
-  restriction[0](2,2) = 0.0;
-  restriction[0](2,3) = 0.0;
   restriction[0](2,4) = 1.0/32.0;
-  restriction[0](2,5) = 0.0;
-  restriction[0](2,6) = 0.0;
-  restriction[0](2,7) = 0.0;
   restriction[0](3,0) = -9.0/64.0;
   restriction[0](3,1) = -3.0/32.0;
-  restriction[0](3,2) = 0.0;
-  restriction[0](3,3) = 0.0;
   restriction[0](3,4) = -3.0/32.0;
   restriction[0](3,5) = -1.0/16.0;
-  restriction[0](3,6) = 0.0;
-  restriction[0](3,7) = 0.0;
   restriction[0](4,0) = -9.0/64.0;
   restriction[0](4,1) = -3.0/32.0;
   restriction[0](4,2) = -1.0/16.0;
   restriction[0](4,3) = -3.0/32.0;
-  restriction[0](4,4) = 0.0;
-  restriction[0](4,5) = 0.0;
-  restriction[0](4,6) = 0.0;
-  restriction[0](4,7) = 0.0;
   restriction[0](5,0) = 3.0/64.0;
-  restriction[0](5,1) = 0.0;
-  restriction[0](5,2) = 0.0;
   restriction[0](5,3) = 1.0/32.0;
-  restriction[0](5,4) = 0.0;
-  restriction[0](5,5) = 0.0;
-  restriction[0](5,6) = 0.0;
-  restriction[0](5,7) = 0.0;
   restriction[0](6,0) = -1.0/64.0;
-  restriction[0](6,1) = 0.0;
-  restriction[0](6,2) = 0.0;
-  restriction[0](6,3) = 0.0;
-  restriction[0](6,4) = 0.0;
-  restriction[0](6,5) = 0.0;
-  restriction[0](6,6) = 0.0;
-  restriction[0](6,7) = 0.0;
   restriction[0](7,0) = 3.0/64.0;
   restriction[0](7,1) = 1.0/32.0;
-  restriction[0](7,2) = 0.0;
-  restriction[0](7,3) = 0.0;
-  restriction[0](7,4) = 0.0;
-  restriction[0](7,5) = 0.0;
-  restriction[0](7,6) = 0.0;
-  restriction[0](7,7) = 0.0;
-  restriction[1](0,0) = 0.0;
   restriction[1](0,1) = -9.0/64.0;
   restriction[1](0,2) = -3.0/32.0;
-  restriction[1](0,3) = 0.0;
-  restriction[1](0,4) = 0.0;
   restriction[1](0,5) = -3.0/32.0;
   restriction[1](0,6) = -1.0/16.0;
-  restriction[1](0,7) = 0.0;
   restriction[1](1,0) = 9.0/32.0;
   restriction[1](1,1) = 27.0/64.0;
   restriction[1](1,2) = 9.0/32.0;
@@ -4280,66 +4218,23 @@ FEDG_Q1<3>::FEDG_Q1():
   restriction[1](1,7) = 1.0/8.0;
   restriction[1](2,0) = -3.0/32.0;
   restriction[1](2,1) = -9.0/64.0;
-  restriction[1](2,2) = 0.0;
-  restriction[1](2,3) = 0.0;
   restriction[1](2,4) = -1.0/16.0;
   restriction[1](2,5) = -3.0/32.0;
-  restriction[1](2,6) = 0.0;
-  restriction[1](2,7) = 0.0;
-  restriction[1](3,0) = 0.0;
   restriction[1](3,1) = 3.0/64.0;
-  restriction[1](3,2) = 0.0;
-  restriction[1](3,3) = 0.0;
-  restriction[1](3,4) = 0.0;
   restriction[1](3,5) = 1.0/32.0;
-  restriction[1](3,6) = 0.0;
-  restriction[1](3,7) = 0.0;
-  restriction[1](4,0) = 0.0;
   restriction[1](4,1) = 3.0/64.0;
   restriction[1](4,2) = 1.0/32.0;
-  restriction[1](4,3) = 0.0;
-  restriction[1](4,4) = 0.0;
-  restriction[1](4,5) = 0.0;
-  restriction[1](4,6) = 0.0;
-  restriction[1](4,7) = 0.0;
   restriction[1](5,0) = -3.0/32.0;
   restriction[1](5,1) = -9.0/64.0;
   restriction[1](5,2) = -3.0/32.0;
   restriction[1](5,3) = -1.0/16.0;
-  restriction[1](5,4) = 0.0;
-  restriction[1](5,5) = 0.0;
-  restriction[1](5,6) = 0.0;
-  restriction[1](5,7) = 0.0;
   restriction[1](6,0) = 1.0/32.0;
   restriction[1](6,1) = 3.0/64.0;
-  restriction[1](6,2) = 0.0;
-  restriction[1](6,3) = 0.0;
-  restriction[1](6,4) = 0.0;
-  restriction[1](6,5) = 0.0;
-  restriction[1](6,6) = 0.0;
-  restriction[1](6,7) = 0.0;
-  restriction[1](7,0) = 0.0;
   restriction[1](7,1) = -1.0/64.0;
-  restriction[1](7,2) = 0.0;
-  restriction[1](7,3) = 0.0;
-  restriction[1](7,4) = 0.0;
-  restriction[1](7,5) = 0.0;
-  restriction[1](7,6) = 0.0;
-  restriction[1](7,7) = 0.0;
-  restriction[2](0,0) = 0.0;
-  restriction[2](0,1) = 0.0;
   restriction[2](0,2) = 3.0/64.0;
-  restriction[2](0,3) = 0.0;
-  restriction[2](0,4) = 0.0;
-  restriction[2](0,5) = 0.0;
   restriction[2](0,6) = 1.0/32.0;
-  restriction[2](0,7) = 0.0;
-  restriction[2](1,0) = 0.0;
-  restriction[2](1,1) = 0.0;
   restriction[2](1,2) = -9.0/64.0;
   restriction[2](1,3) = -3.0/32.0;
-  restriction[2](1,4) = 0.0;
-  restriction[2](1,5) = 0.0;
   restriction[2](1,6) = -3.0/32.0;
   restriction[2](1,7) = -1.0/16.0;
   restriction[2](2,0) = 3.0/16.0;
@@ -4350,69 +4245,28 @@ FEDG_Q1<3>::FEDG_Q1():
   restriction[2](2,5) = 3.0/16.0;
   restriction[2](2,6) = 9.0/32.0;
   restriction[2](2,7) = 3.0/16.0;
-  restriction[2](3,0) = 0.0;
   restriction[2](3,1) = -3.0/32.0;
   restriction[2](3,2) = -9.0/64.0;
-  restriction[2](3,3) = 0.0;
-  restriction[2](3,4) = 0.0;
   restriction[2](3,5) = -1.0/16.0;
   restriction[2](3,6) = -3.0/32.0;
-  restriction[2](3,7) = 0.0;
-  restriction[2](4,0) = 0.0;
-  restriction[2](4,1) = 0.0;
   restriction[2](4,2) = -1.0/64.0;
-  restriction[2](4,3) = 0.0;
-  restriction[2](4,4) = 0.0;
-  restriction[2](4,5) = 0.0;
-  restriction[2](4,6) = 0.0;
-  restriction[2](4,7) = 0.0;
-  restriction[2](5,0) = 0.0;
-  restriction[2](5,1) = 0.0;
   restriction[2](5,2) = 3.0/64.0;
   restriction[2](5,3) = 1.0/32.0;
-  restriction[2](5,4) = 0.0;
-  restriction[2](5,5) = 0.0;
-  restriction[2](5,6) = 0.0;
-  restriction[2](5,7) = 0.0;
   restriction[2](6,0) = -1.0/16.0;
   restriction[2](6,1) = -3.0/32.0;
   restriction[2](6,2) = -9.0/64.0;
   restriction[2](6,3) = -3.0/32.0;
-  restriction[2](6,4) = 0.0;
-  restriction[2](6,5) = 0.0;
-  restriction[2](6,6) = 0.0;
-  restriction[2](6,7) = 0.0;
-  restriction[2](7,0) = 0.0;
   restriction[2](7,1) = 1.0/32.0;
   restriction[2](7,2) = 3.0/64.0;
-  restriction[2](7,3) = 0.0;
-  restriction[2](7,4) = 0.0;
-  restriction[2](7,5) = 0.0;
-  restriction[2](7,6) = 0.0;
-  restriction[2](7,7) = 0.0;
-  restriction[3](0,0) = 0.0;
-  restriction[3](0,1) = 0.0;
   restriction[3](0,2) = -3.0/32.0;
   restriction[3](0,3) = -9.0/64.0;
-  restriction[3](0,4) = 0.0;
-  restriction[3](0,5) = 0.0;
   restriction[3](0,6) = -1.0/16.0;
   restriction[3](0,7) = -3.0/32.0;
-  restriction[3](1,0) = 0.0;
-  restriction[3](1,1) = 0.0;
-  restriction[3](1,2) = 0.0;
   restriction[3](1,3) = 3.0/64.0;
-  restriction[3](1,4) = 0.0;
-  restriction[3](1,5) = 0.0;
-  restriction[3](1,6) = 0.0;
   restriction[3](1,7) = 1.0/32.0;
   restriction[3](2,0) = -3.0/32.0;
-  restriction[3](2,1) = 0.0;
-  restriction[3](2,2) = 0.0;
   restriction[3](2,3) = -9.0/64.0;
   restriction[3](2,4) = -1.0/16.0;
-  restriction[3](2,5) = 0.0;
-  restriction[3](2,6) = 0.0;
   restriction[3](2,7) = -3.0/32.0;
   restriction[3](3,0) = 9.0/32.0;
   restriction[3](3,1) = 3.0/16.0;
@@ -4422,70 +4276,24 @@ FEDG_Q1<3>::FEDG_Q1():
   restriction[3](3,5) = 1.0/8.0;
   restriction[3](3,6) = 3.0/16.0;
   restriction[3](3,7) = 9.0/32.0;
-  restriction[3](4,0) = 0.0;
-  restriction[3](4,1) = 0.0;
   restriction[3](4,2) = 1.0/32.0;
   restriction[3](4,3) = 3.0/64.0;
-  restriction[3](4,4) = 0.0;
-  restriction[3](4,5) = 0.0;
-  restriction[3](4,6) = 0.0;
-  restriction[3](4,7) = 0.0;
-  restriction[3](5,0) = 0.0;
-  restriction[3](5,1) = 0.0;
-  restriction[3](5,2) = 0.0;
   restriction[3](5,3) = -1.0/64.0;
-  restriction[3](5,4) = 0.0;
-  restriction[3](5,5) = 0.0;
-  restriction[3](5,6) = 0.0;
-  restriction[3](5,7) = 0.0;
   restriction[3](6,0) = 1.0/32.0;
-  restriction[3](6,1) = 0.0;
-  restriction[3](6,2) = 0.0;
   restriction[3](6,3) = 3.0/64.0;
-  restriction[3](6,4) = 0.0;
-  restriction[3](6,5) = 0.0;
-  restriction[3](6,6) = 0.0;
-  restriction[3](6,7) = 0.0;
   restriction[3](7,0) = -3.0/32.0;
   restriction[3](7,1) = -1.0/16.0;
   restriction[3](7,2) = -3.0/32.0;
   restriction[3](7,3) = -9.0/64.0;
-  restriction[3](7,4) = 0.0;
-  restriction[3](7,5) = 0.0;
-  restriction[3](7,6) = 0.0;
-  restriction[3](7,7) = 0.0;
-  restriction[4](0,0) = 0.0;
-  restriction[4](0,1) = 0.0;
-  restriction[4](0,2) = 0.0;
-  restriction[4](0,3) = 0.0;
   restriction[4](0,4) = -9.0/64.0;
   restriction[4](0,5) = -3.0/32.0;
   restriction[4](0,6) = -1.0/16.0;
   restriction[4](0,7) = -3.0/32.0;
-  restriction[4](1,0) = 0.0;
-  restriction[4](1,1) = 0.0;
-  restriction[4](1,2) = 0.0;
-  restriction[4](1,3) = 0.0;
   restriction[4](1,4) = 3.0/64.0;
-  restriction[4](1,5) = 0.0;
-  restriction[4](1,6) = 0.0;
   restriction[4](1,7) = 1.0/32.0;
-  restriction[4](2,0) = 0.0;
-  restriction[4](2,1) = 0.0;
-  restriction[4](2,2) = 0.0;
-  restriction[4](2,3) = 0.0;
   restriction[4](2,4) = -1.0/64.0;
-  restriction[4](2,5) = 0.0;
-  restriction[4](2,6) = 0.0;
-  restriction[4](2,7) = 0.0;
-  restriction[4](3,0) = 0.0;
-  restriction[4](3,1) = 0.0;
-  restriction[4](3,2) = 0.0;
-  restriction[4](3,3) = 0.0;
   restriction[4](3,4) = 3.0/64.0;
   restriction[4](3,5) = 1.0/32.0;
-  restriction[4](3,6) = 0.0;
-  restriction[4](3,7) = 0.0;
   restriction[4](4,0) = 9.0/32.0;
   restriction[4](4,1) = 3.0/16.0;
   restriction[4](4,2) = 1.0/8.0;
@@ -4495,69 +4303,28 @@ FEDG_Q1<3>::FEDG_Q1():
   restriction[4](4,6) = 3.0/16.0;
   restriction[4](4,7) = 9.0/32.0;
   restriction[4](5,0) = -3.0/32.0;
-  restriction[4](5,1) = 0.0;
-  restriction[4](5,2) = 0.0;
   restriction[4](5,3) = -1.0/16.0;
   restriction[4](5,4) = -9.0/64.0;
-  restriction[4](5,5) = 0.0;
-  restriction[4](5,6) = 0.0;
   restriction[4](5,7) = -3.0/32.0;
   restriction[4](6,0) = 1.0/32.0;
-  restriction[4](6,1) = 0.0;
-  restriction[4](6,2) = 0.0;
-  restriction[4](6,3) = 0.0;
   restriction[4](6,4) = 3.0/64.0;
-  restriction[4](6,5) = 0.0;
-  restriction[4](6,6) = 0.0;
-  restriction[4](6,7) = 0.0;
   restriction[4](7,0) = -3.0/32.0;
   restriction[4](7,1) = -1.0/16.0;
-  restriction[4](7,2) = 0.0;
-  restriction[4](7,3) = 0.0;
   restriction[4](7,4) = -9.0/64.0;
   restriction[4](7,5) = -3.0/32.0;
-  restriction[4](7,6) = 0.0;
-  restriction[4](7,7) = 0.0;
-  restriction[5](0,0) = 0.0;
-  restriction[5](0,1) = 0.0;
-  restriction[5](0,2) = 0.0;
-  restriction[5](0,3) = 0.0;
-  restriction[5](0,4) = 0.0;
   restriction[5](0,5) = 3.0/64.0;
   restriction[5](0,6) = 1.0/32.0;
-  restriction[5](0,7) = 0.0;
-  restriction[5](1,0) = 0.0;
-  restriction[5](1,1) = 0.0;
-  restriction[5](1,2) = 0.0;
-  restriction[5](1,3) = 0.0;
   restriction[5](1,4) = -3.0/32.0;
   restriction[5](1,5) = -9.0/64.0;
   restriction[5](1,6) = -3.0/32.0;
   restriction[5](1,7) = -1.0/16.0;
-  restriction[5](2,0) = 0.0;
-  restriction[5](2,1) = 0.0;
-  restriction[5](2,2) = 0.0;
-  restriction[5](2,3) = 0.0;
   restriction[5](2,4) = 1.0/32.0;
   restriction[5](2,5) = 3.0/64.0;
-  restriction[5](2,6) = 0.0;
-  restriction[5](2,7) = 0.0;
-  restriction[5](3,0) = 0.0;
-  restriction[5](3,1) = 0.0;
-  restriction[5](3,2) = 0.0;
-  restriction[5](3,3) = 0.0;
-  restriction[5](3,4) = 0.0;
   restriction[5](3,5) = -1.0/64.0;
-  restriction[5](3,6) = 0.0;
-  restriction[5](3,7) = 0.0;
-  restriction[5](4,0) = 0.0;
   restriction[5](4,1) = -3.0/32.0;
   restriction[5](4,2) = -1.0/16.0;
-  restriction[5](4,3) = 0.0;
-  restriction[5](4,4) = 0.0;
   restriction[5](4,5) = -9.0/64.0;
   restriction[5](4,6) = -3.0/32.0;
-  restriction[5](4,7) = 0.0;
   restriction[5](5,0) = 3.0/16.0;
   restriction[5](5,1) = 9.0/32.0;
   restriction[5](5,2) = 3.0/16.0;
@@ -4568,66 +4335,23 @@ FEDG_Q1<3>::FEDG_Q1():
   restriction[5](5,7) = 3.0/16.0;
   restriction[5](6,0) = -1.0/16.0;
   restriction[5](6,1) = -3.0/32.0;
-  restriction[5](6,2) = 0.0;
-  restriction[5](6,3) = 0.0;
   restriction[5](6,4) = -3.0/32.0;
   restriction[5](6,5) = -9.0/64.0;
-  restriction[5](6,6) = 0.0;
-  restriction[5](6,7) = 0.0;
-  restriction[5](7,0) = 0.0;
   restriction[5](7,1) = 1.0/32.0;
-  restriction[5](7,2) = 0.0;
-  restriction[5](7,3) = 0.0;
-  restriction[5](7,4) = 0.0;
   restriction[5](7,5) = 3.0/64.0;
-  restriction[5](7,6) = 0.0;
-  restriction[5](7,7) = 0.0;
-  restriction[6](0,0) = 0.0;
-  restriction[6](0,1) = 0.0;
-  restriction[6](0,2) = 0.0;
-  restriction[6](0,3) = 0.0;
-  restriction[6](0,4) = 0.0;
-  restriction[6](0,5) = 0.0;
   restriction[6](0,6) = -1.0/64.0;
-  restriction[6](0,7) = 0.0;
-  restriction[6](1,0) = 0.0;
-  restriction[6](1,1) = 0.0;
-  restriction[6](1,2) = 0.0;
-  restriction[6](1,3) = 0.0;
-  restriction[6](1,4) = 0.0;
-  restriction[6](1,5) = 0.0;
   restriction[6](1,6) = 3.0/64.0;
   restriction[6](1,7) = 1.0/32.0;
-  restriction[6](2,0) = 0.0;
-  restriction[6](2,1) = 0.0;
-  restriction[6](2,2) = 0.0;
-  restriction[6](2,3) = 0.0;
   restriction[6](2,4) = -1.0/16.0;
   restriction[6](2,5) = -3.0/32.0;
   restriction[6](2,6) = -9.0/64.0;
   restriction[6](2,7) = -3.0/32.0;
-  restriction[6](3,0) = 0.0;
-  restriction[6](3,1) = 0.0;
-  restriction[6](3,2) = 0.0;
-  restriction[6](3,3) = 0.0;
-  restriction[6](3,4) = 0.0;
   restriction[6](3,5) = 1.0/32.0;
   restriction[6](3,6) = 3.0/64.0;
-  restriction[6](3,7) = 0.0;
-  restriction[6](4,0) = 0.0;
-  restriction[6](4,1) = 0.0;
   restriction[6](4,2) = 1.0/32.0;
-  restriction[6](4,3) = 0.0;
-  restriction[6](4,4) = 0.0;
-  restriction[6](4,5) = 0.0;
   restriction[6](4,6) = 3.0/64.0;
-  restriction[6](4,7) = 0.0;
-  restriction[6](5,0) = 0.0;
-  restriction[6](5,1) = 0.0;
   restriction[6](5,2) = -3.0/32.0;
   restriction[6](5,3) = -1.0/16.0;
-  restriction[6](5,4) = 0.0;
-  restriction[6](5,5) = 0.0;
   restriction[6](5,6) = -9.0/64.0;
   restriction[6](5,7) = -3.0/32.0;
   restriction[6](6,0) = 1.0/8.0;
@@ -4638,69 +4362,28 @@ FEDG_Q1<3>::FEDG_Q1():
   restriction[6](6,5) = 9.0/32.0;
   restriction[6](6,6) = 27.0/64.0;
   restriction[6](6,7) = 9.0/32.0;
-  restriction[6](7,0) = 0.0;
   restriction[6](7,1) = -1.0/16.0;
   restriction[6](7,2) = -3.0/32.0;
-  restriction[6](7,3) = 0.0;
-  restriction[6](7,4) = 0.0;
   restriction[6](7,5) = -3.0/32.0;
   restriction[6](7,6) = -9.0/64.0;
-  restriction[6](7,7) = 0.0;
-  restriction[7](0,0) = 0.0;
-  restriction[7](0,1) = 0.0;
-  restriction[7](0,2) = 0.0;
-  restriction[7](0,3) = 0.0;
-  restriction[7](0,4) = 0.0;
-  restriction[7](0,5) = 0.0;
   restriction[7](0,6) = 1.0/32.0;
   restriction[7](0,7) = 3.0/64.0;
-  restriction[7](1,0) = 0.0;
-  restriction[7](1,1) = 0.0;
-  restriction[7](1,2) = 0.0;
-  restriction[7](1,3) = 0.0;
-  restriction[7](1,4) = 0.0;
-  restriction[7](1,5) = 0.0;
-  restriction[7](1,6) = 0.0;
   restriction[7](1,7) = -1.0/64.0;
-  restriction[7](2,0) = 0.0;
-  restriction[7](2,1) = 0.0;
-  restriction[7](2,2) = 0.0;
-  restriction[7](2,3) = 0.0;
   restriction[7](2,4) = 1.0/32.0;
-  restriction[7](2,5) = 0.0;
-  restriction[7](2,6) = 0.0;
   restriction[7](2,7) = 3.0/64.0;
-  restriction[7](3,0) = 0.0;
-  restriction[7](3,1) = 0.0;
-  restriction[7](3,2) = 0.0;
-  restriction[7](3,3) = 0.0;
   restriction[7](3,4) = -3.0/32.0;
   restriction[7](3,5) = -1.0/16.0;
   restriction[7](3,6) = -3.0/32.0;
   restriction[7](3,7) = -9.0/64.0;
-  restriction[7](4,0) = 0.0;
-  restriction[7](4,1) = 0.0;
   restriction[7](4,2) = -1.0/16.0;
   restriction[7](4,3) = -3.0/32.0;
-  restriction[7](4,4) = 0.0;
-  restriction[7](4,5) = 0.0;
   restriction[7](4,6) = -3.0/32.0;
   restriction[7](4,7) = -9.0/64.0;
-  restriction[7](5,0) = 0.0;
-  restriction[7](5,1) = 0.0;
-  restriction[7](5,2) = 0.0;
   restriction[7](5,3) = 1.0/32.0;
-  restriction[7](5,4) = 0.0;
-  restriction[7](5,5) = 0.0;
-  restriction[7](5,6) = 0.0;
   restriction[7](5,7) = 3.0/64.0;
   restriction[7](6,0) = -1.0/16.0;
-  restriction[7](6,1) = 0.0;
-  restriction[7](6,2) = 0.0;
   restriction[7](6,3) = -3.0/32.0;
   restriction[7](6,4) = -3.0/32.0;
-  restriction[7](6,5) = 0.0;
-  restriction[7](6,6) = 0.0;
   restriction[7](6,7) = -9.0/64.0;
   restriction[7](7,0) = 3.0/16.0;
   restriction[7](7,1) = 1.0/8.0;

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.