summaryrefslogtreecommitdiff
path: root/modules/sparse/macros/sp2adj.sci
blob: da86ffdbad8267d208c8548f873cec63e5e5848a (plain)
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
45
46
47
48
49
50
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) XXXX-2008 - INRIA
// Copyright (C) 2010 - DIGITEO - Michael Baudin
//
// This file must be used under the terms of the CeCILL.
// This source file is licensed as described in the file COPYING, which
// you should have received as part of this distribution.  The terms
// are also available at
// http://www.cecill.info/licences/Licence_CeCILL_V2.1-en.txt


function [xadj,iadj,v]=sp2adj(A)
    // sparse conversion to adjacency
    // See sp2adj.xml for help.
    //
    // Check number of arguments
    [lhs, rhs] = argn()
    if ( rhs <> 1 ) then
        lstr = gettext("%s: Wrong number of input arguments: %d to %d expected.\n")
        error ( msprintf(lstr,"sp2adj",1,1))
    end
    if ( and(lhs <> [0 1 2 3]) ) then
        lstr = gettext("%s: Wrong number of output arguments: %d to %d expected.\n")
        error ( msprintf(lstr,"sp2adj",0,3))
    end
    //
    // Check type of arguments
    if ( typeof(A) <> "sparse" ) then
        lstr = gettext("%s: Wrong type for input argument #%d.\n")
        error ( msprintf(lstr,"sp2adj",1))
    end
    //
    // Check size of arguments
    // Nothing to do
    //
    // Check content of arguments
    // Nothing to do
    //
    [ij,v,n]=spget(A')
    N=n(1)
    if ( ij == [] ) then,
        xadj=ones(n(2)+1,1)
        iadj=[]
        v=[]
    else,
        [xadj,la,iadj]=ta2lpd(ij(:,1)',ij(:,2)',N+1,N)
        xadj=xadj(:)
        iadj=iadj(:)
    end
endfunction