tabulated_dynamicalDE_EoS.h

Go to the documentation of this file.

 
#pragma once
 
#include <fstream>
#include <iostream>
#include <string>
#include <vector>
 
#include "../global/global.h"
#include "../utils/error_handling.h"
 
class TabulatedDynamicalDarkEnergyEoS
{
 private:
  std::vector<float> dynamicalDE_table_z;
  std::vector<float> dynamicalDE_table_w;
 
  // dynamical dark energy energy density normalized to the present-date rho_DE(z) / rho_DE(z=0)
  std::vector<float> dynamicalDE_table_density;
 
  void Setup_DynamicalDE_EquationOfState_(std::istream& in, const std::string& fname, bool silent);
 
  void Set_DynamicalDE_Density();
 
  void Setup_Full(std::istream& in, const std::string& fname) {}
 
 public:
  Real Get_DynamicalDE_Density_from_a(Real a);
 
  explicit TabulatedDynamicalDarkEnergyEoS(const std::string& path, std::istream* f = nullptr, bool silent = false)
  {
    std::fstream tmp;
    if (f == nullptr) {
      tmp.open(path, std::ios_base::in);
      CHOLLA_ASSERT(tmp.is_open(), "Unable to open DE equation of state file: %s\n", path.c_str());
      f = dynamic_cast<std::istream*>(&tmp);
    }
    Setup_DynamicalDE_EquationOfState_(*f, path, silent);
    Set_DynamicalDE_Density();
  }
};