diff --git a/lib/mesh/utils/marking/markDoerflerQuickRecursive.m b/lib/mesh/utils/marking/markDoerflerQuickRecursive.m index 77ad15f56227bfb5b702f0e19119b15159b0587e..659354fd0cce08803ef8a140191885619b5322cd 100644 --- a/lib/mesh/utils/marking/markDoerflerQuickRecursive.m +++ b/lib/mesh/utils/marking/markDoerflerQuickRecursive.m @@ -39,7 +39,7 @@ function [permutation, nSelected] = recursiveCall(indicators, permutation, lower % and determine its position position = find(abs(abs(indicators(permutation(lower:upper)) - m) ... - min(abs(indicators(permutation(lower:upper)) - m))) < 1e-12); - [~, idx] = min(x(position)); + [~, idx] = min(indicators(position)); position = position(idx); % Choose pivot element of the overall permuation