From 88a25ae82a72355ac1ee5476cef7924d56c84aa8 Mon Sep 17 00:00:00 2001 From: JordiManyer Date: Tue, 1 Oct 2024 10:13:37 +1000 Subject: [PATCH] Bugfix in local_views of a BlockPMatrix --- src/Algebra.jl | 2 +- src/Geometry.jl | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/src/Algebra.jl b/src/Algebra.jl index 57bcbac..fd570a4 100644 --- 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 index 4c78824..d49e035 100644 --- 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