Esempio n. 1
0
void BLASInterface::Mul(const dMatrix& A,const dMatrix& B,dMatrix& X)
{
  X.resize(A.m,B.n);
  X.setZero();
  if(X.isRowMajor()) std::swap(X.istride,X.jstride);
  Madd(A,B,X);
}
Esempio n. 2
0
void Madd(const dMatrix& A,const dVector& x,dVector& y,double alpha,double beta,bool transpose)
{
  if(A.isRowMajor()) {
    Assert(A.jstride == 1);
    transpose = !transpose;
  }
  else {
    Assert(A.istride == 1);
    //Assert(IsCompliant(A));
  }
  //Assert(IsCompliant(A));
  Assert(x.n == A.n);
  Assert(y.n == A.m);
  integer m=A.m;
  integer n=A.n;
  integer lda=A.m;
  integer xinc=x.stride;
  integer yinc=y.stride;
  char trans = (transpose?'T':'N');
  dgemv_(&trans,&m,&n,&alpha,A.getStart(),&lda,x.getStart(),&xinc,&beta,y.getStart(),&yinc);
}
Esempio n. 3
0
bool IsCompliant(const dMatrix& A)
{
  return (!A.isRowMajor() && A.istride == 1);
}