synaps/linalg/Rank.h

Go to the documentation of this file.

/*********************************************************************
*      This file is part of the source code of SYNAPS kernel.        *
*   Author(s): B. Mourrain, GALAAD, INRIA                            *
**********************************************************************
$Id: Rank.h,v 1.1 2005/07/11 08:33:34 mourrain Exp $
**********************************************************************/
#ifndef SYNAPS_LINALG_RANK_H
#define SYNAPS_LINALG_RANK_H

#include <synaps/init.h>
#include <synaps/linalg/MATRIX.m>

__BEGIN_NAMESPACE_SYNAPS


//----------------------------------------------------------------------
namespace lapack {
  template<class R> struct rep;
}
//--------------------------------------------------------------------
template<class R> inline
typename R::size_type Rank(const R & m) 
{
  using namespace MATRIX;
  return rank(m.rep());
}
//--------------------------------------------------------------------
template<class R,class MTH> inline
typename R::size_type Rank(const R & m, const MTH& mth) 
{
  using MATRIX::rank;
  return rank(mth,m.rep());
}
//----------------------------------------------------------------------
template<class R,class V> 
typename R::size_type Rank(const R & M, 
                           V & I, V & J)
{
  lapack::rep<typename R::value_type> A(M.nbrow(),M.nbcol());
  assign(A,M);
  return rank(A,I.begin(), J.begin());
}
//----------------------------------------------------------------------
template<class R,class I> 
R Extract(R & M, const I & ipvt, const I & jpvt, 
                const typename R::size_type & rank)
{
   R SM(rank,rank);
   for (typename R::size_type i=0; i<rank; i++)
     for (typename R::size_type j = 0; j < rank; j++) 
       SM(i,j)=M(ipvt[i],jpvt[j]);
   return (SM);
}
//--------------------------------------------------------------------
__END_NAMESPACE_SYNAPS
//--------------------------------------------------------------------
#endif // SYNAPS_LINALG_RANK_H