Refactor gridbuilder to use meshbuilder directly with ranges and calculate min/max values for improved clarity

Author: ZenanH

src/utils.jl

@@ -122,12 +122,14 @@ function gridbuilder(xr::AbstractRange, yr::AbstractRange)
     h1 ≈ h2 || error("Grid spacing in x and y directions must be equal")
     h     = T2(h1)
     inv_h = T2(1 / h)
-    x1    = T2(x[1])
-    x2    = T2(x[end])
-    y1    = T2(y[1])
-    y2    = T2(y[end])
-    ξ     = meshbuilder(x1:h:x2, y1:h:y2)
-    return (nx=nx, ny=ny, ni=ni, nc=nc, h=h, x1=x[1], x2=x[end], y1=y[1], y2=y[end], inv_h=inv_h, ξ=ξ)
+    ξ     = meshbuilder(xr, yr)
+    vmin  = minimum(ξ, dims=1)
+    vmax  = maximum(ξ, dims=1)
+    x1    = T2(vmin[1])
+    x2    = T2(vmax[1])
+    y1    = T2(vmin[2])
+    y2    = T2(vmax[2])
+    return (nx=nx, ny=ny, ni=ni, nc=nc, h=h, x1=x1, x2=x2, y1=y1, y2=y2, inv_h=inv_h, ξ=ξ)
 end
 
 function gridbuilder(xr::AbstractRange, yr::AbstractRange, zr::AbstractRange)
@@ -146,13 +148,15 @@ function gridbuilder(xr::AbstractRange, yr::AbstractRange, zr::AbstractRange)
     h1 ≈ h2 ≈ h3 || error("Grid spacing in x, y and z directions must be equal")
     h     = T2(h1)
     inv_h = T2(1 / h)
-    x1    = T2(x[1])
-    x2    = T2(x[end])
-    y1    = T2(y[1])
-    y2    = T2(y[end])
-    z1    = T2(z[1])
-    z2    = T2(z[end])
-    ξ     = meshbuilder(x1:h:x2, y1:h:y2, z1:h:z2)
+    ξ     = meshbuilder(xr, yr, zr)
+    vmin  = minimum(ξ, dims=1)
+    vmax  = maximum(ξ, dims=1)
+    x1    = T2(vmin[1])
+    x2    = T2(vmax[1])
+    y1    = T2(vmin[2])
+    y2    = T2(vmax[2])
+    z1    = T2(vmin[3])
+    z2    = T2(vmax[3])
     return (nx=nx, ny=ny, nz=nz, ni=ni, nc=nc, h=h, x1=x1, x2=x2, y1=y1, y2=y2, z1=z1, 
         z2=z2, inv_h=inv_h, ξ=ξ)
 end