view src/samer/maths/random/GeneralisedLaplacian.java @ 8:5e3cbbf173aa tip

Reorganise some more

/*
 *	Copyright (c) 2000, Samer Abdallah, King's College London.
 *	All rights reserved.
 *
 *	This software is provided AS iS and WITHOUT ANY WARRANTY; 
 *	without even the implied warranty of MERCHANTABILITY or 
 *	FITNESS FOR A PARTICULAR PURPOSE.
 */

package samer.maths.random;
import  samer.core.*;
import  samer.core.types.*;

// This is going to work by the rejection method
// (see Numerical Recipes in C)
// We'll use the cauchy distribution as our
// bounding function

// NB ONLY WORKS FOR alpha<=2!

public class GeneralisedLaplacian extends Cauchy
{
	VDouble	alpha = new VDouble("alpha", 1);

	public void dispose() { alpha.dispose(); }
	public double next()
	{ 
		double x, f, g, a=alpha.value;
		do {
			// first, get a Cauchy distributed RV
			x = 2*super.next();
			f = 4/(4+x*x); // cauchy
			g = Math.exp(-Math.pow(Math.abs(x),a)); // exponential
		} while (f*rnd.nextDouble()>g);
		return x; 
	}
}