FULL Convert TTeMPS_op operator to full array X = FULL(A) converts the TTeMPS_op operator A to an (A.order)-dimensional full array. Use with care! Result can easily exceed available memory.
0001 function Afull = full( A ) 0002 %FULL Convert TTeMPS_op operator to full array 0003 % X = FULL(A) converts the TTeMPS_op operator A to an (A.order)-dimensional full array. 0004 % 0005 % Use with care! Result can easily exceed available memory. 0006 % 0007 0008 % TTeMPS Toolbox. 0009 % Michael Steinlechner, 2013-2016 0010 % Questions and contact: michael.steinlechner@epfl.ch 0011 % BSD 2-clause license, see LICENSE.txt 0012 0013 d = A.order; 0014 X = full( TTeMPS_op_to_TTeMPS(A) ); 0015 permutevec = [1:d; d+1:2*d]; 0016 permutevec = permutevec(:)'; 0017 0018 sizes = [A.size_col, A.size_row]; 0019 sizes = sizes( permutevec ); 0020 X = reshape( X, sizes(permutevec) ); 0021 X = ipermute( X, permutevec ); 0022 Afull = reshape( X, prod(A.size_col), prod(A.size_row) ); 0023 end