dust.h

#ifndef COSMO_DUST_SIM_H
#define COSMO_DUST_SIM_H

#include "sim.h"
#include "../components/dust_fluid/dust.h"
#include "../components/Lambda/lambda.h"
#include "../components/phase_space_sheet/sheets.h"

namespace cosmo
{

class DustSim : public CosmoSim
{
protected:
  Sheet * raySheet;
  Dust * dustSim;
  Lambda * lambda;
  bool take_ray_step;
  idx_t raysheet_flip_step;
public:
  DustSim();
  ~DustSim()
  {
    std::cout << "Cleaning up...";
    delete iodata;
    delete bssnSim;
    delete fourier;
    if(use_bardeen)
    {
      delete bardeen;
    }
    std::cout << "done.\n";
    std::cout << std::flush;
  }

  void init();
  void setICs();
  void initDustStep();
  void outputDustStep();
  void runDustStep();
  void runStep();
};

} /* namespace cosmo */

#endif