File: GB_spec_resize.m

package info (click to toggle)
suitesparse 1%3A5.8.1%2Bdfsg-2
  • links: PTS, VCS
  • area: main
  • in suites: bullseye
  • size: 152,716 kB
  • sloc: ansic: 774,385; cpp: 24,213; makefile: 6,310; fortran: 1,927; java: 1,826; csh: 1,686; ruby: 725; sh: 535; perl: 225; python: 209; sed: 164; awk: 60
file content (44 lines) | stat: -rw-r--r-- 1,580 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
function C = GB_spec_resize (A, nrows_new, ncols_new)
%GB_SPEC_RESIZE a MATLAB mimic of GxB_resize
%
% Usage:
% C = GB_spec_resize (A, nrows_new, ncols_new)

% SuiteSparse:GraphBLAS, Timothy A. Davis, (c) 2017-2020, All Rights Reserved.
% http://suitesparse.com   See GraphBLAS/Doc/License.txt for license.

%-------------------------------------------------------------------------------
% get inputs
%-------------------------------------------------------------------------------

if (nargout > 1 || nargin ~= 3)
    error ('usage: C = GB_spec_resize (A, nrows_new, ncols_new)') ;
end

C = GB_spec_matrix (A) ;
c_class = C.class ;

%-------------------------------------------------------------------------------
% do the work via a clean MATLAB interpretation of the entire GraphBLAS spec
%-------------------------------------------------------------------------------

[nrows_old ncols_old] = size (C.matrix) ;

if (ncols_new < ncols_old)
    C.matrix  = C.matrix  (:, 1:ncols_new) ;
    C.pattern = C.pattern (:, 1:ncols_new) ;
elseif (ncols_new > ncols_old)
    z = GB_spec_zeros ([nrows_old, ncols_new - ncols_old], c_class) ;
    C.matrix  = [C.matrix  z ] ;
    C.pattern = [C.pattern (false (nrows_old, ncols_new - ncols_old))] ;
end

if (nrows_new < nrows_old)
    C.matrix  = C.matrix  (1:nrows_new, :) ;
    C.pattern = C.pattern (1:nrows_new, :) ;
elseif (nrows_new > nrows_old)
    z = GB_spec_zeros ([nrows_new - nrows_old, ncols_new], c_class) ;
    C.matrix  = [C.matrix  ; z ] ;
    C.pattern = [C.pattern ; (false (nrows_new - nrows_old, ncols_new))] ;
end