clatm5.c

Go to the documentation of this file.

/* clatm5.f -- translated by f2c (version 20061008).
   You must link the resulting object file with libf2c:
        on Microsoft Windows system, link with libf2c.lib;
        on Linux or Unix systems, link with .../path/to/libf2c.a -lm
        or, if you install libf2c.a in a standard place, with -lf2c -lm
        -- in that order, at the end of the command line, as in
                cc *.o -lf2c -lm
        Source for libf2c is in /netlib/f2c/libf2c.zip, e.g.,

                http://www.netlib.org/f2c/libf2c.zip
*/

#include "f2c.h"
#include "blaswrap.h"

/* Table of constant values */

static complex c_b1 = {1.f,0.f};
static complex c_b3 = {0.f,0.f};
static complex c_b5 = {20.f,0.f};

/* Subroutine */ int clatm5_(integer *prtype, integer *m, integer *n, complex 
        *a, integer *lda, complex *b, integer *ldb, complex *c__, integer *
        ldc, complex *d__, integer *ldd, complex *e, integer *lde, complex *f, 
         integer *ldf, complex *r__, integer *ldr, complex *l, integer *ldl, 
        real *alpha, integer *qblcka, integer *qblckb)
{
    /* System generated locals */
    integer a_dim1, a_offset, b_dim1, b_offset, c_dim1, c_offset, d_dim1, 
            d_offset, e_dim1, e_offset, f_dim1, f_offset, l_dim1, l_offset, 
            r_dim1, r_offset, i__1, i__2, i__3, i__4;
    doublereal d__1;
    complex q__1, q__2, q__3, q__4, q__5;

    /* Builtin functions */
    void c_sin(complex *, complex *), c_div(complex *, complex *, complex *);

    /* Local variables */
    integer i__, j, k;
    extern /* Subroutine */ int cgemm_(char *, char *, integer *, integer *, 
            integer *, complex *, complex *, integer *, complex *, integer *, 
            complex *, complex *, integer *);
    complex imeps, reeps;


/*  -- LAPACK test routine (version 3.1) -- */
/*     Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. */
/*     November 2006 */

/*     .. Scalar Arguments .. */
/*     .. */
/*     .. Array Arguments .. */
/*     .. */

/*  Purpose */
/*  ======= */

/*  CLATM5 generates matrices involved in the Generalized Sylvester */
/*  equation: */

/*      A * R - L * B = C */
/*      D * R - L * E = F */

/*  They also satisfy (the diagonalization condition) */

/*   [ I -L ] ( [ A  -C ], [ D -F ] ) [ I  R ] = ( [ A    ], [ D    ] ) */
/*   [    I ] ( [     B ]  [    E ] ) [    I ]   ( [    B ]  [    E ] ) */


/*  Arguments */
/*  ========= */

/*  PRTYPE  (input) INTEGER */
/*          "Points" to a certian type of the matrices to generate */
/*          (see futher details). */

/*  M       (input) INTEGER */
/*          Specifies the order of A and D and the number of rows in */
/*          C, F,  R and L. */

/*  N       (input) INTEGER */
/*          Specifies the order of B and E and the number of columns in */
/*          C, F, R and L. */

/*  A       (output) COMPLEX array, dimension (LDA, M). */
/*          On exit A M-by-M is initialized according to PRTYPE. */

/*  LDA     (input) INTEGER */
/*          The leading dimension of A. */

/*  B       (output) COMPLEX array, dimension (LDB, N). */
/*          On exit B N-by-N is initialized according to PRTYPE. */

/*  LDB     (input) INTEGER */
/*          The leading dimension of B. */

/*  C       (output) COMPLEX array, dimension (LDC, N). */
/*          On exit C M-by-N is initialized according to PRTYPE. */

/*  LDC     (input) INTEGER */
/*          The leading dimension of C. */

/*  D       (output) COMPLEX array, dimension (LDD, M). */
/*          On exit D M-by-M is initialized according to PRTYPE. */

/*  LDD     (input) INTEGER */
/*          The leading dimension of D. */

/*  E       (output) COMPLEX array, dimension (LDE, N). */
/*          On exit E N-by-N is initialized according to PRTYPE. */

/*  LDE     (input) INTEGER */
/*          The leading dimension of E. */

/*  F       (output) COMPLEX array, dimension (LDF, N). */
/*          On exit F M-by-N is initialized according to PRTYPE. */

/*  LDF     (input) INTEGER */
/*          The leading dimension of F. */

/*  R       (output) COMPLEX array, dimension (LDR, N). */
/*          On exit R M-by-N is initialized according to PRTYPE. */

/*  LDR     (input) INTEGER */
/*          The leading dimension of R. */

/*  L       (output) COMPLEX array, dimension (LDL, N). */
/*          On exit L M-by-N is initialized according to PRTYPE. */

/*  LDL     (input) INTEGER */
/*          The leading dimension of L. */

/*  ALPHA   (input) REAL */
/*          Parameter used in generating PRTYPE = 1 and 5 matrices. */

/*  QBLCKA  (input) INTEGER */
/*          When PRTYPE = 3, specifies the distance between 2-by-2 */
/*          blocks on the diagonal in A. Otherwise, QBLCKA is not */
/*          referenced. QBLCKA > 1. */

/*  QBLCKB  (input) INTEGER */
/*          When PRTYPE = 3, specifies the distance between 2-by-2 */
/*          blocks on the diagonal in B. Otherwise, QBLCKB is not */
/*          referenced. QBLCKB > 1. */


/*  Further Details */
/*  =============== */

/*  PRTYPE = 1: A and B are Jordan blocks, D and E are identity matrices */

/*             A : if (i == j) then A(i, j) = 1.0 */
/*                 if (j == i + 1) then A(i, j) = -1.0 */
/*                 else A(i, j) = 0.0,            i, j = 1...M */

/*             B : if (i == j) then B(i, j) = 1.0 - ALPHA */
/*                 if (j == i + 1) then B(i, j) = 1.0 */
/*                 else B(i, j) = 0.0,            i, j = 1...N */

/*             D : if (i == j) then D(i, j) = 1.0 */
/*                 else D(i, j) = 0.0,            i, j = 1...M */

/*             E : if (i == j) then E(i, j) = 1.0 */
/*                 else E(i, j) = 0.0,            i, j = 1...N */

/*             L =  R are chosen from [-10...10], */
/*                  which specifies the right hand sides (C, F). */

/*  PRTYPE = 2 or 3: Triangular and/or quasi- triangular. */

/*             A : if (i <= j) then A(i, j) = [-1...1] */
/*                 else A(i, j) = 0.0,             i, j = 1...M */

/*                 if (PRTYPE = 3) then */
/*                    A(k + 1, k + 1) = A(k, k) */
/*                    A(k + 1, k) = [-1...1] */
/*                    sign(A(k, k + 1) = -(sin(A(k + 1, k)) */
/*                        k = 1, M - 1, QBLCKA */

/*             B : if (i <= j) then B(i, j) = [-1...1] */
/*                 else B(i, j) = 0.0,            i, j = 1...N */

/*                 if (PRTYPE = 3) then */
/*                    B(k + 1, k + 1) = B(k, k) */
/*                    B(k + 1, k) = [-1...1] */
/*                    sign(B(k, k + 1) = -(sign(B(k + 1, k)) */
/*                        k = 1, N - 1, QBLCKB */

/*             D : if (i <= j) then D(i, j) = [-1...1]. */
/*                 else D(i, j) = 0.0,            i, j = 1...M */


/*             E : if (i <= j) then D(i, j) = [-1...1] */
/*                 else E(i, j) = 0.0,            i, j = 1...N */

/*                 L, R are chosen from [-10...10], */
/*                 which specifies the right hand sides (C, F). */

/*  PRTYPE = 4 Full */
/*             A(i, j) = [-10...10] */
/*             D(i, j) = [-1...1]    i,j = 1...M */
/*             B(i, j) = [-10...10] */
/*             E(i, j) = [-1...1]    i,j = 1...N */
/*             R(i, j) = [-10...10] */
/*             L(i, j) = [-1...1]    i = 1..M ,j = 1...N */

/*             L, R specifies the right hand sides (C, F). */

/*  PRTYPE = 5 special case common and/or close eigs. */

/*  ===================================================================== */

/*     .. Parameters .. */
/*     .. */
/*     .. Local Scalars .. */
/*     .. */
/*     .. Intrinsic Functions .. */
/*     .. */
/*     .. External Subroutines .. */
/*     .. */
/*     .. Executable Statements .. */

    /* Parameter adjustments */
    a_dim1 = *lda;
    a_offset = 1 + a_dim1;
    a -= a_offset;
    b_dim1 = *ldb;
    b_offset = 1 + b_dim1;
    b -= b_offset;
    c_dim1 = *ldc;
    c_offset = 1 + c_dim1;
    c__ -= c_offset;
    d_dim1 = *ldd;
    d_offset = 1 + d_dim1;
    d__ -= d_offset;
    e_dim1 = *lde;
    e_offset = 1 + e_dim1;
    e -= e_offset;
    f_dim1 = *ldf;
    f_offset = 1 + f_dim1;
    f -= f_offset;
    r_dim1 = *ldr;
    r_offset = 1 + r_dim1;
    r__ -= r_offset;
    l_dim1 = *ldl;
    l_offset = 1 + l_dim1;
    l -= l_offset;

    /* Function Body */
    if (*prtype == 1) {
        i__1 = *m;
        for (i__ = 1; i__ <= i__1; ++i__) {
            i__2 = *m;
            for (j = 1; j <= i__2; ++j) {
                if (i__ == j) {
                    i__3 = i__ + j * a_dim1;
                    a[i__3].r = 1.f, a[i__3].i = 0.f;
                    i__3 = i__ + j * d_dim1;
                    d__[i__3].r = 1.f, d__[i__3].i = 0.f;
                } else if (i__ == j - 1) {
                    i__3 = i__ + j * a_dim1;
                    q__1.r = -1.f, q__1.i = -0.f;
                    a[i__3].r = q__1.r, a[i__3].i = q__1.i;
                    i__3 = i__ + j * d_dim1;
                    d__[i__3].r = 0.f, d__[i__3].i = 0.f;
                } else {
                    i__3 = i__ + j * a_dim1;
                    a[i__3].r = 0.f, a[i__3].i = 0.f;
                    i__3 = i__ + j * d_dim1;
                    d__[i__3].r = 0.f, d__[i__3].i = 0.f;
                }
/* L10: */
            }
/* L20: */
        }

        i__1 = *n;
        for (i__ = 1; i__ <= i__1; ++i__) {
            i__2 = *n;
            for (j = 1; j <= i__2; ++j) {
                if (i__ == j) {
                    i__3 = i__ + j * b_dim1;
                    q__1.r = 1.f - *alpha, q__1.i = 0.f;
                    b[i__3].r = q__1.r, b[i__3].i = q__1.i;
                    i__3 = i__ + j * e_dim1;
                    e[i__3].r = 1.f, e[i__3].i = 0.f;
                } else if (i__ == j - 1) {
                    i__3 = i__ + j * b_dim1;
                    b[i__3].r = 1.f, b[i__3].i = 0.f;
                    i__3 = i__ + j * e_dim1;
                    e[i__3].r = 0.f, e[i__3].i = 0.f;
                } else {
                    i__3 = i__ + j * b_dim1;
                    b[i__3].r = 0.f, b[i__3].i = 0.f;
                    i__3 = i__ + j * e_dim1;
                    e[i__3].r = 0.f, e[i__3].i = 0.f;
                }
/* L30: */
            }
/* L40: */
        }

        i__1 = *m;
        for (i__ = 1; i__ <= i__1; ++i__) {
            i__2 = *n;
            for (j = 1; j <= i__2; ++j) {
                i__3 = i__ + j * r_dim1;
                i__4 = i__ / j;
                q__4.r = (real) i__4, q__4.i = 0.f;
                c_sin(&q__3, &q__4);
                q__2.r = .5f - q__3.r, q__2.i = 0.f - q__3.i;
                q__1.r = q__2.r * 20.f - q__2.i * 0.f, q__1.i = q__2.r * 0.f 
                        + q__2.i * 20.f;
                r__[i__3].r = q__1.r, r__[i__3].i = q__1.i;
                i__3 = i__ + j * l_dim1;
                i__4 = i__ + j * r_dim1;
                l[i__3].r = r__[i__4].r, l[i__3].i = r__[i__4].i;
/* L50: */
            }
/* L60: */
        }

    } else if (*prtype == 2 || *prtype == 3) {
        i__1 = *m;
        for (i__ = 1; i__ <= i__1; ++i__) {
            i__2 = *m;
            for (j = 1; j <= i__2; ++j) {
                if (i__ <= j) {
                    i__3 = i__ + j * a_dim1;
                    q__4.r = (real) i__, q__4.i = 0.f;
                    c_sin(&q__3, &q__4);
                    q__2.r = .5f - q__3.r, q__2.i = 0.f - q__3.i;
                    q__1.r = q__2.r * 2.f - q__2.i * 0.f, q__1.i = q__2.r * 
                            0.f + q__2.i * 2.f;
                    a[i__3].r = q__1.r, a[i__3].i = q__1.i;
                    i__3 = i__ + j * d_dim1;
                    i__4 = i__ * j;
                    q__4.r = (real) i__4, q__4.i = 0.f;
                    c_sin(&q__3, &q__4);
                    q__2.r = .5f - q__3.r, q__2.i = 0.f - q__3.i;
                    q__1.r = q__2.r * 2.f - q__2.i * 0.f, q__1.i = q__2.r * 
                            0.f + q__2.i * 2.f;
                    d__[i__3].r = q__1.r, d__[i__3].i = q__1.i;
                } else {
                    i__3 = i__ + j * a_dim1;
                    a[i__3].r = 0.f, a[i__3].i = 0.f;
                    i__3 = i__ + j * d_dim1;
                    d__[i__3].r = 0.f, d__[i__3].i = 0.f;
                }
/* L70: */
            }
/* L80: */
        }

        i__1 = *n;
        for (i__ = 1; i__ <= i__1; ++i__) {
            i__2 = *n;
            for (j = 1; j <= i__2; ++j) {
                if (i__ <= j) {
                    i__3 = i__ + j * b_dim1;
                    i__4 = i__ + j;
                    q__4.r = (real) i__4, q__4.i = 0.f;
                    c_sin(&q__3, &q__4);
                    q__2.r = .5f - q__3.r, q__2.i = 0.f - q__3.i;
                    q__1.r = q__2.r * 2.f - q__2.i * 0.f, q__1.i = q__2.r * 
                            0.f + q__2.i * 2.f;
                    b[i__3].r = q__1.r, b[i__3].i = q__1.i;
                    i__3 = i__ + j * e_dim1;
                    q__4.r = (real) j, q__4.i = 0.f;
                    c_sin(&q__3, &q__4);
                    q__2.r = .5f - q__3.r, q__2.i = 0.f - q__3.i;
                    q__1.r = q__2.r * 2.f - q__2.i * 0.f, q__1.i = q__2.r * 
                            0.f + q__2.i * 2.f;
                    e[i__3].r = q__1.r, e[i__3].i = q__1.i;
                } else {
                    i__3 = i__ + j * b_dim1;
                    b[i__3].r = 0.f, b[i__3].i = 0.f;
                    i__3 = i__ + j * e_dim1;
                    e[i__3].r = 0.f, e[i__3].i = 0.f;
                }
/* L90: */
            }
/* L100: */
        }

        i__1 = *m;
        for (i__ = 1; i__ <= i__1; ++i__) {
            i__2 = *n;
            for (j = 1; j <= i__2; ++j) {
                i__3 = i__ + j * r_dim1;
                i__4 = i__ * j;
                q__4.r = (real) i__4, q__4.i = 0.f;
                c_sin(&q__3, &q__4);
                q__2.r = .5f - q__3.r, q__2.i = 0.f - q__3.i;
                q__1.r = q__2.r * 20.f - q__2.i * 0.f, q__1.i = q__2.r * 0.f 
                        + q__2.i * 20.f;
                r__[i__3].r = q__1.r, r__[i__3].i = q__1.i;
                i__3 = i__ + j * l_dim1;
                i__4 = i__ + j;
                q__4.r = (real) i__4, q__4.i = 0.f;
                c_sin(&q__3, &q__4);
                q__2.r = .5f - q__3.r, q__2.i = 0.f - q__3.i;
                q__1.r = q__2.r * 20.f - q__2.i * 0.f, q__1.i = q__2.r * 0.f 
                        + q__2.i * 20.f;
                l[i__3].r = q__1.r, l[i__3].i = q__1.i;
/* L110: */
            }
/* L120: */
        }

        if (*prtype == 3) {
            if (*qblcka <= 1) {
                *qblcka = 2;
            }
            i__1 = *m - 1;
            i__2 = *qblcka;
            for (k = 1; i__2 < 0 ? k >= i__1 : k <= i__1; k += i__2) {
                i__3 = k + 1 + (k + 1) * a_dim1;
                i__4 = k + k * a_dim1;
                a[i__3].r = a[i__4].r, a[i__3].i = a[i__4].i;
                i__3 = k + 1 + k * a_dim1;
                c_sin(&q__2, &a[k + (k + 1) * a_dim1]);
                q__1.r = -q__2.r, q__1.i = -q__2.i;
                a[i__3].r = q__1.r, a[i__3].i = q__1.i;
/* L130: */
            }

            if (*qblckb <= 1) {
                *qblckb = 2;
            }
            i__2 = *n - 1;
            i__1 = *qblckb;
            for (k = 1; i__1 < 0 ? k >= i__2 : k <= i__2; k += i__1) {
                i__3 = k + 1 + (k + 1) * b_dim1;
                i__4 = k + k * b_dim1;
                b[i__3].r = b[i__4].r, b[i__3].i = b[i__4].i;
                i__3 = k + 1 + k * b_dim1;
                c_sin(&q__2, &b[k + (k + 1) * b_dim1]);
                q__1.r = -q__2.r, q__1.i = -q__2.i;
                b[i__3].r = q__1.r, b[i__3].i = q__1.i;
/* L140: */
            }
        }

    } else if (*prtype == 4) {
        i__1 = *m;
        for (i__ = 1; i__ <= i__1; ++i__) {
            i__2 = *m;
            for (j = 1; j <= i__2; ++j) {
                i__3 = i__ + j * a_dim1;
                i__4 = i__ * j;
                q__4.r = (real) i__4, q__4.i = 0.f;
                c_sin(&q__3, &q__4);
                q__2.r = .5f - q__3.r, q__2.i = 0.f - q__3.i;
                q__1.r = q__2.r * 20.f - q__2.i * 0.f, q__1.i = q__2.r * 0.f 
                        + q__2.i * 20.f;
                a[i__3].r = q__1.r, a[i__3].i = q__1.i;
                i__3 = i__ + j * d_dim1;
                i__4 = i__ + j;
                q__4.r = (real) i__4, q__4.i = 0.f;
                c_sin(&q__3, &q__4);
                q__2.r = .5f - q__3.r, q__2.i = 0.f - q__3.i;
                q__1.r = q__2.r * 2.f - q__2.i * 0.f, q__1.i = q__2.r * 0.f + 
                        q__2.i * 2.f;
                d__[i__3].r = q__1.r, d__[i__3].i = q__1.i;
/* L150: */
            }
/* L160: */
        }

        i__1 = *n;
        for (i__ = 1; i__ <= i__1; ++i__) {
            i__2 = *n;
            for (j = 1; j <= i__2; ++j) {
                i__3 = i__ + j * b_dim1;
                i__4 = i__ + j;
                q__4.r = (real) i__4, q__4.i = 0.f;
                c_sin(&q__3, &q__4);
                q__2.r = .5f - q__3.r, q__2.i = 0.f - q__3.i;
                q__1.r = q__2.r * 20.f - q__2.i * 0.f, q__1.i = q__2.r * 0.f 
                        + q__2.i * 20.f;
                b[i__3].r = q__1.r, b[i__3].i = q__1.i;
                i__3 = i__ + j * e_dim1;
                i__4 = i__ * j;
                q__4.r = (real) i__4, q__4.i = 0.f;
                c_sin(&q__3, &q__4);
                q__2.r = .5f - q__3.r, q__2.i = 0.f - q__3.i;
                q__1.r = q__2.r * 2.f - q__2.i * 0.f, q__1.i = q__2.r * 0.f + 
                        q__2.i * 2.f;
                e[i__3].r = q__1.r, e[i__3].i = q__1.i;
/* L170: */
            }
/* L180: */
        }

        i__1 = *m;
        for (i__ = 1; i__ <= i__1; ++i__) {
            i__2 = *n;
            for (j = 1; j <= i__2; ++j) {
                i__3 = i__ + j * r_dim1;
                i__4 = j / i__;
                q__4.r = (real) i__4, q__4.i = 0.f;
                c_sin(&q__3, &q__4);
                q__2.r = .5f - q__3.r, q__2.i = 0.f - q__3.i;
                q__1.r = q__2.r * 20.f - q__2.i * 0.f, q__1.i = q__2.r * 0.f 
                        + q__2.i * 20.f;
                r__[i__3].r = q__1.r, r__[i__3].i = q__1.i;
                i__3 = i__ + j * l_dim1;
                i__4 = i__ * j;
                q__4.r = (real) i__4, q__4.i = 0.f;
                c_sin(&q__3, &q__4);
                q__2.r = .5f - q__3.r, q__2.i = 0.f - q__3.i;
                q__1.r = q__2.r * 2.f - q__2.i * 0.f, q__1.i = q__2.r * 0.f + 
                        q__2.i * 2.f;
                l[i__3].r = q__1.r, l[i__3].i = q__1.i;
/* L190: */
            }
/* L200: */
        }

    } else if (*prtype >= 5) {
        q__3.r = 1.f, q__3.i = 0.f;
        q__2.r = q__3.r * 20.f - q__3.i * 0.f, q__2.i = q__3.r * 0.f + q__3.i 
                * 20.f;
        q__1.r = q__2.r / *alpha, q__1.i = q__2.i / *alpha;
        reeps.r = q__1.r, reeps.i = q__1.i;
        q__2.r = -1.5f, q__2.i = 0.f;
        q__1.r = q__2.r / *alpha, q__1.i = q__2.i / *alpha;
        imeps.r = q__1.r, imeps.i = q__1.i;
        i__1 = *m;
        for (i__ = 1; i__ <= i__1; ++i__) {
            i__2 = *n;
            for (j = 1; j <= i__2; ++j) {
                i__3 = i__ + j * r_dim1;
                i__4 = i__ * j;
                q__5.r = (real) i__4, q__5.i = 0.f;
                c_sin(&q__4, &q__5);
                q__3.r = .5f - q__4.r, q__3.i = 0.f - q__4.i;
                q__2.r = *alpha * q__3.r, q__2.i = *alpha * q__3.i;
                c_div(&q__1, &q__2, &c_b5);
                r__[i__3].r = q__1.r, r__[i__3].i = q__1.i;
                i__3 = i__ + j * l_dim1;
                i__4 = i__ + j;
                q__5.r = (real) i__4, q__5.i = 0.f;
                c_sin(&q__4, &q__5);
                q__3.r = .5f - q__4.r, q__3.i = 0.f - q__4.i;
                q__2.r = *alpha * q__3.r, q__2.i = *alpha * q__3.i;
                c_div(&q__1, &q__2, &c_b5);
                l[i__3].r = q__1.r, l[i__3].i = q__1.i;
/* L210: */
            }
/* L220: */
        }

        i__1 = *m;
        for (i__ = 1; i__ <= i__1; ++i__) {
            i__2 = i__ + i__ * d_dim1;
            d__[i__2].r = 1.f, d__[i__2].i = 0.f;
/* L230: */
        }

        i__1 = *m;
        for (i__ = 1; i__ <= i__1; ++i__) {
            if (i__ <= 4) {
                i__2 = i__ + i__ * a_dim1;
                a[i__2].r = 1.f, a[i__2].i = 0.f;
                if (i__ > 2) {
                    i__2 = i__ + i__ * a_dim1;
                    q__1.r = reeps.r + 1.f, q__1.i = reeps.i + 0.f;
                    a[i__2].r = q__1.r, a[i__2].i = q__1.i;
                }
                if (i__ % 2 != 0 && i__ < *m) {
                    i__2 = i__ + (i__ + 1) * a_dim1;
                    a[i__2].r = imeps.r, a[i__2].i = imeps.i;
                } else if (i__ > 1) {
                    i__2 = i__ + (i__ - 1) * a_dim1;
                    q__1.r = -imeps.r, q__1.i = -imeps.i;
                    a[i__2].r = q__1.r, a[i__2].i = q__1.i;
                }
            } else if (i__ <= 8) {
                if (i__ <= 6) {
                    i__2 = i__ + i__ * a_dim1;
                    a[i__2].r = reeps.r, a[i__2].i = reeps.i;
                } else {
                    i__2 = i__ + i__ * a_dim1;
                    q__1.r = -reeps.r, q__1.i = -reeps.i;
                    a[i__2].r = q__1.r, a[i__2].i = q__1.i;
                }
                if (i__ % 2 != 0 && i__ < *m) {
                    i__2 = i__ + (i__ + 1) * a_dim1;
                    a[i__2].r = 1.f, a[i__2].i = 0.f;
                } else if (i__ > 1) {
                    i__2 = i__ + (i__ - 1) * a_dim1;
                    q__1.r = -1.f, q__1.i = -0.f;
                    a[i__2].r = q__1.r, a[i__2].i = q__1.i;
                }
            } else {
                i__2 = i__ + i__ * a_dim1;
                a[i__2].r = 1.f, a[i__2].i = 0.f;
                if (i__ % 2 != 0 && i__ < *m) {
                    i__2 = i__ + (i__ + 1) * a_dim1;
                    d__1 = 2.;
                    q__1.r = d__1 * imeps.r, q__1.i = d__1 * imeps.i;
                    a[i__2].r = q__1.r, a[i__2].i = q__1.i;
                } else if (i__ > 1) {
                    i__2 = i__ + (i__ - 1) * a_dim1;
                    q__2.r = -imeps.r, q__2.i = -imeps.i;
                    d__1 = 2.;
                    q__1.r = d__1 * q__2.r, q__1.i = d__1 * q__2.i;
                    a[i__2].r = q__1.r, a[i__2].i = q__1.i;
                }
            }
/* L240: */
        }

        i__1 = *n;
        for (i__ = 1; i__ <= i__1; ++i__) {
            i__2 = i__ + i__ * e_dim1;
            e[i__2].r = 1.f, e[i__2].i = 0.f;
            if (i__ <= 4) {
                i__2 = i__ + i__ * b_dim1;
                q__1.r = -1.f, q__1.i = -0.f;
                b[i__2].r = q__1.r, b[i__2].i = q__1.i;
                if (i__ > 2) {
                    i__2 = i__ + i__ * b_dim1;
                    q__1.r = 1.f - reeps.r, q__1.i = 0.f - reeps.i;
                    b[i__2].r = q__1.r, b[i__2].i = q__1.i;
                }
                if (i__ % 2 != 0 && i__ < *n) {
                    i__2 = i__ + (i__ + 1) * b_dim1;
                    b[i__2].r = imeps.r, b[i__2].i = imeps.i;
                } else if (i__ > 1) {
                    i__2 = i__ + (i__ - 1) * b_dim1;
                    q__1.r = -imeps.r, q__1.i = -imeps.i;
                    b[i__2].r = q__1.r, b[i__2].i = q__1.i;
                }
            } else if (i__ <= 8) {
                if (i__ <= 6) {
                    i__2 = i__ + i__ * b_dim1;
                    b[i__2].r = reeps.r, b[i__2].i = reeps.i;
                } else {
                    i__2 = i__ + i__ * b_dim1;
                    q__1.r = -reeps.r, q__1.i = -reeps.i;
                    b[i__2].r = q__1.r, b[i__2].i = q__1.i;
                }
                if (i__ % 2 != 0 && i__ < *n) {
                    i__2 = i__ + (i__ + 1) * b_dim1;
                    q__1.r = imeps.r + 1.f, q__1.i = imeps.i + 0.f;
                    b[i__2].r = q__1.r, b[i__2].i = q__1.i;
                } else if (i__ > 1) {
                    i__2 = i__ + (i__ - 1) * b_dim1;
                    q__2.r = -1.f, q__2.i = -0.f;
                    q__1.r = q__2.r - imeps.r, q__1.i = q__2.i - imeps.i;
                    b[i__2].r = q__1.r, b[i__2].i = q__1.i;
                }
            } else {
                i__2 = i__ + i__ * b_dim1;
                q__1.r = 1.f - reeps.r, q__1.i = 0.f - reeps.i;
                b[i__2].r = q__1.r, b[i__2].i = q__1.i;
                if (i__ % 2 != 0 && i__ < *n) {
                    i__2 = i__ + (i__ + 1) * b_dim1;
                    d__1 = 2.;
                    q__1.r = d__1 * imeps.r, q__1.i = d__1 * imeps.i;
                    b[i__2].r = q__1.r, b[i__2].i = q__1.i;
                } else if (i__ > 1) {
                    i__2 = i__ + (i__ - 1) * b_dim1;
                    q__2.r = -imeps.r, q__2.i = -imeps.i;
                    d__1 = 2.;
                    q__1.r = d__1 * q__2.r, q__1.i = d__1 * q__2.i;
                    b[i__2].r = q__1.r, b[i__2].i = q__1.i;
                }
            }
/* L250: */
        }
    }

/*     Compute rhs (C, F) */

    cgemm_("N", "N", m, n, m, &c_b1, &a[a_offset], lda, &r__[r_offset], ldr, &
            c_b3, &c__[c_offset], ldc);
    q__1.r = -1.f, q__1.i = -0.f;
    cgemm_("N", "N", m, n, n, &q__1, &l[l_offset], ldl, &b[b_offset], ldb, &
            c_b1, &c__[c_offset], ldc);
    cgemm_("N", "N", m, n, m, &c_b1, &d__[d_offset], ldd, &r__[r_offset], ldr, 
             &c_b3, &f[f_offset], ldf);
    q__1.r = -1.f, q__1.i = -0.f;
    cgemm_("N", "N", m, n, n, &q__1, &l[l_offset], ldl, &e[e_offset], lde, &
            c_b1, &f[f_offset], ldf);

/*     End of CLATM5 */

    return 0;
} /* clatm5_ */