1function testRB1
2%testRB1: test the RBio toolbox.
3%
4% Example:
5%   testRB1
6%
7% See also UFget, RBread, RBreade, testRB2.
8
9% Copyright 2006, Timothy A. Davis
10
11files = {
12'bcsstk01.rb'
13'farm.rb'
14'lap_25.pse'
15'lap_25.rb'
16'west0479.rb'
17'west0479.rua' } ;
18
19mtypes = {
20'rsa',
21'ira',
22'psa',
23'psa',
24'rua'
25'rua' } ;
26
27for k = 1:length(files)
28    file = files {k} ;
29    % fprintf ('%s : ', file) ;
30    if (file (end) == 'e')
31	[A Z] = RBreade (file) ;
32    else
33	[A Z] = RBread (file) ;
34    end
35    mtype = RBtype (A) ;
36    if (any (mtype ~= mtypes {k}))
37        fprintf ('test failed: %s %s type differs\n', mtype, mtypes {k}) ;
38    end
39    RBwrite ('temp.rb', A, Z) ;
40    [A2 Z2] = RBread ('temp.rb') ;
41    if (~isequal (A, A2))
42	fprintf ('test failed: %s (A differs)\n', file) ;
43	error ('!') ;
44    end
45    if (~isequal (Z, Z2))
46	fprintf ('test failed: %s (Z differs)\n', file) ;
47	error ('!') ;
48    end
49end
50
51load west0479
52C = west0479 ;
53RBwrite ('temp.rb', C, 'WEST0479 chemical eng. problem', 'west0479') ;
54A = RBread ('temp.rb') ;
55if (~isequal (A, C))
56    error ('test failed: west0479 (MATLAB version)') ;
57end
58if (any (mtype ~= 'rua'))
59    fprintf ('test failed: %s %s type differs\n', mtype, 'rua') ;
60end
61
62if (~strcmp (RBtype (A), 'rua'))
63    error ('test failed: RBtype(A)\n') ;
64end
65if (~strcmp (RBtype (spones (A)), 'pua'))
66    error ('test failed: RBtype(spones(A))\n') ;
67end
68if (~strcmp (RBtype (2*spones (A)), 'iua'))
69    error ('test failed: RBtype(2*spones(A))\n') ;
70end
71C = A+A' ;
72if (~strcmp (RBtype (C), 'rsa'))
73    error ('test failed: RBtype(A+A'')\n') ;
74end
75
76delete ('temp.rb') ;
77fprintf ('RB test 1: passed\n') ;
78