package tim.prune.function.estimate.jama;

/**
 * Static helper method, taken from public domain NIST code for JAMA
 */
public abstract class Maths
{
	/**
	 * Work out sqrt(a^2 + b^2)
	 */
	public static double pythag(double a, double b)
	{
		double r;

		if (Math.abs(a) > Math.abs(b))
		{
			r = b/a;
			r = Math.abs(a)*Math.sqrt(1+r*r);
		}
		else if (b != 0)
		{
			r = a/b;
			r = Math.abs(b)*Math.sqrt(1+r*r);
		}
		else
		{
			r = 0.0;
		}
		return r;
	}
}