]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Suppress very small constraints to avoid cluttering the sparsity pattern for no good...
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Fri, 16 Mar 2007 17:16:07 +0000 (17:16 +0000)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Fri, 16 Mar 2007 17:16:07 +0000 (17:16 +0000)
git-svn-id: https://svn.dealii.org/trunk@14586 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/source/dofs/dof_tools.cc
tests/hp/crash_12/cmp/generic
tests/hp/hp_constraints_q_03/cmp/generic
tests/hp/hp_constraints_q_05/cmp/generic
tests/hp/hp_constraints_q_system_03/cmp/generic
tests/hp/hp_constraints_q_system_05/cmp/generic
tests/hp/hp_constraints_q_system_x_02/cmp/generic
tests/hp/hp_constraints_q_system_x_03/cmp/generic
tests/hp/hp_constraints_q_system_x_05/cmp/generic

index 0a4291b07ecdbb1b184400e6c2be0437c62a4751..a47c140e62706a711af58639ea6897c1bb64f748 100644 (file)
@@ -1483,12 +1483,20 @@ namespace internal
 
       
                                        /**
-                                        * This method removes zero
+                                       * Copy constraints into a constraint
+                                       * matrix object.
+                                       *
+                                        * This function removes zero
                                         * constraints and those, which
                                         * constrain a DoF which was
                                         * already eliminated in one of
                                         * the previous steps of the hp
                                         * hanging node procedure.
+                                       *
+                                       * It also suppresses very small
+                                       * entries in the constraint matrix to
+                                       * avoid making the sparsity pattern
+                                       * fuller than necessary.
                                         */
 #ifdef DEAL_II_ANON_NAMESPACE_BUG
       static
@@ -1528,10 +1536,36 @@ namespace internal
                     }
 
               if (constraint_already_satisfied == false)
