distcosmos.h

Go to the documentation of this file.

#ifndef dist_cosmos_H
#define dist_cosmos_H

#include <cmath>
#include <iostream>

#include "basecosmos.h"

#include "spline.h"
#include "anchor.h"
#include "miscmath.h"

class distCosmos : public baseCosmos {
protected:
  static const double flattol; 
  static const int ndist; 
  double zmax; 

  double h_h; 

  double Omega_m; 
  double Omega_v; 
  double Omega_rel; 
  double Omega_k; 

  Miscmath::keyword odetype; 

  Anchor SplineAnchor; 
  Spline* distSpline; 
  Spline* dDdzSpline; 
  Spline* Z2iHSpline; 

  virtual bool calc_isrebounding(); 

  void iEode(double, const double*, double*) const; 

public :
  
  distCosmos(); 
  virtual ~distCosmos() {};

  void reset(); 
  virtual void initSplines(); 

  //Scale factor stuff
  double a_0() const {return 1.0;}  
  
  //Hubble constant stuff
  double h() const { return h_h; } 
  double h2() const { return h_h * h_h; } 

  virtual double E(double z) const; 

  double ZMax() const { return zmax; } 

  //Omegas
  double omega_m() { return Omega_m; } 
  double omega_v() { return Omega_v; } 
  double omega_relativistic() { return Omega_rel; } 
  virtual double omega_0() { return omega_m()+omega_v()+omega_relativistic(); } 

  //Current Densities
  double rho_0() const { return  M_p()*M_p()*h2()* 3.3379e-7;  } ; 
  double rho_m() { return rho_0() * omega_m(); } 
  double rho_relativistic() { return rho_0() * omega_relativistic(); } 
  double rho_v() { return rho_0() * omega_v(); } 

  //Initial Densities -- or, really, density at a
  double initialRho_m(const double a) { return rho_m()*pow(a/a_0(),-3); } 
  double initialRho_relativistic(const double a) { return rho_relativistic()*pow(a/a_0(),-4); } 
  double initialRho_v(const double a) { return rho_v(); } 

  double Z2iH(double z) { return Z2iHSpline->fastY(z); }  

  //Distances
  double propermotionDistance(double z) const;
  double angulardiameterDistance(double z) const;
  double luminosityDistance(double z) const;
  double luminosityDistance(double zhel, double zcmb) const;
  
  //Derivatives with distance
  double propermotionDistanceDeriv(double z) const;
  double angulardiameterDistanceDeriv(double z) const;
  double luminosityDistanceDeriv(double z) const;

  //Set various things
  void seth(double h) { h_h = h; setOmegaH2_rel( 2.471e-5 ); }

  void setOmega_m(double x) { Omega_m = x; Omega_k=1.0-omega_0();} 
  void setOmega_v(double x) { Omega_v = x; Omega_k=1.0-omega_0();} 
  void setOmega_relativistic(double x) { Omega_rel = x; Omega_k=1.0-omega_0(); } 
  void setOmegaH2_rel(double x) { Omega_rel = x / h2(); Omega_k=1.0-omega_0(); } 

  void setZMax(double x) { zmax = x; if (validHistory()) reset(); } 

  virtual void history(bool inform=false); 

};

#endif

---