The R Project SVN R

Rev

Rev 35296 | Rev 36679 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 35296 Rev 36367
Line 413... Line 413...
413 
			    x, &nrx, y, &nry, &zero, z, &nrx);
 413 
			    x, &nrx, y, &nry, &zero, z, &nrx);
414 
    } else /* zero-extent operations should return zeroes */
 414 
    } else /* zero-extent operations should return zeroes */
415 
	for(i = 0; i < nrx*ncy; i++) z[i] = 0;
 415 
	for(i = 0; i < nrx*ncy; i++) z[i] = 0;
416 
}
 416 
}
417 
 
 417 
 
418 
#ifdef HAVE_FORTRAN_DOUBLE_COMPLEX
 - 
 
419 
/* ZGEMM - perform one of the matrix-matrix operations    */
 - 
 
420 
/* C := alpha*op( A )*op( B ) + beta*C */
 - 
 
421 
extern void
 - 
 
422 
F77_NAME(zgemm)(const char *transa, const char *transb, const int *m,
 - 
 
423 
		const int *n, const int *k, const Rcomplex *alpha,
 - 
 
424 
		const Rcomplex *a, const int *lda,
 - 
 
425 
		const Rcomplex *b, const int *ldb,
 - 
 
426 
		const Rcomplex *beta, Rcomplex *c, const int *ldc);
 - 
 
427 
#endif
 - 
 
428 
 
 - 
 
429 
static void cmatprod(Rcomplex *x, int nrx, int ncx,
 418 
static void cmatprod(Rcomplex *x, int nrx, int ncx,
430 
		     Rcomplex *y, int nry, int ncy, Rcomplex *z)
 419 
		     Rcomplex *y, int nry, int ncy, Rcomplex *z)
431 
{
 420 
{
432 
#ifdef HAVE_FORTRAN_DOUBLE_COMPLEX
 421 
#ifdef HAVE_FORTRAN_DOUBLE_COMPLEX
433 
    char *transa = "N", *transb = "N";
 422 
    char *transa = "N", *transb = "N";