-                { 
+                {
+                                                  // add up the absolute
+                                                  // values of all
+                                                  // constraints in this line
+                                                  // to get a measure of
+                                                  // their absolute size
+                 double abs_sum = 0;
+                 for (unsigned int i=0; i<n_master_dofs; ++i)
+                   abs_sum += std::abs (face_constraints(row,i));
+                 
+                                                  // then enter those
+                                                  // constraints that are
+                                                  // larger than
+                                                  // 1e-14*abs_sum. everything
+                                                  // else probably originated
+                                                  // from inexact inversion
+                                                  // of matrices and similar
+                                                  // effects. having those
+                                                  // constraints in here will
+                                                  // only lead to problems
+                                                  // because it makes
+                                                  // sparsity patterns fuller
+                                                  // than necessary without
+                                                  // producing any
+                                                  // significant effect
                   constraints.add_line (slave_dofs[row]);
                   for (unsigned int i=0; i<n_master_dofs; ++i)
-                    if (face_constraints(row,i) != 0)
+                    if ((face_constraints(row,i) != 0)
+                       &&
+                       (std::fabs(face_constraints(row,i)) >= 1e-14*abs_sum))
                       {
 #ifdef WOLFGANG
                         std::cout << "   " << slave_dofs[row]
index 008c422c192fe7591951ce27f785d3babfdda888..1141f98ca016055d87db62af5c9bd02a816b4021 100644 (file)
@@ -256,9 +256,7 @@ DEAL::Testing FE_Q<2>(4) vs. FE_Q<2>(3)
     7 6:  1.50
     7 23:  0.250
     7 51:  0.250
-    25 5:  5.55e-17
-    25 6:  1.000
-    25 23:  -1.39e-17
+    25 6:  1.00
     27 5:  0.988
     27 6:  -0.185
     27 23:  0.185
@@ -2888,21 +2886,12 @@ DEAL::Testing FE_Q<3>(3) vs. FE_Q<3>(2)
     27 26:  1.00
     27 105:  -0.333
     27 165:  0.333
-    33 6:  -2.78e-17
     33 22:  -0.333
     33 26:  0.333
     33 32:  1.00
-    33 60:  6.94e-18
-    33 105:  -1.73e-17
-    33 135:  -6.94e-18
     34 6:  -0.333
     34 14:  0.333
-    34 22:  -2.78e-17
     34 32:  1.00
-    34 60:  6.94e-18
-    34 105:  -3.47e-18
-    34 135:  -1.73e-17
-    34 165:  -1.73e-18
     35 6:  -0.333
     35 14:  0.333
     35 22:  -0.333
@@ -5860,61 +5849,22 @@ DEAL::Testing FE_Q<3>(4) vs. FE_Q<3>(3)
     39 38:  1.50
     39 233:  0.250
     39 401:  0.250
-    51 7:  -1.39e-17
-    51 19:  -1.04e-17
-    51 20:  6.94e-18
     51 31:  0.250
     51 37:  0.250
-    51 38:  6.94e-18
     51 49:  -1.00
     51 50:  1.50
-    51 53:  5.55e-17
-    51 121:  -8.67e-19
-    51 317:  1.73e-18
-    51 401:  8.67e-19
-    54 7:  -1.39e-17
-    54 8:  2.78e-17
-    54 19:  1.39e-17
-    54 31:  -1.39e-17
     54 32:  0.250
-    54 37:  1.39e-17
     54 38:  0.250
-    54 49:  1.11e-16
-    54 50:  -1.11e-16
     54 52:  -1.000
     54 53:  1.50
-    54 121:  -3.47e-18
-    54 233:  -1.73e-18
-    54 317:  1.73e-18
-    54 401:  1.73e-18
     55 7:  0.250
-    55 8:  6.94e-18
     55 19:  0.250
-    55 20:  6.94e-18
-    55 31:  -6.94e-18
-    55 37:  -6.94e-18
-    55 38:  6.94e-18
     55 49:  -1.000
     55 52:  1.50
-    55 53:  5.55e-17
-    55 121:  -8.67e-19
-    55 317:  1.73e-18
-    56 7:  -1.39e-17
     56 8:  0.250
-    56 19:  6.94e-18
     56 20:  0.250
-    56 31:  -1.39e-17
-    56 32:  2.78e-17
-    56 37:  2.08e-17
-    56 38:  -1.39e-17
-    56 49:  1.11e-16
     56 50:  -1.000
-    56 52:  -1.11e-16
     56 53:  1.50
-    56 121:  -6.94e-18
-    56 233:  -1.73e-18
-    56 317:  1.73e-18
-    56 401:  1.73e-18
     57 7:  -0.250
     57 8:  0.375
     57 19:  -0.250
@@ -5931,23 +5881,9 @@ DEAL::Testing FE_Q<3>(4) vs. FE_Q<3>(3)
     57 233:  0.0625
     57 317:  0.0625
     57 401:  0.0625
-    123 7:  5.55e-17
-    123 8:  1.000
-    123 121:  -1.39e-17
-    125 31:  5.55e-17
-    125 32:  1.000
-    125 121:  -1.39e-17
-    127 7:  -6.94e-18
-    127 19:  -6.94e-18
-    127 20:  6.94e-18
-    127 31:  -6.94e-18
-    127 38:  6.94e-18
-    127 49:  1.11e-16
-    127 50:  -5.55e-17
-    127 52:  -5.55e-17
+    123 8:  1.00
+    125 32:  1.00
     127 53:  1.00
-    127 317:  3.47e-18
-    127 401:  -8.67e-19
     129 7:  0.988
     129 8:  -0.185
     129 121:  0.185
@@ -5956,33 +5892,14 @@ DEAL::Testing FE_Q<3>(4) vs. FE_Q<3>(3)
     130 8:  0.296
     130 121:  -0.0741
     130 233:  -0.0123
-    137 7:  6.94e-18
-    137 19:  -1.04e-17
-    137 31:  3.47e-18
     137 32:  0.185
-    137 37:  -6.94e-18
     137 38:  0.0123
-    137 49:  2.78e-17
-    137 50:  2.78e-17
     137 52:  0.988
     137 53:  -0.185
-    137 121:  -5.20e-18
-    137 317:  5.20e-18
-    138 7:  -1.39e-17
-    138 8:  3.47e-18
-    138 19:  -1.39e-17
-    138 20:  3.47e-18
-    138 31:  6.94e-18
     138 32:  -0.0741
-    138 37:  -6.07e-18
     138 38:  -0.0123
-    138 49:  1.11e-16
     138 52:  0.790
     138 53:  0.296
-    138 121:  8.67e-19
-    138 233:  6.51e-19
-    138 317:  2.60e-18
-    138 401:  -5.42e-19
     145 31:  0.988
     145 32:  -0.185
     145 121:  0.185
@@ -5991,33 +5908,14 @@ DEAL::Testing FE_Q<3>(4) vs. FE_Q<3>(3)
     146 32:  0.296
     146 121:  -0.0741
     146 317:  -0.0123
-    149 7:  3.47e-18
     149 8:  0.185
-    149 19:  -3.47e-18
     149 20:  0.0123
-    149 31:  3.47e-18
-    149 32:  6.94e-18
-    149 37:  -1.04e-17
-    149 49:  2.78e-17
     149 50:  0.988
-    149 52:  2.78e-17
     149 53:  -0.185
-    149 121:  -5.20e-18
-    149 233:  5.20e-18
-    149 401:  8.67e-19
     150 8:  -0.0741
-    150 19:  -7.81e-18
     150 20:  -0.0123
-    150 32:  3.47e-18
-    150 37:  -1.39e-17
-    150 38:  3.47e-18
-    150 49:  1.11e-16
     150 50:  0.790
     150 53:  0.296
-    150 121:  8.67e-19
-    150 233:  1.73e-18
-    150 317:  1.08e-18
-    150 401:  -5.42e-19
     153 7:  0.183
     153 8:  -0.0343
     153 19:  0.0122
@@ -6082,9 +5980,7 @@ DEAL::Testing FE_Q<3>(4) vs. FE_Q<3>(3)
     156 233:  0.000914
     156 317:  0.000914
     156 401:  0.000152
-    235 37:  5.55e-17
-    235 38:  1.000
-    235 233:  -1.39e-17
+    235 38:  1.00
     237 7:  -0.790
     237 8:  1.48
     237 121:  0.185
@@ -6093,34 +5989,14 @@ DEAL::Testing FE_Q<3>(4) vs. FE_Q<3>(3)
     238 8:  1.30
     238 121:  0.259
     238 233:  0.432
-    243 7:  -2.78e-17
-    243 8:  2.78e-17
-    243 31:  -2.78e-17
     243 32:  0.185
-    243 37:  6.94e-18
     243 38:  0.123
-    243 49:  1.11e-16
-    243 50:  -2.22e-16
     243 52:  -0.790
     243 53:  1.48
-    243 121:  -3.47e-18
-    243 233:  -1.73e-18
-    243 317:  1.73e-18
-    244 8:  -2.78e-17
-    244 19:  6.94e-18
-    244 20:  -1.39e-17
-    244 31:  -1.39e-17
     244 32:  0.259
-    244 37:  6.94e-18
     244 38:  0.432
-    244 49:  5.55e-17
-    244 50:  -1.11e-16
     244 52:  -0.988
     244 53:  1.30
-    244 121:  -3.47e-18
-    244 233:  -8.67e-18
-    244 317:  1.73e-18
-    244 401:  3.47e-18
     249 37:  0.988
     249 38:  -0.185
     249 233:  0.185
@@ -6193,9 +6069,7 @@ DEAL::Testing FE_Q<3>(4) vs. FE_Q<3>(3)
     256 233:  -0.0320
     256 317:  0.00320
     256 401:  0.00533
-    319 19:  5.55e-17
-    319 20:  1.000
-    319 317:  -1.39e-17
+    319 20:  1.00
     321 19:  0.988
     321 20:  -0.185
     321 317:  0.185
@@ -6212,30 +6086,14 @@ DEAL::Testing FE_Q<3>(4) vs. FE_Q<3>(3)
     330 32:  1.30
     330 121:  0.259
     330 317:  0.432
-    333 7:  -2.78e-17
     333 8:  0.185
     333 20:  0.123
-    333 31:  -2.78e-17
-    333 32:  2.78e-17
-    333 37:  6.94e-18
-    333 49:  1.11e-16
     333 50:  -0.790
-    333 52:  -2.22e-16
     333 53:  1.48
-    333 121:  -3.47e-18
-    333 317:  1.73e-18
-    334 7:  -1.39e-17
     334 8:  0.259
     334 20:  0.432
-    334 37:  1.39e-17
-    334 38:  -1.39e-17
-    334 49:  5.55e-17
     334 50:  -0.988
-    334 52:  -1.11e-16
     334 53:  1.30
-    334 121:  -3.47e-18
-    334 317:  -3.47e-18
-    334 401:  3.47e-18
     337 7:  0.183
     337 8:  -0.0343
     337 19:  0.122
@@ -8994,40 +8852,18 @@ DEAL::Testing FE_Q<3>(5) vs. FE_Q<3>(2)
     51 48:  1.000
     51 257:  -0.600
     51 317:  0.600
-    69 8:  1.11e-16
-    69 24:  -2.78e-17
     69 40:  -0.600
     69 48:  0.1000
     69 68:  1.50
-    69 212:  -1.11e-16
-    69 257:  4.16e-17
-    69 287:  2.78e-17
-    69 317:  -1.04e-17
-    70 8:  1.11e-16
-    70 24:  -2.78e-17
     70 40:  -0.800
     70 48:  0.300
     70 68:  1.50
-    70 212:  -1.11e-16
-    70 257:  1.39e-17
-    70 287:  2.78e-17
-    70 317:  -3.47e-18
     71 40:  -0.600
     71 48:  0.600
     71 68:  1.000
-    71 212:  -2.78e-17
-    71 257:  -2.08e-17
-    71 287:  6.94e-18
-    71 317:  5.20e-18
     72 8:  -0.600
     72 24:  0.1000
-    72 40:  1.11e-16
-    72 48:  -2.78e-17
     72 68:  1.50
-    72 212:  -1.11e-16
-    72 257:  2.78e-17
-    72 287:  4.16e-17
-    72 317:  -1.04e-17
     73 8:  -0.900
     73 24:  0.150
     73 40:  -0.900
@@ -9057,13 +8893,7 @@ DEAL::Testing FE_Q<3>(5) vs. FE_Q<3>(2)
     75 317:  0.0600
     76 8:  -0.800
     76 24:  0.300
-    76 40:  1.11e-16
-    76 48:  -2.78e-17
     76 68:  1.50
-    76 212:  -1.11e-16
-    76 257:  2.78e-17
-    76 287:  1.39e-17
-    76 317:  -3.47e-18
     77 8:  -1.20
     77 24:  0.450
     77 40:  -0.900
@@ -9094,10 +8924,6 @@ DEAL::Testing FE_Q<3>(5) vs. FE_Q<3>(2)
     80 8:  -0.600
     80 24:  0.600
     80 68:  1.000
-    80 212:  -5.55e-17
-    80 257:  1.39e-17
-    80 287:  -1.39e-17
-    80 317:  3.47e-18
     81 8:  -0.900
     81 24:  0.900
     81 40:  -0.600
@@ -9275,94 +9101,30 @@ DEAL::Testing FE_Q<3>(5) vs. FE_Q<3>(3)
     51 49:  2.00
     51 324:  0.600
     51 492:  0.400
-    70 9:  5.55e-17
-    70 24:  -2.78e-17
-    70 25:  2.78e-17
     70 40:  0.400
-    70 41:  1.67e-16
     70 48:  0.100
-    70 49:  4.16e-17
     70 68:  -1.50
     70 69:  2.00
-    70 72:  -2.22e-16
-    70 73:  4.44e-16
-    70 324:  6.94e-18
-    70 408:  3.47e-18
-    70 492:  4.34e-18
-    71 8:  -1.11e-16
-    71 9:  1.67e-16
-    71 24:  -4.16e-17
-    71 25:  4.16e-17
     71 40:  0.600
-    71 41:  2.22e-16
     71 48:  0.400
-    71 49:  8.33e-17
     71 68:  -2.00
     71 69:  2.00
-    71 72:  -4.44e-16
-    71 73:  6.66e-16
-    71 324:  6.94e-18
-    71 492:  3.47e-18
-    74 8:  1.11e-16
-    74 9:  2.78e-17
-    74 24:  6.94e-18
-    74 25:  -4.16e-17
     74 41:  0.400
-    74 48:  1.39e-17
     74 49:  0.100
-    74 68:  4.44e-16
     74 72:  -1.50
     74 73:  2.00
-    74 212:  -1.39e-17
-    74 324:  1.04e-17
-    74 408:  6.94e-18
-    74 492:  -8.67e-19
-    75 8:  1.39e-16
-    75 9:  2.78e-17
-    75 24:  1.39e-17
-    75 25:  -4.16e-17
-    75 40:  5.55e-17
     75 41:  0.600
-    75 48:  2.78e-17
     75 49:  0.400
-    75 68:  1.11e-16
-    75 69:  1.11e-16
     75 72:  -2.00
     75 73:  2.00
-    75 212:  -2.08e-17
-    75 324:  6.94e-18
-    75 408:  3.47e-18
-    75 492:  -2.60e-18
     76 8:  0.400
-    76 9:  1.67e-16
-    76 24:  0.1000
-    76 25:  4.16e-17
-    76 40:  -2.78e-17
-    76 41:  1.11e-16
-    76 48:  -2.78e-17
-    76 49:  2.78e-17
+    76 24:  0.100
     76 68:  -1.50
-    76 69:  -2.22e-16
     76 72:  2.00
-    76 73:  4.44e-16
-    76 324:  3.47e-18
-    76 408:  3.47e-18
-    76 492:  1.73e-18
-    77 8:  8.33e-17
     77 9:  0.400
-    77 24:  1.39e-17
     77 25:  0.100
-    77 40:  2.78e-17
-    77 41:  -5.55e-17
-    77 48:  1.39e-17
-    77 49:  -4.16e-17
-    77 68:  3.33e-16
     77 69:  -1.50
     77 73:  2.00
-    77 212:  -1.39e-17
-    77 324:  3.47e-18
-    77 408:  6.94e-18
-    77 492:  -2.60e-18
     78 8:  -0.600
     78 9:  0.800
     78 24:  -0.150
@@ -9396,37 +9158,13 @@ DEAL::Testing FE_Q<3>(5) vs. FE_Q<3>(3)
     79 408:  0.0600
     79 492:  0.0400
     80 8:  0.600
-    80 9:  2.22e-16
     80 24:  0.400
-    80 25:  9.71e-17
-    80 40:  -1.67e-16
-    80 41:  1.67e-16
-    80 48:  -4.16e-17
-    80 49:  4.16e-17
     80 68:  -2.00
-    80 69:  -4.44e-16
     80 72:  2.00
-    80 73:  6.66e-16
-    80 212:  1.39e-17
-    80 324:  3.47e-18
-    80 408:  1.73e-17
-    80 492:  1.73e-18
-    81 8:  1.11e-16
     81 9:  0.600
-    81 24:  2.78e-17
     81 25:  0.400
-    81 40:  5.55e-17
-    81 41:  -5.55e-17
-    81 48:  1.39e-17
-    81 49:  -4.16e-17
-    81 68:  2.22e-16
     81 69:  -2.00
-    81 72:  1.11e-16
     81 73:  2.00
-    81 212:  -2.08e-17
-    81 324:  6.94e-18
-    81 408:  3.47e-18
-    81 492:  -4.34e-18
     82 8:  -0.900
     82 9:  1.20
     82 24:  -0.600
@@ -9960,139 +9698,36 @@ DEAL::Testing FE_Q<3>(5) vs. FE_Q<3>(4)
     51 50:  2.00
     51 437:  -0.200
     51 797:  0.200
-    71 8:  1.73e-18
-    71 9:  2.78e-17
-    71 10:  -2.08e-17
-    71 24:  3.47e-18
-    71 25:  3.47e-18
-    71 26:  6.94e-18
     71 40:  -0.200
-    71 41:  -2.43e-17
     71 48:  0.200
-    71 49:  2.08e-17
     71 68:  1.000
     71 69:  -2.00
     71 70:  2.00
-    71 72:  5.55e-17
-    71 74:  1.11e-16
-    71 212:  -4.34e-19
-    71 437:  2.17e-18
-    71 617:  8.67e-19
-    71 797:  6.51e-19
-    75 8:  -1.39e-17
-    75 9:  -6.94e-18
-    75 10:  1.39e-17
-    75 24:  -1.73e-18
-    75 25:  2.78e-17
-    75 40:  1.39e-17
     75 41:  -0.200
-    75 48:  -1.73e-18
     75 49:  0.200
-    75 50:  -6.94e-18
-    75 68:  -5.55e-17
-    75 69:  -2.78e-17
-    75 70:  -1.67e-16
     75 72:  1.000
     75 73:  -2.00
     75 74:  2.00
-    75 76:  -5.55e-17
-    75 77:  1.11e-16
-    75 78:  1.11e-16
-    75 212:  8.67e-19
-    75 437:  3.47e-18
-    75 617:  1.30e-18
-    75 797:  -4.34e-19
-    79 8:  2.78e-17
-    79 9:  -4.16e-17
-    79 10:  1.39e-17
-    79 24:  1.39e-17
-    79 25:  6.94e-17
-    79 26:  1.39e-17
-    79 40:  6.94e-18
-    79 41:  -2.78e-17
     79 42:  -0.200
-    79 48:  1.39e-17
-    79 49:  4.16e-17
     79 50:  0.200
-    79 68:  -8.33e-17
-    79 69:  3.33e-16
-    79 72:  2.22e-16
-    79 73:  -6.66e-16
-    79 74:  2.22e-16
     79 76:  1.000
     79 77:  -2.00
     79 78:  2.00
-    79 212:  -3.47e-18
-    79 617:  4.34e-18
-    79 797:  1.73e-18
     80 8:  -0.200
-    80 9:  -2.08e-17
     80 24:  0.200
-    80 25:  1.39e-17
-    80 26:  -3.47e-18
-    80 40:  -1.91e-17
-    80 41:  2.78e-17
-    80 42:  -2.78e-17
-    80 48:  5.20e-18
     80 68:  1.000
-    80 69:  2.78e-17
     80 72:  -2.00
-    80 73:  -1.11e-16
     80 76:  2.00
-    80 77:  1.67e-16
-    80 212:  4.77e-18
-    80 437:  -4.34e-19
-    80 797:  1.52e-18
-    81 8:  3.47e-18
     81 9:  -0.200
-    81 10:  1.39e-17
-    81 24:  -6.94e-18
     81 25:  0.200
-    81 26:  -1.39e-17
-    81 40:  -1.39e-17
-    81 41:  -2.78e-17
-    81 42:  3.47e-17
-    81 48:  -6.94e-18
-    81 49:  2.78e-17
-    81 50:  -1.39e-17
-    81 68:  1.39e-17
     81 69:  1.000
-    81 70:  -1.11e-16
-    81 72:  -2.78e-17
     81 73:  -2.00
-    81 74:  1.11e-16
-    81 76:  -2.22e-16
     81 77:  2.00
-    81 78:  -1.11e-16
-    81 212:  -8.67e-19
-    81 437:  8.67e-19
-    81 617:  4.34e-19
-    81 797:  -1.30e-18
-    82 8:  2.78e-17
-    82 9:  -5.55e-17
     82 10:  -0.200
-    82 24:  1.39e-17
-    82 25:  5.55e-17
     82 26:  0.200
-    82 40:  6.94e-18
-    82 41:  -1.39e-17
-    82 42:  5.55e-17
-    82 48:  6.94e-18
-    82 49:  5.55e-17
-    82 50:  -1.39e-17
-    82 68:  -1.11e-16
-    82 69:  1.11e-16
     82 70:  1.000
-    82 72:  1.11e-16
-    82 73:  -4.44e-16
     82 74:  -2.00
-    82 76:  -3.33e-16
-    82 77:  6.66e-16
     82 78:  2.00
-    82 212:  -5.20e-18
-    82 437:  2.60e-18
-    82 617:  2.60e-18
-    82 797:  8.67e-19
     83 8:  -0.200
     83 9:  0.400
     83 10:  -0.400
index eb23c2225123b7fafcdda212ca931094ce73fa7d..cb06a03730ff98637d7568e497db973737eaba61 100644 (file)
@@ -1032,7 +1032,6 @@ DEAL::n_dofs=379
     41 58:  0.111
     41 87:  -0.222
     41 89:  -0.111
-    41 92:  -2.47e-17
     42 40:  1.00
     42 57:  0.0370
     42 58:  0.111
@@ -1552,12 +1551,7 @@ DEAL::n_dofs=379
     269 58:  -0.0139
     269 87:  0.278
     269 89:  0.139
-    269 92:  -2.47e-17
     270 40:  1.00
-    270 57:  -3.85e-18
-    270 87:  -1.77e-17
-    270 89:  -1.62e-17
-    270 92:  -2.47e-17
     279 57:  0.833
     279 58:  0.167
     280 57:  0.667
@@ -1697,7 +1691,6 @@ DEAL::n_dofs=379
     303 58:  0.111
     303 87:  -0.222
     303 89:  -0.111
-    303 92:  -2.47e-17
     304 40:  0.625
     304 57:  0.417
     304 58:  0.208
index 5006923bb19fd1833c20dab7722f1a1957e3c796..d8eaf7b07a3668d38768b059c2b1c4ebe9e0bcb4 100644 (file)
@@ -87,7 +87,6 @@ DEAL::n_dofs=46
     26 0:  0.417
     26 2:  0.625
     26 14:  -0.0417
-    27 0:  -2.78e-17
     27 2:  1.00
     28 40:  0.667
     28 41:  0.333
index 8aaa007067a007d759cf6b6b38197c25014ae09e..f1c449a939940e2e4e25da361daf4d8acc5cb547 100644 (file)
@@ -1032,7 +1032,6 @@ DEAL::n_dofs=1493
     41 183:  0.111
     41 303:  -0.222
     41 305:  -0.111
-    41 335:  -2.47e-17
     42 40:  1.00
     42 182:  0.0370
     42 183:  0.111
@@ -1552,12 +1551,7 @@ DEAL::n_dofs=1493
     944 183:  -0.0139
     944 303:  0.278
     944 305:  0.139
-    944 335:  -2.47e-17
     945 40:  1.00
-    945 182:  -3.85e-18
-    945 303:  -1.77e-17
-    945 305:  -1.62e-17
-    945 335:  -2.47e-17
     1079 182:  0.833
     1079 183:  0.167
     1080 182:  0.667
@@ -1697,7 +1691,6 @@ DEAL::n_dofs=1493
     1103 183:  0.111
     1103 303:  -0.222
     1103 305:  -0.111
-    1103 335:  -2.47e-17
     1104 40:  0.625
     1104 182:  0.417
     1104 183:  0.208
index 81cd13989996f79da1f86a6b1dcb38dd9ac09317..b5c830512630debeba0c4d0676c45447ecdb91db 100644 (file)
@@ -87,7 +87,6 @@ DEAL::n_dofs=146
     67 0:  0.417
     67 2:  0.625
     67 39:  -0.0417
-    68 0:  -2.78e-17
     68 2:  1.00
     69 122:  0.667
     69 123:  0.333
index 068da7d88e44ce557d70d78e3ca3f2dbf39ba21c..481412f20603833aa770b6e674abf67b5ea64c61 100644 (file)
@@ -534,15 +534,12 @@ DEAL::n_dofs=5800
     454 851:  -0.0988
     454 855:  0.198
     454 862:  0.790
-    456 99:  1.16e-18
     456 102:  0.111
     456 443:  -0.333
     456 455:  1.00
     456 615:  0.222
-    456 3579:  -2.48e-17
     457 99:  0.222
     457 102:  0.111
-    457 443:  -3.12e-17
     457 455:  1.00
     457 615:  -0.111
     457 3579:  -0.222
@@ -3683,7 +3680,6 @@ DEAL::n_dofs=5800
     3771 443:  0.625
     3771 3579:  0.417
     3772 443:  1.00
-    3772 3579:  -2.78e-17
     3773 99:  0.00694
     3773 102:  -0.0122
     3773 443:  0.0312
@@ -3710,12 +3706,10 @@ DEAL::n_dofs=5800
     3777 455:  0.703
     3777 615:  0.260
     3777 3579:  0.208
-    3778 99:  5.78e-18
     3778 102:  0.0417
     3778 443:  -0.250
     3778 455:  1.12
     3778 615:  0.0833
-    3778 3579:  -1.70e-17
     3779 3577:  0.667
     3779 3581:  0.333
     3780 3577:  0.333
@@ -3736,7 +3730,6 @@ DEAL::n_dofs=5800
     3785 3579:  -0.208
     3786 99:  0.0833
     3786 102:  0.0417
-    3786 443:  -1.73e-17
     3786 455:  1.12
     3786 615:  -0.0833
     3786 3579:  -0.167
@@ -3824,22 +3817,14 @@ DEAL::n_dofs=5800
     3799 3579:  0.347
     3800 99:  -0.0278
     3800 102:  -0.0139
-    3800 443:  -1.73e-17
     3800 455:  0.625
     3800 615:  0.139
     3800 3579:  0.278
-    3801 99:  -3.47e-18
     3801 102:  -0.0139
     3801 443:  0.417
     3801 455:  0.625
     3801 615:  -0.0278
-    3801 3579:  -3.92e-17
-    3802 99:  1.16e-18
-    3802 102:  3.90e-18
-    3802 443:  -3.12e-17
     3802 455:  1.00
-    3802 615:  -9.25e-18
-    3802 3579:  -2.48e-17
     3807 99:  0.00231
     3807 102:  0.000579
     3807 443:  0.0174
@@ -3851,12 +3836,12 @@ DEAL::n_dofs=5800
     3807 3978:  0.198
     3807 3979:  -0.0988
     3807 3989:  0.790
-    3808 99:  2.17e-19
+    3808 99:  6.51e-19
     3808 102:  0.00463
     3808 443:  0.0278
     3808 455:  -0.125
     3808 615:  0.00926
-    3808 3579:  1.73e-18
+    3808 3579:  8.67e-19
     3808 3581:  0.0741
     3808 3755:  0.148
     3808 3978:  -0.0988
@@ -3873,12 +3858,12 @@ DEAL::n_dofs=5800
     3809 3978:  0.198
     3809 3979:  -0.0988
     3809 3989:  0.790
-    3810 99:  -1.30e-18
+    3810 99:  -3.90e-18
     3810 102:  0.00926
     3810 443:  -0.0556
     3810 455:  0.250
     3810 615:  0.0185
-    3810 3579:  -3.47e-18
+    3810 3579:  0.00
     3810 3581:  -0.0370
     3810 3755:  -0.0741
     3810 3978:  -0.0988
index 053e2a57057a65daf78b39f4a236ab033498523d..463394325a3b0f1e870376914c53edc83c2218ad 100644 (file)
@@ -2555,17 +2555,11 @@ DEAL::n_dofs=796
     21 218:  0.222
     21 220:  0.444
     31 30:  1.00
-    31 151:  2.74e-18
     31 166:  0.333
     31 206:  -0.222
     31 208:  -0.111
-    31 220:  -2.08e-17
-    31 256:  1.73e-18
     32 30:  1.00
-    32 151:  1.01e-18
-    32 166:  1.73e-18
     32 206:  -0.222
-    32 208:  -2.08e-17
     32 220:  -0.111
     32 256:  0.333
     33 30:  1.00
@@ -2660,22 +2654,13 @@ DEAL::n_dofs=796
     179 218:  0.222
     179 219:  0.222
     179 220:  0.444
-    189 150:  -1.04e-17
-    189 151:  2.74e-18
     189 162:  0.333
-    189 164:  -3.12e-17
-    189 166:  1.73e-18
     189 188:  1.00
     189 219:  -0.222
     189 220:  -0.111
-    190 150:  3.47e-18
-    190 151:  1.01e-18
-    190 162:  1.73e-18
     190 164:  -0.333
     190 166:  0.333
     190 188:  1.00
-    190 219:  -2.48e-17
-    190 220:  -2.08e-17
     191 150:  -0.111
     191 151:  0.111
     191 162:  0.333
@@ -2722,19 +2707,10 @@ DEAL::n_dofs=796
     271 204:  0.222
     271 207:  0.222
     271 208:  0.444
-    273 150:  -1.04e-17
-    273 151:  2.74e-18
     273 162:  0.333
     273 207:  -0.222
     273 208:  -0.111
-    273 254:  -3.12e-17
-    273 256:  1.73e-18
     273 272:  1.00
-    274 150:  3.47e-18
-    274 151:  1.01e-18
-    274 162:  1.73e-18
-    274 207:  -2.48e-17
-    274 208:  -2.08e-17
     274 254:  -0.333
     274 256:  0.333
     274 272:  1.00
@@ -2965,7 +2941,7 @@ DEAL::n_dofs=796
     402 150:  0.125
     402 164:  1.12
     402 219:  -0.250
-    403 150:  -0.0312
+    403 150:  -0.0313
     403 151:  0.0156
     403 162:  0.141
     403 164:  -0.281
@@ -3132,13 +3108,9 @@ DEAL::n_dofs=796
     493 188:  0.703
     493 219:  0.208
     493 220:  0.208
-    494 151:  3.08e-18
-    494 162:  8.67e-18
     494 164:  -0.250
     494 166:  0.125
     494 188:  1.12
-    494 219:  -1.70e-17
-    494 220:  -9.25e-18
     495 30:  0.703
     495 151:  0.00521
     495 166:  0.0781
@@ -3147,17 +3119,13 @@ DEAL::n_dofs=796
     495 220:  0.208
     495 256:  -0.0469
     496 30:  1.13
-    496 151:  3.08e-18
     496 166:  0.125
     496 206:  -0.167
     496 208:  -0.0833
-    496 220:  -9.25e-18
-    496 256:  8.67e-18
     497 151:  -0.0417
     497 166:  0.625
     497 220:  0.417
     498 166:  1.00
-    498 220:  -2.78e-17
     499 30:  -0.0781
     499 151:  0.000579
     499 166:  0.00868
@@ -3185,7 +3153,7 @@ DEAL::n_dofs=796
     500 465:  -0.0988
     500 469:  0.790
     501 30:  -0.125
-    501 151:  -2.71e-19
+    501 151:  1.08e-19
     501 166:  -0.0139
     501 205:  -0.0185
     501 206:  -0.0247
@@ -3198,7 +3166,7 @@ DEAL::n_dofs=796
     501 465:  0.198
     501 469:  0.790
     502 30:  0.250
-    502 151:  4.34e-19
+    502 151:  -3.25e-19
     502 166:  0.0278
     502 205:  0.00926
     502 206:  -0.0710
@@ -3218,14 +3186,9 @@ DEAL::n_dofs=796
     503 188:  0.625
     503 219:  0.139
     503 220:  0.278
-    504 150:  3.47e-18
-    504 151:  1.01e-18
-    504 162:  1.73e-18
     504 164:  -0.333
     504 166:  0.333
     504 188:  1.00
-    504 219:  -2.48e-17
-    504 220:  -2.08e-17
     505 150:  0.0104
     505 151:  -0.0260
     505 162:  -0.0260
@@ -3234,14 +3197,9 @@ DEAL::n_dofs=796
     505 188:  0.391
     505 219:  0.0694
     505 220:  0.347
-    506 150:  -1.73e-18
-    506 151:  -8.89e-19
-    506 162:  4.34e-18
     506 164:  -0.250
     506 166:  0.625
     506 188:  0.625
-    506 219:  -9.17e-18
-    506 220:  -2.54e-17
     507 150:  0.00116
     507 151:  0.000579
     507 162:  0.00521
@@ -3256,8 +3214,8 @@ DEAL::n_dofs=796
     507 418:  0.790
     507 453:  -0.0247
     507 465:  0.198
-    508 150:  -2.17e-19
-    508 151:  -2.71e-19
+    508 150:  6.51e-19
+    508 151:  1.08e-19
     508 162:  8.67e-19
     508 164:  0.0278
     508 166:  -0.0139
@@ -3284,8 +3242,8 @@ DEAL::n_dofs=796
     509 418:  0.790
     509 453:  0.0123
     509 465:  -0.0988
-    510 150:  6.51e-19
-    510 151:  4.34e-19
+    510 150:  -1.08e-18
+    510 151:  -3.25e-19
     510 162:  -2.60e-18
     510 164:  -0.0556
     510 166:  0.0278
@@ -3301,7 +3259,7 @@ DEAL::n_dofs=796
     511 30:  0.625
     511 151:  -0.0139
     511 166:  0.208
-    511 206:  0.00
+    511 206:  -5.55e-17
     511 208:  -0.0556
     511 220:  0.278
     511 256:  -0.0417
@@ -3313,19 +3271,13 @@ DEAL::n_dofs=796
     512 220:  0.347
     512 256:  -0.0260
     513 30:  1.00
-    513 151:  2.74e-18
     513 166:  0.333
     513 206:  -0.222
     513 208:  -0.111
-    513 220:  -2.08e-17
-    513 256:  1.73e-18
     514 30:  0.625
-    514 151:  1.71e-18
     514 166:  0.625
     514 206:  -0.167
     514 208:  -0.0833
-    514 220:  -1.85e-17
-    514 256:  4.34e-18
     515 205:  0.111
     515 206:  0.222
     515 219:  0.222
@@ -3349,7 +3301,7 @@ DEAL::n_dofs=796
     559 254:  1.12
     560 207:  0.500
     560 208:  0.500
-    561 150:  -0.0312
+    561 150:  -0.0313
     561 151:  0.0156
     561 162:  0.141
     561 207:  -0.125
@@ -3417,11 +3369,7 @@ DEAL::n_dofs=796
     646 188:  0.703
     646 219:  -0.208
     646 220:  -0.0417
-    647 150:  -3.47e-18
-    647 151:  3.08e-18
     647 162:  0.125
-    647 164:  -1.73e-17
-    647 166:  8.67e-18
     647 188:  1.13
     647 219:  -0.167
     647 220:  -0.0833
@@ -3463,18 +3411,13 @@ DEAL::n_dofs=796
     652 254:  0.469
     652 256:  -0.0469
     652 272:  0.703
-    653 150:  -3.47e-18
-    653 151:  3.08e-18
     653 162:  0.125
     653 207:  -0.167
     653 208:  -0.0833
-    653 254:  -1.73e-17
-    653 256:  8.67e-18
     653 272:  1.12
     654 150:  0.417
     654 151:  -0.0417
     654 162:  0.625
-    655 150:  -2.78e-17
     655 162:  1.00
     656 150:  0.333
     656 207:  -0.333
@@ -3564,8 +3507,8 @@ DEAL::n_dofs=796
     668 569:  0.198
     668 571:  -0.0988
     668 573:  0.790
-    669 150:  4.34e-19
-    669 151:  -2.71e-19
+    669 150:  1.30e-18
+    669 151:  1.08e-19
     669 162:  -0.0139
     669 207:  0.0185
     669 208:  0.00926
@@ -3592,8 +3535,8 @@ DEAL::n_dofs=796
     670 569:  0.198
     670 571:  -0.0988
     670 573:  0.790
-    671 150:  -8.67e-19
-    671 151:  4.34e-19
+    671 150:  -1.73e-18
+    671 151:  -3.25e-19
     671 162:  0.0278
     671 207:  -0.0370
     671 208:  -0.0185
@@ -3622,19 +3565,11 @@ DEAL::n_dofs=796
     673 188:  0.391
     673 219:  -0.208
     673 220:  -0.0417
-    674 150:  -1.04e-17
-    674 151:  2.74e-18
     674 162:  0.333
-    674 164:  -3.12e-17
-    674 166:  1.73e-18
     674 188:  1.00
     674 219:  -0.222
     674 220:  -0.111
-    675 150:  -1.56e-17
-    675 151:  1.71e-18
     675 162:  0.625
-    675 164:  -8.67e-18
-    675 166:  4.34e-18
     675 188:  0.625
     675 219:  -0.167
     675 220:  -0.0833
@@ -3666,8 +3601,8 @@ DEAL::n_dofs=796
     685 411:  -0.0988
     685 420:  0.790
     685 453:  0.0123
-    686 150:  4.34e-19
-    686 151:  -2.71e-19
+    686 150:  1.30e-18
+    686 151:  1.08e-19
     686 162:  -0.0139
     686 164:  3.47e-18
     686 166:  8.67e-19
@@ -3680,8 +3615,8 @@ DEAL::n_dofs=796
     686 411:  0.198
     686 420:  0.790
     686 453:  0.0494
-    687 150:  -8.67e-19
-    687 151:  4.34e-19
+    687 150:  -1.73e-18
+    687 151:  -3.25e-19
     687 162:  0.0278
     687 164:  0.00
     687 166:  -2.60e-18
@@ -3710,21 +3645,13 @@ DEAL::n_dofs=796
     689 254:  0.260
     689 256:  -0.0260
     689 272:  0.391
-    690 150:  -1.04e-17
-    690 151:  2.74e-18
     690 162:  0.333
     690 207:  -0.222
     690 208:  -0.111
-    690 254:  -3.12e-17
-    690 256:  1.73e-18
     690 272:  1.00
-    691 150:  -1.56e-17
-    691 151:  1.71e-18
     691 162:  0.625
     691 207:  -0.167
     691 208:  -0.0833
-    691 254:  -8.67e-18
-    691 256:  4.34e-18
     691 272:  0.625
     701 151:  0.125
     701 208:  -0.250
@@ -3775,10 +3702,7 @@ DEAL::n_dofs=796
     708 220:  -0.0417
     708 256:  0.0781
     709 30:  1.13
-    709 151:  3.08e-18
-    709 166:  8.67e-18
     709 206:  -0.167
-    709 208:  -9.25e-18
     709 220:  -0.0833
     709 256:  0.125
     710 207:  0.333
@@ -3792,7 +3716,7 @@ DEAL::n_dofs=796
     712 208:  -0.167
     712 254:  -0.375
     712 256:  0.375
-    712 272:  1.12
+    712 272:  1.13
     713 150:  -0.0313
     713 151:  0.0781
     713 162:  0.0781
@@ -3809,17 +3733,12 @@ DEAL::n_dofs=796
     714 254:  -0.156
     714 256:  0.0781
     714 272:  0.703
-    715 151:  3.08e-18
-    715 162:  8.67e-18
-    715 207:  -1.70e-17
-    715 208:  -9.25e-18
     715 254:  -0.250
     715 256:  0.125
     715 272:  1.12
     716 151:  -0.0417
     716 208:  0.417
     716 256:  0.625
-    717 208:  -2.78e-17
     717 256:  1.00
     718 205:  0.167
     718 206:  0.167
@@ -3833,7 +3752,7 @@ DEAL::n_dofs=796
     720 208:  0.667
     721 206:  0.167
     721 208:  0.833
-    722 30:  1.12
+    722 30:  1.13
     722 151:  0.0417
     722 166:  0.125
     722 206:  -0.333
@@ -3891,8 +3810,8 @@ DEAL::n_dofs=796
     733 462:  -0.0988
     733 571:  -0.0988
     733 575:  0.790
-    734 150:  -2.17e-19
-    734 151:  -2.71e-19
+    734 150:  6.51e-19
+    734 151:  1.08e-19
     734 162:  8.67e-19
     734 205:  -0.0185
     734 206:  -0.0185
@@ -3905,8 +3824,8 @@ DEAL::n_dofs=796
     734 462:  0.198
     734 571:  0.198
     734 575:  0.790
-    735 150:  6.51e-19
-    735 151:  4.34e-19
+    735 150:  -1.08e-18
+    735 151:  -3.25e-19
     735 162:  -2.60e-18
     735 205:  0.00926
     735 206:  0.00926
@@ -3922,15 +3841,12 @@ DEAL::n_dofs=796
     736 30:  0.625
     736 151:  -0.0139
     736 166:  -0.0417
-    736 206:  0.00
+    736 206:  -5.55e-17
     736 208:  0.278
     736 220:  -0.0556
     736 256:  0.208
     737 30:  1.00
-    737 151:  1.01e-18
-    737 166:  1.73e-18
     737 206:  -0.222
-    737 208:  -2.08e-17
     737 220:  -0.111
     737 256:  0.333
     738 30:  0.391
@@ -3941,10 +3857,7 @@ DEAL::n_dofs=796
     738 220:  -0.0417
     738 256:  0.391
     739 30:  0.625
-    739 151:  -8.89e-19
-    739 166:  4.34e-18
     739 206:  -0.167
-    739 208:  -2.54e-17
     739 220:  -0.0833
     739 256:  0.625
     740 30:  -0.0781
@@ -3961,7 +3874,7 @@ DEAL::n_dofs=796
     740 462:  0.198
     740 473:  0.790
     741 30:  -0.125
-    741 151:  -2.71e-19
+    741 151:  1.08e-19
     741 166:  8.67e-19
     741 205:  -0.0185
     741 206:  -0.0247
@@ -3987,7 +3900,7 @@ DEAL::n_dofs=796
     742 462:  -0.0988
     742 473:  0.790
     743 30:  0.250
-    743 151:  4.34e-19
+    743 151:  -3.25e-19
     743 166:  -2.60e-18
     743 205:  0.00926
     743 206:  -0.0710
@@ -4007,11 +3920,6 @@ DEAL::n_dofs=796
     744 254:  -0.208
     744 256:  0.208
     744 272:  0.625
-    745 150:  3.47e-18
-    745 151:  1.01e-18
-    745 162:  1.73e-18
-    745 207:  -2.48e-17
-    745 208:  -2.08e-17
     745 254:  -0.333
     745 256:  0.333
     745 272:  1.00
@@ -4023,11 +3931,6 @@ DEAL::n_dofs=796
     746 254:  -0.156
     746 256:  0.391
     746 272:  0.391
-    747 150:  -1.73e-18
-    747 151:  -8.89e-19
-    747 162:  4.34e-18
-    747 207:  -9.17e-18
-    747 208:  -2.54e-17
     747 254:  -0.250
     747 256:  0.625
     747 272:  0.625
@@ -4039,7 +3942,7 @@ DEAL::n_dofs=796
     758 162:  0.125
     758 164:  -0.375
     758 166:  0.375
-    758 188:  1.12
+    758 188:  1.13
     758 219:  -0.0833
     758 220:  -0.167
     759 150:  -0.0313
@@ -4050,7 +3953,7 @@ DEAL::n_dofs=796
     759 188:  0.703
     759 219:  -0.0417
     759 220:  -0.208
-    760 30:  1.12
+    760 30:  1.13
     760 151:  0.0417
     760 166:  0.375
     760 206:  -0.333
@@ -5885,15 +5788,12 @@ DEAL::n_dofs=645
     165 290:  0.790
     167 67:  0.222
     167 69:  0.111
-    167 144:  -3.12e-17
     167 166:  1.00
     167 244:  -0.222
     167 245:  -0.111
-    168 67:  4.63e-18
     168 69:  0.111
     168 144:  -0.333
     168 166:  1.00
-    168 244:  -2.48e-17
     168 245:  0.222
     169 67:  0.111
     169 69:  0.333
@@ -6265,7 +6165,6 @@ DEAL::n_dofs=645
     536 245:  -0.0729
     537 67:  0.0833
     537 69:  0.0417
-    537 144:  -1.73e-17
     537 166:  1.13
     537 244:  -0.167
     537 245:  -0.0833
@@ -6336,7 +6235,7 @@ DEAL::n_dofs=645
     554 60:  -0.0216
     554 64:  -0.0216
     554 67:  0.00926
-    554 69:  -3.47e-17
+    554 69:  -4.16e-17
     554 144:  -0.0556
     554 166:  0.250
     554 244:  -0.0370
@@ -6429,13 +6328,11 @@ DEAL::n_dofs=645
     565 245:  -0.0590
     566 67:  0.222
     566 69:  0.111
-    566 144:  -3.12e-17
     566 166:  1.00
     566 244:  -0.222
     566 245:  -0.111
     567 67:  0.417
     567 69:  0.208
-    567 144:  -8.67e-18
     567 166:  0.625
     567 244:  -0.167
     567 245:  -0.0833
index a5870ef8662c3460e425fe2eb2ed998257c8b040..4824575fc71cebf4aa0e807447c127e0edc198c6 100644 (file)
@@ -13,7 +13,6 @@ DEAL::n_dofs=119
     47 46:  0.417
     47 89:  -0.0417
     48 29:  1.00
-    48 46:  -2.78e-17
     49 63:  0.222
     49 66:  0.889
     49 88:  -0.111

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.