219 lines
7.4 KiB
Matlab
219 lines
7.4 KiB
Matlab
function failed = trysuperlu(A);
|
||
% TRYSUPERLU : Test the Matlab interface to superlu.
|
||
%
|
||
% failed = trysuperlu;
|
||
% This runs several tests on superlu, using a matrix with the
|
||
% structure of "smallmesh". It returns a list of failed tests.
|
||
%
|
||
% failed = trysuperlu(A);
|
||
% This just times the factorization of A.
|
||
%
|
||
% Copyright (c) 1995 by Xerox Corporation. All rights reserved.
|
||
% HELP COPYRIGHT for complete copyright and licensing notice.
|
||
|
||
disp(' ');
|
||
disp(['Testing SUPERLU on ' time]);
|
||
disp(' ');
|
||
|
||
resetrandoms;
|
||
if nargin == 0
|
||
A = hbo('smallmesh');
|
||
end;
|
||
[n,m] = size(A);
|
||
if n~=m, error('matrix must be square'), end;
|
||
|
||
failed = [];
|
||
ntest = 0;
|
||
tol = 100 * eps * n^2;
|
||
|
||
if nargin == 1 % Just try to factor the given input matrix.
|
||
ntest=ntest+1;
|
||
fprintf('Matrix size: %d by %d with %d nonzeros.\n',n,m,nnz(A));
|
||
r = trytime(A);
|
||
if r > tol, disp('*** FAILED ***'); failed = [failed ntest]; end;
|
||
return;
|
||
end;
|
||
|
||
% With no input argument, try several tests on the small mesh.
|
||
|
||
PivPostA = sprandn(A);
|
||
A = (n+1)*speye(n) - spones(A); % Diagonally dominant
|
||
[t,post] = etree(A'*A);
|
||
|
||
if post == [1:n], disp('note: column etree already in postorder'), end;
|
||
|
||
NoPivPostA = A(post,:);
|
||
NoPivNoPostA = A(post,post);
|
||
PivNoPostA = sprandn(NoPivNoPostA);
|
||
p = randperm(n);
|
||
pinv = zeros(1,n);
|
||
pinv(p) = 1:n;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: Input perm 0, no pivot, no postorder, 4 outputs.\n',ntest);
|
||
r = try4(NoPivNoPostA,0);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: Input perm 0, no pivot, no postorder, 3 outputs.\n',ntest);
|
||
r = try3(NoPivNoPostA,0);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: Input perm 0, no pivot, no postorder, 2 outputs.\n',ntest);
|
||
r = try2(NoPivNoPostA,0);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: Input perm 0, pivot, no postorder, 4 outputs.\n',ntest);
|
||
r = try4(PivNoPostA,0);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: Input perm 0, pivot, no postorder, 3 outputs.\n',ntest);
|
||
r = try3(PivNoPostA,0);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: Input perm 0, pivot, no postorder, 2 outputs.\n',ntest);
|
||
r = try2(PivNoPostA,0);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: Input perm 0, no pivot, postorder, 4 outputs.\n',ntest);
|
||
r = try4(NoPivPostA,0);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: Input perm 0, no pivot, postorder, 3 outputs.\n',ntest);
|
||
r = try3(NoPivPostA,0);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: Input perm 0, no pivot, postorder, 2 outputs.\n',ntest);
|
||
r = try2(NoPivPostA,0);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: Input perm 0, pivot, postorder, 4 outputs.\n',ntest);
|
||
r = try4(PivPostA,0);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: Input perm 0, pivot, postorder, 3 outputs.\n',ntest);
|
||
r = try3(PivPostA,0);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: Input perm 0, pivot, postorder, 2 outputs.\n',ntest);
|
||
r = try2(PivPostA,0);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: Input perm given, no pivot, no postorder, 4 outputs.\n',ntest);
|
||
r = try4(NoPivNoPostA(:,pinv),p);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: Input perm given, no pivot, no postorder, 3 outputs.\n',ntest);
|
||
r = try3(NoPivNoPostA(:,pinv),p);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: Input perm given, no pivot, no postorder, 2 outputs.\n',ntest);
|
||
r = try2(NoPivNoPostA(:,pinv),p);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: Input perm given, pivot, no postorder, 4 outputs.\n',ntest);
|
||
r = try4(PivNoPostA(:,pinv),p);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: Input perm given, pivot, no postorder, 3 outputs.\n',ntest);
|
||
r = try3(PivNoPostA(:,pinv),p);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: Input perm given, pivot, no postorder, 2 outputs.\n',ntest);
|
||
r = try2(PivNoPostA(:,pinv),p);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: Input perm given, no pivot, postorder, 4 outputs.\n',ntest);
|
||
r = try4(NoPivPostA(:,pinv),p);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: Input perm given, no pivot, postorder, 3 outputs.\n',ntest);
|
||
r = try3(NoPivPostA(:,pinv),p);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: Input perm given, no pivot, postorder, 2 outputs.\n',ntest);
|
||
r = try2(NoPivPostA(:,pinv),p);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: Input perm given, pivot, postorder, 4 outputs.\n',ntest);
|
||
r = try4(PivPostA(:,pinv),p);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: Input perm given, pivot, postorder, 3 outputs.\n',ntest);
|
||
r = try3(PivPostA(:,pinv),p);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: Input perm given, pivot, postorder, 2 outputs.\n',ntest);
|
||
r = try2(PivPostA(:,pinv),p);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: No input perm (colamd done internally), pivot, postorder, 4 outputs.\n',ntest);
|
||
r = try4(PivPostA);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: No input perm, pivot, postorder, 3 outputs.\n',ntest);
|
||
r = try3(PivPostA);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: No input perm, pivot, postorder, 2 outputs.\n',ntest);
|
||
r = try3(PivPostA);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: Input perm given as matrix, pivot, postorder, 4 outputs.\n',ntest);
|
||
P = speye(n);
|
||
P = P(p,:);
|
||
r = try4(PivPostA*P,P);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: Empty matrix.\n',ntest);
|
||
[L,U,PROW,PCOL] = superlu([]);
|
||
disp(' ');
|
||
if max([size(L) size(U) size(PROW) size(PCOL)])
|
||
L, U, PROW, PCOL
|
||
disp('*** FAILED ***^G'), disp(' ');
|
||
failed = [failed ntest];
|
||
end;
|
||
|
||
ntest=ntest+1;
|
||
fprintf('Test %d: Timing versus LU on the 3-element airfoil matrix.\n',ntest);
|
||
A = hbo('airfoil2');
|
||
n = max(size(A));
|
||
A = sprandn(A);
|
||
pmmd = colamd(A);
|
||
A = A(:,pmmd);
|
||
r = trytime(A);
|
||
if r > tol, disp('*** FAILED ***'), disp(' '), failed = [failed ntest]; end;
|
||
|
||
nfailed = length(failed);
|
||
if nfailed
|
||
fprintf('\n%d tests failed.\n\n',nfailed);
|
||
else
|
||
fprintf('\nAll tests passed.\n\n',nfailed);
|
||
end;
|