From 8ed0fec364a011b9f10d42c748940b0cf0b5c27a Mon Sep 17 00:00:00 2001
From: Stuart Buchanan <stuart_d_buchanan@yahoo.co.uk>
Date: Fri, 13 Jul 2012 20:00:06 +0100
Subject: [PATCH] Fix degenerate random building placement case which resulted
 in buildings being placed outside of the triangle bounds.

---
 simgear/scene/tgdb/obj.cxx | 15 +++++++++++----
 1 file changed, 11 insertions(+), 4 deletions(-)

diff --git a/simgear/scene/tgdb/obj.cxx b/simgear/scene/tgdb/obj.cxx
index 558d98c1..4263060e 100644
--- a/simgear/scene/tgdb/obj.cxx
+++ b/simgear/scene/tgdb/obj.cxx
@@ -553,12 +553,12 @@ struct SGTileGeometryBin {
         float b = mt_rand(&seed) * stepv1;
 
         // Place an object each unit of area
-        while ( num > 1.0 ) {
+        while (num > 1.0) {
 
           // Set the next location to place a building          
           a += stepv0;
           
-          if ( a + b > 1.0f ) {
+          if ((a + b) > 1.0f) {
             // Reached the end of the scan-line on v0. Reset and increment
             // scan-line on v1
             a = mt_rand(&seed) * stepv0;
@@ -567,8 +567,15 @@ struct SGTileGeometryBin {
           
           if (b > 1.0f) {
             // In a degenerate case of a single point, we might be outside the 
-            // scanline.
-            b = mt_rand(&seed) * stepv1;
+            // scanline.  Note that we need to still ensure that a+b < 1.
+            b = mt_rand(&seed) * stepv1 * (1.0f - a);
+          }
+          
+          if ((a + b) > 1.0f ) {
+            // Truly degenerate case - simply choose a random point guaranteed
+            // to fulfil the constraing of a+b < 1.
+            a = mt_rand(&seed);
+            b = mt_rand(&seed) * (1.0f - a);
           }
           
           SGVec3f randomPoint = vorigin + a*v0 + b*v1;
-- 
2.39.5