]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Better debuggable now.
authorwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 24 Jul 2002 13:36:55 +0000 (13:36 +0000)
committerwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 24 Jul 2002 13:36:55 +0000 (13:36 +0000)
git-svn-id: https://svn.dealii.org/trunk@6271 0785d39b-7218-0410-832d-ea1e28bc413d

tests/fe/shapes.cc

index 44560287d4a6d5322bce9cbc5bb98942d93588a6..b15b9e23ac1204d0515a18604a4f2f44f1c37715 100644 (file)
@@ -193,15 +193,19 @@ plot_face_shape_functions(Mapping<dim>& mapping,
                          for (unsigned int c=0; c<fe.get_fe().n_components(); ++c)
                            {
                              if (fe.get_fe().system_to_component_index(i).first == c)
-                               Assert ((sub.shape_value(i,k) ==
-                                        sub.shape_value_component(i,k,c))
-                                       &&
-                                       (sub.shape_grad(i,k) ==
-                                        sub.shape_grad_component(i,k,c))
-                                       &&
-                                       (sub.shape_2nd_derivative(i,k) ==
-                                        sub.shape_2nd_derivative_component(i,k,c)),
-                                       ExcInternalError())
+                               {
+                                 const double v1 = sub.shape_value(i,k),
+                                              v2 = sub.shape_value_component(i,k,c);
+                                 Assert (v1 == v2, ExcInternalError());
+
+                                 const Tensor<1,dim> g1 = sub.shape_grad(i,k),
+                                                     g2 = sub.shape_grad_component(i,k,c);
+                                 Assert (g1 == g2, ExcInternalError());
+
+                                 const Tensor<2,dim> s1 = sub.shape_2nd_derivative(i,k),
+                                                     s2 = sub.shape_2nd_derivative_component(i,k,c);
+                                 Assert (s1 == s2, ExcInternalError());
+                               }       
                              else
                                Assert ((sub.shape_value_component(i,k,c) == 0) &&
                                        (sub.shape_grad_component(i,k,c) == Tensor<1,dim>()) &&

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.