]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Fix else branch in SphericalManifold::normals_at_vertices
authorDaniel Arndt <daniel.arndt@iwr.uni-heidelberg.de>
Tue, 30 Jan 2018 21:43:10 +0000 (22:43 +0100)
committerDaniel Arndt <daniel.arndt@iwr.uni-heidelberg.de>
Tue, 30 Jan 2018 21:43:45 +0000 (22:43 +0100)
source/grid/manifold_lib.cc

index 78ab035496bbc1ade1937e754ae6aa64bb3c7e19..c17484212ba11cd79b01d08131796e711eb0d86e 100644 (file)
@@ -462,10 +462,12 @@ get_normals_at_vertices (const typename Triangulation<dim,spacedim>::face_iterat
                                                              distances_to_first_vertex.end());
 
   if (*minmax_distance.second-*minmax_distance.first<1.e-10 * *min_distance_to_first_vertex)
-    for (unsigned int vertex=0; vertex<n_vertices; ++vertex)
-      face_vertex_normals[vertex] = face->vertex(vertex)-center;
-
-  Manifold<dim, spacedim>::get_normals_at_vertices(face, face_vertex_normals);
+    {
+      for (unsigned int vertex=0; vertex<n_vertices; ++vertex)
+        face_vertex_normals[vertex] = face->vertex(vertex)-center;
+    }
+  else
+    Manifold<dim, spacedim>::get_normals_at_vertices(face, face_vertex_normals);
 }
 
 

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.