computeUpVector

 THIS IS NO USER FUNCTION

• computePartitionIndexDec determine the partition of a given hC-vector
• mapToValidityCone map a hC vector to the cone of valid polytopes

0001 function upVector = computeUpVector(pbeDef, hcMatrix, zerotol)
0002 %
0003 %
0004 % THIS IS NO USER FUNCTION
0005 
0006 % The elk-library: convex geometry applied to crystallization modeling.
0007 %   Copyright (C) 2012 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 if pbeDef.dec.isProper
0022     upVector = computePartitionIndexDec(pbeDef.dec, hcMatrix, zerotol, 1);
0023 elseif pbeDef.dec.isComplete && ~pbeDef.dec.isProper
0024     if pbeDef.optionStruct.useEmbeddingPartitions
0025         hcMatrix = mapToValidityCone(...
0026             pbeDef.dec.properData.mappingNewToProper*hcMatrix, ...
0027             pbeDef.dec.properData, zerotol);
0028         upVector = computePartitionIndexDec(pbeDef.dec.properData, ...
0029             hcMatrix, zerotol, 1);
0030     else
0031         upVector = computePartitionIndexDec(pbeDef.dec, hcMatrix, zerotol, 1);
0032     end
0033 elseif ~pbeDef.dec.isComplete
0034     upVector(1, end + (1:nNew)) = 1;
0035 else
0036     error('elk:pbeSolver:inputError', ['The given representation is neither ' ...
0037         'proper, nor improper but complete, nor incomplete. This does not ' ...
0038         'makes sense. Please check your decomposition, dude.']);
0039 end