gemma: pkg/octree/tree.go annotate

annotate pkg/octree/tree.go @ 4650:f5fce22184da stree-experiment

Added a deserializer from STRTrees.

author	Sascha L. Teichmann <sascha.teichmann@intevation.de>
date	Mon, 14 Oct 2019 01:28:18 +0200
parents	8f745c353784
children	3eda5a7215ab

rev	line source
1017 a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	1 // This is Free Software under GNU Affero General Public License v >= 3.0
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	2 // without warranty, see README.md and license for details.
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	3 //
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	4 // SPDX-License-Identifier: AGPL-3.0-or-later
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	5 // License-Filename: LICENSES/AGPL-3.0.txt
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	6 //
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	7 // Copyright (C) 2018 by via donau
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	8 // – Österreichische Wasserstraßen-Gesellschaft mbH
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	9 // Software engineering by Intevation GmbH
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	10 //
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	11 // Author(s):
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	12 // * Sascha L. Teichmann <sascha.teichmann@intevation.de>
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	13
727 41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	14 package octree
707 9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	15
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	16 import (
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	17 "math"
3655 a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	18 "runtime"
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	19 "sync"
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	20
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	21 "gemma.intevation.de/gemma/pkg/common"
707 9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	22 )
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	23
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1017 diff changeset	24 // Tree is an Octree holding triangles.
727 41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	25 type Tree struct {
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1017 diff changeset	26 // EPSG is the projection.
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 707 diff changeset	27 EPSG uint32
707 9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	28
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 707 diff changeset	29 vertices []Vertex
707 9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	30 triangles [][]int32
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	31 index []int32
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	32
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1017 diff changeset	33 // Min is the lower left corner of the bbox.
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 707 diff changeset	34 Min Vertex
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1017 diff changeset	35 // Max is the upper right corner of the bbox.
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 707 diff changeset	36 Max Vertex
707 9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	37 }
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	38
2488 cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2479 diff changeset	39 type boxFrame struct {
cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2479 diff changeset	40 pos int32
cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2479 diff changeset	41 Box2D
cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2479 diff changeset	42 }
cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2479 diff changeset	43
2466 a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	44 func (ot *Tree) Vertices() []Vertex {
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	45 return ot.vertices
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	46 }
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	47
759 46fe2ae761e8 Check bounding boxes against plane, too. WIP Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 758 diff changeset	48 var scale = [4][4]float64{
46fe2ae761e8 Check bounding boxes against plane, too. WIP Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 758 diff changeset	49 {0.0, 0.0, 0.5, 0.5},
46fe2ae761e8 Check bounding boxes against plane, too. WIP Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 758 diff changeset	50 {0.5, 0.0, 1.0, 0.5},
46fe2ae761e8 Check bounding boxes against plane, too. WIP Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 758 diff changeset	51 {0.0, 0.5, 0.5, 1.0},
46fe2ae761e8 Check bounding boxes against plane, too. WIP Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 758 diff changeset	52 {0.5, 0.5, 1.0, 1.0},
46fe2ae761e8 Check bounding boxes against plane, too. WIP Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 758 diff changeset	53 }
46fe2ae761e8 Check bounding boxes against plane, too. WIP Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 758 diff changeset	54
4649 8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	55 func (ot Tree) Tin() Tin {
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	56 return &Tin{
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	57 EPSG: ot.EPSG,
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	58 Vertices: ot.vertices,
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	59 Triangles: ot.triangles,
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	60 Min: ot.Min,
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	61 Max: ot.Max,
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	62 }
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	63 }
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	64
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	65 func (ot *Tree) FindUnused() map[int32]struct{} {
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	66
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	67 used := make(map[int32]struct{}, len(ot.triangles))
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	68 for i := int32(0); i < int32(len(ot.triangles)); i++ {
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	69 used[i] = struct{}{}
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	70 }
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	71
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	72 stack := []int32{1}
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	73
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	74 for len(stack) > 0 {
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	75 top := stack[len(stack)-1]
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	76 stack = stack[:len(stack)-1]
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	77
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	78 if top > 0 { // node
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	79 if index := ot.index[top:]; len(index) > 7 {
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	80 for _, idx := range index[:8] {
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	81 if idx != 0 {
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	82 stack = append(stack, idx)
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	83 }
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	84 }
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	85 }
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	86 } else { // leaf
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	87 pos := -top - 1
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	88 n := ot.index[pos]
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	89 indices := ot.index[pos+1 : pos+1+n]
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	90 for _, idx := range indices {
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	91 delete(used, idx)
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	92 }
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	93 }
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	94 }
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	95
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	96 return used
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	97 }
8f745c353784 Finished first version of conversion tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4151 diff changeset	98
2466 a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	99 func (ot *Tree) Value(x, y float64) (float64, bool) {
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	100
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	101 // out of bounding box
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	102 if x < ot.Min.X \|\| ot.Max.X < x \|\|
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	103 y < ot.Min.Y \|\| ot.Max.Y < y {
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	104 return 0, false
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	105 }
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	106
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	107 all := Box2D{ot.Min.X, ot.Min.Y, ot.Max.X, ot.Max.Y}
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	108
2488 cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2479 diff changeset	109 stack := []boxFrame{{1, all}}
2466 a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	110
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	111 for len(stack) > 0 {
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	112 top := stack[len(stack)-1]
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	113 stack = stack[:len(stack)-1]
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	114
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	115 if top.pos > 0 { // node
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	116 if index := ot.index[top.pos:]; len(index) > 7 {
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	117 for i := 0; i < 4; i++ {
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	118 a := index[i]
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	119 b := index[i+4]
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	120 if a == 0 && b == 0 {
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	121 continue
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	122 }
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	123 dx := top.X2 - top.X1
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	124 dy := top.Y2 - top.Y1
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	125 nbox := Box2D{
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	126 dx*scale[i][0] + top.X1,
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	127 dy*scale[i][1] + top.Y1,
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	128 dx*scale[i][2] + top.X1,
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	129 dy*scale[i][3] + top.Y1,
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	130 }
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	131 if nbox.Contains(x, y) {
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	132 if a != 0 {
2488 cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2479 diff changeset	133 stack = append(stack, boxFrame{a, nbox})
2466 a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	134 }
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	135 if b != 0 {
2488 cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2479 diff changeset	136 stack = append(stack, boxFrame{b, nbox})
2466 a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	137 }
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	138 break
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	139 }
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	140 }
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	141 }
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	142 } else { // leaf
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	143 pos := -top.pos - 1
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	144 n := ot.index[pos]
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	145 indices := ot.index[pos+1 : pos+1+n]
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	146
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	147 for _, idx := range indices {
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	148 tri := ot.triangles[idx]
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	149 t := Triangle{
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	150 ot.vertices[tri[0]],
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	151 ot.vertices[tri[1]],
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	152 ot.vertices[tri[2]],
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	153 }
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	154 if t.Contains(x, y) {
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	155 return t.Plane3D().Z(x, y), true
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	156 }
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	157 }
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	158 }
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	159 }
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	160
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	161 return 0, false
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	162 }
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1790 diff changeset	163
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1017 diff changeset	164 // Vertical does a vertical cross cut from (x1, y1) to (x2, y2).
755 1a1a8b5f2d02 Vertical traversal of octree for cross sections. WIP. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 728 diff changeset	165 func (ot Tree) Vertical(x1, y1, x2, y2 float64, fn func(Triangle)) {
1a1a8b5f2d02 Vertical traversal of octree for cross sections. WIP. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 728 diff changeset	166
774 c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 761 diff changeset	167 box := Box2D{
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 761 diff changeset	168 X1: math.Min(x1, x2),
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 761 diff changeset	169 Y1: math.Min(y1, y2),
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 761 diff changeset	170 X2: math.Max(x1, x2),
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 761 diff changeset	171 Y2: math.Max(y1, y2),
755 1a1a8b5f2d02 Vertical traversal of octree for cross sections. WIP. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 728 diff changeset	172 }
1a1a8b5f2d02 Vertical traversal of octree for cross sections. WIP. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 728 diff changeset	173
1a1a8b5f2d02 Vertical traversal of octree for cross sections. WIP. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 728 diff changeset	174 // out of bounding box
774 c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 761 diff changeset	175 if box.X2 < ot.Min.X \|\| ot.Max.X < box.X1 \|\|
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 761 diff changeset	176 box.Y2 < ot.Min.Y \|\| ot.Max.Y < box.Y1 {
755 1a1a8b5f2d02 Vertical traversal of octree for cross sections. WIP. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 728 diff changeset	177 return
1a1a8b5f2d02 Vertical traversal of octree for cross sections. WIP. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 728 diff changeset	178 }
1a1a8b5f2d02 Vertical traversal of octree for cross sections. WIP. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 728 diff changeset	179
774 c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 761 diff changeset	180 line := NewPlane2D(x1, y1, x2, y2)
756 5e14000829d1 Complete vertical octree traversal. Seems not very selective. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 755 diff changeset	181
5e14000829d1 Complete vertical octree traversal. Seems not very selective. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 755 diff changeset	182 dupes := map[int32]struct{}{}
5e14000829d1 Complete vertical octree traversal. Seems not very selective. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 755 diff changeset	183
774 c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 761 diff changeset	184 all := Box2D{ot.Min.X, ot.Min.Y, ot.Max.X, ot.Max.Y}
761 033975d49c90 Removed a bit debug output on vertical traversal of octree. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 760 diff changeset	185 //log.Printf("area: %f\n", (box.x2-box.x1)*(box.y2-box.y1))
033975d49c90 Removed a bit debug output on vertical traversal of octree. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 760 diff changeset	186 //log.Printf("all: %f\n", (all.x2-all.x1)*(all.y2-all.y1))
758 0f3ba8bfa641 Simplified vertical traversal of octree. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 757 diff changeset	187
2488 cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2479 diff changeset	188 stack := []boxFrame{{1, all}}
755 1a1a8b5f2d02 Vertical traversal of octree for cross sections. WIP. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 728 diff changeset	189
1a1a8b5f2d02 Vertical traversal of octree for cross sections. WIP. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 728 diff changeset	190 for len(stack) > 0 {
1a1a8b5f2d02 Vertical traversal of octree for cross sections. WIP. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 728 diff changeset	191 top := stack[len(stack)-1]
1a1a8b5f2d02 Vertical traversal of octree for cross sections. WIP. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 728 diff changeset	192 stack = stack[:len(stack)-1]
1a1a8b5f2d02 Vertical traversal of octree for cross sections. WIP. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 728 diff changeset	193
1a1a8b5f2d02 Vertical traversal of octree for cross sections. WIP. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 728 diff changeset	194 if top.pos > 0 { // node
1790 fe1aa62195c2 Octree: Same emtry tree fix for horizontal as for vertical traversal. Re-worked both to be BCE friendly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1788 diff changeset	195 if index := ot.index[top.pos:]; len(index) > 7 {
fe1aa62195c2 Octree: Same emtry tree fix for horizontal as for vertical traversal. Re-worked both to be BCE friendly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1788 diff changeset	196 for i := 0; i < 4; i++ {
fe1aa62195c2 Octree: Same emtry tree fix for horizontal as for vertical traversal. Re-worked both to be BCE friendly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1788 diff changeset	197 a := index[i]
fe1aa62195c2 Octree: Same emtry tree fix for horizontal as for vertical traversal. Re-worked both to be BCE friendly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1788 diff changeset	198 b := index[i+4]
fe1aa62195c2 Octree: Same emtry tree fix for horizontal as for vertical traversal. Re-worked both to be BCE friendly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1788 diff changeset	199 if a == 0 && b == 0 {
fe1aa62195c2 Octree: Same emtry tree fix for horizontal as for vertical traversal. Re-worked both to be BCE friendly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1788 diff changeset	200 continue
758 0f3ba8bfa641 Simplified vertical traversal of octree. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 757 diff changeset	201 }
1790 fe1aa62195c2 Octree: Same emtry tree fix for horizontal as for vertical traversal. Re-worked both to be BCE friendly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1788 diff changeset	202 dx := top.X2 - top.X1
fe1aa62195c2 Octree: Same emtry tree fix for horizontal as for vertical traversal. Re-worked both to be BCE friendly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1788 diff changeset	203 dy := top.Y2 - top.Y1
fe1aa62195c2 Octree: Same emtry tree fix for horizontal as for vertical traversal. Re-worked both to be BCE friendly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1788 diff changeset	204 nbox := Box2D{
fe1aa62195c2 Octree: Same emtry tree fix for horizontal as for vertical traversal. Re-worked both to be BCE friendly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1788 diff changeset	205 dx*scale[i][0] + top.X1,
fe1aa62195c2 Octree: Same emtry tree fix for horizontal as for vertical traversal. Re-worked both to be BCE friendly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1788 diff changeset	206 dy*scale[i][1] + top.Y1,
fe1aa62195c2 Octree: Same emtry tree fix for horizontal as for vertical traversal. Re-worked both to be BCE friendly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1788 diff changeset	207 dx*scale[i][2] + top.X1,
fe1aa62195c2 Octree: Same emtry tree fix for horizontal as for vertical traversal. Re-worked both to be BCE friendly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1788 diff changeset	208 dy*scale[i][3] + top.Y1,
fe1aa62195c2 Octree: Same emtry tree fix for horizontal as for vertical traversal. Re-worked both to be BCE friendly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1788 diff changeset	209 }
fe1aa62195c2 Octree: Same emtry tree fix for horizontal as for vertical traversal. Re-worked both to be BCE friendly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1788 diff changeset	210 if nbox.Intersects(box) && nbox.IntersectsPlane(line) {
fe1aa62195c2 Octree: Same emtry tree fix for horizontal as for vertical traversal. Re-worked both to be BCE friendly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1788 diff changeset	211 if a != 0 {
2488 cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2479 diff changeset	212 stack = append(stack, boxFrame{a, nbox})
1790 fe1aa62195c2 Octree: Same emtry tree fix for horizontal as for vertical traversal. Re-worked both to be BCE friendly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1788 diff changeset	213 }
fe1aa62195c2 Octree: Same emtry tree fix for horizontal as for vertical traversal. Re-worked both to be BCE friendly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1788 diff changeset	214 if b != 0 {
2488 cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2479 diff changeset	215 stack = append(stack, boxFrame{b, nbox})
1790 fe1aa62195c2 Octree: Same emtry tree fix for horizontal as for vertical traversal. Re-worked both to be BCE friendly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1788 diff changeset	216 }
755 1a1a8b5f2d02 Vertical traversal of octree for cross sections. WIP. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 728 diff changeset	217 }
1a1a8b5f2d02 Vertical traversal of octree for cross sections. WIP. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 728 diff changeset	218 }
1a1a8b5f2d02 Vertical traversal of octree for cross sections. WIP. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 728 diff changeset	219 }
756 5e14000829d1 Complete vertical octree traversal. Seems not very selective. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 755 diff changeset	220 } else { // leaf
5e14000829d1 Complete vertical octree traversal. Seems not very selective. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 755 diff changeset	221 pos := -top.pos - 1
5e14000829d1 Complete vertical octree traversal. Seems not very selective. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 755 diff changeset	222 n := ot.index[pos]
5e14000829d1 Complete vertical octree traversal. Seems not very selective. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 755 diff changeset	223 indices := ot.index[pos+1 : pos+1+n]
755 1a1a8b5f2d02 Vertical traversal of octree for cross sections. WIP. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 728 diff changeset	224
756 5e14000829d1 Complete vertical octree traversal. Seems not very selective. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 755 diff changeset	225 for _, idx := range indices {
5e14000829d1 Complete vertical octree traversal. Seems not very selective. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 755 diff changeset	226 if _, found := dupes[idx]; found {
5e14000829d1 Complete vertical octree traversal. Seems not very selective. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 755 diff changeset	227 continue
5e14000829d1 Complete vertical octree traversal. Seems not very selective. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 755 diff changeset	228 }
5e14000829d1 Complete vertical octree traversal. Seems not very selective. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 755 diff changeset	229 tri := ot.triangles[idx]
5e14000829d1 Complete vertical octree traversal. Seems not very selective. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 755 diff changeset	230 t := Triangle{
5e14000829d1 Complete vertical octree traversal. Seems not very selective. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 755 diff changeset	231 ot.vertices[tri[0]],
5e14000829d1 Complete vertical octree traversal. Seems not very selective. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 755 diff changeset	232 ot.vertices[tri[1]],
5e14000829d1 Complete vertical octree traversal. Seems not very selective. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 755 diff changeset	233 ot.vertices[tri[2]],
5e14000829d1 Complete vertical octree traversal. Seems not very selective. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 755 diff changeset	234 }
5e14000829d1 Complete vertical octree traversal. Seems not very selective. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 755 diff changeset	235
774 c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 761 diff changeset	236 v0 := line.Eval(t[0].X, t[0].Y)
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 761 diff changeset	237 v1 := line.Eval(t[1].X, t[1].Y)
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 761 diff changeset	238 v2 := line.Eval(t[2].X, t[2].Y)
756 5e14000829d1 Complete vertical octree traversal. Seems not very selective. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 755 diff changeset	239
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	240 if onPlane(v0) \|\| onPlane(v1) \|\| onPlane(v2) \|\|
759 46fe2ae761e8 Check bounding boxes against plane, too. WIP Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 758 diff changeset	241 sides(sides(sides(0, v0), v1), v2) == 3 {
756 5e14000829d1 Complete vertical octree traversal. Seems not very selective. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 755 diff changeset	242 fn(&t)
5e14000829d1 Complete vertical octree traversal. Seems not very selective. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 755 diff changeset	243 }
758 0f3ba8bfa641 Simplified vertical traversal of octree. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 757 diff changeset	244 dupes[idx] = struct{}{}
756 5e14000829d1 Complete vertical octree traversal. Seems not very selective. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 755 diff changeset	245 }
755 1a1a8b5f2d02 Vertical traversal of octree for cross sections. WIP. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 728 diff changeset	246 }
1a1a8b5f2d02 Vertical traversal of octree for cross sections. WIP. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 728 diff changeset	247 }
1a1a8b5f2d02 Vertical traversal of octree for cross sections. WIP. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 728 diff changeset	248 }
1a1a8b5f2d02 Vertical traversal of octree for cross sections. WIP. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 728 diff changeset	249
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1017 diff changeset	250 // Horizontal does a horizontal cross cut.
727 41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	251 func (ot Tree) Horizontal(h float64, fn func(Triangle)) {
707 9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	252
755 1a1a8b5f2d02 Vertical traversal of octree for cross sections. WIP. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 728 diff changeset	253 if h < ot.Min.Z \|\| ot.Max.Z < h {
1a1a8b5f2d02 Vertical traversal of octree for cross sections. WIP. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 728 diff changeset	254 return
1a1a8b5f2d02 Vertical traversal of octree for cross sections. WIP. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 728 diff changeset	255 }
1a1a8b5f2d02 Vertical traversal of octree for cross sections. WIP. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 728 diff changeset	256
707 9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	257 type frame struct {
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	258 pos int32
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	259 min float64
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	260 max float64
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	261 }
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	262
756 5e14000829d1 Complete vertical octree traversal. Seems not very selective. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 755 diff changeset	263 dupes := map[int32]struct{}{}
707 9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	264
756 5e14000829d1 Complete vertical octree traversal. Seems not very selective. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 755 diff changeset	265 stack := []frame{{1, ot.Min.Z, ot.Max.Z}}
707 9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	266
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	267 for len(stack) > 0 {
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	268 top := stack[len(stack)-1]
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	269 stack = stack[:len(stack)-1]
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	270
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	271 pos := top.pos
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	272 if pos == 0 {
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	273 continue
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	274 }
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	275 min, max := top.min, top.max
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	276
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	277 if pos > 0 { // node
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	278 if mid := (max-min)*0.5 + min; h >= mid {
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	279 pos += 4 // nodes with z-bit set
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	280 min = mid
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	281 } else {
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	282 max = mid
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	283 }
2479 c85b16db8a02 Calculate better triangulation and store it into database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2466 diff changeset	284 if pos < int32(len(ot.index)) {
c85b16db8a02 Calculate better triangulation and store it into database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2466 diff changeset	285 if index := ot.index[pos:]; len(index) > 3 {
c85b16db8a02 Calculate better triangulation and store it into database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2466 diff changeset	286 stack = append(stack,
c85b16db8a02 Calculate better triangulation and store it into database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2466 diff changeset	287 frame{index[0], min, max},
c85b16db8a02 Calculate better triangulation and store it into database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2466 diff changeset	288 frame{index[1], min, max},
c85b16db8a02 Calculate better triangulation and store it into database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2466 diff changeset	289 frame{index[2], min, max},
c85b16db8a02 Calculate better triangulation and store it into database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2466 diff changeset	290 frame{index[3], min, max})
c85b16db8a02 Calculate better triangulation and store it into database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2466 diff changeset	291 }
1788 00f34a00e6d5 Octree: Don't crash in horizontal traversal if the tree is empty. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	292 }
707 9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	293 } else { // leaf
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	294 pos = -pos - 1
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	295 n := ot.index[pos]
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	296 //log.Printf("%d %d %d\n", pos, n, len(ot.index))
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	297 indices := ot.index[pos+1 : pos+1+n]
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	298
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	299 for _, idx := range indices {
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	300 if _, found := dupes[idx]; found {
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	301 continue
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	302 }
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	303 tri := ot.triangles[idx]
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 707 diff changeset	304 t := Triangle{
707 9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	305 ot.vertices[tri[0]],
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	306 ot.vertices[tri[1]],
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	307 ot.vertices[tri[2]],
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	308 }
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	309
727 41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	310 if !(math.Min(t[0].Z, math.Min(t[1].Z, t[2].Z)) > h \|\|
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	311 math.Max(t[0].Z, math.Max(t[1].Z, t[2].Z)) < h) {
707 9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	312 dupes[idx] = struct{}{}
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	313 fn(&t)
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	314 }
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	315 }
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	316 }
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	317 }
9db4ae29ded9 octree: Moved octree traveral code to own file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	318 }
2572 7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	319
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	320 func (ot Tree) Diff(other Tree) PointMap {
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	321
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	322 firstVs, secondVs := ot.Vertices(), other.Vertices()
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	323
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	324 result := make(PointMap, len(firstVs)+len(secondVs))
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	325
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	326 sliceWork(
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	327 firstVs,
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	328 result,
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	329 func(slice []Vertex, turn func([]Vertex) []Vertex) {
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	330 p := turn(nil)
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	331 for i := range slice {
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	332 v := &slice[i]
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	333 if z, found := other.Value(v.X, v.Y); found {
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	334 p = append(p, Vertex{v.X, v.Y, v.Z - z})
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	335 if len(p) == cap(p) {
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	336 p = turn(p)
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	337 }
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	338 }
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	339 }
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	340 if len(p) > 0 {
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	341 turn(p)
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	342 }
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	343 })
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	344
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	345 sliceWork(
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	346 secondVs,
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	347 result,
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	348 func(
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	349 slice []Vertex, turn func([]Vertex) []Vertex) {
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	350 p := turn(nil)
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	351 for i := range slice {
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	352 v := &slice[i]
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	353 if z, found := ot.Value(v.X, v.Y); found {
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	354 p = append(p, Vertex{v.X, v.Y, z - v.Z})
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	355 if len(p) == cap(p) {
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	356 p = turn(p)
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	357 }
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	358 }
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	359 }
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	360 if len(p) > 0 {
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	361 turn(p)
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	362 }
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	363 })
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	364
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	365 return result
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2516 diff changeset	366 }
3655 a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	367
4151 d12c2f4d3483 Made 'golint' more happy with the octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3655 diff changeset	368 func (ot *Tree) GenerateRandomVertices(n int, callback func([]Vertex)) {
3655 a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	369 var wg sync.WaitGroup
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	370
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	371 jobs := make(chan int)
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	372 out := make(chan []Vertex)
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	373 done := make(chan struct{})
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	374
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	375 cpus := runtime.NumCPU()
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	376
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	377 free := make(chan []Vertex, cpus)
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	378
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	379 for i := 0; i < cpus; i++ {
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	380 wg.Add(1)
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	381 go func() {
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	382 defer wg.Done()
4151 d12c2f4d3483 Made 'golint' more happy with the octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3655 diff changeset	383 xRange := common.Random(ot.Min.X, ot.Max.X)
d12c2f4d3483 Made 'golint' more happy with the octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3655 diff changeset	384 yRange := common.Random(ot.Min.Y, ot.Max.Y)
3655 a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	385
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	386 for size := range jobs {
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	387 var vertices []Vertex
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	388 select {
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	389 case vertices = <-free:
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	390 default:
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	391 vertices = make([]Vertex, 0, 1000)
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	392 }
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	393 for len(vertices) < size {
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	394 x, y := xRange(), yRange()
4151 d12c2f4d3483 Made 'golint' more happy with the octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3655 diff changeset	395 if z, ok := ot.Value(x, y); ok {
3655 a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	396 vertices = append(vertices, Vertex{X: x, Y: y, Z: z})
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	397 }
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	398 }
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	399 out <- vertices
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	400 }
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	401 }()
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	402 }
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	403
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	404 go func() {
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	405 defer close(done)
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	406 for vertices := range out {
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	407 callback(vertices)
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	408 select {
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	409 case free <- vertices[:0]:
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	410 default:
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	411 }
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	412 }
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	413 }()
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	414
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	415 for remain := n; remain > 0; {
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	416 if remain > 1000 {
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	417 jobs <- 1000
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	418 remain -= 1000
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	419 } else {
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	420 jobs <- remain
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	421 remain = 0
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	422 }
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	423 }
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	424 close(jobs)
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	425 wg.Wait()
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	426 close(out)
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	427 <-done
a6c671abbc35 Started with cleaning up the single beam import. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	428 }

Mercurial > gemma

annotate pkg/octree/tree.go @ 4650:f5fce22184da stree-experiment