#ifndef _GAMMALIKE_H_
#define _GAMMALIKE_H_
// gamma deviated mutation rate among loci
//
//
// code for the calculation of the gamma likelihood
// and for its derivatives
#include "migration.h"

extern MYREAL gamma_loci_like (helper_fmt * helper, MYREAL *oparam,
                                   MYREAL *olparam, MYREAL denom);
extern MYREAL gamma_locus_like (nr_fmt * nr, MYREAL *oparam, MYREAL *olparam,
                                    MYREAL denom, long locus);
extern void gamma_loci_derivative (helper_fmt * helper);
extern void gamma_locus_derivative (helper_fmt * helper, long locus);
extern void gamma_loci_difference (helper_fmt * helper);
#endif