Exemple #1
0
void MatrixTemplate<T>::maddTranspose(const VectorT& a, VectorT& b) const
{
  if(m != a.n)
  {
    RaiseErrorFmt(WHERE_AM_I,MatrixError_ArgIncompatibleDimensions);
  }
  if(b.n == 0)
  {
    b.resize(n);
  }
  else if(b.n != n)
  {
    RaiseErrorFmt(WHERE_AM_I,MatrixError_DestIncompatibleDimensions);
  }

  gen_array2d_vector_madd_transpose(b.getStart(),b.stride,
    getStart(),istride,jstride, 
    a.getStart(),a.stride, 
    m, n);
}