FFastTrig::FFastTrig()
{
const double pimul = M_PI * 2 / TBLPERIOD;
for (int i = 0; i < 2049; i++)
{
sinetable[i] = (float)c_sin(i*pimul);
}
}
/* 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 */
static integer i__, j, k;
extern /* Subroutine */ int cgemm_(char *, char *, integer *, integer *,
integer *, complex *, complex *, integer *, complex *, integer *,
complex *, complex *, integer *);
static complex imeps, reeps;
#define a_subscr(a_1,a_2) (a_2)*a_dim1 + a_1
#define a_ref(a_1,a_2) a[a_subscr(a_1,a_2)]
#define b_subscr(a_1,a_2) (a_2)*b_dim1 + a_1
#define b_ref(a_1,a_2) b[b_subscr(a_1,a_2)]
#define d___subscr(a_1,a_2) (a_2)*d_dim1 + a_1
#define d___ref(a_1,a_2) d__[d___subscr(a_1,a_2)]
#define e_subscr(a_1,a_2) (a_2)*e_dim1 + a_1
#define e_ref(a_1,a_2) e[e_subscr(a_1,a_2)]
#define l_subscr(a_1,a_2) (a_2)*l_dim1 + a_1
#define l_ref(a_1,a_2) l[l_subscr(a_1,a_2)]
#define r___subscr(a_1,a_2) (a_2)*r_dim1 + a_1
#define r___ref(a_1,a_2) r__[r___subscr(a_1,a_2)]
/* -- LAPACK test routine (version 3.0) --
Univ. of Tennessee, Univ. of California Berkeley, NAG Ltd.,
Courant Institute, Argonne National Lab, and Rice University
June 30, 1999
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.
=====================================================================
Parameter adjustments */
a_dim1 = *lda;
a_offset = 1 + a_dim1 * 1;
a -= a_offset;
b_dim1 = *ldb;
b_offset = 1 + b_dim1 * 1;
b -= b_offset;
c_dim1 = *ldc;
c_offset = 1 + c_dim1 * 1;
c__ -= c_offset;
d_dim1 = *ldd;
d_offset = 1 + d_dim1 * 1;
d__ -= d_offset;
e_dim1 = *lde;
e_offset = 1 + e_dim1 * 1;
e -= e_offset;
f_dim1 = *ldf;
f_offset = 1 + f_dim1 * 1;
f -= f_offset;
r_dim1 = *ldr;
r_offset = 1 + r_dim1 * 1;
r__ -= r_offset;
l_dim1 = *ldl;
l_offset = 1 + l_dim1 * 1;
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 = a_subscr(i__, j);
a[i__3].r = 1.f, a[i__3].i = 0.f;
i__3 = d___subscr(i__, j);
d__[i__3].r = 1.f, d__[i__3].i = 0.f;
} else if (i__ == j - 1) {
i__3 = a_subscr(i__, j);
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 = d___subscr(i__, j);
d__[i__3].r = 0.f, d__[i__3].i = 0.f;
} else {
i__3 = a_subscr(i__, j);
a[i__3].r = 0.f, a[i__3].i = 0.f;
i__3 = d___subscr(i__, j);
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 = b_subscr(i__, j);
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 = e_subscr(i__, j);
e[i__3].r = 1.f, e[i__3].i = 0.f;
} else if (i__ == j - 1) {
i__3 = b_subscr(i__, j);
b[i__3].r = 1.f, b[i__3].i = 0.f;
i__3 = e_subscr(i__, j);
e[i__3].r = 0.f, e[i__3].i = 0.f;
} else {
i__3 = b_subscr(i__, j);
b[i__3].r = 0.f, b[i__3].i = 0.f;
i__3 = e_subscr(i__, j);
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 = r___subscr(i__, j);
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 = l_subscr(i__, j);
i__4 = r___subscr(i__, j);
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 = a_subscr(i__, j);
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 = d___subscr(i__, j);
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 = a_subscr(i__, j);
a[i__3].r = 0.f, a[i__3].i = 0.f;
i__3 = d___subscr(i__, j);
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 = b_subscr(i__, j);
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 = e_subscr(i__, j);
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 = b_subscr(i__, j);
b[i__3].r = 0.f, b[i__3].i = 0.f;
i__3 = e_subscr(i__, j);
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 = r___subscr(i__, j);
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 = l_subscr(i__, j);
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 = a_subscr(k + 1, k + 1);
i__4 = a_subscr(k, k);
a[i__3].r = a[i__4].r, a[i__3].i = a[i__4].i;
i__3 = a_subscr(k + 1, k);
c_sin(&q__2, &a_ref(k, k + 1));
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 = b_subscr(k + 1, k + 1);
i__4 = b_subscr(k, k);
b[i__3].r = b[i__4].r, b[i__3].i = b[i__4].i;
i__3 = b_subscr(k + 1, k);
c_sin(&q__2, &b_ref(k, k + 1));
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 = a_subscr(i__, j);
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 = d___subscr(i__, j);
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 = b_subscr(i__, j);
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 = e_subscr(i__, j);
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 = r___subscr(i__, j);
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 = l_subscr(i__, j);
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 = r___subscr(i__, j);
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 = l_subscr(i__, j);
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 = d___subscr(i__, i__);
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 = a_subscr(i__, i__);
a[i__2].r = 1.f, a[i__2].i = 0.f;
if (i__ > 2) {
i__2 = a_subscr(i__, i__);
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 = a_subscr(i__, i__ + 1);
a[i__2].r = imeps.r, a[i__2].i = imeps.i;
} else if (i__ > 1) {
i__2 = a_subscr(i__, i__ - 1);
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 = a_subscr(i__, i__);
a[i__2].r = reeps.r, a[i__2].i = reeps.i;
} else {
i__2 = a_subscr(i__, i__);
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 = a_subscr(i__, i__ + 1);
a[i__2].r = 1.f, a[i__2].i = 0.f;
} else if (i__ > 1) {
i__2 = a_subscr(i__, i__ - 1);
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 = a_subscr(i__, i__);
a[i__2].r = 1.f, a[i__2].i = 0.f;
if (i__ % 2 != 0 && i__ < *m) {
i__2 = a_subscr(i__, i__ + 1);
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 = a_subscr(i__, i__ - 1);
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 = e_subscr(i__, i__);
e[i__2].r = 1.f, e[i__2].i = 0.f;
if (i__ <= 4) {
i__2 = b_subscr(i__, i__);
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 = b_subscr(i__, i__);
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 = b_subscr(i__, i__ + 1);
b[i__2].r = imeps.r, b[i__2].i = imeps.i;
} else if (i__ > 1) {
i__2 = b_subscr(i__, i__ - 1);
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 = b_subscr(i__, i__);
b[i__2].r = reeps.r, b[i__2].i = reeps.i;
} else {
i__2 = b_subscr(i__, i__);
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 = b_subscr(i__, i__ + 1);
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 = b_subscr(i__, i__ - 1);
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 = b_subscr(i__, i__);
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 = b_subscr(i__, i__ + 1);
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 = b_subscr(i__, i__ - 1);
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_ */
