
index = UFget ;
[ignore f] = sort (max (index.nrows, index.ncols)) ;
f = f (1:100) ;

for i = f
    Prob = UFget (i)
    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)
	if (err > 1e-14)
	    err
	    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_triplet (i,j,x) ;
    % [i j x]
    F = cs_triplet2 (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

