Merge pull request #192 from gridap/fix_1d_cartesian_grid

Small bugfix in 1d CartesianGrid

Author: fverdugo

src/Geometry/CartesianGrids.jl

@@ -98,6 +98,12 @@ end
   getindex!(cache,a,Tuple(ci)...)
 end
 
+@inline function getindex!(v,a::CartesianCellNodes{1},i::Integer)
+  v[1] = i
+  v[2] = i+1
+  v
+end
+
 @inline function getindex!(v,a::CartesianCellNodes{D},i::Integer...) where D
   nodes = LinearIndices(size(a).+1)
   lnodes = CartesianIndices(tfill(2,Val{D}()))

test/GeometryTests/CartesianGridsTests.jl

@@ -10,6 +10,12 @@ using Gridap.Fields
 using Gridap.Geometry
 using Gridap.ReferenceFEs
 
+domain = (0,1)
+partition = (1,)
+grid = CartesianGrid(domain,partition)
+test_grid(grid)
+test_grid(UnstructuredGrid(grid))
+
 domain = (0.0,1.0,-1.0,2.0)
 partition = (3,4)