Merge pull request #267 from gridap/fix_bug_restrict_array_triangulation_portion_boundary_function

Bug fixed `restrict(::AbstractArray,::TriangulationPortion)` for portions of BoundaryTriangulation(s)

Author: fverdugo

NEWS.md

@@ -16,6 +16,7 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0
 ### Fixed
  - Fixed some methods of the `sparsecsr` generic function. Since PR [#262](https://github.com/gridap/Gridap.jl/pull/262).
  - Fixed BUG in `findnz` function for `SparseMatrixCSR`. Since PR [#264](https://github.com/gridap/Gridap.jl/pull/264).
+ - Fixed `restrict(::AbstractArray,::TriangulationPortion)` for portions of triangulations extending `BoundaryTriangulation`. Since PR [#267](https://github.com/gridap/Gridap.jl/pull/267).
 
 ## [0.10.2] - 2020-5-21
 

src/Geometry/BoundaryTriangulations.jl

@@ -131,6 +131,3 @@ end
 #  D = num_cell_dims(model)
 #  face_to_mask = get_face_mask(labeling,tags,D-1)
 #end
-
-
-

src/Geometry/TriangulationPortions.jl

@@ -45,7 +45,7 @@ function get_cell_id(trian::TriangulationPortion)
 end
 
 function restrict(f::AbstractArray,trian::TriangulationPortion)
-  reindex(f,trian)
+  reindex(restrict(f,trian.oldtrian),trian.cell_to_oldcell)
 end
 
 function get_normal_vector(trian::TriangulationPortion)

test/GeometryTests/TriangulationPortionsTests.jl

@@ -30,11 +30,17 @@ test_triangulation(strian)
 ns = get_normal_vector(strian)
 @test length(ns) == num_cells(strian)
 
+v = [i for i=1:num_cells(trian)]
+cell_to_oldcell = [2,9,7]
+trian_portion = TriangulationPortion(oldtrian,cell_to_oldcell)
+trian_portion_portion = TriangulationPortion(trian_portion,[3,1])
+@test get_cell_id(trian_portion_portion)==[7,2]
+@test restrict(v, trian_portion_portion)==[7,2]
+
 #using Gridap.Visualization
 #writevtk(oldtrian,"oldtrian")
 #writevtk(btrian,"btrian",cellfields=["normal"=>nb],celldata=["oldcell"=>get_cell_id(btrian)])
 #writevtk(strian,"strian",cellfields=["normal"=>ns],
 #  celldata=["oldcell_left"=>get_cell_id(strian).left,"oldcell_right"=>get_cell_id(strian).right])
 
 end # module
-