comparison 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;
+			}
+		};
+	}
+}
+