<div dir="ltr"> Shewchuk has a good number of writings in this topic  including this random one i found. page 9 appears to be the releavant one?<div><a href="http://www.cs.berkeley.edu/~jrs/meshpapers/robnotes.pdf">http://www.cs.berkeley.edu/~jrs/meshpapers/robnotes.pdf</a><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Apr 25, 2015 at 10:20 AM, Ertugrul Söylemez <span dir="ltr"><<a href="mailto:ertesx@gmx.de" target="_blank">ertesx@gmx.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">>> My real problem is that I've got a list of points in R3 and want to<br>
>> decide if they determine a plane, meaning they are coplanar but not<br>
>> colinear.  Similarly, given a list of points in R2, I want to verify<br>
>> that they aren't colinear. Both of these can be done by converting the<br>
>> list of points to a matrix and finding the rank of the matrix, but I<br>
>> only use the rank function in the definitions of colinear and<br>
>> coplanar.<br>
>><br>
>> Maybe there's an easier way to tackle the underlying problems. Anyone<br>
>> got a suggestion for such?<br>
><br>
> I have written an experimental [implementation] of a Gauss-Jordan solver<br>
> and matrix inverter.  You might find some use for it.  It does work and<br>
> is reasonably fast, though not as fast as hmatrix.  One advantage is<br>
> that you can feed the points incrementally, and it will tell you<br>
> immediately when there is no solution.  It will also quickly reject<br>
> redundant points, even in the presence of rounding errors.<br>
<br>
</span>I should note:  The `solve` function isn't yet written, but it also<br>
doesn't do much.  Once you have fed enough relations, the matrix will<br>
already have been reduced to the identity, so you can simply extract the<br>
solutions from the relations.<br>
<br>
<br>
Greets,<br>
Ertugrul<br>
<br>_______________________________________________<br>
Haskell-Cafe mailing list<br>
<a href="mailto:Haskell-Cafe@haskell.org">Haskell-Cafe@haskell.org</a><br>
<a href="http://mail.haskell.org/cgi-bin/mailman/listinfo/haskell-cafe" target="_blank">http://mail.haskell.org/cgi-bin/mailman/listinfo/haskell-cafe</a><br>
<br></blockquote></div><br></div>