From 33cd18d64ea3a67b90aa87b225b0d18447f01138 Mon Sep 17 00:00:00 2001 From: Eduard Braun Date: Tue, 2 Oct 2018 01:21:40 +0200 Subject: 2geom: update to c709d6b6780098d3e53363d925f7aee3c2640478 Update README and remove obsolete files --- src/2geom/sbasis-roots.cpp | 14 +++++++------- 1 file changed, 7 insertions(+), 7 deletions(-) (limited to 'src/2geom/sbasis-roots.cpp') diff --git a/src/2geom/sbasis-roots.cpp b/src/2geom/sbasis-roots.cpp index e3e5e4441..244b7efe4 100644 --- a/src/2geom/sbasis-roots.cpp +++ b/src/2geom/sbasis-roots.cpp @@ -46,11 +46,11 @@ * the gsl poly roots finder is faster than bernstein too, but we don't use it for 3 reasons: - a) it requires convertion to poly, which is numerically unstable + a) it requires conversion to poly, which is numerically unstable b) it requires gsl (which is currently not a dependency, and would bring in a whole slew of unrelated stuff) - c) it finds all roots, even complex ones. We don't want to accidently treat a nearly real root as a real root + c) it finds all roots, even complex ones. We don't want to accidentally treat a nearly real root as a real root From memory gsl poly roots was about 10 times faster than bernstein in the case where all the roots are in [0,1] for polys of order 5. I spent some time working out whether eigenvalue root finding @@ -85,7 +85,7 @@ namespace Geom{ /** Find the smallest interval that bounds a \param a sbasis function - \returns inteval + \returns interval */ @@ -113,7 +113,7 @@ OptInterval bounds_exact(SBasis const &a) { /** Find a small interval that bounds a \param a sbasis function - \returns inteval + \returns interval */ // I have no idea how this works, some clever bounding argument by jfb. @@ -231,7 +231,7 @@ static void multi_roots_internal(SBasis const &f, } return; } -////usefull? +////useful? // if (f.size()==1){ // int idxa=upper_level(levels,fa); // int idxb=upper_level(levels,fb); @@ -390,8 +390,8 @@ static std::vector fuseContiguous(std::vector const &sets, d -compute f at both ends of the given segment [a,b]. -compute bounds m