remove the prior closure of the cone THIS FUNCTION IS NO USER FUNCTION.
0001 function data = openCone(data, zerotol) 0002 % remove the prior closure of the cone 0003 % 0004 % THIS FUNCTION IS NO USER FUNCTION. 0005 0006 % The elk-library: convex geometry applied to crystallization modeling. 0007 % Copyright (C) 2013 Alexander Reinhold 0008 % 0009 % This program is free software: you can redistribute it and/or modify it 0010 % under the terms of the GNU General Public License as published by the 0011 % Free Software Foundation, either version 3 of the License, or (at your 0012 % option) any later version. 0013 % 0014 % This program is distributed in the hope that it will be useful, but 0015 % WITHOUT ANY WARRANTY; without even the implied warranty of 0016 % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 0017 % General Public License for more details. 0018 % 0019 % You should have received a copy of the GNU General Public License along 0020 % with this program. If not, see <http://www.gnu.org/licenses/>. 0021 0022 if isstruct(data) && isfield(data, 'A') 0023 %% hrep 0024 removeFilter = data.h == max(data.h); 0025 data.A(removeFilter, :) = []; 0026 data.h(removeFilter) = []; 0027 elseif isstruct(data) && isfield(data, 'V') 0028 %% vrep 0029 distanceVector = sum(data.V.^2, 2); 0030 removeFilter = distanceVector < zerotol*max(distanceVector); 0031 data.V(removeFilter, :) = []; 0032 elseif isnumeric(data) 0033 %% rayMatrix (vrep without structure) 0034 distanceVector = sum(data.^2, 2); 0035 removeFilter = distanceVector < zerotol*max(distanceVector); 0036 data(removeFilter, :) = []; 0037 else 0038 error('elk:decompoisition:internalError', 'This should not happen'); 0039 end 0040