Bugfix in local_views of a BlockPMatrix

gridap · Oct 1, 2024 · 88a25ae · 88a25ae
1 parent a3309ee
commit 88a25ae
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/src/Algebra.jl b/src/Algebra.jl
@@ -373,7 +373,7 @@ end
 
 function local_views(a::BlockPMatrix,new_rows::BlockPRange,new_cols::BlockPRange)
   vals = map(CartesianIndices(blocksize(a))) do I
-    local_views(blocks(a)[I],blocks(new_rows)[I],blocks(new_cols)[I])
+    local_views(blocks(a)[I],blocks(new_rows)[I[1]],blocks(new_cols)[I[2]])
   end |> to_parray_of_arrays
   return map(mortar,vals)
 end

diff --git a/src/Geometry.jl b/src/Geometry.jl
@@ -482,7 +482,7 @@ end
 
 function Geometry.Triangulation(
   model::DistributedDiscreteModel;kwargs...)
-  D=num_cell_dims(model)
+  D = num_cell_dims(model)
   Triangulation(no_ghost,ReferenceFE{D},model;kwargs...)
 end