StokesCosine tested

git-svn-id: https://svn.dealii.org/trunk@14943 0785d39b-7218-0410-832d-ea1e28bc413d

diff --git a/tests/base/functions_01.cc b/tests/base/functions_01.cc
index 4cf49787d0f8dd5fd5009252031fbff9230bd69d..e60e68373459d96fd7fc08c07e9ec9d78966e9e9 100644 (file)
--- a/tests/base/functions_01.cc
+++ b/tests/base/functions_01.cc
@@ -143,8 +143,8 @@ check_function(const Functions::FlowFunction<dim>& f,
                                             // bit generous
            if (d2.norm() < 12.* d1.norm())
              {
-               deallog << "Gradient error " << k << " " << i
-                       << " " << d1.norm() << " " << d2.norm()
+               deallog << "Gradient error: point " << k << " comp " << i
+                       << " norms " << d1.norm() << " " << d2.norm()
                        << std::endl;
 //             for (unsigned int i=0;i<f.n_components;++i)
                  for (unsigned int d=0;d<dim;++d)
@@ -209,7 +209,14 @@ int main()
     {
       deallog << " Functions::StokesCosine<2>" << std::endl;
       Functions::StokesCosine<2> f(1.);
-      check_function(f, 10, logfile);
+      check_function(f, 4, logfile);
+    }
+
+  if (true)
+    {
+      deallog << " Functions::StokesCosine<3>" << std::endl;
+      Functions::StokesCosine<3> f(1.);
+      check_function(f, 4, logfile);
     }
 
   if (false)

diff --git a/tests/base/functions_01/cmp/generic b/tests/base/functions_01/cmp/generic
index 3c61dd99d661d7a872c350fee15b47a5cff610ea..5cd8b6eca2d57749e3d45975cfe1fdbf0517f364 100644 (file)
--- a/tests/base/functions_01/cmp/generic
+++ b/tests/base/functions_01/cmp/generic
@@ -1,154 +1,370 @@
-object "vertices" class array type float rank 1 shape 2 items 25 data follows
--1     -1      
--0.5   -1      
-0      -1      
-0.5    -1      
-1      -1      
--1     -0.5    
--0.5   -0.5    
-0      -0.5    
-0.5    -0.5    
-1      -0.5    
--1     0       
--0.5   0       
-0      0       
-0.5    0       
-1      0       
--1     0.5     
--0.5   0.5     
-0      0.5     
-0.5    0.5     
-1      0.5     
--1     1       
--0.5   1       
-0      1       
-0.5    1       
-1      1       
-object "cells" class array type int rank 1 shape 4 items 16 data follows
-0      5       1       6
-1      6       2       7
-2      7       3       8
-3      8       4       9
-5      10      6       11
-6      11      7       12
-7      12      8       13
-8      13      9       14
-10     15      11      16
-11     16      12      17
-12     17      13      18
-13     18      14      19
-15     20      16      21
-16     21      17      22
-17     22      18      23
-18     23      19      24
 
-attribute "element type" string "quads"
-attribute "ref" string "positions"
-object "data" class array type float rank 1 shape 3 items 25 data follows
-       4.26453 -4.26453        6.95821e-14
-       4.04692 -3.18452        -0.633437
-       3.64683 -1.9805 -1.59808
-       3.03204 -1.05036        -2.33684
-       2.47239 -0.566216       -2.55657
-       3.18452 -4.04692        0.633437
-       2.92392 -2.92392        9.54161e-14
-       2.5004  -1.3579 -2.1914
-       1.69516 -0.388218       -3.50576
-       1.2366  -0.154216       -3.29301
-       1.9805  -3.64683        1.59808
-       1.3579  -2.5004 2.1914
-       0       0       0
-       0       0       0
-       0       0       0
-       1.05036 -3.03204        2.33684
-       0.388218        -1.69516        3.50576
-       0       0       0
-       0       0       0
-       0       0       0
-       0.566216        -2.47239        2.55657
-       0.154216        -1.2366 3.29301
-       0       0       0
-       0       0       0
-       0       0       0
-attribute "dep" string "positions"
-object "deal data" class field
-component "positions" value "vertices"
-component "connections" value "cells"
-component "data" value "data"
-attribute 
-end
-object "vertices" class array type float rank 1 shape 2 items 25 data follows
--1     -1      
--0.5   -1      
-0      -1      
-0.5    -1      
-1      -1      
--1     -0.5    
--0.5   -0.5    
-0      -0.5    
-0.5    -0.5    
-1      -0.5    
--1     0       
--0.5   0       
-0      0       
-0.5    0       
-1      0       
--1     0.5     
--0.5   0.5     
-0      0.5     
-0.5    0.5     
-1      0.5     
--1     1       
--0.5   1       
-0      1       
-0.5    1       
-1      1       
-object "cells" class array type int rank 1 shape 4 items 16 data follows
-0      5       1       6
-1      6       2       7
-2      7       3       8
-3      8       4       9
-5      10      6       11
-6      11      7       12
-7      12      8       13
-8      13      9       14
-10     15      11      16
-11     16      12      17
-12     17      13      18
-13     18      14      19
-15     20      16      21
-16     21      17      22
-17     22      18      23
-18     23      19      24
+DEAL:: Functions::StokesCosine<2>
+DEAL::Gradients tested
+DEAL::Divergence tested
+DEAL::Laplacians 0.719980
+# This file was generated by the deal.II library.
+
+
+#
+# For a description of the GNUPLOT format see the GNUPLOT manual.
+#
+# <x> <y> <comp> <comp> <comp> 
+-1.00000 -1.00000 0.00000 0.00000 3.74940e-33 
+-0.500000 -1.00000 -3.06162e-17 1.87470e-33 3.06162e-17 
+0.00000 -1.00000 -6.12323e-17 0.00000 0.00000 
+0.500000 -1.00000 -3.06162e-17 -1.87470e-33 -3.06162e-17 
+1.00000 -1.00000 0.00000 0.00000 -3.74940e-33 
+
+-1.00000 -0.500000 -1.87470e-33 3.06162e-17 3.06162e-17 
+-0.500000 -0.500000 -0.250000 0.250000 0.250000 
+0.00000 -0.500000 -0.500000 0.00000 0.00000 
+0.500000 -0.500000 -0.250000 -0.250000 -0.250000 
+1.00000 -0.500000 -1.87470e-33 -3.06162e-17 -3.06162e-17 
+
+-1.00000 0.00000 0.00000 6.12323e-17 0.00000 
+-0.500000 0.00000 0.00000 0.500000 0.00000 
+0.00000 0.00000 0.00000 0.00000 0.00000 
+0.500000 0.00000 0.00000 -0.500000 0.00000 
+1.00000 0.00000 0.00000 -6.12323e-17 0.00000 
+
+-1.00000 0.500000 1.87470e-33 3.06162e-17 -3.06162e-17 
+-0.500000 0.500000 0.250000 0.250000 -0.250000 
+0.00000 0.500000 0.500000 0.00000 0.00000 
+0.500000 0.500000 0.250000 -0.250000 0.250000 
+1.00000 0.500000 1.87470e-33 -3.06162e-17 3.06162e-17 
+
+-1.00000 1.00000 0.00000 0.00000 -3.74940e-33 
+-0.500000 1.00000 3.06162e-17 1.87470e-33 -3.06162e-17 
+0.00000 1.00000 6.12323e-17 0.00000 0.00000 
+0.500000 1.00000 3.06162e-17 -1.87470e-33 3.06162e-17 
+1.00000 1.00000 0.00000 0.00000 3.74940e-33 
+
+
+DEAL:: Functions::StokesCosine<3>
+DEAL::Gradients tested
+DEAL::Divergence tested
+DEAL::Laplacians 0.119270
+object "vertices" class array type float rank 1 shape 3 items 125 data follows
+-1.00000       -1.00000        -1.00000        
+-0.500000      -1.00000        -1.00000        
+0.00000        -1.00000        -1.00000        
+0.500000       -1.00000        -1.00000        
+1.00000        -1.00000        -1.00000        
+-1.00000       -0.500000       -1.00000        
+-0.500000      -0.500000       -1.00000        
+0.00000        -0.500000       -1.00000        
+0.500000       -0.500000       -1.00000        
+1.00000        -0.500000       -1.00000        
+-1.00000       0.00000 -1.00000        
+-0.500000      0.00000 -1.00000        
+0.00000        0.00000 -1.00000        
+0.500000       0.00000 -1.00000        
+1.00000        0.00000 -1.00000        
+-1.00000       0.500000        -1.00000        
+-0.500000      0.500000        -1.00000        
+0.00000        0.500000        -1.00000        
+0.500000       0.500000        -1.00000        
+1.00000        0.500000        -1.00000        
+-1.00000       1.00000 -1.00000        
+-0.500000      1.00000 -1.00000        
+0.00000        1.00000 -1.00000        
+0.500000       1.00000 -1.00000        
+1.00000        1.00000 -1.00000        
+-1.00000       -1.00000        -0.500000       
+-0.500000      -1.00000        -0.500000       
+0.00000        -1.00000        -0.500000       
+0.500000       -1.00000        -0.500000       
+1.00000        -1.00000        -0.500000       
+-1.00000       -0.500000       -0.500000       
+-0.500000      -0.500000       -0.500000       
+0.00000        -0.500000       -0.500000       
+0.500000       -0.500000       -0.500000       
+1.00000        -0.500000       -0.500000       
+-1.00000       0.00000 -0.500000       
+-0.500000      0.00000 -0.500000       
+0.00000        0.00000 -0.500000       
+0.500000       0.00000 -0.500000       
+1.00000        0.00000 -0.500000       
+-1.00000       0.500000        -0.500000       
+-0.500000      0.500000        -0.500000       
+0.00000        0.500000        -0.500000       
+0.500000       0.500000        -0.500000       
+1.00000        0.500000        -0.500000       
+-1.00000       1.00000 -0.500000       
+-0.500000      1.00000 -0.500000       
+0.00000        1.00000 -0.500000       
+0.500000       1.00000 -0.500000       
+1.00000        1.00000 -0.500000       
+-1.00000       -1.00000        0.00000 
+-0.500000      -1.00000        0.00000 
+0.00000        -1.00000        0.00000 
+0.500000       -1.00000        0.00000 
+1.00000        -1.00000        0.00000 
+-1.00000       -0.500000       0.00000 
+-0.500000      -0.500000       0.00000 
+0.00000        -0.500000       0.00000 
+0.500000       -0.500000       0.00000 
+1.00000        -0.500000       0.00000 
+-1.00000       0.00000 0.00000 
+-0.500000      0.00000 0.00000 
+0.00000        0.00000 0.00000 
+0.500000       0.00000 0.00000 
+1.00000        0.00000 0.00000 
+-1.00000       0.500000        0.00000 
+-0.500000      0.500000        0.00000 
+0.00000        0.500000        0.00000 
+0.500000       0.500000        0.00000 
+1.00000        0.500000        0.00000 
+-1.00000       1.00000 0.00000 
+-0.500000      1.00000 0.00000 
+0.00000        1.00000 0.00000 
+0.500000       1.00000 0.00000 
+1.00000        1.00000 0.00000 
+-1.00000       -1.00000        0.500000        
+-0.500000      -1.00000        0.500000        
+0.00000        -1.00000        0.500000        
+0.500000       -1.00000        0.500000        
+1.00000        -1.00000        0.500000        
+-1.00000       -0.500000       0.500000        
+-0.500000      -0.500000       0.500000        
+0.00000        -0.500000       0.500000        
+0.500000       -0.500000       0.500000        
+1.00000        -0.500000       0.500000        
+-1.00000       0.00000 0.500000        
+-0.500000      0.00000 0.500000        
+0.00000        0.00000 0.500000        
+0.500000       0.00000 0.500000        
+1.00000        0.00000 0.500000        
+-1.00000       0.500000        0.500000        
+-0.500000      0.500000        0.500000        
+0.00000        0.500000        0.500000        
+0.500000       0.500000        0.500000        
+1.00000        0.500000        0.500000        
+-1.00000       1.00000 0.500000        
+-0.500000      1.00000 0.500000        
+0.00000        1.00000 0.500000        
+0.500000       1.00000 0.500000        
+1.00000        1.00000 0.500000        
+-1.00000       -1.00000        1.00000 
+-0.500000      -1.00000        1.00000 
+0.00000        -1.00000        1.00000 
+0.500000       -1.00000        1.00000 
+1.00000        -1.00000        1.00000 
+-1.00000       -0.500000       1.00000 
+-0.500000      -0.500000       1.00000 
+0.00000        -0.500000       1.00000 
+0.500000       -0.500000       1.00000 
+1.00000        -0.500000       1.00000 
+-1.00000       0.00000 1.00000 
+-0.500000      0.00000 1.00000 
+0.00000        0.00000 1.00000 
+0.500000       0.00000 1.00000 
+1.00000        0.00000 1.00000 
+-1.00000       0.500000        1.00000 
+-0.500000      0.500000        1.00000 
+0.00000        0.500000        1.00000 
+0.500000       0.500000        1.00000 
+1.00000        0.500000        1.00000 
+-1.00000       1.00000 1.00000 
+-0.500000      1.00000 1.00000 
+0.00000        1.00000 1.00000 
+0.500000       1.00000 1.00000 
+1.00000        1.00000 1.00000 
+object "cells" class array type int rank 1 shape 8 items 64 data follows
+0      25      5       30      1       26      6       31
+1      26      6       31      2       27      7       32
+2      27      7       32      3       28      8       33
+3      28      8       33      4       29      9       34
+5      30      10      35      6       31      11      36
+6      31      11      36      7       32      12      37
+7      32      12      37      8       33      13      38
+8      33      13      38      9       34      14      39
+10     35      15      40      11      36      16      41
+11     36      16      41      12      37      17      42
+12     37      17      42      13      38      18      43
+13     38      18      43      14      39      19      44
+15     40      20      45      16      41      21      46
+16     41      21      46      17      42      22      47
+17     42      22      47      18      43      23      48
+18     43      23      48      19      44      24      49
+25     50      30      55      26      51      31      56
+26     51      31      56      27      52      32      57
+27     52      32      57      28      53      33      58
+28     53      33      58      29      54      34      59
+30     55      35      60      31      56      36      61
+31     56      36      61      32      57      37      62
+32     57      37      62      33      58      38      63
+33     58      38      63      34      59      39      64
+35     60      40      65      36      61      41      66
+36     61      41      66      37      62      42      67
+37     62      42      67      38      63      43      68
+38     63      43      68      39      64      44      69
+40     65      45      70      41      66      46      71
+41     66      46      71      42      67      47      72
+42     67      47      72      43      68      48      73
+43     68      48      73      44      69      49      74
+50     75      55      80      51      76      56      81
+51     76      56      81      52      77      57      82
+52     77      57      82      53      78      58      83
+53     78      58      83      54      79      59      84
+55     80      60      85      56      81      61      86
+56     81      61      86      57      82      62      87
+57     82      62      87      58      83      63      88
+58     83      63      88      59      84      64      89
+60     85      65      90      61      86      66      91
+61     86      66      91      62      87      67      92
+62     87      67      92      63      88      68      93
+63     88      68      93      64      89      69      94
+65     90      70      95      66      91      71      96
+66     91      71      96      67      92      72      97
+67     92      72      97      68      93      73      98
+68     93      73      98      69      94      74      99
+75     100     80      105     76      101     81      106
+76     101     81      106     77      102     82      107
+77     102     82      107     78      103     83      108
+78     103     83      108     79      104     84      109
+80     105     85      110     81      106     86      111
+81     106     86      111     82      107     87      112
+82     107     87      112     83      108     88      113
+83     108     88      113     84      109     89      114
+85     110     90      115     86      111     91      116
+86     111     91      116     87      112     92      117
+87     112     92      117     88      113     93      118
+88     113     93      118     89      114     94      119
+90     115     95      120     91      116     96      121
+91     116     96      121     92      117     97      122
+92     117     97      122     93      118     98      123
+93     118     98      123     94      119     99      124
 
-attribute "element type" string "quads"
+attribute "element type" string "cubes"
 attribute "ref" string "positions"
-object "data" class array type float rank 1 shape 3 items 25 data follows
-       -0.5625 0       0.3125
-       -0.5625 0       0.15625
-       -0.5625 0       0
-       -0.5625 0       -0.15625
-       -0.5625 0       -0.3125
-       0.609375        0       0.3125
-       0.609375        0       0.15625
-       0.609375        0       0
-       0.609375        0       -0.15625
-       0.609375        0       -0.3125
-       1       0       0.3125
-       1       0       0.15625
-       1       0       0
-       1       0       -0.15625
-       1       0       -0.3125
-       0.609375        0       0.3125
-       0.609375        0       0.15625
-       0.609375        0       0
-       0.609375        0       -0.15625
-       0.609375        0       -0.3125
-       -0.5625 0       0.3125
-       -0.5625 0       0.15625
-       -0.5625 0       0
-       -0.5625 0       -0.15625
-       -0.5625 0       -0.3125
+object "data" class array type float rank 1 shape 4 items 125 data follows
+       0.00000 0.00000 0.00000 0.00000
+       1.87470e-33     0.00000 0.00000 -1.87470e-33
+       3.74940e-33     0.00000 0.00000 0.00000
+       1.87470e-33     0.00000 0.00000 1.87470e-33
+       0.00000 0.00000 0.00000 0.00000
+       0.00000 1.87470e-33     0.00000 -1.87470e-33
+       1.53081e-17     1.53081e-17     -1.87470e-33    -1.53081e-17
+       3.06162e-17     0.00000 0.00000 0.00000
+       1.53081e-17     -1.53081e-17    1.87470e-33     1.53081e-17
+       0.00000 -1.87470e-33    0.00000 1.87470e-33
+       0.00000 3.74940e-33     0.00000 0.00000
+       0.00000 3.06162e-17     0.00000 0.00000
+       0.00000 0.00000 0.00000 0.00000
+       0.00000 -3.06162e-17    0.00000 0.00000
+       0.00000 -3.74940e-33    0.00000 0.00000
+       0.00000 1.87470e-33     0.00000 1.87470e-33
+       -1.53081e-17    1.53081e-17     1.87470e-33     1.53081e-17
+       -3.06162e-17    0.00000 0.00000 0.00000
+       -1.53081e-17    -1.53081e-17    -1.87470e-33    -1.53081e-17
+       0.00000 -1.87470e-33    0.00000 -1.87470e-33
+       0.00000 0.00000 0.00000 0.00000
+       -1.87470e-33    0.00000 0.00000 1.87470e-33
+       -3.74940e-33    0.00000 0.00000 0.00000
+       -1.87470e-33    0.00000 0.00000 -1.87470e-33
+       0.00000 0.00000 0.00000 0.00000
+       0.00000 0.00000 -3.74940e-33    -1.87470e-33
+       1.53081e-17     9.37350e-34     -3.06162e-17    -1.53081e-17
+       3.06162e-17     0.00000 0.00000 0.00000
+       1.53081e-17     -9.37350e-34    3.06162e-17     1.53081e-17
+       0.00000 0.00000 3.74940e-33     1.87470e-33
+       9.37350e-34     1.53081e-17     -3.06162e-17    -1.53081e-17
+       0.125000        0.125000        -0.250000       -0.125000
+       0.250000        0.00000 0.00000 0.00000
+       0.125000        -0.125000       0.250000        0.125000
+       9.37350e-34     -1.53081e-17    3.06162e-17     1.53081e-17
+       0.00000 3.06162e-17     0.00000 0.00000
+       0.00000 0.250000        0.00000 0.00000
+       0.00000 0.00000 0.00000 0.00000
+       0.00000 -0.250000       0.00000 0.00000
+       0.00000 -3.06162e-17    0.00000 0.00000
+       -9.37350e-34    1.53081e-17     3.06162e-17     1.53081e-17
+       -0.125000       0.125000        0.250000        0.125000
+       -0.250000       0.00000 0.00000 0.00000
+       -0.125000       -0.125000       -0.250000       -0.125000
+       -9.37350e-34    -1.53081e-17    -3.06162e-17    -1.53081e-17
+       0.00000 0.00000 3.74940e-33     1.87470e-33
+       -1.53081e-17    9.37350e-34     3.06162e-17     1.53081e-17
+       -3.06162e-17    0.00000 0.00000 0.00000
+       -1.53081e-17    -9.37350e-34    -3.06162e-17    -1.53081e-17
+       0.00000 0.00000 -3.74940e-33    -1.87470e-33
+       0.00000 0.00000 -7.49880e-33    0.00000
+       0.00000 0.00000 -6.12323e-17    0.00000
+       0.00000 0.00000 0.00000 0.00000
+       0.00000 0.00000 6.12323e-17     0.00000
+       0.00000 0.00000 7.49880e-33     0.00000
+       0.00000 0.00000 -6.12323e-17    0.00000
+       0.00000 0.00000 -0.500000       0.00000
+       0.00000 0.00000 0.00000 0.00000
+       0.00000 0.00000 0.500000        0.00000
+       0.00000 0.00000 6.12323e-17     0.00000
+       0.00000 0.00000 0.00000 0.00000
+       0.00000 0.00000 0.00000 0.00000
+       0.00000 0.00000 0.00000 0.00000
+       0.00000 0.00000 0.00000 0.00000
+       0.00000 0.00000 0.00000 0.00000
+       0.00000 0.00000 6.12323e-17     0.00000
+       0.00000 0.00000 0.500000        0.00000
+       0.00000 0.00000 0.00000 0.00000
+       0.00000 0.00000 -0.500000       0.00000
+       0.00000 0.00000 -6.12323e-17    0.00000
+       0.00000 0.00000 7.49880e-33     0.00000
+       0.00000 0.00000 6.12323e-17     0.00000
+       0.00000 0.00000 0.00000 0.00000
+       0.00000 0.00000 -6.12323e-17    0.00000
+       0.00000 0.00000 -7.49880e-33    0.00000
+       0.00000 0.00000 -3.74940e-33    1.87470e-33
+       -1.53081e-17    -9.37350e-34    -3.06162e-17    1.53081e-17
+       -3.06162e-17    0.00000 0.00000 0.00000
+       -1.53081e-17    9.37350e-34     3.06162e-17     -1.53081e-17
+       0.00000 0.00000 3.74940e-33     -1.87470e-33
+       -9.37350e-34    -1.53081e-17    -3.06162e-17    1.53081e-17
+       -0.125000       -0.125000       -0.250000       0.125000
+       -0.250000       0.00000 0.00000 0.00000
+       -0.125000       0.125000        0.250000        -0.125000
+       -9.37350e-34    1.53081e-17     3.06162e-17     -1.53081e-17
+       0.00000 -3.06162e-17    0.00000 0.00000
+       0.00000 -0.250000       0.00000 0.00000
+       0.00000 0.00000 0.00000 0.00000
+       0.00000 0.250000        0.00000 0.00000
+       0.00000 3.06162e-17     0.00000 0.00000
+       9.37350e-34     -1.53081e-17    3.06162e-17     -1.53081e-17
+       0.125000        -0.125000       0.250000        -0.125000
+       0.250000        0.00000 0.00000 0.00000
+       0.125000        0.125000        -0.250000       0.125000
+       9.37350e-34     1.53081e-17     -3.06162e-17    1.53081e-17
+       0.00000 0.00000 3.74940e-33     -1.87470e-33
+       1.53081e-17     -9.37350e-34    3.06162e-17     -1.53081e-17
+       3.06162e-17     0.00000 0.00000 0.00000
+       1.53081e-17     9.37350e-34     -3.06162e-17    1.53081e-17
+       0.00000 0.00000 -3.74940e-33    1.87470e-33
+       0.00000 0.00000 0.00000 0.00000
+       -1.87470e-33    0.00000 0.00000 1.87470e-33
+       -3.74940e-33    0.00000 0.00000 0.00000
+       -1.87470e-33    0.00000 0.00000 -1.87470e-33
+       0.00000 0.00000 0.00000 0.00000
+       0.00000 -1.87470e-33    0.00000 1.87470e-33
+       -1.53081e-17    -1.53081e-17    -1.87470e-33    1.53081e-17
+       -3.06162e-17    0.00000 0.00000 0.00000
+       -1.53081e-17    1.53081e-17     1.87470e-33     -1.53081e-17
+       0.00000 1.87470e-33     0.00000 -1.87470e-33
+       0.00000 -3.74940e-33    0.00000 0.00000
+       0.00000 -3.06162e-17    0.00000 0.00000
+       0.00000 0.00000 0.00000 0.00000
+       0.00000 3.06162e-17     0.00000 0.00000
+       0.00000 3.74940e-33     0.00000 0.00000
+       0.00000 -1.87470e-33    0.00000 -1.87470e-33
+       1.53081e-17     -1.53081e-17    1.87470e-33     -1.53081e-17
+       3.06162e-17     0.00000 0.00000 0.00000
+       1.53081e-17     1.53081e-17     -1.87470e-33    1.53081e-17
+       0.00000 1.87470e-33     0.00000 1.87470e-33
+       0.00000 0.00000 0.00000 0.00000
+       1.87470e-33     0.00000 0.00000 -1.87470e-33
+       3.74940e-33     0.00000 0.00000 0.00000
+       1.87470e-33     0.00000 0.00000 1.87470e-33
+       0.00000 0.00000 0.00000 0.00000
 attribute "dep" string "positions"
 object "deal data" class field
 component "positions" value "vertices"
@@ -156,132 +372,178 @@ component "connections" value "cells"
 component "data" value "data"
 attribute 
 end
+DEAL::Functions::PoisseuilleFlow<2>
+DEAL::Gradients tested
+DEAL::Divergence tested
+DEAL::Laplacians 0
+# This file was generated by the deal.II library.
+
+
+#
+# For a description of the GNUPLOT format see the GNUPLOT manual.
+#
+# <x> <y> <comp> <comp> <comp> 
+-1.00000 -1.00000 -0.562500 0.00000 0.312500 
+-0.500000 -1.00000 -0.562500 0.00000 0.156250 
+0.00000 -1.00000 -0.562500 0.00000 0.00000 
+0.500000 -1.00000 -0.562500 0.00000 -0.156250 
+1.00000 -1.00000 -0.562500 0.00000 -0.312500 
+
+-1.00000 -0.500000 0.609375 0.00000 0.312500 
+-0.500000 -0.500000 0.609375 0.00000 0.156250 
+0.00000 -0.500000 0.609375 0.00000 0.00000 
+0.500000 -0.500000 0.609375 0.00000 -0.156250 
+1.00000 -0.500000 0.609375 0.00000 -0.312500 
+
+-1.00000 0.00000 1.00000 0.00000 0.312500 
+-0.500000 0.00000 1.00000 0.00000 0.156250 
+0.00000 0.00000 1.00000 0.00000 0.00000 
+0.500000 0.00000 1.00000 0.00000 -0.156250 
+1.00000 0.00000 1.00000 0.00000 -0.312500 
+
+-1.00000 0.500000 0.609375 0.00000 0.312500 
+-0.500000 0.500000 0.609375 0.00000 0.156250 
+0.00000 0.500000 0.609375 0.00000 0.00000 
+0.500000 0.500000 0.609375 0.00000 -0.156250 
+1.00000 0.500000 0.609375 0.00000 -0.312500 
+
+-1.00000 1.00000 -0.562500 0.00000 0.312500 
+-0.500000 1.00000 -0.562500 0.00000 0.156250 
+0.00000 1.00000 -0.562500 0.00000 0.00000 
+0.500000 1.00000 -0.562500 0.00000 -0.156250 
+1.00000 1.00000 -0.562500 0.00000 -0.312500 
+
+
+DEAL::Functions::PoisseuilleFlow<3>
+DEAL::Gradients tested
+DEAL::Divergence tested
+DEAL::Laplacians 0
 object "vertices" class array type float rank 1 shape 3 items 125 data follows
--1     -1      -1      
--0.5   -1      -1      
-0      -1      -1      
-0.5    -1      -1      
-1      -1      -1      
--1     -0.5    -1      
--0.5   -0.5    -1      
-0      -0.5    -1      
-0.5    -0.5    -1      
-1      -0.5    -1      
--1     0       -1      
--0.5   0       -1      
-0      0       -1      
-0.5    0       -1      
-1      0       -1      
--1     0.5     -1      
--0.5   0.5     -1      
-0      0.5     -1      
-0.5    0.5     -1      
-1      0.5     -1      
--1     1       -1      
--0.5   1       -1      
-0      1       -1      
-0.5    1       -1      
-1      1       -1      
--1     -1      -0.5    
--0.5   -1      -0.5    
-0      -1      -0.5    
-0.5    -1      -0.5    
-1      -1      -0.5    
--1     -0.5    -0.5    
--0.5   -0.5    -0.5    
-0      -0.5    -0.5    
-0.5    -0.5    -0.5    
-1      -0.5    -0.5    
--1     0       -0.5    
--0.5   0       -0.5    
-0      0       -0.5    
-0.5    0       -0.5    
-1      0       -0.5    
--1     0.5     -0.5    
--0.5   0.5     -0.5    
-0      0.5     -0.5    
-0.5    0.5     -0.5    
-1      0.5     -0.5    
--1     1       -0.5    
--0.5   1       -0.5    
-0      1       -0.5    
-0.5    1       -0.5    
-1      1       -0.5    
--1     -1      0       
--0.5   -1      0       
-0      -1      0       
-0.5    -1      0       
-1      -1      0       
--1     -0.5    0       
--0.5   -0.5    0       
-0      -0.5    0       
-0.5    -0.5    0       
-1      -0.5    0       
--1     0       0       
--0.5   0       0       
-0      0       0       
-0.5    0       0       
-1      0       0       
--1     0.5     0       
--0.5   0.5     0       
-0      0.5     0       
-0.5    0.5     0       
-1      0.5     0       
--1     1       0       
--0.5   1       0       
-0      1       0       
-0.5    1       0       
-1      1       0       
--1     -1      0.5     
--0.5   -1      0.5     
-0      -1      0.5     
-0.5    -1      0.5     
-1      -1      0.5     
--1     -0.5    0.5     
--0.5   -0.5    0.5     
-0      -0.5    0.5     
-0.5    -0.5    0.5     
-1      -0.5    0.5     
--1     0       0.5     
--0.5   0       0.5     
-0      0       0.5     
-0.5    0       0.5     
-1      0       0.5     
--1     0.5     0.5     
--0.5   0.5     0.5     
-0      0.5     0.5     
-0.5    0.5     0.5     
-1      0.5     0.5     
--1     1       0.5     
--0.5   1       0.5     
-0      1       0.5     
-0.5    1       0.5     
-1      1       0.5     
--1     -1      1       
--0.5   -1      1       
-0      -1      1       
-0.5    -1      1       
-1      -1      1       
--1     -0.5    1       
--0.5   -0.5    1       
-0      -0.5    1       
-0.5    -0.5    1       
-1      -0.5    1       
--1     0       1       
--0.5   0       1       
-0      0       1       
-0.5    0       1       
-1      0       1       
--1     0.5     1       
--0.5   0.5     1       
-0      0.5     1       
-0.5    0.5     1       
-1      0.5     1       
--1     1       1       
--0.5   1       1       
-0      1       1       
-0.5    1       1       
-1      1       1       
+-1.00000       -1.00000        -1.00000        
+-0.500000      -1.00000        -1.00000        
+0.00000        -1.00000        -1.00000        
+0.500000       -1.00000        -1.00000        
+1.00000        -1.00000        -1.00000        
+-1.00000       -0.500000       -1.00000        
+-0.500000      -0.500000       -1.00000        
+0.00000        -0.500000       -1.00000        
+0.500000       -0.500000       -1.00000        
+1.00000        -0.500000       -1.00000        
+-1.00000       0.00000 -1.00000        
+-0.500000      0.00000 -1.00000        
+0.00000        0.00000 -1.00000        
+0.500000       0.00000 -1.00000        
+1.00000        0.00000 -1.00000        
+-1.00000       0.500000        -1.00000        
+-0.500000      0.500000        -1.00000        
+0.00000        0.500000        -1.00000        
+0.500000       0.500000        -1.00000        
+1.00000        0.500000        -1.00000        
+-1.00000       1.00000 -1.00000        
+-0.500000      1.00000 -1.00000        
+0.00000        1.00000 -1.00000        
+0.500000       1.00000 -1.00000        
+1.00000        1.00000 -1.00000        
+-1.00000       -1.00000        -0.500000       
+-0.500000      -1.00000        -0.500000       
+0.00000        -1.00000        -0.500000       
+0.500000       -1.00000        -0.500000       
+1.00000        -1.00000        -0.500000       
+-1.00000       -0.500000       -0.500000       
+-0.500000      -0.500000       -0.500000       
+0.00000        -0.500000       -0.500000       
+0.500000       -0.500000       -0.500000       
+1.00000        -0.500000       -0.500000       
+-1.00000       0.00000 -0.500000       
+-0.500000      0.00000 -0.500000       
+0.00000        0.00000 -0.500000       
+0.500000       0.00000 -0.500000       
+1.00000        0.00000 -0.500000       
+-1.00000       0.500000        -0.500000       
+-0.500000      0.500000        -0.500000       
+0.00000        0.500000        -0.500000       
+0.500000       0.500000        -0.500000       
+1.00000        0.500000        -0.500000       
+-1.00000       1.00000 -0.500000       
+-0.500000      1.00000 -0.500000       
+0.00000        1.00000 -0.500000       
+0.500000       1.00000 -0.500000       
+1.00000        1.00000 -0.500000       
+-1.00000       -1.00000        0.00000 
+-0.500000      -1.00000        0.00000 
+0.00000        -1.00000        0.00000 
+0.500000       -1.00000        0.00000 
+1.00000        -1.00000        0.00000 
+-1.00000       -0.500000       0.00000 
+-0.500000      -0.500000       0.00000 
+0.00000        -0.500000       0.00000 
+0.500000       -0.500000       0.00000 
+1.00000        -0.500000       0.00000 
+-1.00000       0.00000 0.00000 
+-0.500000      0.00000 0.00000 
+0.00000        0.00000 0.00000 
+0.500000       0.00000 0.00000 
+1.00000        0.00000 0.00000 
+-1.00000       0.500000        0.00000 
+-0.500000      0.500000        0.00000 
+0.00000        0.500000        0.00000 
+0.500000       0.500000        0.00000 
+1.00000        0.500000        0.00000 
+-1.00000       1.00000 0.00000 
+-0.500000      1.00000 0.00000 
+0.00000        1.00000 0.00000 
+0.500000       1.00000 0.00000 
+1.00000        1.00000 0.00000 
+-1.00000       -1.00000        0.500000        
+-0.500000      -1.00000        0.500000        
+0.00000        -1.00000        0.500000        
+0.500000       -1.00000        0.500000        
+1.00000        -1.00000        0.500000        
+-1.00000       -0.500000       0.500000        
+-0.500000      -0.500000       0.500000        
+0.00000        -0.500000       0.500000        
+0.500000       -0.500000       0.500000        
+1.00000        -0.500000       0.500000        
+-1.00000       0.00000 0.500000        
+-0.500000      0.00000 0.500000        
+0.00000        0.00000 0.500000        
+0.500000       0.00000 0.500000        
+1.00000        0.00000 0.500000        
+-1.00000       0.500000        0.500000        
+-0.500000      0.500000        0.500000        
+0.00000        0.500000        0.500000        
+0.500000       0.500000        0.500000        
+1.00000        0.500000        0.500000        
+-1.00000       1.00000 0.500000        
+-0.500000      1.00000 0.500000        
+0.00000        1.00000 0.500000        
+0.500000       1.00000 0.500000        
+1.00000        1.00000 0.500000        
+-1.00000       -1.00000        1.00000 
+-0.500000      -1.00000        1.00000 
+0.00000        -1.00000        1.00000 
+0.500000       -1.00000        1.00000 
+1.00000        -1.00000        1.00000 
+-1.00000       -0.500000       1.00000 
+-0.500000      -0.500000       1.00000 
+0.00000        -0.500000       1.00000 
+0.500000       -0.500000       1.00000 
+1.00000        -0.500000       1.00000 
+-1.00000       0.00000 1.00000 
+-0.500000      0.00000 1.00000 
+0.00000        0.00000 1.00000 
+0.500000       0.00000 1.00000 
+1.00000        0.00000 1.00000 
+-1.00000       0.500000        1.00000 
+-0.500000      0.500000        1.00000 
+0.00000        0.500000        1.00000 
+0.500000       0.500000        1.00000 
+1.00000        0.500000        1.00000 
+-1.00000       1.00000 1.00000 
+-0.500000      1.00000 1.00000 
+0.00000        1.00000 1.00000 
+0.500000       1.00000 1.00000 
+1.00000        1.00000 1.00000 
 object "cells" class array type int rank 1 shape 8 items 64 data follows
 0      25      5       30      1       26      6       31
 1      26      6       31      2       27      7       32
@@ -351,131 +613,131 @@ object "cells" class array type int rank 1 shape 8 items 64 data follows
 attribute "element type" string "cubes"
 attribute "ref" string "positions"
 object "data" class array type float rank 1 shape 4 items 125 data follows
-       -2.125  0       0       0.625
-       -2.125  0       0       0.3125
-       -2.125  0       0       0
-       -2.125  0       0       -0.3125
-       -2.125  0       0       -0.625
-       -0.953125       0       0       0.625
-       -0.953125       0       0       0.3125
-       -0.953125       0       0       0
-       -0.953125       0       0       -0.3125
-       -0.953125       0       0       -0.625
-       -0.5625 0       0       0.625
-       -0.5625 0       0       0.3125
-       -0.5625 0       0       0
-       -0.5625 0       0       -0.3125
-       -0.5625 0       0       -0.625
-       -0.953125       0       0       0.625
-       -0.953125       0       0       0.3125
-       -0.953125       0       0       0
-       -0.953125       0       0       -0.3125
-       -0.953125       0       0       -0.625
-       -2.125  0       0       0.625
-       -2.125  0       0       0.3125
-       -2.125  0       0       0
-       -2.125  0       0       -0.3125
-       -2.125  0       0       -0.625
-       -0.953125       0       0       0.625
-       -0.953125       0       0       0.3125
-       -0.953125       0       0       0
-       -0.953125       0       0       -0.3125
-       -0.953125       0       0       -0.625
-       0.21875 0       0       0.625
-       0.21875 0       0       0.3125
-       0.21875 0       0       0
-       0.21875 0       0       -0.3125
-       0.21875 0       0       -0.625
-       0.609375        0       0       0.625
-       0.609375        0       0       0.3125
-       0.609375        0       0       0
-       0.609375        0       0       -0.3125
-       0.609375        0       0       -0.625
-       0.21875 0       0       0.625
-       0.21875 0       0       0.3125
-       0.21875 0       0       0
-       0.21875 0       0       -0.3125
-       0.21875 0       0       -0.625
-       -0.953125       0       0       0.625
-       -0.953125       0       0       0.3125
-       -0.953125       0       0       0
-       -0.953125       0       0       -0.3125
-       -0.953125       0       0       -0.625
-       -0.5625 0       0       0.625
-       -0.5625 0       0       0.3125
-       -0.5625 0       0       0
-       -0.5625 0       0       -0.3125
-       -0.5625 0       0       -0.625
-       0.609375        0       0       0.625
-       0.609375        0       0       0.3125
-       0.609375        0       0       0
-       0.609375        0       0       -0.3125
-       0.609375        0       0       -0.625
-       1       0       0       0.625
-       1       0       0       0.3125
-       1       0       0       0
-       1       0       0       -0.3125
-       1       0       0       -0.625
-       0.609375        0       0       0.625
-       0.609375        0       0       0.3125
-       0.609375        0       0       0
-       0.609375        0       0       -0.3125
-       0.609375        0       0       -0.625
-       -0.5625 0       0       0.625
-       -0.5625 0       0       0.3125
-       -0.5625 0       0       0
-       -0.5625 0       0       -0.3125
-       -0.5625 0       0       -0.625
-       -0.953125       0       0       0.625
-       -0.953125       0       0       0.3125
-       -0.953125       0       0       0
-       -0.953125       0       0       -0.3125
-       -0.953125       0       0       -0.625
-       0.21875 0       0       0.625
-       0.21875 0       0       0.3125
-       0.21875 0       0       0
-       0.21875 0       0       -0.3125
-       0.21875 0       0       -0.625
-       0.609375        0       0       0.625
-       0.609375        0       0       0.3125
-       0.609375        0       0       0
-       0.609375        0       0       -0.3125
-       0.609375        0       0       -0.625
-       0.21875 0       0       0.625
-       0.21875 0       0       0.3125
-       0.21875 0       0       0
-       0.21875 0       0       -0.3125
-       0.21875 0       0       -0.625
-       -0.953125       0       0       0.625
-       -0.953125       0       0       0.3125
-       -0.953125       0       0       0
-       -0.953125       0       0       -0.3125
-       -0.953125       0       0       -0.625
-       -2.125  0       0       0.625
-       -2.125  0       0       0.3125
-       -2.125  0       0       0
-       -2.125  0       0       -0.3125
-       -2.125  0       0       -0.625
-       -0.953125       0       0       0.625
-       -0.953125       0       0       0.3125
-       -0.953125       0       0       0
-       -0.953125       0       0       -0.3125
-       -0.953125       0       0       -0.625
-       -0.5625 0       0       0.625
-       -0.5625 0       0       0.3125
-       -0.5625 0       0       0
-       -0.5625 0       0       -0.3125
-       -0.5625 0       0       -0.625
-       -0.953125       0       0       0.625
-       -0.953125       0       0       0.3125
-       -0.953125       0       0       0
-       -0.953125       0       0       -0.3125
-       -0.953125       0       0       -0.625
-       -2.125  0       0       0.625
-       -2.125  0       0       0.3125
-       -2.125  0       0       0
-       -2.125  0       0       -0.3125
-       -2.125  0       0       -0.625
+       -2.12500        0.00000 0.00000 0.625000
+       -2.12500        0.00000 0.00000 0.312500
+       -2.12500        0.00000 0.00000 0.00000
+       -2.12500        0.00000 0.00000 -0.312500
+       -2.12500        0.00000 0.00000 -0.625000
+       -0.953125       0.00000 0.00000 0.625000
+       -0.953125       0.00000 0.00000 0.312500
+       -0.953125       0.00000 0.00000 0.00000
+       -0.953125       0.00000 0.00000 -0.312500
+       -0.953125       0.00000 0.00000 -0.625000
+       -0.562500       0.00000 0.00000 0.625000
+       -0.562500       0.00000 0.00000 0.312500
+       -0.562500       0.00000 0.00000 0.00000
+       -0.562500       0.00000 0.00000 -0.312500
+       -0.562500       0.00000 0.00000 -0.625000
+       -0.953125       0.00000 0.00000 0.625000
+       -0.953125       0.00000 0.00000 0.312500
+       -0.953125       0.00000 0.00000 0.00000
+       -0.953125       0.00000 0.00000 -0.312500
+       -0.953125       0.00000 0.00000 -0.625000
+       -2.12500        0.00000 0.00000 0.625000
+       -2.12500        0.00000 0.00000 0.312500
+       -2.12500        0.00000 0.00000 0.00000
+       -2.12500        0.00000 0.00000 -0.312500
+       -2.12500        0.00000 0.00000 -0.625000
+       -0.953125       0.00000 0.00000 0.625000
+       -0.953125       0.00000 0.00000 0.312500
+       -0.953125       0.00000 0.00000 0.00000
+       -0.953125       0.00000 0.00000 -0.312500
+       -0.953125       0.00000 0.00000 -0.625000
+       0.218750        0.00000 0.00000 0.625000
+       0.218750        0.00000 0.00000 0.312500
+       0.218750        0.00000 0.00000 0.00000
+       0.218750        0.00000 0.00000 -0.312500
+       0.218750        0.00000 0.00000 -0.625000
+       0.609375        0.00000 0.00000 0.625000
+       0.609375        0.00000 0.00000 0.312500
+       0.609375        0.00000 0.00000 0.00000
+       0.609375        0.00000 0.00000 -0.312500
+       0.609375        0.00000 0.00000 -0.625000
+       0.218750        0.00000 0.00000 0.625000
+       0.218750        0.00000 0.00000 0.312500
+       0.218750        0.00000 0.00000 0.00000
+       0.218750        0.00000 0.00000 -0.312500
+       0.218750        0.00000 0.00000 -0.625000
+       -0.953125       0.00000 0.00000 0.625000
+       -0.953125       0.00000 0.00000 0.312500
+       -0.953125       0.00000 0.00000 0.00000
+       -0.953125       0.00000 0.00000 -0.312500
+       -0.953125       0.00000 0.00000 -0.625000
+       -0.562500       0.00000 0.00000 0.625000
+       -0.562500       0.00000 0.00000 0.312500
+       -0.562500       0.00000 0.00000 0.00000
+       -0.562500       0.00000 0.00000 -0.312500
+       -0.562500       0.00000 0.00000 -0.625000
+       0.609375        0.00000 0.00000 0.625000
+       0.609375        0.00000 0.00000 0.312500
+       0.609375        0.00000 0.00000 0.00000
+       0.609375        0.00000 0.00000 -0.312500
+       0.609375        0.00000 0.00000 -0.625000
+       1.00000 0.00000 0.00000 0.625000
+       1.00000 0.00000 0.00000 0.312500
+       1.00000 0.00000 0.00000 0.00000
+       1.00000 0.00000 0.00000 -0.312500
+       1.00000 0.00000 0.00000 -0.625000
+       0.609375        0.00000 0.00000 0.625000
+       0.609375        0.00000 0.00000 0.312500
+       0.609375        0.00000 0.00000 0.00000
+       0.609375        0.00000 0.00000 -0.312500
+       0.609375        0.00000 0.00000 -0.625000
+       -0.562500       0.00000 0.00000 0.625000
+       -0.562500       0.00000 0.00000 0.312500
+       -0.562500       0.00000 0.00000 0.00000
+       -0.562500       0.00000 0.00000 -0.312500
+       -0.562500       0.00000 0.00000 -0.625000
+       -0.953125       0.00000 0.00000 0.625000
+       -0.953125       0.00000 0.00000 0.312500
+       -0.953125       0.00000 0.00000 0.00000
+       -0.953125       0.00000 0.00000 -0.312500
+       -0.953125       0.00000 0.00000 -0.625000
+       0.218750        0.00000 0.00000 0.625000
+       0.218750        0.00000 0.00000 0.312500
+       0.218750        0.00000 0.00000 0.00000
+       0.218750        0.00000 0.00000 -0.312500
+       0.218750        0.00000 0.00000 -0.625000
+       0.609375        0.00000 0.00000 0.625000
+       0.609375        0.00000 0.00000 0.312500
+       0.609375        0.00000 0.00000 0.00000
+       0.609375        0.00000 0.00000 -0.312500
+       0.609375        0.00000 0.00000 -0.625000
+       0.218750        0.00000 0.00000 0.625000
+       0.218750        0.00000 0.00000 0.312500
+       0.218750        0.00000 0.00000 0.00000
+       0.218750        0.00000 0.00000 -0.312500
+       0.218750        0.00000 0.00000 -0.625000
+       -0.953125       0.00000 0.00000 0.625000
+       -0.953125       0.00000 0.00000 0.312500
+       -0.953125       0.00000 0.00000 0.00000
+       -0.953125       0.00000 0.00000 -0.312500
+       -0.953125       0.00000 0.00000 -0.625000
+       -2.12500        0.00000 0.00000 0.625000
+       -2.12500        0.00000 0.00000 0.312500
+       -2.12500        0.00000 0.00000 0.00000
+       -2.12500        0.00000 0.00000 -0.312500
+       -2.12500        0.00000 0.00000 -0.625000
+       -0.953125       0.00000 0.00000 0.625000
+       -0.953125       0.00000 0.00000 0.312500
+       -0.953125       0.00000 0.00000 0.00000
+       -0.953125       0.00000 0.00000 -0.312500
+       -0.953125       0.00000 0.00000 -0.625000
+       -0.562500       0.00000 0.00000 0.625000
+       -0.562500       0.00000 0.00000 0.312500
+       -0.562500       0.00000 0.00000 0.00000
+       -0.562500       0.00000 0.00000 -0.312500
+       -0.562500       0.00000 0.00000 -0.625000
+       -0.953125       0.00000 0.00000 0.625000
+       -0.953125       0.00000 0.00000 0.312500
+       -0.953125       0.00000 0.00000 0.00000
+       -0.953125       0.00000 0.00000 -0.312500
+       -0.953125       0.00000 0.00000 -0.625000
+       -2.12500        0.00000 0.00000 0.625000
+       -2.12500        0.00000 0.00000 0.312500
+       -2.12500        0.00000 0.00000 0.00000
+       -2.12500        0.00000 0.00000 -0.312500
+       -2.12500        0.00000 0.00000 -0.625000
 attribute "dep" string "positions"
 object "deal data" class field
 component "positions" value "vertices"