Mercurial > repos > lsong10 > psiclass

comparison PsiCLASS-1.0.2/stats.hpp @ 0:903fc43d6227 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author lsong10
date Fri, 26 Mar 2021 16:52:45 +0000
parents
children
comparison
equal deleted inserted replaced
-1:000000000000 0:903fc43d6227
1 #ifndef _LSONG_STATS
2 #define _LSONG_STATS
3
4 #include <stdio.h>
5 #include <math.h>
6
7 #ifndef M_PIl
8 /** The constant Pi in high precision */
9 #define M_PIl 3.1415926535897932384626433832795029L
10 #endif
11 #ifndef M_GAMMAl
12 /** Euler's constant in high precision */
13 #define M_GAMMAl 0.5772156649015328606065120900824024L
14 #endif
15 #ifndef M_LN2l
16 /** the natural logarithm of 2 in high precision */
17 #define M_LN2l 0.6931471805599453094172321214581766L
18 #endif
19
20 /** The digamma function in long double precision.
21 * @param x the real value of the argument
22 * @return the value of the digamma (psi) function at that point
23 * @author Richard J. Mathar
24 * @since 2005-11-24
25 */
26 long double digammal(long double x) ;
27
28 //https://people.sc.fsu.edu/~jburkardt/cpp_src/asa121/asa121.hpp
29 double trigamma ( double x, int *ifault );
30
31
32 double LogGammaDensity( double x, double k, double theta ) ;
33 double MixtureGammaAssignment( double x, double pi, double* k, double *theta ) ;
34
35
36 // http://people.sc.fsu.edu/~jburkardt/c_src/asa091/asa091.hpp
37 double alnorm ( double x, bool upper );
38 double gammad ( double x, double p, int *ifault );
39 double r8_min ( double x, double y );
40
41 double chicdf( double x, double df ) ;
42
43 #endif