dlaqz1.f - Man Page

subroutine dlaqz1 (a, lda, b, ldb, sr1, sr2, si, beta1, beta2, v)
DLAQZ1

DLAQZ1  

Purpose:

      Given a 3-by-3 matrix pencil (A,B), DLAQZ1 sets v to a
      scalar multiple of the first column of the product

      (*)  K = (A - (beta2*sr2 - i*si)*B)*B^(-1)*(beta1*A - (sr2 + i*si2)*B)*B^(-1).

      It is assumed that either

              1) sr1 = sr2
          or
              2) si = 0.

      This is useful for starting double implicit shift bulges
      in the QZ algorithm.

Parameters

A

          A is DOUBLE PRECISION array, dimension (LDA,N)
              The 3-by-3 matrix A in (*).

LDA

          LDA is INTEGER
              The leading dimension of A as declared in
              the calling procedure.

B

          B is DOUBLE PRECISION array, dimension (LDB,N)
              The 3-by-3 matrix B in (*).

LDB

          LDB is INTEGER
              The leading dimension of B as declared in
              the calling procedure.

SR1

          SR1 is DOUBLE PRECISION

SR2

          SR2 is DOUBLE PRECISION

SI

          SI is DOUBLE PRECISION

BETA1

          BETA1 is DOUBLE PRECISION

BETA2

          BETA2 is DOUBLE PRECISION

V

          V is DOUBLE PRECISION array, dimension (N)
              A scalar multiple of the first column of the
              matrix K in (*).

Author

Thijs Steel, KU Leuven

Date

May 2020

Definition at line 125 of file dlaqz1.f.