#ifndef RADIX_RADIXIO_F71STREAM_HH_
#define RADIX_RADIXIO_F71STREAM_HH_

#include <algorithm>
#include <cmath>
#include <iomanip>
#include <memory>
#include <sstream>
#include <string>
#include <vector>

#include "radixio/eafstream.hh"

#include "radixcore/visibility.hh"

namespace radix
{
class RADIX_PUBLIC F71Stream
{
 private:
  std::string mFile;
  eafstream mStream;

 public:
  F71Stream() = delete;
  F71Stream(const std::string& mFile);

  /**
   * @brief read_case Reads a single case from the stream
   * @param f71case_type& inventory
   *
   * Requires the f71case_type to have the following methods:
   *
   * setLightElementNuclides(const std::vector<int>& zaids)
   * setActinideNuclides(const std::vector<int>& zaids)
   * setFissionProductNuclides(const std::vector<int>& zaids)
   *
   * setLightElementNuclideAbundances(const std::vector<float>& abs)
   * setActinideNuclideAbundances(const std::vector<float>& abs)
   * setFissionProductNuclideAbundances(const std::vector<float>& abs)
   *
   * setTitle(const std::string& title)
   * setBasis(const std::string& basis)
   *
   * setLightElementGammaEnergyBounds(const std::vector<float>& bounds)
   * setActinideGammaEnergyBounds(const std::vector<float>& bounds)
   * setFissionProductGammaEnergyBounds(const std::vector<float>& bounds)
   *
   * setLightElementSpectrum(const std::vector<float>& spectrum)
   * setActinideSpectrum(const std::vector<float>& spectrum)
   * setFissionProductSpectrum(const std::vector<float>& spectrum)
   *
   * setGammaEnergyBounds(const std::vector<float>& bounds)
   * setPhotonPerSSpectrum(const std::vector<float>& spectrum)
   * setMevPerSSpectrum(const std::vector<float>& spectrum)
   *
   * setPhotonSpectrumTitle(const std::string& title)
   * setSpontFissionNeutronSource(const std::vector<float>& source)
   * setActinideAlphaNSource(const std::vector<float>& source)
   * setTotalNeutronSpectrum(const std::vector<float>& spectrum)
   * setAlphaNSpectrum(const std::vector<float>& spectrum)
   * setSpontaneousFissionSpectrum(const std::vector<float>& spectrum)
   * setNeutronEnergyBounds(const std::vector<float>& bounds)
   *
   * setTime(double seconds)
   * setPower(float power)
   * setFlux(float flux)
   * @return bool on if inventory was populated
   */
  template <typename f71case_type>
  bool read_case(f71case_type& inventory);

};  // class F71Stream

class RADIX_PUBLIC F71Case
{
  // forward declare private implementation
  class PImpl;
  // unique pointer to private implmentation
  std::unique_ptr<PImpl, void (*)(PImpl*)> p;

 public:
  F71Case();
  F71Case(const F71Case& orig);

  const std::vector<int>& lightElementNuclides() const;
  void setLightElementNuclides(const std::vector<int>& zaids);

  const std::vector<int>& actinideNuclides() const;
  void setActinideNuclides(const std::vector<int>& zaids);

  const std::vector<int>& fissionProductNuclides() const;
  void setFissionProductNuclides(const std::vector<int>& zaids);

  const std::vector<float>& lightElementNuclideAbundances() const;
  void setLightElementNuclideAbundances(const std::vector<float>& abs);

  const std::vector<float>& actinideNuclideAbundances() const;
  void setActinideNuclideAbundances(const std::vector<float>& abs);

  const std::vector<float>& fissionProductNuclideAbundances() const;
  void setFissionProductNuclideAbundances(const std::vector<float>& abs);

  std::string title() const;
  void setTitle(const std::string& title);

  std::string basis() const;
  void setBasis(const std::string& basis);

  const std::vector<float>& lightElementGammaEnergyBounds() const;
  void setLightElementGammaEnergyBounds(const std::vector<float>& bounds);

  const std::vector<float>& actinideGammaEnergyBounds() const;
  void setActinideGammaEnergyBounds(const std::vector<float>& bounds);

  const std::vector<float>& fissionProductGammaEnergyBounds() const;
  void setFissionProductGammaEnergyBounds(const std::vector<float>& bounds);

  const std::vector<float>& lightElementSpectrum() const;
  void setLightElementSpectrum(const std::vector<float>& spectrum);

  const std::vector<float>& actinitdeSpectrum() const;
  void setActinideSpectrum(const std::vector<float>& spectrum);

  const std::vector<float>& fissionProductSpectrum() const;
  void setFissionProductSpectrum(const std::vector<float>& spectrum);

  const std::vector<float>& gammaEnergyBounds() const;
  void setGammaEnergyBounds(const std::vector<float>& bounds);

  const std::vector<float>& photonPerSSpectrum() const;
  void setPhotonPerSSpectrum(const std::vector<float>& spectrum);

  const std::vector<float>& mevPerSSpectrum() const;
  void setMevPerSSpectrum(const std::vector<float>& spectrum);

  std::string photonSpectrumTitle() const;
  void setPhotonSpectrumTitle(const std::string& title);

  const std::vector<float>& spontaneousFissionNeutronSource() const;
  void setSpontaneousFissionNeutronSource(const std::vector<float>& source);

  const std::vector<float>& actinideAlphaNSource() const;
  void setActinideAlphaNSource(const std::vector<float>& source);

  const std::vector<float>& totalNeutronSpectrum() const;
  void setTotalNeutronSpectrum(const std::vector<float>& spectrum);

  const std::vector<float>& alphaNSpectrum() const;
  void setAlphaNSpectrum(const std::vector<float>& spectrum);

  const std::vector<float>& spontaneousFissionSpectrum() const;
  void setSpontaneousFissionSpectrum(const std::vector<float>& spectrum);

  const std::vector<float>& neutronEnergyBounds() const;
  void setNeutronEnergyBounds(const std::vector<float>& bounds);

  double time() const;
  void setTime(double seconds);

  float power() const;
  void setPower(float power);

  float flux() const;
  void setFlux(float flux);
};  // F71Case
}  // namespace radix

#include "radixio/f71stream.i.hh"
#endif /** RADIX_RADIXIO_F71STREAM_HH_ */