DistrArrayDisk.h

#ifndef LINEARALGEBRA_SRC_MOLPRO_LINALG_ARRAY_DISTRARRAYDISK_H
#define LINEARALGEBRA_SRC_MOLPRO_LINALG_ARRAY_DISTRARRAYDISK_H
 
#include <future>
#include <molpro/Profiler.h>
#include <molpro/linalg/array/DistrArray.h>
#include <molpro/linalg/array/Span.h>
 
namespace molpro::linalg::array {
class DistrArrayDisk : public DistrArray {
public:
  using disk_array = void; 
protected:
  bool m_allocated = false;                     
  std::unique_ptr<Distribution> m_distribution; 
  size_t m_buffer_size = 8192;                  
  using DistrArray::DistrArray;
 
  DistrArrayDisk(std::unique_ptr<Distribution> distr, MPI_Comm commun);
  DistrArrayDisk();
  DistrArrayDisk(const DistrArray& source);
  DistrArrayDisk(DistrArrayDisk&& source) noexcept;
  ~DistrArrayDisk() override;
 
public:
  virtual void erase() = 0;
  [[nodiscard]] const Distribution& distribution() const override;
  [[nodiscard]] value_type dot(const DistrArrayDisk& y) const { return DistrArray::dot(y); }
  using DistrArray::dot;
  void set_buffer_size(size_t buffer_size) { m_buffer_size = buffer_size; }
  [[nodiscard]] size_t disk_page_size() const override { return m_buffer_size; }
 
  [[nodiscard]] std::unique_ptr<LocalBuffer> local_buffer() override;
  [[nodiscard]] std::unique_ptr<const LocalBuffer> local_buffer() const override;
};
 
double dot(const DistrArrayDisk& x, const DistrArrayDisk& y);
double dot(const DistrArrayDisk& x, const DistrArray& y);
double dot(const DistrArray& x, const DistrArrayDisk& y);
 
} // namespace molpro::linalg::array
 
#endif // LINEARALGEBRA_SRC_MOLPRO_LINALG_ARRAY_DISTRARRAYDISK_H