LA05AS(3)

SYNOPSIS
SUBROUTINE LA05AS (A, IND, NZ, IA, N, IP, IW, W, G, U)

NAME
la05as - Subsidiary to SPLP

LIBRARY
SLATEC

TYPE
SINGLE PRECISION (LA05AS-S, LA05AD-D)

AUTHOR
(UNKNOWN)

DESCRIPTION
THIS SUBPROGRAM IS A SLIGHT MODIFICATION OF A SUBPROGRAM FROM THE C. 1979 AERE HARWELL LIBRARY. THE NAME OF THE CORRESPONDING HARWELL CODE CAN BE OBTAINED BY DELETING THE FINAL LETTER =S= IN THE NAMES USED HERE. REVISIONS MADE BY R J HANSON, SNLA, AUGUST, 1979. REVISED SEP. 13, 1979. ROYALTIES HAVE BEEN PAID TO AERE-UK FOR USE OF THEIR CODES IN THE PACKAGE GIVEN HERE. ANY PRIMARY USAGE OF THE HARWELL SUBROUTINES REQUIRES A ROYALTY AGREEMENT AND PAYMENT BETWEEN THE USER AND AERE-UK. ANY USAGE OF THE SANDIA WRITTEN CODES SPLP( ) (WHICH USES THE HARWELL SUBROUTINES) IS PERMITTED. IP(I,1),IP(I,2) POINT TO THE START OF ROW/COL I. IW(I,1),IW(I,2) HOLD THE NUMBER OF NON-ZEROS IN ROW/COL I. DURING THE MAIN BODY OF THIS SUBROUTINE THE VECTORS IW(.,3),IW(.,5), IW(.,7) ARE USED TO HOLD DOUBLY LINKED LISTS OF ROWS THAT HAVE NOT BEEN PIVOTAL AND HAVE EQUAL NUMBERS OF NON-ZEROS. IW(.,4),IW(.,6),IW(.,8) HOLD SIMILAR LISTS FOR THE COLUMNS. IW(I,3),IW(I,4) HOLD FIRST ROW/COLUMN TO HAVE I NON-ZEROS OR ZERO IF THERE ARE NONE. IW(I,5), IW(I,6) HOLD ROW/COL NUMBER OF ROW/COL PRIOR TO ROW/COL I IN ITS LIST, OR ZERO IF NONE. IW(I,7), IW(I,8) HOLD ROW/COL NUMBER OF ROW/COL AFTER ROW/COL I IN ITS LIST, OR ZERO IF NONE. FOR ROWS/COLS THAT HAVE BEEN PIVOTAL IW(I,5),IW(I,6) HOLD NEGATION OF POSITION OF ROW/COL I IN THE PIVOTAL ORDERING.

SEE ALSO
SPLP

ROUTINES CALLED
LA05ES, MC20AS, R1MACH, XERMSG, XSETUN

COMMON BLOCKS
LA05DS

REVISION HISTORY (YYMMDD)
811215 DATE WRITTEN
890531 Changed all specific intrinsics to generic. (WRB)
890605 Corrected references to XERRWV. (WRB)
890831 Modified array declarations. (WRB)
891214 Prologue converted to Version 4.0 format. (BAB)
900315 CALLs to XERROR changed to CALLs to XERMSG. (THJ)
900402 Added TYPE section. (WRB)
900510 Convert XERRWV calls to XERMSG calls. (RWC)