Added Graham Scan in Haskell

chapters/computational_geometry/gift_wrapping/graham_scan/code/haskell/grahamScan.hs

@@ -0,0 +1,31 @@
+import Data.List (sortOn, minimumBy)
+import Data.Function (on)
+
+type Point = (Double, Double)
+
+angle :: Point -> Point -> Point -> Double
+angle a@(xa, ya) b@(xb, yb) c@(xc, yc)
+  | a==b || c==b = 0
+  | theta<0      = theta+2*pi
+  | otherwise    = theta
+  where thetaA = atan2 (ya-yb) (xa-xb)
+        thetaC = atan2 (yc-yb) (xc-xb)
+        theta = thetaC - thetaA
+
+ccw :: Point -> Point -> Point -> Double
+ccw (xa, ya) (xb, yb) (xc, yc) = (xb - xa) * (yc - ya) - (yb - ya) * (xc - xa)
+
+grahamScan :: [Point] -> [Point]
+grahamScan [] = []
+grahamScan pts = wrap sortedPts [p0]
+  where p0@(x, y)= minimumBy (compare `on` snd) pts
+        sortedPts = init $ sortOn (negate . angle (x, y-1)  p0) pts
+        wrap [] p = p
+        wrap (s:ss) [p] = wrap ss [s, p]
+        wrap (s:ss) (p1:p2:ps)
+          | ccw s p1 p2 < 0 = wrap (s:ss) (p2:ps)
+          | otherwise       = wrap ss (s:p1:p2:ps)
+
+main = do
+  let pts = [(x,y) | x<-[-5..5], y<-[-5..5], x^2+y^2<=5^2]
+  print $ grahamScan pts

chapters/computational_geometry/gift_wrapping/graham_scan/graham_scan.md

@@ -1,10 +1,10 @@
 # Graham Scan
 
 At around the same time of the [Jarvis March](jarvis_march.md), R. L. Graham was also developing an algorithm to find the convex hull of a random set of points {{ "gs1972" | cite }}.
-Unlike the Jarvis March, which is an $$\mathcal{O}(nh)$$ operation, the Graham Scan is $$\mathcal{O}(n\log(n))$$, where $$n$$ is the number of points and $$h$$ is the size fo the hull. 
+Unlike the Jarvis March, which is an $$\mathcal{O}(nh)$$ operation, the Graham Scan is $$\mathcal{O}(n\log(n))$$, where $$n$$ is the number of points and $$h$$ is the size fo the hull.
 This means that the complexity of the Graham Scan is not output-sensitive; moreover, there are some cases where the Jarvis March is more optimal, depending on the size of the hull and the number of points to wrap.
 
-Rather than starting at the leftmost point like the Jarvis March, the Graham scan starts at the bottom. 
+Rather than starting at the leftmost point like the Jarvis March, the Graham scan starts at the bottom.
 We then sort the distribution of points based on the angle between the bottom-most point, the origin, and each other point.
 After sorting, we go through point-by-point, searching for points that are on the convex hull and throwing out any other points.
 We do this by looking for counter-clockwise rotations.
@@ -14,6 +14,8 @@ We can find whether a rotation is counter-clockwise with trigonometric functions
 {% method %}
 {% sample lang="jl" %}
 [import:30-32, lang:"julia"](code/julia/graham.jl)
+{% sample lang="hs" %}
+[import:15-16, lang:"haskell"](code/haskell/grahamScan.hs)
 {% endmethod %}
 
 If the output of this function is 0, the points are collinear.
@@ -24,12 +26,14 @@ We basically do not want clockwise rotations, because this means we are at an in
 <!---ADD FIGURE--->
 
 To save memory and expensive `append()` operations, we ultimately look for points that should be on the hull and swap them with the first elements in the array.
-If there are $$M$$ elements on the hull, then the first $$M$$ elements in our output random distribution of points will be the hull. 
+If there are $$M$$ elements on the hull, then the first $$M$$ elements in our output random distribution of points will be the hull.
 In the end, the code should look something like this:
 
 {% method %}
 {% sample lang="jl" %}
 [import:34-69, lang:"julia"](code/julia/graham.jl)
+{% sample lang="hs" %}
+[import:18-27, lang:"haskell"](code/haskell/grahamScan.hs)
 {% endmethod %}
 
 ### Bibliography
@@ -40,13 +44,17 @@ In the end, the code should look something like this:
 
 {% method %}
 {% sample lang="jl" %}
+### Julia
 [import, lang:"julia"](code/julia/graham.jl)
+{% sample lang="hs" %}
+### Haskell
+[import, lang:"haskell"](code/haskell/grahamScan.hs)
 {% endmethod %}
 
 <script>
 MathJax.Hub.Queue(["Typeset",MathJax.Hub]);
 </script>
-$$ 
+$$
 \newcommand{\d}{\mathrm{d}}
 \newcommand{\bff}{\boldsymbol{f}}
 \newcommand{\bfg}{\boldsymbol{g}}
@@ -65,4 +73,3 @@ $$
 \newcommand{\bfomega}{\boldsymbol{\omega}}
 \newcommand{\bftau}{\boldsymbol{\tau}}
 $$
-