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