]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
test for existence of restriction matrices added
authorguido <guido@0785d39b-7218-0410-832d-ea1e28bc413d>
Sat, 16 Apr 2005 07:11:11 +0000 (07:11 +0000)
committerguido <guido@0785d39b-7218-0410-832d-ea1e28bc413d>
Sat, 16 Apr 2005 07:11:11 +0000 (07:11 +0000)
git-svn-id: https://svn.dealii.org/trunk@10508 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/include/fe/fe_base.h
deal.II/deal.II/source/fe/fe.cc
tests/bits/fe_tools_08.cc
tests/results/i686-pc-linux-gnu+gcc3.2/bits/fe_tools_08.output

index 2183916fac7af086a8f7a38fbc8c3bfbcdc0a9e4..2f266802d8c87e4decfb4388f2b1b18468edbcbe 100644 (file)
@@ -771,6 +771,34 @@ class FiniteElementBase : public Subscriptor,
                                       */
     bool prolongation_is_implemented () const;
 
+                                     /**
+                                      * Return whether this element implements
+                                      * its restriction matrices. The return
+                                      * value also indicates whether a call to
+                                      * the @p get_restriction_matrix
+                                      * function will generate an error or
+                                      * not.
+                                      *
+                                      * This function is mostly here in order
+                                      * to allow us to write more efficient
+                                      * test programs which we run on all
+                                      * kinds of weird elements, and for which
+                                      * we simply need to exclude certain
+                                      * tests in case something is not
+                                      * implemented. It will in general
+                                      * probably not be a great help in
+                                      * applications, since there is not much
+                                      * one can do if one needs these features
+                                      * and they are not implemented. This
+                                      * function could be used to check
+                                      * whether a call to
+                                      * <tt>get_restriction_matrix()</tt> will
+                                      * succeed; however, one then still needs
+                                      * to cope with the lack of information
+                                      * this just expresses.
+                                      */
+    bool restriction_is_implemented () const;
+
                                     /**
                                      * Return a readonly reference to
                                      * the matrix which describes the
index 008ecdbb2a51e0d41df5637b913a59b37d64dce2..6dd0c2d3eb42e7eba64adaadd5aa79b52a1355e7 100644 (file)
@@ -278,6 +278,27 @@ FiniteElementBase<dim>::prolongation_is_implemented () const
 
 
 
+template <int dim>
+bool
+FiniteElementBase<dim>::restriction_is_implemented () const
+{
+  for (unsigned int c=0; c<GeometryInfo<dim>::children_per_cell; ++c)
+    {
+      Assert ((restriction[c].m() == this->dofs_per_cell) ||
+              (restriction[c].m() == 0),
+              ExcInternalError());
+      Assert ((restriction[c].n() == this->dofs_per_cell) ||
+              (restriction[c].n() == 0),
+              ExcInternalError());
+      if ((restriction[c].m() == 0) ||
+          (restriction[c].n() == 0))
+        return false;
+    }
+  return true;
+}
+
+
+
 template <int dim>
 bool
 FiniteElementBase<dim>::constraints_are_implemented () const
index f52ed12b6c41fdb98bebfa8f02a67f0b2b58305d..302c142f88bfc40253fcea4f0c2a5534dd1d3fdf 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2003, 2004 by the deal.II authors
+//    Copyright (C) 2003, 2004, 2005 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -43,7 +43,7 @@ check_this (const FiniteElement<dim> &fe1,
     return;
                                    // we need prolongation matrices in
                                    // fe2
-  if (!fe2.prolongation_is_implemented())
+  if (!fe2.restriction_is_implemented())
     return;
   
   std::auto_ptr<Triangulation<dim> > tria(make_tria<dim>());
index c726acf2202617bf6e45d7d0950cff8a4784c42e..9258d77de295b2247da0df50570dbdf38e958137 100644 (file)
@@ -3,86 +3,86 @@ DEAL::Checking Q1 against Q1 in 1d:
 DEAL::Checking Q1 against Q1 in 1d:
 DEAL::Checking Q1 against Q1 in 2d:
 DEAL::99. 29. 15.
-DEAL::0.0 1.0 2.0 3.5 4.0 5.0 6.0 8.0 5.0 5.0 10. 5.0 2.0 12. 6.5 15. 8.0 0.50 
+DEAL::0 1.0 2.0 3.5 4.0 5.0 6.0 8.0 5.0 5.0 10. 5.0 2.0 12. 6.5 15. 8.0 0.50 
 DEAL::Checking Q1 against Q1 in 2d:
 DEAL::99. 29. 15.
-DEAL::0.0 1.0 2.0 3.5 4.0 5.0 6.0 8.0 5.0 5.0 10. 5.0 2.0 12. 6.5 15. 8.0 0.50 
+DEAL::0 1.0 2.0 3.5 4.0 5.0 6.0 8.0 5.0 5.0 10. 5.0 2.0 12. 6.5 15. 8.0 0.50 
 DEAL::Checking Q1 against Q1 in 3d:
 DEAL::8.8e+02 1.4e+02 49.
-DEAL::0.0 5.5 12. 9.0 12. 16. 18. 21. 24. 16. 13. 16. 4.0 9.0 11. 40. 15. 34. 3.0 10. 
+DEAL::0 5.5 12. 9.0 12. 16. 18. 21. 24. 16. 13. 16. 4.0 9.0 11. 40. 15. 34. 3.0 10. 
 DEAL::Checking Q1 against Q1 in 3d:
 DEAL::8.8e+02 1.4e+02 49.
-DEAL::0.0 5.5 12. 9.0 12. 16. 18. 21. 24. 16. 13. 16. 4.0 9.0 11. 40. 15. 34. 3.0 10. 
+DEAL::0 5.5 12. 9.0 12. 16. 18. 21. 24. 16. 13. 16. 4.0 9.0 11. 40. 15. 34. 3.0 10. 
 DEAL::Checking Q1 against Q2 in 1d:
 DEAL::Checking Q2 against Q1 in 1d:
 DEAL::Checking Q1 against Q2 in 2d:
 DEAL::4.1e+02 66. 18.
-DEAL::0.0 2.0 -0.12 2.4 -0.75 5.0 3.9 1.5 7.0 9.0 8.4 8.5 10. 8.2 9.0 1.4 0.53 6.1 7.5 13. 12. 13. 15. 14. 15. -0.12 0.28 9.9 8.1 
+DEAL::0 2.0 -0.12 2.4 -0.75 5.0 3.9 1.5 7.0 9.0 8.4 8.5 10. 8.2 9.0 1.4 0.53 6.1 7.5 13. 12. 13. 15. 14. 15. -0.12 0.28 9.9 8.1 
 DEAL::Checking Q2 against Q1 in 2d:
 DEAL::2.5e+02 76. 44.
-DEAL::0.0 1.0 2.0 6.0 9.0 10. 15. 20. 12. 12. 24. 12. 4.5 34. 17. 44. 22. 0.50 
+DEAL::0 1.0 2.0 6.0 9.0 10. 15. 20. 12. 12. 24. 12. 4.5 34. 17. 44. 22. 0.50 
 DEAL::Checking Q1 against Q2 in 3d:
 DEAL::6.4e+03 4.5e+02 62.
-DEAL::0.0 1.0 5.4 15. 11. 17. 10. 31. 34. 32. 3.1 24. 23. 4.8 45. 46. 54. 33. 35. 36. 20. 
+DEAL::0 1.0 5.4 15. 11. 17. 10. 31. 34. 32. 3.1 24. 23. 4.8 45. 46. 54. 33. 35. 36. 20. 
 DEAL::Checking Q2 against Q1 in 3d:
 DEAL::3.5e+03 5.8e+02 2.4e+02
-DEAL::0.0 15. 44. 28. 45. 58. 69. 85. 1.0e+02 66. 48. 66. 14. 28. 37. 1.9e+02 62. 1.6e+02 3.0 31. 
+DEAL::0 15. 44. 28. 45. 58. 69. 85. 1.0e+02 66. 48. 66. 14. 28. 37. 1.9e+02 62. 1.6e+02 3.0 31. 
 DEAL::Checking Q1 against Q3 in 1d:
 DEAL::Checking Q3 against Q1 in 1d:
 DEAL::Checking Q1 against Q3 in 2d:
 DEAL::9.1e+02 95. 16.
-DEAL::0.0 6.8 1.5 -1.4 4.3 2.1 8.9 8.1 6.9 9.5 9.2 9.5 3.7 6.3 6.0 13. 14. 12. 14. 13. 11. 9.4 4.2 8.6 
+DEAL::0 6.8 1.5 -1.4 4.3 2.1 8.9 8.1 6.9 9.5 9.2 9.5 3.7 6.3 6.0 13. 14. 12. 14. 13. 11. 9.4 4.2 8.6 
 DEAL::Checking Q3 against Q1 in 2d:
 DEAL::4.6e+02 1.5e+02 87.
-DEAL::0.0 1.0 2.0 9.5 16. 17. 28. 36. 22. 22. 44. 22. 8.0 66. 33. 87. 44. 0.50 
+DEAL::0 1.0 2.0 9.5 16. 17. 28. 36. 22. 22. 44. 22. 8.0 66. 33. 87. 44. 0.50 
 DEAL::Checking Q1 against Q3 in 3d:
 DEAL::2.3e+04 8.6e+02 56.
-DEAL::0.0 4.0 6.4 8.5 18. 11. 14. 26. 34. 35. 24. 24. 30. 27. 50. 49. 43. 40. 33. 20. 41. 
+DEAL::0 4.0 6.4 8.5 18. 11. 14. 26. 34. 35. 24. 24. 30. 27. 50. 49. 43. 40. 33. 20. 41. 
 DEAL::Checking Q3 against Q1 in 3d:
 DEAL::8.9e+03 1.5e+03 6.6e+02
-DEAL::0.0 34. 1.1e+02 65. 1.1e+02 1.4e+02 1.7e+02 2.2e+02 2.8e+02 1.7e+02 1.2e+02 1.7e+02 32. 65. 89. 5.0e+02 1.6e+02 4.2e+02 3.0 74. 
+DEAL::0 34. 1.1e+02 65. 1.1e+02 1.4e+02 1.7e+02 2.2e+02 2.8e+02 1.7e+02 1.2e+02 1.7e+02 32. 65. 89. 5.0e+02 1.6e+02 4.2e+02 3.0 74. 
 DEAL::Checking Q2 against Q2 in 1d:
 DEAL::Checking Q2 against Q2 in 1d:
 DEAL::Checking Q2 against Q2 in 2d:
 DEAL::1.0e+03 1.7e+02 50.
-DEAL::0.0 2.0 -0.75 1.8 -4.0 10. 5.8 2.3 16. 18. 17. 18. 24. 19. 20. 3.2 1.3 15. 15. 38. 32. 33. 44. 42. 43. -0.75 -0.19 28. 22. 
+DEAL::0 2.0 -0.75 1.8 -4.0 10. 5.8 2.3 16. 18. 17. 18. 24. 19. 20. 3.2 1.3 15. 15. 38. 32. 33. 44. 42. 43. -0.75 -0.19 28. 22. 
 DEAL::Checking Q2 against Q2 in 2d:
 DEAL::1.0e+03 1.7e+02 50.
-DEAL::0.0 2.0 -0.75 1.8 -4.0 10. 5.8 2.3 16. 18. 17. 18. 24. 19. 20. 3.2 1.3 15. 15. 38. 32. 33. 44. 42. 43. -0.75 -0.19 28. 22. 
+DEAL::0 2.0 -0.75 1.8 -4.0 10. 5.8 2.3 16. 18. 17. 18. 24. 19. 20. 3.2 1.3 15. 15. 38. 32. 33. 44. 42. 43. -0.75 -0.19 28. 22. 
 DEAL::Checking Q2 against Q2 in 3d:
 DEAL::2.6e+04 2.0e+03 3.0e+02
-DEAL::0.0 -3.1 12. 48. 35. 53. 25. 1.1e+02 1.3e+02 1.3e+02 8.0 87. 91. 7.8 2.2e+02 2.0e+02 2.6e+02 1.6e+02 1.7e+02 1.8e+02 98. 
+DEAL::0 -3.1 12. 48. 35. 53. 25. 1.1e+02 1.3e+02 1.3e+02 8.0 87. 91. 7.8 2.2e+02 2.0e+02 2.6e+02 1.6e+02 1.7e+02 1.8e+02 98. 
 DEAL::Checking Q2 against Q2 in 3d:
 DEAL::2.6e+04 2.0e+03 3.0e+02
-DEAL::0.0 -3.1 12. 48. 35. 53. 25. 1.1e+02 1.3e+02 1.3e+02 8.0 87. 91. 7.8 2.2e+02 2.0e+02 2.6e+02 1.6e+02 1.7e+02 1.8e+02 98. 
+DEAL::0 -3.1 12. 48. 35. 53. 25. 1.1e+02 1.3e+02 1.3e+02 8.0 87. 91. 7.8 2.2e+02 2.0e+02 2.6e+02 1.6e+02 1.7e+02 1.8e+02 98. 
 DEAL::Checking Q2 against Q3 in 1d:
 DEAL::Checking Q3 against Q2 in 1d:
 DEAL::Checking Q2 against Q3 in 2d:
 DEAL::3.2e+03 3.4e+02 55.
-DEAL::0.0 34. 9.8 -10. 11. 12. 30. 25. 15. 32. 31. 28. 15. 26. 30. 44. 47. 38. 44. 42. 55. 40. 19. 41. 
+DEAL::0 34. 9.8 -10. 11. 12. 30. 25. 15. 32. 31. 28. 15. 26. 30. 44. 47. 38. 44. 42. 55. 40. 19. 41. 
 DEAL::Checking Q3 against Q2 in 2d:
 DEAL::1.9e+03 3.4e+02 98.
-DEAL::0.0 2.0 -1.6 0.88 -8.7 17. 8.4 3.5 29. 31. 31. 32. 44. 35. 36. 5.9 2.4 27. 26. 75. 62. 63. 87. 83. 84. -1.6 -0.84 54. 41. 
+DEAL::0 2.0 -1.6 0.88 -8.7 17. 8.4 3.5 29. 31. 31. 32. 44. 35. 36. 5.9 2.4 27. 26. 75. 62. 63. 87. 83. 84. -1.6 -0.84 54. 41. 
 DEAL::Checking Q2 against Q3 in 3d:
 DEAL::1.3e+05 5.0e+03 3.3e+02
-DEAL::0.0 22. 48. 22. 72. 72. 87. 1.2e+02 1.1e+02 1.9e+02 1.6e+02 1.3e+02 1.8e+02 1.6e+02 2.9e+02 2.8e+02 2.8e+02 2.8e+02 2.3e+02 1.4e+02 2.8e+02 
+DEAL::0 22. 48. 22. 72. 72. 87. 1.2e+02 1.1e+02 1.9e+02 1.6e+02 1.3e+02 1.8e+02 1.6e+02 2.9e+02 2.8e+02 2.8e+02 2.8e+02 2.3e+02 1.4e+02 2.8e+02 
 DEAL::Checking Q3 against Q2 in 3d:
 DEAL::6.9e+04 5.3e+03 8.2e+02
-DEAL::0.0 -12. 26. 1.2e+02 84. 1.3e+02 56. 2.9e+02 3.3e+02 3.4e+02 17. 2.2e+02 2.3e+02 14. 6.2e+02 5.5e+02 7.3e+02 4.4e+02 4.5e+02 5.0e+02 2.6e+02 
+DEAL::0 -12. 26. 1.2e+02 84. 1.3e+02 56. 2.9e+02 3.3e+02 3.4e+02 17. 2.2e+02 2.3e+02 14. 6.2e+02 5.5e+02 7.3e+02 4.4e+02 4.5e+02 5.0e+02 2.6e+02 
 DEAL::Checking Q3 against Q3 in 1d:
 DEAL::Checking Q3 against Q3 in 1d:
 DEAL::Checking Q3 against Q3 in 2d:
 DEAL::6.3e+03 6.7e+02 1.1e+02
-DEAL::0.0 67. 19. -25. 23. 25. 62. 49. 23. 57. 62. 55. 31. 51. 59. 86. 94. 71. 89. 86. 1.0e+02 84. 40. 85. 
+DEAL::0 67. 19. -25. 23. 25. 62. 49. 23. 57. 62. 55. 31. 51. 59. 86. 94. 71. 89. 86. 1.0e+02 84. 40. 85. 
 DEAL::Checking Q3 against Q3 in 2d:
 DEAL::6.3e+03 6.7e+02 1.1e+02
-DEAL::0.0 67. 19. -25. 23. 25. 62. 49. 23. 57. 62. 55. 31. 51. 59. 86. 94. 71. 89. 86. 1.0e+02 84. 40. 85. 
+DEAL::0 67. 19. -25. 23. 25. 62. 49. 23. 57. 62. 55. 31. 51. 59. 86. 94. 71. 89. 86. 1.0e+02 84. 40. 85. 
 DEAL::Checking Q3 against Q3 in 3d:
 DEAL::3.7e+05 1.4e+04 9.4e+02
-DEAL::0.0 45. 1.6e+02 33. 1.9e+02 2.2e+02 2.9e+02 3.2e+02 2.9e+02 5.5e+02 4.6e+02 3.8e+02 5.0e+02 4.8e+02 8.2e+02 8.1e+02 7.7e+02 7.6e+02 6.5e+02 4.2e+02 8.0e+02 
+DEAL::0 45. 1.6e+02 33. 1.9e+02 2.2e+02 2.9e+02 3.2e+02 2.9e+02 5.5e+02 4.6e+02 3.8e+02 5.0e+02 4.8e+02 8.2e+02 8.1e+02 7.7e+02 7.6e+02 6.5e+02 4.2e+02 8.0e+02 
 DEAL::Checking Q3 against Q3 in 3d:
 DEAL::3.7e+05 1.4e+04 9.4e+02
-DEAL::0.0 45. 1.6e+02 33. 1.9e+02 2.2e+02 2.9e+02 3.2e+02 2.9e+02 5.5e+02 4.6e+02 3.8e+02 5.0e+02 4.8e+02 8.2e+02 8.1e+02 7.7e+02 7.6e+02 6.5e+02 4.2e+02 8.0e+02 
+DEAL::0 45. 1.6e+02 33. 1.9e+02 2.2e+02 2.9e+02 3.2e+02 2.9e+02 5.5e+02 4.6e+02 3.8e+02 5.0e+02 4.8e+02 8.2e+02 8.1e+02 7.7e+02 7.6e+02 6.5e+02 4.2e+02 8.0e+02 
 DEAL::Checking DGQ0 against DGQ0 in 1d:
 DEAL::6.0 3.6 2.5
 DEAL::0.50 0.50 2.5 2.5 
@@ -97,10 +97,10 @@ DEAL::50. 17. 7.5
 DEAL::3.0 3.0 3.5 3.5 3.5 3.5 7.5 7.5 7.5 7.5 
 DEAL::Checking DGQ0 against DGQ0 in 3d:
 DEAL::2.5e+02 58. 18.
-DEAL::5.2 5.2 5.2 5.2 5.2 5.2 9.5 9.5 9.5 9.5 9.5 9.5 9.5 9.5 18. 18. 18. 18. 18. 18. 18. 18. 
+DEAL::5.3 5.3 5.3 5.3 5.3 5.3 9.5 9.5 9.5 9.5 9.5 9.5 9.5 9.5 18. 18. 18. 18. 18. 18. 18. 18. 
 DEAL::Checking DGQ0 against DGQ0 in 3d:
 DEAL::2.5e+02 58. 18.
-DEAL::5.2 5.2 5.2 5.2 5.2 5.2 9.5 9.5 9.5 9.5 9.5 9.5 9.5 9.5 18. 18. 18. 18. 18. 18. 18. 18. 
+DEAL::5.3 5.3 5.3 5.3 5.3 5.3 9.5 9.5 9.5 9.5 9.5 9.5 9.5 9.5 18. 18. 18. 18. 18. 18. 18. 18. 
 DEAL::Checking DGQ0 against DGQ1 in 1d:
 DEAL::12. 5.3 3.2
 DEAL::-0.25 0.50 0.50 1.2 1.8 2.5 2.5 3.2 
@@ -109,19 +109,19 @@ DEAL::14. 8.1 5.5
 DEAL::1.5 1.5 5.5 5.5 
 DEAL::Checking DGQ0 against DGQ1 in 2d:
 DEAL::1.9e+02 34. 10.
-DEAL::3.0 -0.38 1.9 1.5 0.88 3.5 2.0 3.5 3.5 5.0 3.5 6.1 4.9 7.5 6.0 7.5 7.5 9.0 7.5 10. 
+DEAL::3.0 -0.38 1.9 1.5 0.87 3.5 2.0 3.5 3.5 5.0 3.5 6.1 4.9 7.5 6.0 7.5 7.5 9.0 7.5 10. 
 DEAL::Checking DGQ1 against DGQ0 in 2d:
 DEAL::2.2e+02 73. 32.
 DEAL::14. 14. 16. 16. 16. 16. 32. 32. 32. 32. 
 DEAL::Checking DGQ0 against DGQ1 in 3d:
 DEAL::1.9e+03 1.7e+02 23.
-DEAL::7.5 4.2 4.7 9.2 5.2 3.9 3.9 5.0 6.5 6.5 6.9 8.0 9.5 9.5 12. 13. 14. 14. 15. 16. 18. 18
+DEAL::7.5 4.2 4.7 9.2 5.2 3.9 3.9 5.0 6.5 6.5 6.9 8.0 9.5 9.5 12. 13. 14. 14. 15. 16. 17. 17
 DEAL::Checking DGQ1 against DGQ0 in 3d:
 DEAL::2.1e+03 4.8e+02 1.4e+02
 DEAL::46. 46. 46. 46. 46. 46. 80. 80. 80. 80. 80. 80. 80. 80. 1.4e+02 1.4e+02 1.4e+02 1.4e+02 1.4e+02 1.4e+02 1.4e+02 1.4e+02 
 DEAL::Checking DGQ0 against DGQ2 in 1d:
-DEAL::18. 6.5 3.2
-DEAL::-0.25 0.12 0.50 0.50 0.88 1.2 1.8 2.1 2.5 2.5 2.9 3.2 
+DEAL::18. 6.5 3.3
+DEAL::-0.25 0.12 0.50 0.50 0.87 1.3 1.8 2.1 2.5 2.5 2.9 3.3 
 DEAL::Checking DGQ2 against DGQ0 in 1d:
 DEAL::22. 13. 8.5
 DEAL::2.5 2.5 8.5 8.5 
@@ -150,6 +150,8 @@ DEAL::Checking DGQ4 against DGQ0 in 2d:
 DEAL::1.4e+03 4.6e+02 2.0e+02
 DEAL::87. 87. 1.e+02 1.e+02 1.e+02 1.e+02 2.0e+02 2.0e+02 2.0e+02 2.0e+02 
 DEAL::Checking DGQ0 against DGQ4 in 3d:
+DEAL::3.0e+04 6.5e+02 22.
+DEAL::5.5 5.8 5.1 1.5 3.9 5.3 7.7 6.7 10. 8.0 11. 11. 12. 13. 16. 16. 18. 17. 19. 18. 21. 
 DEAL::Checking DGQ4 against DGQ0 in 3d:
 DEAL::3.2e+04 7.5e+03 2.2e+03
 DEAL::7.2e+02 7.2e+02 7.2e+02 7.2e+02 7.2e+02 7.2e+02 1.2e+03 1.2e+03 1.2e+03 1.2e+03 1.2e+03 1.2e+03 1.2e+03 1.2e+03 2.2e+03 2.2e+03 2.2e+03 2.2e+03 2.2e+03 2.2e+03 2.2e+03 2.2e+03 
@@ -161,28 +163,28 @@ DEAL::28. 12. 7.2
 DEAL::-0.25 1.5 1.5 3.2 3.8 5.5 5.5 7.2 
 DEAL::Checking DGQ1 against DGQ1 in 2d:
 DEAL::8.2e+02 1.5e+02 42.
-DEAL::14. 0.38 8.8 7.4 4.2 15. 9.0 16. 16. 22. 15. 26. 20. 31. 25. 32. 32. 38. 31. 42. 
+DEAL::13. 0.37 8.8 7.4 4.2 15. 9.0 16. 15. 22. 15. 26. 20. 31. 25. 32. 31. 38. 31. 42. 
 DEAL::Checking DGQ1 against DGQ1 in 2d:
 DEAL::8.2e+02 1.5e+02 42.
-DEAL::14. 0.38 8.8 7.4 4.2 15. 9.0 16. 16. 22. 15. 26. 20. 31. 25. 32. 32. 38. 31. 42. 
+DEAL::13. 0.37 8.8 7.4 4.2 15. 9.0 16. 15. 22. 15. 26. 20. 31. 25. 32. 31. 38. 31. 42. 
 DEAL::Checking DGQ1 against DGQ1 in 3d:
 DEAL::1.6e+04 1.4e+03 1.9e+02
-DEAL::63. 37. 40. 76. 46. 35. 33. 42. 55. 55. 57. 66. 80. 79. 97. 1.1e+02 1.2e+02 1.2e+02 1.2e+02 1.3e+02 1.4e+02 1.4e+02 
+DEAL::63. 37. 40. 76. 45. 35. 33. 42. 55. 55. 57. 66. 79. 79. 97. 1.1e+02 1.2e+02 1.2e+02 1.2e+02 1.3e+02 1.4e+02 1.4e+02 
 DEAL::Checking DGQ1 against DGQ1 in 3d:
 DEAL::1.6e+04 1.4e+03 1.9e+02
-DEAL::63. 37. 40. 76. 46. 35. 33. 42. 55. 55. 57. 66. 80. 79. 97. 1.1e+02 1.2e+02 1.2e+02 1.2e+02 1.3e+02 1.4e+02 1.4e+02 
+DEAL::63. 37. 40. 76. 45. 35. 33. 42. 55. 55. 57. 66. 79. 79. 97. 1.1e+02 1.2e+02 1.2e+02 1.2e+02 1.3e+02 1.4e+02 1.4e+02 
 DEAL::Checking DGQ1 against DGQ3 in 1d:
 DEAL::56. 17. 6.8
 DEAL::0.19 0.22 0.74 1.5 1.5 2.3 2.8 2.8 4.2 4.2 4.7 5.5 5.5 6.3 6.8 6.8 
 DEAL::Checking DGQ3 against DGQ1 in 1d:
 DEAL::60. 25. 15.
-DEAL::-0.25 3.5 3.5 7.2 7.8 12. 12. 15. 
+DEAL::-0.25 3.5 3.5 7.2 7.8 11. 12. 15. 
 DEAL::Checking DGQ1 against DGQ3 in 2d:
 DEAL::3.3e+03 2.9e+02 39.
 DEAL::15. 1.1 15. 2.1 10. 11. 12. 12. 16. 21. 15. 21. 26. 27. 28. 28. 32. 37. 31. 37. 
 DEAL::Checking DGQ3 against DGQ1 in 2d:
 DEAL::3.3e+03 5.9e+02 1.7e+02
-DEAL::56. 3.8 37. 32. 18. 63. 36. 64. 64. 90. 63. 1.1e+02 82. 1.3e+02 1.0e+02 1.3e+02 1.3e+02 1.5e+02 1.3e+02 1.7e+02 
+DEAL::55. 3.7 37. 32. 18. 63. 36. 64. 63. 91. 63. 1.1e+02 82. 1.3e+02 1.0e+02 1.3e+02 1.3e+02 1.5e+02 1.3e+02 1.7e+02 
 DEAL::Checking DGQ1 against DGQ3 in 3d:
 DEAL::1.3e+05 3.8e+03 1.8e+02
 DEAL::41. 51. 28. 1.0e+02 41. 38. 50. 64. 76. 70. 90. 80. 99. 1.3e+02 1.1e+02 1.3e+02 1.4e+02 1.5e+02 1.5e+02 1.5e+02 1.7e+02 
@@ -190,6 +192,8 @@ DEAL::Checking DGQ3 against DGQ1 in 3d:
 DEAL::1.3e+05 1.1e+04 1.5e+03
 DEAL::5.1e+02 3.0e+02 3.2e+02 6.1e+02 3.7e+02 2.8e+02 2.6e+02 3.4e+02 4.4e+02 4.4e+02 4.6e+02 5.3e+02 6.4e+02 6.4e+02 7.8e+02 8.5e+02 9.6e+02 9.6e+02 9.7e+02 1.0e+03 1.2e+03 1.2e+03 
 DEAL::Checking DGQ1 against DGQ5 in 1d:
+DEAL::84. 20. 7.2
+DEAL::-0.16 0.27 0.34 0.51 0.92 1.5 1.5 2.1 2.5 2.7 2.7 3.2 3.8 4.3 4.3 4.5 4.9 5.5 5.5 6.1 6.5 6.7 6.7 7.2 
 DEAL::Checking DGQ5 against DGQ1 in 1d:
 DEAL::93. 38. 23.
 DEAL::-0.25 5.5 5.5 11. 12. 18. 18. 23. 
@@ -204,11 +208,11 @@ DEAL::Checking DGQ5 against DGQ1 in 3d:
 DEAL::4.4e+05 3.7e+04 5.2e+03
 DEAL::1.7e+03 1.0e+03 1.1e+03 2.1e+03 1.2e+03 9.6e+02 8.9e+02 1.1e+03 1.5e+03 1.5e+03 1.5e+03 1.8e+03 2.2e+03 2.2e+03 2.6e+03 2.9e+03 3.2e+03 3.2e+03 3.3e+03 3.5e+03 3.9e+03 3.9e+03 
 DEAL::Checking DGQ2 against DGQ2 in 1d:
-DEAL::66. 23. 11.
-DEAL::-0.25 1.1 2.5 2.5 3.9 5.2 5.8 7.1 8.5 8.5 9.9 11. 
+DEAL::67. 23. 11.
+DEAL::-0.25 1.1 2.5 2.5 3.9 5.3 5.8 7.1 8.5 8.5 9.9 11. 
 DEAL::Checking DGQ2 against DGQ2 in 1d:
-DEAL::66. 23. 11.
-DEAL::-0.25 1.1 2.5 2.5 3.9 5.2 5.8 7.1 8.5 8.5 9.9 11. 
+DEAL::67. 23. 11.
+DEAL::-0.25 1.1 2.5 2.5 3.9 5.3 5.8 7.1 8.5 8.5 9.9 11. 
 DEAL::Checking DGQ2 against DGQ2 in 2d:
 DEAL::4.2e+03 5.0e+02 96.
 DEAL::35. 16. -18. 20. 0.61 21. 35. 26. 34. 35. 41. 41. 48. 55. 57. 71. 62. 70. 71. 77. 77. 84. 91. 
@@ -222,11 +226,11 @@ DEAL::Checking DGQ2 against DGQ2 in 3d:
 DEAL::1.8e+05 8.5e+03 6.4e+02
 DEAL::2.0e+02 2.0e+02 1.0e+02 -65. 52. 1.1e+02 1.9e+02 2.2e+02 2.8e+02 2.6e+02 2.7e+02 3.1e+02 4.0e+02 4.0e+02 3.7e+02 4.4e+02 4.4e+02 5.1e+02 4.6e+02 4.9e+02 5.6e+02 
 DEAL::Checking DGQ2 against DGQ2 in 1d:
-DEAL::66. 23. 11.
-DEAL::-0.25 1.1 2.5 2.5 3.9 5.2 5.8 7.1 8.5 8.5 9.9 11. 
+DEAL::67. 23. 11.
+DEAL::-0.25 1.1 2.5 2.5 3.9 5.3 5.8 7.1 8.5 8.5 9.9 11. 
 DEAL::Checking DGQ2 against DGQ2 in 1d:
-DEAL::66. 23. 11.
-DEAL::-0.25 1.1 2.5 2.5 3.9 5.2 5.8 7.1 8.5 8.5 9.9 11. 
+DEAL::67. 23. 11.
+DEAL::-0.25 1.1 2.5 2.5 3.9 5.3 5.8 7.1 8.5 8.5 9.9 11. 
 DEAL::Checking DGQ2 against DGQ2 in 2d:
 DEAL::4.2e+03 5.0e+02 96.
 DEAL::35. 16. -18. 20. 0.61 21. 35. 26. 34. 35. 41. 41. 48. 55. 57. 71. 62. 70. 71. 77. 77. 84. 91. 
@@ -252,6 +256,8 @@ DEAL::Checking DGQ4 against DGQ2 in 2d:
 DEAL::1.2e+04 1.4e+03 2.7e+02
 DEAL::97. 45. -50. 58. 4.0 57. 99. 72. 93. 99. 1.1e+02 1.1e+02 1.3e+02 1.6e+02 1.6e+02 2.0e+02 1.7e+02 1.9e+02 2.0e+02 2.1e+02 2.1e+02 2.3e+02 2.6e+02 
 DEAL::Checking DGQ2 against DGQ4 in 3d:
+DEAL::8.4e+05 1.8e+04 6.0e+02
+DEAL::1.6e+02 1.7e+02 1.4e+02 51. 1.2e+02 1.5e+02 2.2e+02 1.9e+02 2.9e+02 2.4e+02 3.1e+02 3.0e+02 3.3e+02 3.7e+02 4.3e+02 4.5e+02 5.0e+02 4.7e+02 5.4e+02 5.1e+02 5.8e+02 
 DEAL::Checking DGQ4 against DGQ2 in 3d:
 DEAL::8.5e+05 3.9e+04 3.0e+03
 DEAL::9.4e+02 9.4e+02 4.7e+02 -3.1e+02 2.4e+02 5.2e+02 8.8e+02 1.0e+03 1.3e+03 1.2e+03 1.2e+03 1.4e+03 1.8e+03 1.8e+03 1.7e+03 2.1e+03 2.1e+03 2.3e+03 2.1e+03 2.3e+03 2.6e+03 
@@ -274,6 +280,8 @@ DEAL::Checking DGQ3 against DGQ3 in 3d:
 DEAL::1.0e+06 3.1e+04 1.4e+03
 DEAL::3.3e+02 4.1e+02 2.1e+02 8.4e+02 3.3e+02 3.0e+02 4.1e+02 5.1e+02 6.2e+02 5.7e+02 7.3e+02 6.4e+02 8.0e+02 1.1e+03 9.1e+02 1.1e+03 1.2e+03 1.2e+03 1.2e+03 1.2e+03 1.4e+03 
 DEAL::Checking DGQ3 against DGQ5 in 1d:
+DEAL::1.8e+02 43. 15.
+DEAL::-0.16 0.67 1.1 1.7 2.5 3.5 3.5 4.5 5.3 5.9 6.3 7.2 7.8 8.7 9.1 9.7 11. 12. 12. 12. 13. 14. 14. 15. 
 DEAL::Checking DGQ5 against DGQ3 in 1d:
 DEAL::1.8e+02 54. 23.
 DEAL::0.19 1.6 3.4 5.5 5.5 7.6 9.4 11. 12. 14. 15. 18. 18. 20. 21. 23. 
@@ -288,6 +296,8 @@ DEAL::Checking DGQ5 against DGQ3 in 3d:
 DEAL::3.5e+06 1.0e+05 4.9e+03
 DEAL::1.1e+03 1.4e+03 7.1e+02 2.8e+03 1.1e+03 1.0e+03 1.4e+03 1.7e+03 2.1e+03 1.9e+03 2.4e+03 2.2e+03 2.7e+03 3.6e+03 3.1e+03 3.6e+03 4.0e+03 4.0e+03 4.0e+03 4.1e+03 4.7e+03 
 DEAL::Checking DGQ4 against DGQ5 in 1d:
+DEAL::2.3e+02 54. 19.
+DEAL::-0.16 0.87 1.5 2.3 3.3 4.5 4.5 5.7 6.7 7.5 8.1 9.2 9.8 11. 12. 12. 13. 15. 15. 16. 17. 17. 18. 19. 
 DEAL::Checking DGQ5 against DGQ4 in 1d:
 DEAL::2.3e+02 60. 23.
 DEAL::0.19 1.2 2.4 3.9 5.5 5.5 7.1 8.6 9.8 11. 12. 13. 14. 16. 17. 17. 19. 21. 22. 23. 
@@ -299,6 +309,8 @@ DEAL::4.7e+04 3.3e+03 3.6e+02
 DEAL::1.5e+02 -9.2 88. 45. 1.0e+02 89. 1.4e+02 1.1e+02 1.4e+02 1.5e+02 1.9e+02 1.8e+02 1.9e+02 2.1e+02 2.6e+02 2.4e+02 2.8e+02 2.9e+02 3.2e+02 2.9e+02 3.4e+02 
 DEAL::Checking DGQ4 against DGQ5 in 3d:
 DEAL::Checking DGQ5 against DGQ4 in 3d:
+DEAL::6.7e+06 1.4e+05 4.8e+03
+DEAL::1.3e+03 1.3e+03 1.1e+03 4.1e+02 9.4e+02 1.2e+03 1.8e+03 1.6e+03 2.4e+03 1.9e+03 2.4e+03 2.4e+03 2.7e+03 2.9e+03 3.5e+03 3.6e+03 4.0e+03 3.8e+03 4.3e+03 4.1e+03 4.7e+03 
 DEAL::Checking DGP0 against DGP0 in 1d:
 DEAL::Checking DGP0 against DGP0 in 1d:
 DEAL::Checking DGP0 against DGP0 in 2d:
@@ -385,18 +397,18 @@ DEAL::Checking FE_Q<1>(1)3 against FE_Q<1>(2)3 in 1d:
 DEAL::Checking FE_Q<1>(2)3 against FE_Q<1>(1)3 in 1d:
 DEAL::Checking FE_DGQ<1>(2)2 against FE_DGQ<1>(3)2 in 1d:
 DEAL::3.7e+02 76. 23.
-DEAL::0.75 0.21 1.6 4.0 3.8 3.2 4.6 7.0 4.0 6.4 7.8 7.2 7.0 9.4 11. 10. 13. 12. 14. 16. 16. 15. 17. 19. 16. 18. 20. 19. 19. 21. 23. 22. 
+DEAL::0.75 0.21 1.6 4.0 3.7 3.2 4.6 7.0 4.0 6.4 7.8 7.3 7.0 9.4 11. 10. 13. 12. 14. 16. 16. 15. 17. 19. 16. 18. 20. 19. 19. 21. 23. 22. 
 DEAL::Checking FE_DGQ<1>(3)2 against FE_DGQ<1>(2)2 in 1d:
 DEAL::3.7e+02 89. 32.
-DEAL::-1.2 2.1 5.5 2.8 6.1 9.5 5.5 8.9 12. 9.5 13. 16. 15. 18. 22. 19. 22. 26. 22. 25. 28. 26. 29. 32. 
+DEAL::-1.3 2.1 5.5 2.8 6.1 9.5 5.5 8.9 12. 9.5 13. 16. 15. 18. 22. 19. 22. 26. 22. 25. 28. 26. 29. 32. 
 DEAL::Checking FE_DGP<1>(3)1 against FE_DGP<1>(1)1 in 1d:
 DEAL::Checking FE_DGP<1>(1)1 against FE_DGP<1>(3)1 in 1d:
 DEAL::Checking FE_Q<2>(1)3 against FE_Q<2>(2)3 in 2d:
 DEAL::3.8e+03 3.6e+02 54.
-DEAL::0.0 8.0 -0.12 -2.2 17. 13. 21. 29. 25. 30. 27. 30. 1.6 20. 20. 35. 40. 49. 45. 1.6 28. 24. 
+DEAL::0 8.0 -0.12 -2.2 17. 13. 21. 29. 25. 30. 27. 30. 1.6 20. 20. 35. 40. 49. 45. 1.6 28. 24. 
 DEAL::Checking FE_Q<2>(2)3 against FE_Q<2>(1)3 in 2d:
 DEAL::2.3e+03 4.0e+02 1.3e+02
-DEAL::0.0 2.0 4.0 6.0 8.0 19. 27. 29. 31. 45. 47. 60. 36. 38. 37. 72. 74. 37. 14. 16. 1.0e+02 52. 54. 1.3e+02 68. 70. 2.5 
+DEAL::0 2.0 4.0 6.0 8.0 19. 27. 29. 31. 45. 47. 60. 36. 38. 37. 72. 74. 37. 14. 16. 1.0e+02 52. 54. 1.3e+02 68. 70. 2.5 
 DEAL::Checking FE_DGQ<2>(2)2 against FE_DGQ<2>(3)2 in 2d:
 DEAL::3.0e+04 1.9e+03 1.8e+02
 DEAL::66. 75. 65. 74. 47. 56. 52. 61. 67. 76. 66. 75. 1.2e+02 1.3e+02 1.2e+02 1.3e+02 1.4e+02 1.5e+02 1.4e+02 1.5e+02 
@@ -407,10 +419,10 @@ DEAL::Checking FE_DGP<2>(3)1 against FE_DGP<2>(1)1 in 2d:
 DEAL::Checking FE_DGP<2>(1)1 against FE_DGP<2>(3)1 in 2d:
 DEAL::Checking FE_Q<3>(1)3 against FE_Q<3>(2)3 in 3d:
 DEAL::5.8e+04 2.4e+03 1.9e+02
-DEAL::0.0 3.0 16. 45. 33. 50. 30. 93. 1.0e+02 95. 9.2 71. 69. 15. 1.4e+02 1.4e+02 1.6e+02 99. 1.1e+02 1.1e+02 60. 
+DEAL::0 3.0 16. 45. 33. 50. 30. 93. 1.0e+02 95. 9.2 71. 69. 15. 1.4e+02 1.4e+02 1.6e+02 99. 1.1e+02 1.1e+02 60. 
 DEAL::Checking FE_Q<3>(2)3 against FE_Q<3>(1)3 in 3d:
 DEAL::3.1e+04 3.0e+03 7.3e+02
-DEAL::0.0 45. 1.3e+02 84. 1.4e+02 1.8e+02 2.1e+02 2.6e+02 3.2e+02 2.0e+02 1.5e+02 2.0e+02 40. 84. 1.1e+02 5.6e+02 1.9e+02 4.7e+02 9.0 93. 
+DEAL::0 45. 1.3e+02 84. 1.4e+02 1.8e+02 2.1e+02 2.6e+02 3.2e+02 2.0e+02 1.5e+02 2.0e+02 40. 84. 1.1e+02 5.6e+02 1.9e+02 4.7e+02 9.0 93. 
 DEAL::Checking FE_DGQ<3>(2)2 against FE_DGQ<3>(3)2 in 3d:
 DEAL::1.7e+06 3.7e+04 1.2e+03
 DEAL::2.7e+02 3.2e+02 46. 2.2e+02 2.3e+02 2.5e+02 4.0e+02 3.5e+02 5.5e+02 6.4e+02 6.8e+02 6.2e+02 5.8e+02 8.7e+02 8.2e+02 8.7e+02 9.1e+02 1.0e+03 1.1e+03 1.1e+03 1.0e+03 

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.