File: GB_spec_resize.m

package info (click to toggle)
suitesparse-graphblas 7.4.0%2Bdfsg-2
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 67,112 kB
  • sloc: ansic: 1,072,243; cpp: 8,081; sh: 512; makefile: 506; asm: 369; python: 125; awk: 10
file content (44 lines) | stat: -rw-r--r-- 1,538 bytes parent folder | download | duplicates (3)
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 mimic of GxB_resize
%
% Usage:
% C = GB_spec_resize (A, nrows_new, ncols_new)

% SuiteSparse:GraphBLAS, Timothy A. Davis, (c) 2017-2022, All Rights Reserved.
% SPDX-License-Identifier: Apache-2.0

%-------------------------------------------------------------------------------
% 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 *.m 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