Return to ko-perturbation.c CVS log

Up to [local] / OpenXM / src / hgm / fisher-bingham / src

A manual and a library for nk_fb_ten_c.
nk_fb_gen_c-ja.texi generated from nk_fb_gen_c.oxg gives an explanation
how to make a maximal likelihood estimate (MLE) for the Fisher-Bingham distribution
with this package.

(commit by proxy of T.Koyama, H.Nakayama, K.Nishiyama).

#include <stdio.h>
#include <stdlib.h>
#include <gsl/gsl_rng.h>
#include <gsl/gsl_randist.h>

static int init=1;
static const gsl_rng_type *T;
static gsl_rng * r;
static int length;
static double eps;

double
*perturbation(const double v[], double err[])
{
  if(init){
    fprintf(stderr, "perturbation: erorr! Forget perturbation_init?\n");
    exit(1);
  }

  int i;
  double *vp = (double *) malloc(sizeof(double) * length);
  for(i=0; i<length; i++){
    err[i]= gsl_ran_gaussian (r, eps);
    vp[i] = v[i]+err[i];
  }
  return vp;
}

void 
perturbation_init(int l,  double e)
{
  if(init){
    gsl_rng_env_setup();     
    T = gsl_rng_default;
    r = gsl_rng_alloc (T);
    init = 0;
  } 
  length = l;
  eps = e;
}

void
perturbation_free(void)
{
  gsl_rng_free (r);
}