view xthresh.c @ 584:e3790284fd4a

Merged through librdf storage hooks and apache2 module.

Squashed commit of the following:

commit a6cfca8f04036e12e7d7fcd55c47224e802582f0
Author: Michael Jewell <mjewell@harrison.(none)>
Date:   Fri Jul 31 15:23:32 2009 +0100

    Removed leftover bits and bobs.

commit f1f0dd074d0767de3e24ba636779fd8701d73d9e
Author: Michael Jewell <mjewell@harrison.(none)>
Date:   Fri Jul 31 15:07:20 2009 +0100

    Simple test of database creation via librdf.

commit 90e6350538e004d8785137e5ff2ac878c22a5d42
Author: Michael Jewell <mjewell@harrison.(none)>
Date:   Fri Jul 31 15:05:10 2009 +0100

    Added the apache2 module which will hook into the librdf storage module

commit c75bf53763b7078c83ae97fcf247da2576baa79a
Author: Michael Jewell <mjewell@harrison.(none)>
Date:   Fri Jul 31 15:04:53 2009 +0100

    Added sparql librdf source - requires the librdf sources to compile.

commit 0646f0190112a73ddb2533537e2cc9832c066b52
Author: Michael Jewell <mjewell@harrison.(none)>
Date:   Mon Jul 27 12:12:26 2009 +0100

    Adding execution to mod_audiodb

commit 8f83f27ba4d917278bca0c7cb665d930e28c86df
Author: Michael Jewell <mjewell@harrison.(none)>
Date:   Wed Jul 22 12:15:57 2009 +0100

    Some initial returns for the sparql handler.

commit dc639aed11943a5b0c379eb47cf293f76908b1b7
Author: Michael Jewell <mjewell@harrison.(none)>
Date:   Wed Jul 22 12:06:20 2009 +0100

    Added a little setup.sh script to do libtoolize/autoconf etc.

commit 3a679da499db647fc82cf2797daeb5cc44ed7655
Author: Michael Jewell <mjewell@harrison.(none)>
Date:   Wed Jul 22 12:03:42 2009 +0100

    Adding initial bits for apache mod

#include <gsl/gsl_sf.h>
#include <stdio.h>
#include <stdlib.h>
#include <math.h>

int main(int argc, char *argv[]) {
  if(argc != 4) {
    fprintf(stderr, "Wrong number of arguments: %d\n", argc);
    exit(1);
  }

  long int meanN = strtol(argv[1], NULL, 10);

  double d = strtod(argv[2], NULL);
  double sigma2 = strtod(argv[3], NULL);

  double logw = (2 / d) * gsl_sf_log(-gsl_sf_log(0.99));
  double logxthresh = gsl_sf_log(sigma2) + logw
    - (2 / d) * gsl_sf_log(meanN)
    - gsl_sf_log(d/2)
    - (2 / d) * gsl_sf_log(2 / d)
    + (2 / d) * gsl_sf_lngamma(d / 2);

  printf("w: %f\n", exp(logw));
  printf("x_thresh: %f\n", exp(logxthresh));
  exit(0);
}