From 0ae1414ac3dcd77c41a3324c8fcfdd2d5268a24c Mon Sep 17 00:00:00 2001
From: Thomas Fischer <thomas.fischer@ufz.de>
Date: Tue, 13 Oct 2015 07:58:26 +0200
Subject: [PATCH] Init closest_geo_point with -1 instead of 0.
Saves writing it later.
42.7 seconds ( +- 2.54% )
# Conflicts:
# MeshGeoToolsLib/GeoMapper.cpp
---
MeshGeoToolsLib/GeoMapper.cpp | 5 ++---
1 file changed, 2 insertions(+), 3 deletions(-)
diff --git a/MeshGeoToolsLib/GeoMapper.cpp b/MeshGeoToolsLib/GeoMapper.cpp
index d5931f964d4..3980275ae9f 100644
--- a/MeshGeoToolsLib/GeoMapper.cpp
+++ b/MeshGeoToolsLib/GeoMapper.cpp
@@ -276,12 +276,11 @@ void GeoMapper::advancedMapOnMesh(
const unsigned nMeshNodes ( mesh->getNNodes() );
// index of closest geo point for each mesh node in (x,y)-plane
- std::vector<int> closest_geo_point(nMeshNodes);
+ std::vector<int> closest_geo_point(nMeshNodes, -1);
// distance between geo points and mesh nodes in (x,y)-plane
std::vector<double> dist(nMeshNodes);
auto zero_coords = GeoLib::Point{}; // All coordinates zero.
- for (std::size_t i=0; i<nMeshNodes; ++i)
- {
+ for (std::size_t i=0; i<nMeshNodes; ++i) {
zero_coords[0] = (*mesh->getNode(i))[0];
zero_coords[1] = (*mesh->getNode(i))[1];
GeoLib::Point* pnt = grid.getNearestPoint(zero_coords);
--
GitLab