The R Project SVN – R – Diff – /trunk/src/nmath/pgamma.c



# define LOG logl
#else
# define EXP exp
# define LOG log
#endif
double logspace_sum (const double* logx, int n)
{
    if(n == 0) return ML_NEGINF; // = log( sum(<empty>) )
    if(n == 1) return logx[0];
    if(n == 2) return logspace_add(logx[0], logx[1]);
    // else (n >= 3) :

Rev 68947	Rev 69008
Line 257...	Line 257...
257	`# define LOG logl`	257	`# define LOG logl`
258	`#else`	258	`#else`
259	`# define EXP exp`	259	`# define EXP exp`
260	`# define LOG log`	260	`# define LOG log`
261	`#endif`	261	`#endif`
262	`double logspace_sum (double* logx, int n)`	262	`double logspace_sum (const double* logx, int n)`
263	`{`	263	`{`
264	`if(n == 0) return ML_NEGINF; // = log( sum(<empty>) )`	264	`if(n == 0) return ML_NEGINF; // = log( sum(<empty>) )`
265	`if(n == 1) return logx[0];`	265	`if(n == 1) return logx[0];`
266	`if(n == 2) return logspace_add(logx[0], logx[1]);`	266	`if(n == 2) return logspace_add(logx[0], logx[1]);`
267	`// else (n >= 3) :`	267	`// else (n >= 3) :`

The R Project SVN R

(root)/trunk/src/nmath/pgamma.c – Rev 68947 → 69008