00001 00008 #include "../Matrix/Matrix.h" 00009 #include "variousConstants.h" 00010 00012 namespace VF { 00013 00014 double G(double x); 00015 //-------------------------------------------------------------------------------------------------------- 00016 // Computation of dipole magnetic field 00017 //-------------------------------------------------------------------------------------------------------- 00018 double B (double Lparam); 00019 double Df (double L, double Kp); 00020 //-------------------------------------------------------------------------------------------------------- 00021 // Computation of moumentum from Kinetic energy 00022 //-------------------------------------------------------------------------------------------------------- 00023 double pfunc(double K);//, double mc2); 00024 //-------------------------------------------------------------------------------------------------------- 00025 // Computation of Kinetic energy from given momentum 00026 //-------------------------------------------------------------------------------------------------------- 00027 double Kfunc(double pc);// , double mc2 = 0.511); 00028 double bounce_time(double E, double L); 00029 00030 //-------------------------------------------------------------------------------------------------------- 00031 // Specifying outer boundary 00032 //-------------------------------------------------------------------------------------------------------- 00033 // double J_L7 (double K, double x1, double y1, double x2, double y2); 00034 double J_L7_corrected (double K); 00035 double J_L7 (double K, double x1, double y1, double x2, double y2); 00036 00037 //-------------------------------------------------------------------------------------------------------- 00038 // Converting mu to pc 00039 //-------------------------------------------------------------------------------------------------------- 00040 double mu2pc (double L, double mu, double alpha); 00041 double pc2mu (double L, double pc, double alpha); 00042 00043 double alc(double L); 00044 00045 double f_interp(double E, double f1, double E1, double f2, double E2); 00046 00047 00048 double Jc2Alpha (double L, double pc, double Jc); 00049 double Alpha2Jc (double L, double pc, double Alpha); 00050 00051 double Y(double alpha); 00052 00053 double mu_calc(double L, double pc, double Alpha); 00054 double Jc_calc(double L, double pc, double Alpha); 00055 00056 //double density(double L); 00057 //Matrix3D<double> density(Matrix3D<double> &L); 00058 00059 double max(double v1, double v2); 00060 00061 string dtostr(double n); 00062 }