From: wolf Date: Wed, 24 May 2000 12:47:41 +0000 (+0000) Subject: Eliminate stores of zeroes. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=a7d9fb48f763d03fef450461b31c6fde1ef0b665;p=dealii-svn.git Eliminate stores of zeroes. git-svn-id: https://svn.dealii.org/trunk@2937 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/deal.II/source/fe/fe_lib.dg.cc b/deal.II/deal.II/source/fe/fe_lib.dg.cc index 8849e66748..0a3386f340 100644 --- a/deal.II/deal.II/source/fe/fe_lib.dg.cc +++ b/deal.II/deal.II/source/fe/fe_lib.dg.cc @@ -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::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;