.TH "TESTING/EIG/sort01.f" 3 "Version 3.12.0" "LAPACK" \" -*- nroff -*- .ad l .nh .SH NAME TESTING/EIG/sort01.f .SH SYNOPSIS .br .PP .SS "Functions/Subroutines" .in +1c .ti -1c .RI "subroutine \fBsort01\fP (rowcol, m, n, u, ldu, work, lwork, resid)" .br .RI "\fBSORT01\fP " .in -1c .SH "Function/Subroutine Documentation" .PP .SS "subroutine sort01 (character rowcol, integer m, integer n, real, dimension( ldu, * ) u, integer ldu, real, dimension( * ) work, integer lwork, real resid)" .PP \fBSORT01\fP .PP \fBPurpose:\fP .RS 4 .PP .nf !> !> SORT01 checks that the matrix U is orthogonal by computing the ratio !> !> RESID = norm( I - U*U' ) / ( n * EPS ), if ROWCOL = 'R', !> or !> RESID = norm( I - U'*U ) / ( m * EPS ), if ROWCOL = 'C'\&. !> !> Alternatively, if there isn't sufficient workspace to form !> I - U*U' or I - U'*U, the ratio is computed as !> !> RESID = abs( I - U*U' ) / ( n * EPS ), if ROWCOL = 'R', !> or !> RESID = abs( I - U'*U ) / ( m * EPS ), if ROWCOL = 'C'\&. !> !> where EPS is the machine precision\&. ROWCOL is used only if m = n; !> if m > n, ROWCOL is assumed to be 'C', and if m < n, ROWCOL is !> assumed to be 'R'\&. !> .fi .PP .RE .PP \fBParameters\fP .RS 4 \fIROWCOL\fP .PP .nf !> ROWCOL is CHARACTER !> Specifies whether the rows or columns of U should be checked !> for orthogonality\&. Used only if M = N\&. !> = 'R': Check for orthogonal rows of U !> = 'C': Check for orthogonal columns of U !> .fi .PP .br \fIM\fP .PP .nf !> M is INTEGER !> The number of rows of the matrix U\&. !> .fi .PP .br \fIN\fP .PP .nf !> N is INTEGER !> The number of columns of the matrix U\&. !> .fi .PP .br \fIU\fP .PP .nf !> U is REAL array, dimension (LDU,N) !> The orthogonal matrix U\&. U is checked for orthogonal columns !> if m > n or if m = n and ROWCOL = 'C'\&. U is checked for !> orthogonal rows if m < n or if m = n and ROWCOL = 'R'\&. !> .fi .PP .br \fILDU\fP .PP .nf !> LDU is INTEGER !> The leading dimension of the array U\&. LDU >= max(1,M)\&. !> .fi .PP .br \fIWORK\fP .PP .nf !> WORK is REAL array, dimension (LWORK) !> .fi .PP .br \fILWORK\fP .PP .nf !> LWORK is INTEGER !> The length of the array WORK\&. For best performance, LWORK !> should be at least N*(N+1) if ROWCOL = 'C' or M*(M+1) if !> ROWCOL = 'R', but the test will be done even if LWORK is 0\&. !> .fi .PP .br \fIRESID\fP .PP .nf !> RESID is REAL !> RESID = norm( I - U * U' ) / ( n * EPS ), if ROWCOL = 'R', or !> RESID = norm( I - U' * U ) / ( m * EPS ), if ROWCOL = 'C'\&. !> .fi .PP .RE .PP \fBAuthor\fP .RS 4 Univ\&. of Tennessee .PP Univ\&. of California Berkeley .PP Univ\&. of Colorado Denver .PP NAG Ltd\&. .RE .PP .PP Definition at line \fB115\fP of file \fBsort01\&.f\fP\&. .SH "Author" .PP Generated automatically by Doxygen for LAPACK from the source code\&.