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 )
 
#endif