From: Daniel Arndt <daniel.arndt@iwr.uni-heidelberg.de>
Date: Tue, 30 Jan 2018 21:43:10 +0000 (+0100)
Subject: Fix else branch in SphericalManifold::normals_at_vertices
X-Git-Tag: v9.0.0-rc1~494^2~2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=7e3de63856cdf335c6560884cba98ba3121b9d98;p=dealii.git

Fix else branch in SphericalManifold::normals_at_vertices
---

diff --git a/source/grid/manifold_lib.cc b/source/grid/manifold_lib.cc
index 78ab035496..c17484212b 100644
--- a/source/grid/manifold_lib.cc
+++ b/source/grid/manifold_lib.cc
@@ -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);
 }