functions.h

Go to the documentation of this file.

#ifndef FUNCTIONS_H
#define FUNCTIONS_H
 
#include "warpxm/warpxm_config.h"
 
#include <cmath>
#include <vector>
 
inline real calculate_collision_coefficient(real nu_p_norm,
                                            real coulomb_log_a,
                                            real A_a,
                                            real Z_a,
                                            real A_b,
                                            real Z_b)
{
 
    // K = (nu_p_tau)*sqrt(2)*Z_a*Z_a*Z_b*Z_b * coulomb_log * (A_a+A_b)/(A_a*A_a*A_b*A_b)
 
    return nu_p_norm * sqrt(2) * Z_a * Z_a * Z_b * Z_b * coulomb_log_a * (A_a + A_b) /
           (A_a * A_a * A_b * A_b);
}
 
inline real calculate_interspecies_collision_frequency(real collision_coefficient,
                                                       real minimum,
                                                       real A_a,
                                                       real A_b,
                                                       real theta_a,
                                                       real rho_b,
                                                       real theta_b)
{
    // theta_a = T_a/A_a = pressure_a/rho_a
    // nu_a_b = K * rho_b/(theta_a+theta_b)^(3/2)
    return fmax(minimum, collision_coefficient * rho_b / pow(theta_a + theta_b, 1.5));
}
 
inline real calculate_intraspecies_collision_frequency(real collision_coefficient,
                                                       real minimum,
                                                       real rho,
                                                       real theta)
{
    // here a=b
    // same calculation coefficient (K)
    // rho_b=rho
    // theta_a=theta_b=theta
    return fmax(minimum, collision_coefficient * rho / pow(2.0 * theta, 1.5));
}
 
void mat_mult(const std::vector<std::vector<real>>& a,
              const std::vector<std::vector<real>>& b,
              std::vector<std::vector<real>>& c);
 
 
real heaviside(real x);
 
#endif // FUNCTIONS_H