indirectArray.hpp

Go to the documentation of this file.

#ifndef INDIRECTARRAY_HPP_
#define INDIRECTARRAY_HPP_

#include <iostream>

#include "boost/numeric/ublas/operation.hpp"
#include "boost/numeric/ublas/io.hpp"

#include "jmath/jblas.hpp"

namespace std {
  template<class A>
  ostream& operator <<(ostream & s, const boost::numeric::ublas::indirect_array<A> & ia_) {
    s << "[" << ia_.size() << "]{";
    for (size_t i = 0; i < ia_.size(); i++) {
      s << ia_(i);
      if (i + 1 < ia_.size())
        s << ",";
    }
    s << "}";
    return s;
  }
}

namespace jafar {
  namespace jmath {
    namespace ublasExtra {

      jblas::ind_array ia_set(const jblas::vecb & vb);

      jblas::ind_array ia_set(const jblas::veci & vi);

      jblas::ind_array ia_set(const ublas::range & r);

      jblas::ind_array ia_set(size_t begin, size_t end);

      jblas::ind_array ia_set(const ublas::slice & s);

      jblas::ind_array ia_complement(const jblas::ind_array & ifull, const jblas::ind_array & ipartial);

      jblas::ind_array ia_union(const jblas::ind_array & ia1, const jblas::ind_array & ia2);

      jblas::ind_array ia_concat(const jblas::ind_array & ia1, const jblas::ind_array & ia2);

      jblas::ind_array ia_head(const jblas::ind_array & ia, const size_t N);

      jblas::ind_array ia_head(const jblas::vecb & vb, const size_t N);

      jblas::ind_array ia_popfront(jblas::vecb & vb, const size_t N);

      jblas::ind_array ia_pushfront(jblas::vecb & vb, const size_t N);

      jblas::ind_array operator+(const jblas::ind_array &ia, int n);
      inline jblas::ind_array operator-(const jblas::ind_array &ia, int n) { return ia+(-n); }
      
    }
  }
}

#endif /* INDIRECTARRAY_HPP_ */