more/math/xtmath.h Source File

00001 #ifndef XMATH_H_INC
00002 #define XMATH_H_INC
00003 
00004 /*   C o n s t a n t s   */
00005 
00006 #define c_e             2.7182818284590452354
00007 #define c_e_log2        1.4426950408889634074
00008 #define c_e_log10       0.43429448190325182765
00009 #define c_2_ln          0.69314718055994530942
00010 #define c_10_ln         2.30258509299404568402
00011 #define c_2pi           6.28318530717958647693
00012 #define c_pi            3.14159265358979323846
00013 #define c_pi_2_f        1.57079632679489661923
00014 #define c_pi_4_f        0.78539816339744830962
00015 #define c_1_pi_f        0.31830988618379067154
00016 #define c_2_pi_f        0.63661977236758134308
00017 #define c_2_pi_sqrt_f   1.12837916709551257390
00018 #define c_2_sqrt        1.41421356237309504880
00019 #define c_1_2_sqrt_f    0.70710678118654752440
00020 #define c_pi_sqrt       1.77245385090551602729817
00021 
00022 
00023 /*   F a c t o r i a l s   */
00024 
00025 #define MAX_FACTAB 100
00026 /* fac.c */
00027 /* #define f1c(n) (((n)<80)? __fac[n] : __f1c(n)) */
00028 double f1c(int n);  /* fac.c  n! */
00029 double fhc(int n);  /* fac.c  (n+1/2)! */
00030 double f2c(int n);  /* fac.c  n!! */
00031 double f1h(int n);  /* fac.c  (n/2)! */
00032 /* double __fac[], __oddfac[]; */
00033 
00034 
00035 /*   S p e c i a l   F u n c t i o n s   */
00036 
00037 #if 0
00038 double fulaglh(double x, int n, int l);
00039 double fusphosc(double x, int n, int l);
00040 double fubesjim(double x, int l);  /* j_l(ix)/i^l   */
00041 double fubeshim(double x, int l);  /* -h_l(ix)/i^l */
00042 double fulegp(double x, int l);    /* P[l](x)       */
00043 #endif
00044 
00045 double besjr(int, double);
00046 double besIr(double, double);
00047 double besir(int, double);
00048 double lagLc(int, int);
00049 int lagLmesh(int, double *, double *);
00050 
00051 
00052 /*   A n g u l a r   M o m e n t u m   A l g e b r a   */
00053 
00054 /* The Racha 9j and 6j expressions. */
00055 double sixj_2j(int, int, int, int, int, int);                  /* ang.c */
00056 double ninej_2j(int, int, int, int, int, int, int, int, int);  /* ang.c */
00057 /* The Clebsch-Gordan coefficients and their 3jm symmertic form */
00058 double clebsch_2j(int, int, int, int, int, int);                    /* clebsch.c */
00059 double threejm_2j(int, int, int, int, int, int);            /* clebsch.c */
00060 
00061 #endif /* XMATH_H_INC */