From c10b63457775c205ece1ad73e0706774a3c03124 Mon Sep 17 00:00:00 2001
From: Thomas Fischer <thomas.fischer@ufz.de>
Date: Mon, 19 May 2014 10:56:45 +0200
Subject: [PATCH] [FileIO] Using station mesh density at station points.

---
 FileIO/GmshIO/GMSHPolygonTree.cpp | 14 ++++++++------
 1 file changed, 8 insertions(+), 6 deletions(-)

diff --git a/FileIO/GmshIO/GMSHPolygonTree.cpp b/FileIO/GmshIO/GMSHPolygonTree.cpp
index e8899e9a46e..a766651823a 100644
--- a/FileIO/GmshIO/GMSHPolygonTree.cpp
+++ b/FileIO/GmshIO/GMSHPolygonTree.cpp
@@ -255,13 +255,15 @@ void GMSHPolygonTree::writeSubPolygonsAsLineConstraints(std::size_t &line_offset
 
 void GMSHPolygonTree::writeStations(std::size_t & pnt_id_offset, std::size_t sfc_number, std::ostream& out) const
 {
-	const std::size_t n_stations(_stations.size());
-	for (std::size_t k(0); k<n_stations; k++) {
-		out << "Point(" << pnt_id_offset + k << ") = {" << (*(_stations[k]))[0] << "," << (*(_stations[k]))[1] << ", 0.0, ";
-		out << _mesh_density_strategy->getMeshDensityAtPoint(_stations[k]) << "};\n";
-		out << "Point { " << pnt_id_offset + k << " } In Surface { " << sfc_number << " };\n";
+	for (auto const* station : _stations) {
+		out << "Point(" << pnt_id_offset << ") = {" << (*station)[0] << ", "
+		    << (*station)[1] << ", 0.0, "
+		    << _mesh_density_strategy->getMeshDensityAtStation(station)
+		    << "}; // Station "
+		    << static_cast<GeoLib::Station const*>(station)->getName() << " \n";
+		out << "Point { " << pnt_id_offset << " } In Surface { " << sfc_number << " };\n";
+		++pnt_id_offset;
 	}
-	pnt_id_offset += n_stations;
 }
 
 void GMSHPolygonTree::writeAdditionalPointData(std::size_t & pnt_id_offset, std::size_t sfc_number, std::ostream& out) const
-- 
GitLab