Save only used data. No id. Update memory not realloc.

44.4 seconds on 1/8 unstrut.vtu ~ 1100769 cells

info: Read geometry from file "ArcGisRivers.gml".
info: first part of advancedMapOnMesh(): 21.645936 seconds
info: Extracting mesh surface...
info: mapOnMesh(): 0.202910 seconds

Base line is: 44.55 seconds
info: Read geometry from file "ArcGisRivers.gml".
info: first part of advancedMapOnMesh(): 21.936575 seconds
info: Extracting mesh surface...
info: mapOnMesh(): 0.200034 seconds






Previous

MeshGeoToolsLib/GeoMapper.cpp

@@ -279,10 +279,11 @@ void GeoMapper::advancedMapOnMesh(
 	std::vector<int> closest_geo_point(nMeshNodes);
 	// 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)
 	{
-		auto const zero_coords = GeoLib::Point((*mesh->getNode(i))[0],
-			(*mesh->getNode(i))[1], 0.0, mesh->getNode(i)->getID());
+		zero_coords[0] = (*mesh->getNode(i))[0];
+		zero_coords[1] = (*mesh->getNode(i))[1];
 		GeoLib::Point* pnt = grid.getNearestPoint(zero_coords);
 		dist[i] = MathLib::sqrDist(*pnt, zero_coords);
 		closest_geo_point[i] = (dist[i]<=max_segment_length) ? pnt->getID() : -1;