Fix of ind2sub issue

rafavzqz · Nov 18, 2024 · 1d8546a · 1d8546a
1 parent 7266cf4
commit 1d8546a
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/geopdes/inst/msh/@msh_cartesian/msh_eval_boundary_side.m b/geopdes/inst/msh/@msh_cartesian/msh_eval_boundary_side.m
@@ -154,7 +154,7 @@
     % normal direction (the boundary manifold is not enough).
     element_list = element_list(:)';
     indices = cell (msh_side.ndim, 1);
-    [indices{:}] = ind2sub (msh_side.nel_dir, element_list);
+    [indices{:}] = ind2sub ([msh_side.nel_dir, 1], element_list);
     indices = cell2mat (indices);
 
     qn_elems = arrayfun(@(ii) {msh.boundary(iside).qn{ii}(:,indices(ii,:))}, 1:msh_side.ndim);