CPPLapack Code

//=============================================================================
template<long m, long n>
inline dgematrix_small<n,m> t(const dgematrix_small<m,n>& A)
{
  dgematrix_small<n,m> X;
  for(long i=0; i<m; i++){
    for(long j=0; j<n; j++){
      X(j,i) =A(i,j);
    }
  }
  return X;
}