From 467ff666cc84a5e6caf20bcce92bd82ade41bdc2 Mon Sep 17 00:00:00 2001
From: "SND\\wo80_cp" <SND\wo80_cp@0e2699bc-83d4-4a8f-98e7-55e24ab8c7a5>
Date: Mon, 1 Jul 2013 14:54:33 +0000
Subject: [PATCH] Fix index out of range exception in NewLocation code

git-svn-id: https://triangle.svn.codeplex.com/svn@73245 0e2699bc-83d4-4a8f-98e7-55e24ab8c7a5
---
 Triangle.NET/Triangle/NewLocation.cs | 40 ++++++++++++++++++++++------
 1 file changed, 32 insertions(+), 8 deletions(-)

diff --git a/Triangle.NET/Triangle/NewLocation.cs b/Triangle.NET/Triangle/NewLocation.cs
index 25ccd86..cd4036e 100644
--- a/Triangle.NET/Triangle/NewLocation.cs
+++ b/Triangle.NET/Triangle/NewLocation.cs
@@ -25,6 +25,12 @@ namespace TriangleNet
         Mesh mesh;
         Behavior behavior;
 
+        // Work arrays for wegde intersection
+        double[] petalx = new double[20];
+        double[] petaly = new double[20];
+        double[] petalr = new double[20];
+        double[] wedges = new double[500];
+
         public NewLocation(Mesh mesh)
         {
             this.mesh = mesh;
@@ -2500,11 +2506,20 @@ namespace TriangleNet
 
             //double p[5];
 
-            double[] petalx = new double[2 * numpoints];
-            double[] petaly = new double[2 * numpoints];
-            double[] petalr = new double[2 * numpoints];
+            // Resize work arrays
+            if (2 * numpoints > petalx.Length)
+            {
+                petalx = new double[2 * numpoints];
+                petaly = new double[2 * numpoints];
+                petalr = new double[2 * numpoints];
+                wedges = new double[2 * numpoints * 16 + 36];
+
+                if (Behavior.Verbose)
+                {
+                    Log.SimpleLog.Instance.Info("NewLocation: resized work arrays (" + 2 * numpoints + ")");
+                }
+            }
 
-            double[] wedges = new double[2000];
             double xmid, ymid, dist, x3, y3;
             double x_1, y_1, x_2, y_2, x_3, y_3, x_4, y_4, tempx, tempy;
             double ux, uy;
@@ -2758,11 +2773,20 @@ namespace TriangleNet
 
             //double p[5];
 
-            double[] petalx = new double[2 * numpoints];
-            double[] petaly = new double[2 * numpoints];
-            double[] petalr = new double[2 * numpoints];
+            // Resize work arrays
+            if (2 * numpoints > petalx.Length)
+            {
+                petalx = new double[2 * numpoints];
+                petaly = new double[2 * numpoints];
+                petalr = new double[2 * numpoints];
+                wedges = new double[2 * numpoints * 20 + 40];
+
+                if (Behavior.Verbose)
+                {
+                    Log.SimpleLog.Instance.Info("NewLocation: resized work arrays (" + 2 * numpoints + ")");
+                }
+            }
 
-            double[] wedges = new double[2000];
             double xmid, ymid, dist, x3, y3;
             double x_1, y_1, x_2, y_2, x_3, y_3, x_4, y_4, tempx, tempy, x_5, y_5, x_6, y_6;
             double ux, uy;