hack_array_seq.icc

Go to the documentation of this file.

# include "rheolef/hack_array.h"
namespace rheolef {
 
// ===============================================================
// allocators
// ===============================================================
template<class T, class A>
hack_array_seq_rep<T,A>::hack_array_seq_rep (const A& alloc) 
 : base (alloc),
   _ownership  (),
   _parameter  (),
   _value_size (),
   _data_size  ()
{
}
template<class T, class A>
hack_array_seq_rep<T,A>::hack_array_seq_rep (const distributor& ownership, const parameter_type& param, const A& alloc) 
 : base (alloc),
   _ownership  (ownership),
   _parameter  (param),
   _value_size (T::_size_of  (_parameter)),
   _data_size  (T::_data_size(_parameter))
{
  _init (_ownership, _parameter);
}
template<class T, class A>
hack_array_seq_rep<T,A>::hack_array_seq_rep (size_type n, const parameter_type& param, const A& alloc) 
 : base (alloc),
   _ownership  (distributor::decide, communicator(), n),
   _parameter  (param),
   _value_size (T::_size_of  (_parameter)),
   _data_size  (T::_data_size(_parameter))
{
  _init (_ownership, _parameter);
}
template<class T, class A>
void
hack_array_seq_rep<T,A>::resize (const distributor& ownership, const parameter_type& param) 
{
  _ownership  = ownership;
  _parameter  = param;
  _value_size = T::_size_of  (_parameter);
  _data_size  = T::_data_size(_parameter);
  _init (_ownership, _parameter);
}
template<class T, class A>
void
hack_array_seq_rep<T,A>::resize (size_type n, const parameter_type& param)
{
  _ownership  = distributor  (distributor::decide, communicator(), n);
  _parameter  = param;
  _value_size = T::_size_of  (_parameter);
  _data_size  = T::_data_size(_parameter);
  _init (_ownership, _parameter);
} 
template<class T, class A>
void
hack_array_seq_rep<T,A>::_init (const distributor& ownership, const parameter_type& param) 
{
  base::_ownership = distributor (ownership.dis_size()*_value_size, ownership.comm(), ownership.size()*_value_size),
  base::resize (base::_ownership, std::numeric_limits<raw_type>::max());
  size_type *p = base::begin().operator->();
  for (size_type i = 0, n = ownership.size(); i < n; i++) {
    new (p) T;
    T* q = (T*)p;
    q -> _set_parameter(param);
    p += _value_size;
  }
}
// ===============================================================
// put & get
// ===============================================================
template <class T, class A>
template <class PutFunction>
odiststream&
hack_array_seq_rep<T,A>::put_values (odiststream& ops, PutFunction put_element) const
{
  std::ostream& os = ops.os();
  for (size_type i = 0; i < size(); i++) {
    put_element (os, operator[](i));
    os << std::endl;
  }
  return ops;
}
template <class T, class A>
odiststream&
hack_array_seq_rep<T,A>::put_values (odiststream& ops) const
{
  return put_values (ops, _disarray_put_element_type<generic_value_type>());
}
template <class T, class A>
template <class GetFunction>
idiststream&
hack_array_seq_rep<T,A>::get_values (idiststream& ips, GetFunction get_element) {
  check_macro (load_chunk (ips.is(), begin(), end(), get_element), "read failed on input stream.");
  return ips;
}
template <class T, class A>
idiststream&
hack_array_seq_rep<T,A>::get_values (idiststream& ips) 
{
  return get_values (ips, _disarray_get_element_type<generic_value_type>());
}
 
} // namespace rheolef