function test1 %TEST1 test cs_transpose % % Example: % test1 % See also: testall % Copyright 2006, Timothy A. Davis. % http://www.cise.ufl.edu/research/sparse index = UFget ; [ignore f] = sort (max (index.nrows, index.ncols)) ; f = f (1:100) ; for ii = f Prob = UFget (ii) ; disp (Prob) ; A = Prob.A ; if (~isreal (A)) continue end B = A' ; C = cs_transpose (A) ; if (nnz (B-C) ~= 0) error ('!') end [m n] = size (A) ; % if (m == n) x = rand (n,1) ; y = rand (m,1) ; z = y+A*x ; q = cs_gaxpy (A,x,y) ; err = norm (z-q,1) / norm (z,1) ; disp (err) ; if (err > 1e-14) error ('!') end % end [i j x] = find (A) ; p = randperm (length (i)) ; i = i (p) ; j = j (p) ; x = x (p) ; D = sparse (i,j,x) ; E = cs_sparse (i,j,x) ; % [i j x] F = cs_sparse2 (i,j,x) ; if (nnz (D-E) ~= 0) error ('!') end if (nnz (F-E) ~= 0) error ('!') end clear A B C D E F % pause end