matroid_permuted.hpp

Go to the documentation of this file.

#pragma once
 
#include "matroid.hpp"
 
namespace tu
{
 
  template <typename MatroidType>
  class matroid_permuted
  {
  public:
    typedef MatroidType matroid_type;
    typedef typename MatroidType::size_type size_type;
    typedef typename MatroidType::name_type name_type;
    typedef typename boost::mpl::if_ <boost::is_const <MatroidType>, typename MatroidType::const_reference_type, typename MatroidType::reference_type>::type
        reference_type;
    typedef const name_type& const_reference_type;
    typedef matroid_permuted <name_type> self_type;
    typedef permutation permutation_type;
 
    matroid_permuted(MatroidType& matroid) :
      _matroid(matroid), _perm1(matroid.size1()), _perm2(matroid.size2())
    {
 
    }
 
    inline size_type size1() const
    {
      return _perm1.size();
    }
 
    inline size_type size2() const
    {
      return _perm2.size();
    }
 
    inline const permutation_type& perm1() const
    {
      return _perm1;
    }
 
    inline permutation_type& perm1()
    {
      return _perm1;
    }
 
    inline const permutation_type& perm2() const
    {
      return _perm2;
    }
 
    inline permutation_type& perm2()
    {
      return _perm2;
    }
 
    inline reference_type name1(size_type index)
    {
      return _matroid.name1(_perm1(index));
    }
 
    inline const_reference_type name1(size_type index) const
    {
      return _matroid.name1(_perm1(index));
    }
 
    inline reference_type name2(size_type index)
    {
      return _matroid.name2(_perm2(index));
    }
 
    inline const_reference_type name2(size_type index) const
    {
      return _matroid.name2(_perm2(index));
    }
 
    inline matroid_type& data()
    {
      return _matroid;
    }
 
  private:
    MatroidType& _matroid;
    permutation_type _perm1;
    permutation_type _perm2;
 
  };
 
  template <typename MatroidType>
  inline void matroid_permute1(matroid_permuted <MatroidType>& matroid, size_t index1, size_t index2)
  {
    matroid.perm1().swap(index1, index2);
  }
 
  template <typename MatroidType>
  inline void matroid_permute2(matroid_permuted <MatroidType>& matroid, size_t index1, size_t index2)
  {
    matroid.perm2().swap(index1, index2);
  }
 
  template <typename MatroidType>
  void matroid_binary_pivot(matroid_permuted <MatroidType>& matroid, size_t i, size_t j)
  {
    matroid_binary_pivot(matroid.data(), matroid.perm1()(i), matroid.perm2()(j));
  }
 
  template <typename MatroidType>
  inline void matroid_set_perm1(matroid_permuted <MatroidType>& matroid, const permutation& permutation)
  {
    matroid.perm1() = permutation;
  }
 
  template <typename MatroidType>
  inline void matroid_set_perm2(matroid_permuted <MatroidType>& matroid, const permutation& permutation)
  {
    matroid.perm2() = permutation;
  }
 
  template <typename MatroidType>
  inline void matroid_set_perm1(matroid_transposed <MatroidType>& matroid, const permutation& permutation)
  {
    matroid_set_perm2(matroid.data(), permutation);
  }
 
  template <typename MatroidType>
  inline void matroid_set_perm2(matroid_transposed <MatroidType>& matroid, const permutation& permutation)
  {
    matroid_set_perm1(matroid.data(), permutation);
  }
 
} /* namespace tu */