fix IsOnLine and add unit test

modules/geom/CMakeLists.txt

 add_subdirectory(src)
+add_subdirectory(tests)
 add_subdirectory(pymod)
\ No newline at end of file

modules/geom/src/composite3_op.cc

@@ -92,8 +92,9 @@ Real Distance(const Plane& p, const Vec3& v)
 
 bool IsOnLine(const Line3& l, const Vec3& v, Real ephilon)
 {
-  Vec3 tmp=CompDivide(v-l.GetOrigin(),l.GetDirection());
-  return std::fabs(tmp[0]-tmp[1])<ephilon && std::fabs(tmp[0]-tmp[2])<ephilon;
+  Vec3 s=v-l.GetOrigin();
+  Vec3 s_on_line=Dot(s, l.GetDirection())*l.GetDirection();
+  return Length2(s-s_on_line)<(ephilon*ephilon);
 }
 bool IsInPlane(const Plane& p,  const Line3& l,Real ephilon)
 {

modules/geom/tests/CMakeLists.txt

+set(GEOM_UNITTESTS
+  test_composite3.cc
+  tests.cc
+)
+
+ost_unittest(geom "${GEOM_UNITTESTS}")

modules/geom/tests/test_composite3.cc

+//------------------------------------------------------------------------------
+// This file is part of the OpenStructure project <www.openstructure.org>
+//
+// Copyright (C) 2008-2011 by the OpenStructure authors
+//
+// This library is free software; you can redistribute it and/or modify it under
+// the terms of the GNU Lesser General Public License as published by the Free
+// Software Foundation; either version 3.0 of the License, or (at your option)
+// any later version.
+// This library is distributed in the hope that it will be useful, but WITHOUT
+// ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
+// FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public License for more
+// details.
+//
+// You should have received a copy of the GNU Lesser General Public License
+// along with this library; if not, write to the Free Software Foundation, Inc.,
+// 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
+//------------------------------------------------------------------------------
+
+#include <ost/geom/geom.hh>
+
+#define BOOST_TEST_DYN_LINK
+#include <boost/test/unit_test.hpp>
+
+using namespace geom;
+
+BOOST_AUTO_TEST_SUITE( geom )
+
+
+BOOST_AUTO_TEST_CASE(line_init)
+{
+  Line3 line(geom::Vec3(0,0,0), geom::Vec3(2,0,0));
+  BOOST_CHECK_EQUAL(geom::Length(line.GetDirection()), 1.0);
+}
+
+BOOST_AUTO_TEST_CASE(is_on_line)
+{
+  Line3 line(geom::Vec3(0,0,0), geom::Vec3(1,0,0));
+  BOOST_CHECK(IsOnLine(line, geom::Vec3(0.5,0.0,0.0)));
+  BOOST_CHECK(IsOnLine(line, geom::Vec3(1.0,0.0,0.0)));
+  BOOST_CHECK(IsOnLine(line, geom::Vec3(0.0,0.0,0.0)));
+  BOOST_CHECK(IsOnLine(line, geom::Vec3(-5,0.0,0.0)));
+  BOOST_CHECK(IsOnLine(line, geom::Vec3(10.0,0.0,0.0)));
+  BOOST_CHECK(!IsOnLine(line, geom::Vec3(0.5,0.1,0.0)));
+  BOOST_CHECK(!IsOnLine(line, geom::Vec3(1.0,0.0,0.1)));
+  
+  line=Line3(geom::Vec3(1,0,0), geom::Vec3(1,1,1));
+  for (int i=-10; i<10; ++i) {
+    BOOST_CHECK(IsOnLine(line, line.At(i), 1e-6));
+  }
+  BOOST_CHECK(!IsOnLine(line, geom::Vec3(1,2,2.1), 1e-6));
+}
+
+BOOST_AUTO_TEST_SUITE_END()
\ No newline at end of file

modules/geom/tests/tests.cc

+//------------------------------------------------------------------------------
+// This file is part of the OpenStructure project <www.openstructure.org>
+//
+// Copyright (C) 2008-2011 by the OpenStructure authors
+//
+// This library is free software; you can redistribute it and/or modify it under
+// the terms of the GNU Lesser General Public License as published by the Free
+// Software Foundation; either version 3.0 of the License, or (at your option)
+// any later version.
+// This library is distributed in the hope that it will be useful, but WITHOUT
+// ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
+// FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public License for more
+// details.
+//
+// You should have received a copy of the GNU Lesser General Public License
+// along with this library; if not, write to the Free Software Foundation, Inc.,
+// 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
+//------------------------------------------------------------------------------
+#define BOOST_TEST_DYN_LINK
+#define BOOST_TEST_MODULE geom
+#include <boost/test/unit_test.hpp>
+