optimize scaled polygon.

MFraters · MFraters · commit a9125775663a · 2025-08-26T02:33:24.000+02:00
diff --git a/source/world_builder/parameters.cc b/source/world_builder/parameters.cc
@@ -796,9 +796,9 @@ namespace WorldBuilder
                     layer_type = LayerType::ICE_TOP;
                   }
                 WBAssertThrow(layer_type != LayerType::MAX_LAYERTYPE, "Could not find litho1.0 layer type " << rest_of_string);
+                const std::vector<Point<2>> scaled_polygon = Utilities::get_scaled_polygon(addition_points,1.5);
                 for (size_t index = 0; index < n_nodes; ++index )
                   {
-                    const std::vector<Point<2>> scaled_polygon = Utilities::get_scaled_polygon(addition_points,1.5);
                     // only add the point if it is actually in the feature
                     if (Utilities::polygon_contains_point(scaled_polygon, Point<2>(Datasets::LITHO1_0::coordinates_lat_long[index*2+1],Datasets::LITHO1_0::coordinates_lat_long[index*2],CoordinateSystem::spherical)))
                       {

Original file line number	Diff line number	Diff line change
`@@ -796,9 +796,9 @@ namespace WorldBuilder`
`796`	`796`	`layer_type = LayerType::ICE_TOP;`
`797`	`797`	`}`
`798`	`798`	`WBAssertThrow(layer_type != LayerType::MAX_LAYERTYPE, "Could not find litho1.0 layer type " << rest_of_string);`
	`799`	`+ const std::vector<Point<2>> scaled_polygon = Utilities::get_scaled_polygon(addition_points,1.5);`
`799`	`800`	`for (size_t index = 0; index < n_nodes; ++index )`
`800`	`801`	`{`
`801`		`- const std::vector<Point<2>> scaled_polygon = Utilities::get_scaled_polygon(addition_points,1.5);`
`802`	`802`	`// only add the point if it is actually in the feature`
`803`	`803`	`if (Utilities::polygon_contains_point(scaled_polygon, Point<2>(Datasets::LITHO1_0::coordinates_lat_long[index2+1],Datasets::LITHO1_0::coordinates_lat_long[index2],CoordinateSystem::spherical)))`
`804`	`804`	`{`