annotate src/samer/functions/BiLaplacian.java @ 0:bf79fb79ee13

Initial Mercurial check in.

/*
 *	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.functions;
import  samer.maths.*;

/**
Sparsified Laplacian, v1
	Log prior for an laplacian distribution with
	an extra mass at zero, eg a mixture of a laplacian
	with another much narrower laplacian. The log prior is
	basically an abs function with with a extra-deep pit at
	zero.
	Paramters control the width and depth of the pit.
  */

public class BiLaplacian extends Function
{
	double width=1, A=1, slope=2;

	public void setWidth(double b) { width=b; slope=(width+A)/width; }
	public void setStrength(double a) { A=a; slope=(width+A)/width; }
	public double getJump() { return slope; }

	public String format(String t) { return "BiLaplace("+t+")"; }

	public final double apply(double t) { 
		t = Math.abs(t);
		if  (t<width) return -A+slope*t;
		else return t;
	}

	public Function derivative() { 
		return new Function() {
			public String format(String t) { return "dBiLaplace("+t+")"; }
			public final double apply(double t) { 
				if (t>0) {
					if  (t<width) return slope;
					else return 1;
				} else if (t<0) {
					if  (t>-width) return -slope;
					else return -1;
				} else return 0;
			}
		};
	}
}