hasufell/cga
1
1
Fork 0
Code Issues Pull Requests Releases Wiki Activity

ALGO: fix bug when a pointer is between multiple quads

Browse Source 
We have to set a priority, otherwise a point might belong to
multiple quads (max 4 if it's in the very middle).
This commit is contained in:
hasufell 2014-11-14 22:45:12 +01:00
parent 60dc8ae535
commit 068ea04d56
No known key found for this signature in database
GPG Key ID: 220CD1C5BDEED020
1 changed files with 10 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
Algorithms/RangeSearch/Core.hs
View File

@ -18,6 +18,7 @@ module Algorithms.RangeSearch.Core
import Algebra.VectorTypes
import Algebra.Vector
import Data.Foldable (foldlM)
import Data.List (partition)
import Data.Maybe (fromJust)
import Diagrams.TwoD.Types
@ -83,12 +84,16 @@ isSEchild _ = False
quadTree :: [PT] -- ^ the points to divide
-> Square -- ^ the initial square around the points
-> QuadTree PT -- ^ the quad tree
quadTree pts' sq' = go (flip filter pts' . inRange $ sq') sq'
quadTree [] _ = TNil
quadTree [pt] _ = TLeaf pt
quadTree pts sq = TNode (quadTree nWPT . nwSq $ sq) (quadTree nEPT . neSq $ sq)
(quadTree sWPT . swSq $ sq) (quadTree sEPT . seSq $ sq)
where
go [] _ = TNil
go [pt] _ = TLeaf pt
go pts sq = TNode (quadTree pts . nwSq $ sq) (quadTree pts . neSq $ sq)
(quadTree pts . swSq $ sq) (quadTree pts . seSq $ sq)
-- this sets the priority in case a point is between multiple quads
(sWPT, sWO) = flip partition pts . inRange . swSq $ sq
(nWPT, nWO) = flip partition sWO . inRange . nwSq $ sq
(nEPT, nEO) = flip partition nWO . inRange . neSq $ sq
sEPT = flip filter nEO . inRange . seSq $ sq
-- |Get all squares of a quad tree.