input_module.h

Go to the documentation of this file.

 
#ifndef __INPUT__
#define __INPUT__
 
#include <memory>
#include <string>
#include <vector>
 
#include "../species/base_species.h"
#include "../species/species_build_context.h"
#include "../species/species_collection.h"
#include "background.h"
#include "common.h"
#include "lensing.h"
#include "nonlinear.h"
#include "output.h"
#include "parser.h"
#include "perturbations.h"
#include "primordial.h"
#include "quadrature.h"
#include "spectra.h"
#include "thermodynamics.h"
#include "transfer.h"
 
class InputModule {
 public:
  InputModule(FileContent& fc);
  static void file_content_from_arguments(int argc, char** argv, FileContent& fc);
 
  static InputModulePtr DoShooting(InputModulePtr input_module);
 
  FileContent& file_content_;
  precision precision_;
  background background_;
  thermo thermodynamics_;
  perturbs perturbations_;
  transfers transfers_;
  primordial primordial_;
  spectra spectra_;
  nonlinear nonlinear_;
  lensing lensing_;
  output output_;
  SpeciesCollection all_species_;
  SpeciesOmegaBudget omega_budget_;
  CoupledClusterInputs coupled_inputs_;
 
 private:
  // Hook-based shooting (used by DoShooting). theta_s is identified by
  // target_name == "100*theta_s" / unknown_param == "h"; all targets are scalar.
  struct ShootingWorkspace {
    explicit ShootingWorkspace(FileContent& fc_ref) : fc(fc_ref) {}
    FileContent& fc;
    std::vector<ShootingTarget> targets;  // theta_s first (if present), then species in lex order
    // Parallel to targets: the all_species_ key owning each target ("" for the module-level
    // theta_s). Lets ShootingResidual route each slot to its species' ComputeShootingResidual
    // with the authoritative target (so the species never re-derives it from the file content).
    std::vector<std::string> target_species_keys;
  };
  static void ShootingResidual(double* x, int x_size, void* pworkspace, double* output);
 
  void ConstructSpecies();
 
  void ReadCoupledCluster();
 
  void ReadContext();          // phase i: inputs needed to build species
  void ReadDerived();          // phase iii: everything else + species-dependent reads
  void WriteParameterFiles();  // read/unread parameter dump (runs after ReadDerived)
  void input_read_precisions();
};
 
#endif