]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
added test for Hierarchical. Also changed test (val != values[k]) to allow
authorbrian <brian@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 14 Feb 2003 18:11:46 +0000 (18:11 +0000)
committerbrian <brian@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 14 Feb 2003 18:11:46 +0000 (18:11 +0000)
for some difference in the last decimal place.  Results for Hierarchical have been
spot checked.

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

tests/base/polynomial_test.cc

index 0b65e3a74cfffb3e1791a8018b3525b90fef826d..758ed8a77a214e227d99279d713e3ef199acf89e 100644 (file)
@@ -44,11 +44,12 @@ void check_poly(const Point<dim>& x,
     {
                                       // Check if compute_value is ok
       double val = p.compute_value(k,x);
-      if (val != values[k])
+      //if (val != values[k])
+      if (fabs(val - values[k]) > 5.0E-16)
        deallog << 'P' << k << ": values differ " << val << " != "
                << values[k] << std::endl;
 
-                                      // Check if compute_grad is ok
+                                       // Check if compute_grad is ok
       Tensor<1,dim> grad = p.compute_grad(k,x);
       if (grad != gradients[k])
        deallog << 'P' << k << ": gradients differ " << grad << " != "
@@ -59,7 +60,7 @@ void check_poly(const Point<dim>& x,
       if (grad2 != second[k])
        deallog << 'P' << k << ": second derivatives differ " << grad2 << " != "
                << second[k] << std::endl;
-      
+
       values[k] *= pow(10, dim);
       gradients[k] *= pow(10, dim);
       
@@ -141,5 +142,14 @@ int main()
 
   check_dimensions(p);
 
+  deallog.pop();
+  deallog.push("Hierarchical");
+
+  p.clear ();
+  for (unsigned int i=0;i<3;++i)
+    p.push_back (Hierarchical<double>(i));
+
+  check_dimensions(p);
+
   deallog.pop();
 }

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.