AX + XB = C paper 2

AX + XB = C paper 2 - SU BR O U TIN E PIVO T ( N ) C...

Info iconThis preview shows pages 1–2. Sign up to view the full content.

View Full Document Right Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: SU BR O U TIN E PIVO T ( N ) C PURPOSE - TO PERFORM tHE P IV O I OPEIxATION BY U PD ATIN G ]H E C IN VER SE OF THE BASIS AND 0 VECTOR. C COMMON AM, O*LI , B, NLI , NL2, A*NEI sNE2, I ~JMBASI S, ~, Z D IM EN SIO N AM C SO *SO ), O C bO ), BC50*50) , AC bO ) D IM EN SIO N W (50)J Z ( S O ) * M BASIS(IO O ) DO I I = I , N I B ( I R , I ) = B ( I R * I ) / A ( I R ) Q ( I R ) = Q ( I R ) / A ( I R ) D0 3 I = I ~ N I F { I . E O . I R } GO 7 0 3 Q ( 1 ) : Q ( 1 ) - O ( I R ) * A ( 1 ) DO 2 J = I * N B ( I ~ J ) = B ( I , J ) - B ( I R * J ) * A ( 1 ) 0 C O N T I N U E 3 C O N T IN U E C U PD ATE THE IN D IC ATO R VECTOR OF BASIC VAR IABLES N L I = M B A S I S ( I K ) L = N + IR N L 2 = M B A S IS (L ) M B A S I S I I R ) = N E I M B A S IS (L )= N E 2 LI f LI I RETURN END SU BR O U TIN E PPR IN T ( N ) C PURPOSE - TO PR IN T THE CURRENT SO LU TIO N TO COMPLEMENTARY C PROBLEM AND THE ITER ATIO N NUM BER. O COMMON A M * O n L I , B * N L I * N L 2 , A s N E I * N E 2 * I I ( ~ M B A S I S * ~ * Z D IM EN SIO N AM (SO ~bO )s O (SO )* B ( 5 0 , 5 O ) * A ( 5 0 ) D IM EN SIO N ~ ( 5 0 } * Z ( 5 0 ) * N B A S IS ( IO O ) ~RI TE( 6*I ) LI 1 FORMAT ( 10X*I 3RI TERATI ON N O . ~ I 4 ) I = N + I d = l 2 K I = M B A S I S ( 1 ) KO =M BASIS(J) I F ( O ( J ) * G E - O . O ) GO TO 3 O ( J ) = O . O 3 I F ( R 2 . E O . l ) GO TO 5 WRI TE( 6~4) KI , O( J) A FORMAT ( I O X t ~ H Z ( J I 4 , 2 H ) = J F I 5 5 ) GO TO 7 5 WRI TE( 6*6' KI *Q( J) 6 FORMAT ( I O X t 2 H W ( ~ I 4 ~ O H ) = ~ F I S o S ) 7 I =I +l J=J+l I F ( J * L E . N ) GO TO 2 RETURN END Editor's note: Algorithm 432 described here is available on magnetic tape from the Department of Computer Science, University of Colorado, Boulder, CO 80302. The cost for the tape is $16.00 (U.S. and Canada) or $18.00 (elsewhere). If the user sends a small tape (wt. less than 1 lb.) the algorithm will be copied on it and returned to him at a charge of $10.O0 ( U.S. only). All orders are to be prepaid with checks payable to ACM Algorithms. The algorithm is re corded as one file of BCD 80 character card images at 556 B.P.I., even parity, on seven ~rack tape. We will supply the algorithm at a density of 800 B.P.I. if requested. The cards for the algorithm are sequenced starting at 10 and incremented by 10. The sequence number is right justified in colums 80. Although we will make every attempt to insure that the algorithm conforms to the description printed here, we cannot guarantee it, nor can we guarantee that the algorithm is correct.--L.D.F. Algorithm 432 Solution of the Matrix Equation AX + XB = C [F4] R.H. Bartels and G.W. Stewart [Recd. 21 Oct. 1970 and 7 March 1971] Center for Numerical Analysis, The University of Texas at Austin, Austin, TX 78712 Key Words and Phrases: linear algebra, matrices, linear equa- tions CR Categories: 5.14 Language: Fortran Descdption The following programs are a collection of Fortran IV sub- routines to solve the matrix equation AX -.}- XB = C (1) where A, B, and C are real matrices of dimensions m X m, n X n,...
View Full Document

This note was uploaded on 01/14/2011 for the course ECE 222 taught by Professor Qgsdxjhf during the Fall '07 term at UCSC.

Page1 / 7

AX + XB = C paper 2 - SU BR O U TIN E PIVO T ( N ) C...

This preview shows document pages 1 - 2. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online