test against tol for lmismatch in geomspace

WIAS-PDELib · Nov 13, 2021 · a3a1e57 · a3a1e57 · j-fu · Nov 13, 2021
1 parent 37416fb
commit a3a1e57
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 4 deletions.
diff --git a/Project.toml b/Project.toml
@@ -1,7 +1,7 @@
 name = "ExtendableGrids"
 uuid = "cfc395e8-590f-11e8-1f13-43a2532b2fa8"
 authors = ["Juergen Fuhrmann <[email protected]>, Christian Merdon  <[email protected]>"]
-version = "0.8.9"
+version = "0.8.10"
 
 [deps]
 AbstractTrees = "1520ce14-60c1-5f80-bbc7-55ef81b5835c"

diff --git a/src/arraytools.jl b/src/arraytools.jl
@@ -32,7 +32,7 @@ function geomspace(a, b, ha, hb ; tol=1.0e-10, maxiterations=100)
         @assert (h0>0.0)
         @assert (hl>=h0)
         @assert((hl+h0)<l)
-        
+
         #  We need to  adjust two things:
 
         # The sum of the geometric progression must
@@ -88,7 +88,6 @@ function geomspace(a, b, ha, hb ; tol=1.0e-10, maxiterations=100)
                 q*=0.99
                 ns+=1
             end
-
             ql=q*0.9
             qr=q*1.1
             @assert ns<maxiterations "Unable to determine geomspace data after $(maxiterations) iterations"
@@ -99,7 +98,7 @@ function geomspace(a, b, ha, hb ; tol=1.0e-10, maxiterations=100)
             while abs(qr-ql)>tol && ns<maxiterations
                 ns+=1
                 mmm=lmismatch(qm,n)
-                if abs(mmm)<0
+                if abs(mmm)<tol
                     break
                 elseif lmismatch(ql,n)*mmm<0
                     qr=qm
@@ -111,6 +110,7 @@ function geomspace(a, b, ha, hb ; tol=1.0e-10, maxiterations=100)
             # increase q slightly to increase probability
             # for last interval to be <=hl
             q=qm*(1.0+tol)
+
             @assert ns<maxiterations "Unable to determine geomspace data after $(maxiterations) iterations"
             hmiss=hmismatch(q,n)
             if hmiss>1.0+tol