diff pkg/mesh/triangulation.go @ 4837:3b0ec9b56d7c

Merged
author Sascha Wilde <wilde@intevation.de>
date Fri, 08 Nov 2019 13:10:56 +0100
parents 3f0382e9f302
children 5f47eeea988d
line wrap: on
line diff
--- a/pkg/mesh/triangulation.go	Fri Nov 08 13:09:59 2019 +0100
+++ b/pkg/mesh/triangulation.go	Fri Nov 08 13:10:56 2019 +0100
@@ -20,7 +20,6 @@
 package mesh
 
 import (
-	"fmt"
 	"log"
 	"math"
 
@@ -269,50 +268,3 @@
 		Max:       max,
 	}
 }
-
-func (t *Triangulation) area() float64 {
-	var result float64
-	points := t.Points
-	ts := t.Triangles
-	for i := 0; i < len(ts); i += 3 {
-		p0 := points[ts[i+0]]
-		p1 := points[ts[i+1]]
-		p2 := points[ts[i+2]]
-		result += area(p0, p1, p2)
-	}
-	return result / 2
-}
-
-// Validate performs several sanity checks on the Triangulation to check for
-// potential errors. Returns nil if no issues were found. You normally
-// shouldn't need to call this function but it can be useful for debugging.
-func (t *Triangulation) Validate() error {
-	// verify halfedges
-	for i1, i2 := range t.Halfedges {
-		if i1 != -1 && t.Halfedges[i1] != i2 {
-			return fmt.Errorf("invalid halfedge connection")
-		}
-		if i2 != -1 && t.Halfedges[i2] != int32(i1) {
-			return fmt.Errorf("invalid halfedge connection")
-		}
-	}
-
-	// verify convex hull area vs sum of triangle areas
-	hull1 := t.ConvexHull
-	hull2 := ConvexHull(t.Points)
-	area1 := polygonArea(hull1)
-	area2 := polygonArea(hull2)
-	area3 := t.area()
-	if math.Abs(area1-area2) > 1e-9 || math.Abs(area1-area3) > 1e-9 {
-		return fmt.Errorf("hull areas disagree: %f, %f, %f", area1, area2, area3)
-	}
-
-	// verify convex hull perimeter
-	perimeter1 := polygonPerimeter(hull1)
-	perimeter2 := polygonPerimeter(hull2)
-	if math.Abs(perimeter1-perimeter2) > 1e-9 {
-		return fmt.Errorf("hull perimeters disagree: %f, %f", perimeter1, perimeter2)
-	}
-
-	return nil
-}