From: Rene Gassmoeller Date: Sat, 3 Dec 2016 00:44:29 +0000 (-0700) Subject: Optimize SphericalManifold::get_new_point() X-Git-Tag: v8.5.0-rc1~347^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=a18c40988142b727dda27851363e1efb2d991869;p=dealii.git Optimize SphericalManifold::get_new_point() --- diff --git a/source/grid/manifold_lib.cc b/source/grid/manifold_lib.cc index 431cbe8d66..ef7e48ec69 100644 --- a/source/grid/manifold_lib.cc +++ b/source/grid/manifold_lib.cc @@ -275,13 +275,17 @@ SphericalManifold:: get_new_point (const std::vector > &vertices, const std::vector &weights) const { + const unsigned int n_points = vertices.size(); + double rho = 0.0; Tensor<1,spacedim> candidate; - for (unsigned int i = 0; i direction (vertices[i]-center); + rho += direction.norm()*weights[i]; + candidate += direction*weights[i]; } + // Unit norm direction. candidate /= candidate.norm();