diff --git a/BGL/include/CGAL/boost/graph/IO/STL.h b/BGL/include/CGAL/boost/graph/IO/STL.h
index 8c558d2a6a0..d9b51a11176 100644
--- a/BGL/include/CGAL/boost/graph/IO/STL.h
+++ b/BGL/include/CGAL/boost/graph/IO/STL.h
@@ -133,6 +133,11 @@ bool read_STL(const char* fname, FaceGraph& g, const CGAL_BGL_NP_CLASS& np,
               bool verbose = true)
 {
   std::ifstream is(fname);
+  if(!is)
+  {
+    std::cerr<<"File doesn't exist."<<std::endl;
+    return false;
+  }
   return read_STL(is, g, np, verbose);
 }
 
@@ -229,7 +234,7 @@ bool write_STL(std::ostream& os,
   }
   else
   {
-    os << "solid\n";
+    os << "solid"<<std::endl;;
     for(const face_descriptor f : faces(g))
     {
       halfedge_descriptor h = halfedge(f, g);
@@ -238,13 +243,13 @@ bool write_STL(std::ostream& os,
       Point_ref r = get(vpm, source(h, g));
       Vector n = collinear(p, q, r) ? Vector(1, 0, 0) : unit_normal(p, q, r);
 
-      os << "facet normal " << n << "\nouter loop\n";
-      os << "vertex " << p << "\n";
-      os << "vertex " << q << "\n";
-      os << "vertex " << r << "\n";
-      os << "endloop\nendfacet\n";
+      os << "facet normal " << n << "\nouter loop"<<std::endl;;
+      os << "vertex " << p <<std::endl;
+      os << "vertex " << q <<std::endl;
+      os << "vertex " << r <<std::endl;
+      os << "endloop\nendfacet"<<std::endl;
     }
-    os << "endsolid\n";
+    os << "endsolid"<<std::endl;;
   }
 
   return os.good();
diff --git a/BGL/test/BGL/test_bgl_read_write.cpp b/BGL/test/BGL/test_bgl_read_write.cpp
index 149d61f3118..ddcfcd77db8 100644
--- a/BGL/test/BGL/test_bgl_read_write.cpp
+++ b/BGL/test/BGL/test_bgl_read_write.cpp
@@ -22,6 +22,7 @@
 #include <iostream>
 #include <fstream>
 
+
 typedef CGAL::Simple_cartesian<double>                                                   Kernel;
 typedef Kernel::Point_2                                                                  Point_2;
 typedef Kernel::Point_3                                                                  Point;
@@ -42,6 +43,19 @@ typedef OpenMesh::PolyMesh_ArrayKernelT</* MyTraits*/>
 
 #endif
 
+const double epsilon=1e-6;
+template < typename Type >
+bool approx_equal_nt(const Type &t1, const Type &t2)
+{
+      if (t1 == t2)
+              return true;
+      if (CGAL::abs(t1 - t2) / (CGAL::max)(CGAL::abs(t1), CGAL::abs(t2)) < std::abs(t1)*epsilon)
+              return true;
+      std::cout << " Approximate comparison failed between : " << t1 << "  and  " << t2 << "\n";
+      return false;
+}
+
+
 template <typename Mesh, typename VPM1, typename VPM2, typename VIM1, typename VIM2>
 bool are_equal_meshes(const Mesh& fg1, const VPM1 vpm1, const Mesh& fg2, const VPM2 vpm2,
                       const VIM1 vim1, const VIM2 vim2)
@@ -53,20 +67,23 @@ bool are_equal_meshes(const Mesh& fg1, const VPM1 vpm1, const Mesh& fg2, const V
      num_edges(fg1) != num_edges(fg2) ||
      num_faces(fg1) != num_faces(fg2))
     return false;
-
   std::vector<P> fg1_points, fg2_points;
   for(auto v : vertices(fg1))
+  {
     fg1_points.push_back(get(vpm1, v));
+  }
 
   for(auto v : vertices(fg2))
     fg2_points.push_back(get(vpm2, v));
 
   for(std::size_t id = 0; id < fg1_points.size(); ++id)
   {
-    if(CGAL::squared_distance(fg1_points[id], fg2_points[id]) > 1e-6)
-    {
+    P p1 = fg1_points[id];
+    P p2 = fg2_points[id];
+    if (! (approx_equal_nt(p1.x(), p2.x())
+           && approx_equal_nt(p1.y(), p2.y())
+           && approx_equal_nt(p1.z(), p2.z())))
       return false;
-    }
   }
   fg1_points.clear();
   fg1_points.shrink_to_fit();
@@ -394,7 +411,7 @@ void test_bgl_OBJ(const std::string filename)
 
   ok = CGAL::read_OBJ("data/sphere.obj", fg, CGAL::parameters::vertex_normal_map(vnm));
   assert(ok);
-  assert(num_vertices(fg) == 434 && num_faces(fg) == 864);
+  assert(num_vertices(fg) == 162 && num_faces(fg) == 320);
 
   for(const auto v : vertices(fg))
     assert(get(vnm, v) != CGAL::NULL_VECTOR);
@@ -525,7 +542,7 @@ void test_bgl_PLY(const std::string filename,
     std::ifstream is_rpm("tmp.ply");
     if(binary)
       CGAL::set_mode(is_rpm, CGAL::IO::BINARY);
-    ok = CGAL::read_PLY(is_c, fg, CGAL::parameters::vertex_color_map(vcm2)
+    ok = CGAL::read_PLY(is_rpm, fg2, CGAL::parameters::vertex_color_map(vcm2)
                                                    .face_color_map(fcm2));
     assert(ok);
     assert(are_equal_meshes(fg, fg2));
@@ -635,7 +652,7 @@ void test_bgl_STL(const std::string filename)
   {
     ok = CGAL::write_polygon_mesh("tmp.stl", fg, CGAL::parameters::vertex_point_map(cvpm));
     assert(ok);
-
+    cpoints.clear();
     Mesh fg2;
     ok = CGAL::read_polygon_mesh("tmp.stl", fg2, CGAL::parameters::vertex_point_map(cvpm));
     assert(ok);
@@ -654,6 +671,8 @@ void test_bgl_GOCAD(const char* filename)
   assert(ok);
   assert(num_vertices(fg) != 0 && num_faces(fg) != 0);
 
+  is.seekg(0);
+  fg.clear();
   clear(fg);
   std::pair<std::string, std::string> name_and_color;
   ok = CGAL::read_GOCAD(is, name_and_color, fg);
@@ -705,8 +724,8 @@ void test_bgl_GOCAD(const char* filename)
     ok = CGAL::read_GOCAD(is, cnn, fg2, CGAL::parameters::vertex_point_map(vpm2));
     assert(ok);
     assert(cnn.second.empty());
-    assert(num_vertices(fg2) == 4);
-    assert(num_faces(fg2) == 4);
+    assert(num_vertices(fg2) == 12491);
+    assert(num_faces(fg2) == 24191);
   }
 
   // @todo test wrong inputs (see tests of other formats)
@@ -768,7 +787,7 @@ int main(int argc, char** argv)
 #ifdef CGAL_USE_OPENMESH
   test_bgl_OFF<OMesh>(off_file);
 #endif
-*/
+
   // OBJ
   const char* obj_file = (argc > 2) ? argv[2] : "data/sphere.obj";
   test_bgl_OBJ<Polyhedron>(obj_file);
@@ -795,12 +814,12 @@ int main(int argc, char** argv)
 #ifdef CGAL_USE_OPENMESH
   test_bgl_STL<OMesh>(stl_file);
 #endif
-
+*/
   // GOCAD
   const char* gocad_file = (argc > 5) ? argv[5] : "data/2016206_MHT_surface.ts";
   test_bgl_GOCAD<Polyhedron>(gocad_file);
   test_bgl_GOCAD<SM>(gocad_file);
-  test_bgl_GOCAD<LCC>(gocad_file);
+ test_bgl_GOCAD<LCC>(gocad_file);
 #ifdef CGAL_USE_OPENMESH
   test_bgl_GOCAD<OMesh>(gocad_file);
 #endif
diff --git a/Distance_2/test/Distance_2/test_distance_2.cpp b/Distance_2/test/Distance_2/test_distance_2.cpp
index 0268e606295..29423a018b5 100644
--- a/Distance_2/test/Distance_2/test_distance_2.cpp
+++ b/Distance_2/test/Distance_2/test_distance_2.cpp
@@ -73,7 +73,7 @@ struct Test {
   void check_squared_distance(const O1& o1, const O2& o2, const FT& result)
   {
         assert(approx_equal_nt(CGAL::squared_distance(o1, o2), result));
-        assert(approx_equal_nt(CGAL::squared_distance(o2, o1), result));
+        assert(approx_equal_nt(CGAL::squared_distance(o2, o1), rseult));
   }
 
 
diff --git a/Linear_cell_complex/include/CGAL/Linear_cell_complex_for_combinatorial_map.h b/Linear_cell_complex/include/CGAL/Linear_cell_complex_for_combinatorial_map.h
index 63a37abde72..c7e55403115 100644
--- a/Linear_cell_complex/include/CGAL/Linear_cell_complex_for_combinatorial_map.h
+++ b/Linear_cell_complex/include/CGAL/Linear_cell_complex_for_combinatorial_map.h
@@ -205,7 +205,7 @@ namespace CGAL {
         import_from_halfedge_graph(heg, default_point_converter,
                                    origin_to_copy, copy_to_origin);
       }
-
+      using Base::clear;
     };
 
 } // namespace CGAL
diff --git a/Stream_support/include/CGAL/IO/GOCAD.h b/Stream_support/include/CGAL/IO/GOCAD.h
index 98c1d20cc50..8e0baba1e68 100644
--- a/Stream_support/include/CGAL/IO/GOCAD.h
+++ b/Stream_support/include/CGAL/IO/GOCAD.h
@@ -43,6 +43,11 @@ bool read_GOCAD(std::istream& is,
                 const CGAL_BGL_NP_CLASS&,
                 bool verbose = true)
 {
+  if(!is)
+  {
+    std::cerr<<"File doesn't exist."<<std::endl;
+    return false;
+  }
   typedef typename boost::range_value<PointRange>::type     Point;
   typedef typename boost::range_value<PolygonRange>::type   Poly;