/* zlacrm.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 doublereal c_b6 = 1.;
static doublereal c_b7 = 0.;
/* Subroutine */ int zlacrm_(integer *m, integer *n, doublecomplex *a,
integer *lda, doublereal *b, integer *ldb, doublecomplex *c__,
integer *ldc, doublereal *rwork)
{
/* System generated locals */
integer b_dim1, b_offset, a_dim1, a_offset, c_dim1, c_offset, i__1, i__2,
i__3, i__4, i__5;
doublereal d__1;
doublecomplex z__1;
/* Builtin functions */
double d_imag(doublecomplex *);
/* Local variables */
integer i__, j, l;
extern /* Subroutine */ int dgemm_(char *, char *, integer *, integer *,
integer *, doublereal *, doublereal *, integer *, doublereal *,
integer *, doublereal *, doublereal *, integer *);
/* -- LAPACK auxiliary routine (version 3.2) -- */
/* Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. */
/* November 2006 */
/* .. Scalar Arguments .. */
/* .. */
/* .. Array Arguments .. */
/* .. */
/* Purpose */
/* ======= */
/* ZLACRM performs a very simple matrix-matrix multiplication: */
/* C := A * B, */
/* where A is M by N and complex; B is N by N and real; */
/* C is M by N and complex. */
/* Arguments */
/* ========= */
/* M (input) INTEGER */
/* The number of rows of the matrix A and of the matrix C. */
/* M >= 0. */
/* N (input) INTEGER */
/* The number of columns and rows of the matrix B and */
/* the number of columns of the matrix C. */
/* N >= 0. */
/* A (input) COMPLEX*16 array, dimension (LDA, N) */
/* A contains the M by N matrix A. */
/* LDA (input) INTEGER */
/* The leading dimension of the array A. LDA >=max(1,M). */
/* B (input) DOUBLE PRECISION array, dimension (LDB, N) */
/* B contains the N by N matrix B. */
/* LDB (input) INTEGER */
/* The leading dimension of the array B. LDB >=max(1,N). */
/* C (input) COMPLEX*16 array, dimension (LDC, N) */
/* C contains the M by N matrix C. */
/* LDC (input) INTEGER */
/* The leading dimension of the array C. LDC >=max(1,N). */
/* RWORK (workspace) DOUBLE PRECISION array, dimension (2*M*N) */
/* ===================================================================== */
/* .. Parameters .. */
/* .. */
/* .. Local Scalars .. */
/* .. */
/* .. Intrinsic Functions .. */
/* .. */
/* .. External Subroutines .. */
/* .. */
/* .. Executable Statements .. */
/* Quick return if possible. */
/* 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;
--rwork;
/* Function Body */
if (*m == 0 || *n == 0) {
return 0;
}
i__1 = *n;
for (j = 1; j <= i__1; ++j) {
i__2 = *m;
for (i__ = 1; i__ <= i__2; ++i__) {
i__3 = i__ + j * a_dim1;
rwork[(j - 1) * *m + i__] = a[i__3].r;
/* L10: */
}
/* L20: */
}
l = *m * *n + 1;
dgemm_("N", "N", m, n, n, &c_b6, &rwork[1], m, &b[b_offset], ldb, &c_b7, &
rwork[l], m);
i__1 = *n;
for (j = 1; j <= i__1; ++j) {
i__2 = *m;
for (i__ = 1; i__ <= i__2; ++i__) {
i__3 = i__ + j * c_dim1;
i__4 = l + (j - 1) * *m + i__ - 1;
c__[i__3].r = rwork[i__4], c__[i__3].i = 0.;
/* L30: */
}
/* L40: */
}
i__1 = *n;
for (j = 1; j <= i__1; ++j) {
i__2 = *m;
for (i__ = 1; i__ <= i__2; ++i__) {
rwork[(j - 1) * *m + i__] = d_imag(&a[i__ + j * a_dim1]);
/* L50: */
}
/* L60: */
}
dgemm_("N", "N", m, n, n, &c_b6, &rwork[1], m, &b[b_offset], ldb, &c_b7, &
rwork[l], m);
i__1 = *n;
for (j = 1; j <= i__1; ++j) {
i__2 = *m;
for (i__ = 1; i__ <= i__2; ++i__) {
i__3 = i__ + j * c_dim1;
i__4 = i__ + j * c_dim1;
d__1 = c__[i__4].r;
i__5 = l + (j - 1) * *m + i__ - 1;
z__1.r = d__1, z__1.i = rwork[i__5];
c__[i__3].r = z__1.r, c__[i__3].i = z__1.i;
/* L70: */
}
/* L80: */
}
return 0;
/* End of ZLACRM */
} /* zlacrm_ */