physlib.h

Go to the documentation of this file.

#ifndef PHYSLIB_H
#define PHYSLIB_H
 
#include <math.h>
#include "math.h"
#include "consts.h"

#define physlib_gamma_vnorm(v) (                                        \
        sqrt( 1.0 / ( (1.0 - v / CONST_C) * (1.0 + v / CONST_C) ) ) )

#define physlib_vnorm_gamma(gamma) (                                   \
        sqrt( 1.0 - 1.0 / ( gamma * gamma ) ) * CONST_C )

#define physlib_gamma_pnorm(m, p) (                  \
        sqrt(1.0 + ( p * p ) / ( m * m * CONST_C2 ) ) )

#define physlib_gamma_vpar(m, mu, vpar, B) (                            \
        sqrt( ( 1.0 + (2.0 * mu * B) / ( m * CONST_C2 ) ) /             \
              ( (1.0 - vpar / CONST_C) * (1.0 + vpar / CONST_C) ) ) )

#define physlib_gamma_ppar(m, mu, ppar, B) (             \
        sqrt( 1.0 + 2 * mu * B / ( m * CONST_C2 ) +      \
              ppar * ppar / ( m * m * CONST_C2 ) ) )

#define physlib_Ekin_gamma(m, gamma) ( ( gamma - 1.0 ) * m * CONST_C2 )

#define physlib_gamma_Ekin(m, ekin) ( ekin / ( m * CONST_C2 ) + 1.0 )

#define physlib_Ekin_pnorm(m, p) (                              \
        ( physlib_gamma_pnorm(m, p) - 1.0 ) * m * CONST_C2 )

#define physlib_Ekin_ppar(m, mu, ppar, B) (                             \
        ( physlib_gamma_ppar(m, mu, ppar, B) - 1.0 ) * m * CONST_C2 )

#define physlib_vnorm_pnorm(m, p) (                  \
        p / sqrt(m * m + ( p * p ) / CONST_C2 ) )

#define physlib_pnorm_vnorm(m, v) ( m * v * physlib_gamma_vnorm(v) )

#define physlib_gc_ppar(p, xi) ( p * xi )

#define physlib_gc_mu(m, p, xi, B) (                    \
        p * p * ( 1.0 - xi * xi ) / ( 2 * B * m ) )

#define physlib_gc_p(m, mu, ppar, B) ( \
        m * CONST_C * sqrt( pow( physlib_gamma_ppar(m, mu, ppar, B), 2) - 1 ) )

#define physlib_gc_xi(m, mu, ppar, B) (                                 \
        ppar / physlib_gc_p(m, mu, ppar, B) )

#define physlib_gyrolength_p(q, p, B) (                              \
        math_dot(p, B) / ( fabs(q) * math_dot(B, B) ) )

#define phys_gyrolength_ppar(m, q, mu, ppar, B) (                   \
        sqrt( 2 * m * mu *                                          \
              physlib_gamma_ppar(m, mu, ppar, B) / B ) / fabs(q) )

#define phys_gyrofreq_pnorm(m, q, p, B) (               \
        fabs(q) * B / ( m * physlib_gamma_pnorm(m, p) ) )

#define phys_gyrofreq_ppar(m, q, mu, ppar, B) (                 \
        fabs(q) * B / ( m * physlib_gamma_ppar(m, mu, ppar, B) ) )

#define phys_ptoroid_fo(q, R, pphi, psi) ( \
        R * pphi  +  q * psi )

#define phys_ptoroid_gc(q, R, ppar, psi, B, Bphi) (             \
        ppar * R * (Bphi / B)  +  q * psi )

#define phys_vperp_gc(v,vpar) sqrt(pow(v,2) - pow(vpar,2))

#define phys_pperp_gc(p,ppar) sqrt(pow(p,2) - pow(ppar,2))

#define phys_ppar_Ekin(m, ekin, mu, B) (sqrt((physlib_gamma_Ekin(m, ekin)*\
        physlib_gamma_Ekin(m, ekin) - 1.0)*m*m*CONST_C2 - 2.0*m*mu*B))

#define phys_vperp_mu(m, mu, B) sqrt(2*mu*B/m)

#define phys_ppar_mu(m, mu, B, q, gyrof) (m*CONST_C*sqrt((q*B/(gyrof*m)) * \
        (q*B/(gyrof*m)) - 2*mu*B/(m*CONST_C2) - 1))
 
 
// ISSUE: SEEMS TO CONTAIN A BUG!
#define phys_ppar_pphi(B, R, B_phi, p_phi, q, psi) (B/(R*B_phi)*(p_phi - q*psi))

#define phys_ald_force_chartime(q, m, B, gamma) ( (q*q*q*q) * (B*B) / \
        (6.0 * CONST_PI * CONST_E0 * gamma * (m*m*m) * CONST_C3 ) )
 
#endif