/*CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C                                                                             C
C  COMPUTE MERTENS FUNCTION (using Deleglise and Rivat's algorithm)           C
C  10/30/15 (DKC)							      C
C                                                                             C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC*/
#include <math.h>
double newlocl(unsigned long long x, unsigned int u, char *mob, int *M);
double newhicl(unsigned long long x, unsigned int u, char *mob, int *M);
int newrivl(unsigned long long x, unsigned int u, char *mobb, int *M) {
double f1,f2,f3;
f1=newlocl(x,u,mobb,M);
f2=newhicl(x,u,mobb,M);
f3=(double)M[u-1]-f1-f2;
return((int)f3);
}