]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Output tensors element by element; this way, the are filtered for very small values...
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Tue, 25 Mar 2008 16:12:11 +0000 (16:12 +0000)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Tue, 25 Mar 2008 16:12:11 +0000 (16:12 +0000)
git-svn-id: https://svn.dealii.org/trunk@15933 0785d39b-7218-0410-832d-ea1e28bc413d

tests/deal.II/fe_values_view_01.cc
tests/deal.II/fe_values_view_02.cc
tests/deal.II/fe_values_view_03.cc
tests/deal.II/fe_values_view_04.cc
tests/deal.II/rt_covariant.cc
tests/deal.II/rt_hessian.cc
tests/deal.II/rt_hessian/cmp/generic

index 9a6b189685b18a6789ef87bf6ceeb8972ffaafd9..d27c97acdf55fa7853e0ea5d471de5ac25549222 100644 (file)
@@ -61,8 +61,10 @@ void test (const Triangulation<dim>& tr,
                deallog << "i=" << i << ", q=" << q << std::endl;
                deallog << "   "
                        << fe_values[single_component].value (i,q) << ' '
-                       << fe_values[single_component].gradient (i,q) << ' '
-                       << fe_values[single_component].hessian (i,q) << std::endl;
+                       << fe_values[single_component].gradient (i,q) << std::endl;
+               for (unsigned int k=0; k<dim; ++k)
+                 for (unsigned int l=0; l<dim; ++l)
+                   deallog << fe_values[single_component].hessian (i,q)[k][l] << std::endl;
 
                Assert (fe_values[single_component].value (i,q)
                        ==
@@ -108,7 +110,7 @@ int main()
 
   deallog.attach(logfile);
   deallog.depth_console (0);
-  deallog.threshold_double(1.e-12);
+  deallog.threshold_double(1.e-10);
 
   test_hyper_sphere<2>();
   test_hyper_sphere<3>();
index 23ca64d614c5ae840d2fea30fc6727fc55e6d016..a7d48f646a839650e205a32c5d42238e1695786c 100644 (file)
@@ -68,9 +68,11 @@ void test (const Triangulation<dim>& tr,
                          << fe_values[vec_components].value (i,q) << ' '
                          << fe_values[vec_components].gradient (i,q) << ' '
                          << fe_values[vec_components].divergence (i,q) << ' '
-                         << fe_values[vec_components].symmetric_gradient (i,q) << ' '
-                         << fe_values[vec_components].hessian (i,q)
-                         << std::endl;
+                         << fe_values[vec_components].symmetric_gradient (i,q) << std::endl;
+                 for (unsigned int k=0; k<dim; ++k)
+                   for (unsigned int l=0; l<dim; ++l)
+                     deallog << fe_values[vec_components].hessian (i,q)[k][l]
+                             << std::endl;
                  
                  for (unsigned int d=0; d<dim; ++d)
                    {
@@ -130,7 +132,7 @@ int main()
 
   deallog.attach(logfile);
   deallog.depth_console (0);
-  deallog.threshold_double(1.e-12);
+  deallog.threshold_double(1.e-10);
 
   test_hyper_sphere<2>();
   test_hyper_sphere<3>();
index cbb9c620ba2d6b64c3dd3ba21682775e5400580c..bb49131bc6d5c1ac9af4b79aa68a170e7ec39c85 100644 (file)
@@ -63,8 +63,10 @@ void test (const Triangulation<dim>& tr,
                deallog << "i=" << i << ", q=" << q << std::endl;
                deallog << "   "
                        << fe_values[single_component].value (i,q) << ' '
-                       << fe_values[single_component].gradient (i,q) << ' '
-                       << fe_values[single_component].hessian (i,q) << std::endl;
+                       << fe_values[single_component].gradient (i,q) << std::endl;
+               for (unsigned int k=0; k<dim; ++k)
+                 for (unsigned int l=0; l<dim; ++l)
+                   deallog << fe_values[single_component].hessian (i,q)[k][l] << std::endl;
 
                Assert (fe_values[single_component].value (i,q)
                        ==
@@ -110,7 +112,7 @@ int main()
 
   deallog.attach(logfile);
   deallog.depth_console (0);
-  deallog.threshold_double(1.e-12);
+  deallog.threshold_double(1.e-10);
 
   test_hyper_sphere<2>();
   test_hyper_sphere<3>();
index 9368301ba2400e19a55c5738d2472aa2d221389d..1b4dba201876f2a7baa19493d4edb0b9eb825ac9 100644 (file)
@@ -68,10 +68,12 @@ void test (const Triangulation<dim>& tr,
                  deallog << "i=" << i << ", q=" << q << std::endl;
                  deallog << "   "
                          << fe_values[vec_components].value (i,q) << ' '
-                         << fe_values[vec_components].gradient (i,q) << ' '
-                         << fe_values[vec_components].symmetric_gradient (i,q)
-                         << fe_values[vec_components].hessian (i,q)
-                         << std::endl;
+                         << fe_values[vec_components].gradient (i,q) << std::endl;
+                 for (unsigned int k=0; k<dim; ++k)
+                   for (unsigned int l=0; l<dim; ++l)
+                     deallog << fe_values[vec_components].symmetric_gradient (i,q)[k][l]
+                             << fe_values[vec_components].hessian (i,q)[k][l]
+                             << std::endl;
 
                  for (unsigned int d=0; d<dim; ++d)
                    {
index b78c4379ee7bd162ab3d0930aa0bc4fe20c39e77..685d4872a35d07199bc5a00872a77da13ee71836 100644 (file)
@@ -73,7 +73,7 @@ int main()
 
   deallog.attach(logfile);
   deallog.depth_console (0);
-  deallog.threshold_double(1.e-12);
+  deallog.threshold_double(1.e-10);
 
   test_hyper_sphere<2>();
   test_hyper_sphere<3>();
index c3acdc9ceec4e43075ec13af54815ec558533bb1..fad8e999b43e5ef1996f1fd4ba81b8d28117a4b5 100644 (file)
@@ -46,8 +46,10 @@ void test (const Triangulation<dim>& tr,
   FEValues<dim> fe_values (fe, quadrature, update_covariant_transformation | update_hessians);
 
   fe_values.reinit (dof.begin_active());
-  
-  deallog << fe_values.shape_hessian_component (0,0,0) << std::endl;
+
+  for (unsigned int i=0; i<dim; ++i)
+    for (unsigned int j=0; j<dim; ++j)
+      deallog << fe_values.shape_hessian_component (0,0,0)[i][j] << std::endl;
   
                                   // compare the hessian with
                                   // itself. this fails if the
@@ -79,7 +81,7 @@ int main()
 
   deallog.attach(logfile);
   deallog.depth_console (0);
-  deallog.threshold_double(1.e-12);
+  deallog.threshold_double(1.e-10);
 
   test_hyper_sphere<2>();
   test_hyper_sphere<3>();
index 7c4f8099803fae612f94a3cc983a4186b000669e..d8325b412c2859423f9b14a36ce43e39a9f56d77 100644 (file)
@@ -1,3 +1,14 @@
 
-DEAL::6.0 2.2e-10 -2.8e-11 0.0
-DEAL::6.0 4.4e-10 6.7e-10 0.0 0.0 0.0 -2.6e-11 1.9e-11 -4.6e-27
+DEAL::6.0
+DEAL::4.4e-10
+DEAL::0
+DEAL::0
+DEAL::6.0
+DEAL::6.7e-10
+DEAL::4.4e-10
+DEAL::0
+DEAL::0
+DEAL::0
+DEAL::0
+DEAL::0
+DEAL::0

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.