jiri/LA_library
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

*** empty log message ***

Browse Source
This commit is contained in:
jiri 2006-09-19 15:59:49 +00:00
parent 4f55fe4fc4
commit 6d7b245653
1 changed files with 2 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
mat.cc
View File

@ -696,6 +696,7 @@ const NRMat<double> NRMat<double>::operator*(const NRMat<double> &rhs) const
{ {
#ifdef DEBUG #ifdef DEBUG
if (mm != rhs.nn) laerror("product of incompatible matrices"); if (mm != rhs.nn) laerror("product of incompatible matrices");
if (rhs.mm <=0) laerror("illegal matrix dimension in gemm");
#endif #endif
NRMat<double> result(nn, rhs.mm); NRMat<double> result(nn, rhs.mm);
cblas_dgemm(CblasRowMajor, CblasNoTrans, CblasNoTrans, nn, rhs.mm, mm, 1.0, cblas_dgemm(CblasRowMajor, CblasNoTrans, CblasNoTrans, nn, rhs.mm, mm, 1.0,
@ -868,6 +869,7 @@ void NRMat<double>::gemm(const double &beta, const NRMat<double> &a,
int kk(transb=='n'?b.nn:b.mm); int kk(transb=='n'?b.nn:b.mm);
int ll(transb=='n'?b.mm:b.nn); int ll(transb=='n'?b.mm:b.nn);
if (l!=nn || ll!=mm || k!=kk) laerror("incompatible matrices in Mat:gemm()"); if (l!=nn || ll!=mm || k!=kk) laerror("incompatible matrices in Mat:gemm()");
if(b.mm <=0 || mm<=0) laerror("illegal matrix dimension in gemm");
#endif #endif
if (alpha==0.0 && beta==1.0) return; if (alpha==0.0 && beta==1.0) return;