polygon :: [Coordinate] -> Region
polygon [] = Empty
polygon [x] = Empty
polygon (x:xs) = halfPlanes x (xs ++ [x])
where halfPlanes x1 (x2:xs) = HalfPlane x1 x2 `Intersect` halfPlanes x2 xs
halfPlanes xlast [] = Complement Empty |

polygon :: [Coordinate] -> Region
polygon [] = Empty
polygon [x] = Empty
polygon (x:xs) = halfPlanes x (xs ++ [x])
where halfPlanes x1 (x2:xs) = HalfPlane x1 x2 `Intersect` halfPlanes x2 xs
halfPlanes xlast [] = Complement Empty

This entry was posted on Sunday, July 29th, 2007 at 11:55 pm and is filed under Uncategorized. You can follow any responses to this entry through the RSS 2.0 feed.
You can leave a response, or trackback from your own site.

I surrender. Correct text is at

http://www.mroweczka.pl/exerc806.hs

or if you don't like Empty constructor

[Edited by admin to add code from the above URL]

I’ve removed your false starts and added the code. Try using <pre> tags around the code – that should preserve it.

Hi Wieslaw, I think this is no correct

polygon ver =

foldl Intersect Empty [HalfPlane a b| (a,b)<-zip ver (tail ver++[head ver])]

it should be

polygon ver =

foldl Intersect (Complement Empty) [HalfPlane a b| (a,b)<-zip ver (tail ver++[head ver])]

beacause “anything `Intersect` Empty” is “Empty”. You should Intersect it with all the space which is “Complement Empty”. I wrote the same solution with “Complement Empty” and it is working.

No. There is no such thing as a polygon with one or two sides. There needs to be a check, on the length of the list of coordinates, and if less than 3 it should be Empty. Then proceed as folding the Intersection into these planes.

Fero,Jon you are both right. Thanks.

Jon – assuming there is no polygon with less than 3 sides then the check should return error and not Empty.

The solution given does not guard against the case of a 2-point list. The region will then become the extended infinite line created from the 2 points. I suggest to modify it as follows:

polygon :: [Coordinate] -> Region

polygon x

| (length x) < 3 = error "polygon: requires at least 3 points"

| otherwise = halfPlane (head x) (tail x ++ [head x])

where halfPlane x1 (x2:xs) = HalfPlane x1 x2 `Intersect` halfPlane x2 xs

halfPlane _ [] = Complement Empty