Hi,
In the CMKL 8.1 manual, page 153, description of ?herk routine.
One can read
---
a is an n-by-k matrix in the first case and a n-by-k matrix in the second case.
---
on the next page
---
a
COMPLEX for cherk
DOUBLE COMPLEX for zherk
Array, DIMENSION (lda, ka), where ka is k when trans = 'N' or 'n', and is n otherwise. Before entry with trans = 'N' or 'n', the leading n-by-k part of the array a must contain the matrix a, otherwise the leading n-by-k part of the array a must contain the matrix a.
---
However in MKL 8.0.1 documentation (page 151) the same things are explained as
---
a is an n-by-k matrix in the first case and a k-by-n matrix in the second case.
---
a
COMPLEX for cherk
DOUBLE COMPLEX for zherk
Array, DIMENSION (lda, ka), where ka is k when trans = 'N' or 'n', and is n otherwise. Before entry with trans = 'N' or 'n', the leading n-by-k part of the array a must contain the matrix a, otherwise the leading k-by-n part of the array a must contain the matrix a.
---
MKL documentation agrees with the description of zherk from
http://www.netlib.org/blas/zherk.f .
Alin
Message Edited by alin.elena@qub.ac.uk on 05-15-200611:02 AM
Message Edited by alin.elena@qub.ac.uk on 05-15-200611:04 AM
