core_app.h

Go to the documentation of this file.

#pragma once
#include "stdafx.h"
#include "app_var.hpp"
#include "core_power.h"
#include "params.hpp"
#include "precision.hpp"
#include "class_particles.hpp"

void set_unpert_pos(const Sim_Param &sim, std::vector<Particle_x<FTYPE_t>>& particles);
void set_unpert_pos_w_vel(const Sim_Param &sim, std::vector<Particle_v<FTYPE_t>>& particles, const std::vector< Mesh> &vel_field);
void set_pert_pos(const Sim_Param &sim, const FTYPE_t a, std::vector<Particle_x<FTYPE_t>>& particles, const std::vector< Mesh> &vel_field);
void set_pert_pos(const Sim_Param &sim, const FTYPE_t a, std::vector<Particle_v<FTYPE_t>>& particles, const std::vector< Mesh> &vel_field);

void gen_rho_dist_k(const Sim_Param &sim, Mesh& rho, const FFTW_PLAN_TYPE &p_F);
void gen_pot_k(const Mesh& rho_k, Mesh& pot_k);
void gen_pot_k(Mesh& rho_k);
void gen_displ_k(std::vector<Mesh>& vel_field, const Mesh& pot_k);
void gen_displ_k_cic(std::vector<Mesh>& vel_field, const Mesh& pot_k);
void gen_displ_k_S2(std::vector<Mesh>& vel_field, const Mesh& pot_k, const FTYPE_t a);

template <class T>
void get_rho_from_par(const std::vector<T>& particles, Mesh& rho, const Sim_Param &sim);
bool get_vel_from_par(const std::vector<Particle_v<FTYPE_t>>& particles, std::vector<Mesh>& vel_field, const Sim_Param &sim);
bool get_vel_from_par(const std::vector<Particle_x<FTYPE_t>>& particles, std::vector<Mesh>& vel_field, const Sim_Param &sim);

void pwr_spec_k(const Mesh &rho_k, Mesh& power_aux);
void pwr_spec_k_init(const Mesh &rho_k, Mesh& power_aux);
void vel_pwr_spec_k(const std::vector<Mesh> &vel_field, Mesh& power_aux);
void gen_pow_spec_binned(const Sim_Param &sim, const Mesh &power_aux, Data_Vec<FTYPE_t, 2>& pwr_spec_binned);
void gen_pow_spec_binned_init(const Sim_Param &sim, const Mesh &power_aux, const size_t half_length, Data_Vec<FTYPE_t, 2>& pwr_spec_binned);
template<class P, typename T, size_t N> // P = everything callable P_k(k), T = float-type, N = number
void gen_pow_spec_binned_from_extrap(const Sim_Param &sim, const P &P_k, Data_Vec<T, N>& pwr_spec_binned);
void gen_dens_binned(const Mesh& rho, std::vector<size_t> &dens_binned, const Sim_Param &sim);