gaia.cu9.tools.parallax.util

Class CdfIntegration

java.lang.Object

gaia.cu9.tools.parallax.util.CdfIntegration

public class CdfIntegration
extends Object

Constructor Summary

Constructors

Constructor and Description
CdfIntegration()

Method Summary

Modifier and Type	Method and Description
double[][]	getCdf(PDF pdf, boolean logAxis, int nPoints, double rMin, double rMax, boolean normalisePdf, boolean integrateToInfinite) Calculate the CDF by integrating the PDF in a finite interval.
double	getDeltaThreshold()
int	getIterationLimit()
double[]	getPercentiles(double[][] cdf, double... percentiles)
protected double	integrateBooleRule(PDF pdf, double min, double max) Integrate the PDF over a finite interval using Boole's rule
protected double	integrateBooleRule(PDF pdf, double min, double max, double pdfAtMin, double pdfAtMax) Integrate the PDF over a finite interval using Boole's rule This implementation uses the value of the PDF at the boundaries that has already been calculated externally, thus avoiding unnecesary calls to PDF.getUnnormalizedProbabilityAt(double)
protected double	integrateToInfinitum(PDF pdf, boolean logAxis, double minR, double step)
void	setDeltaThreshold(double deltaThreshold)
void	setIterationLimit(int iterationLimit)

Methods inherited from class java.lang.Object

clone, equals, finalize, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait

Constructor Detail

CdfIntegration

public CdfIntegration()

Method Detail

getCdf

public double[][] getCdf(PDF pdf,
                         boolean logAxis,
                         int nPoints,
                         double rMin,
                         double rMax,
                         boolean normalisePdf,
                         boolean integrateToInfinite)

Calculate the CDF by integrating the PDF in a finite interval. It is assumed the the PDF below the minimum value is 0. The algorithm will continue integrating towards infinity until the increment added is lower than a threshold or a limit of iterations is reached so that the CDF can be normalized To set these thresholds please see setDeltaThreshold(double) and setIterationLimit(int)

Parameters:

pdf - Probability Distribution Function

logAxis - True if the points should be distributed logaritmically, false otherwise

nPoints - Number of points to be calculated in the range

rMin - Minimum value of the range, in parsecs

rMax - Maximum value of the range, in parsecs

Returns:

CDF

integrateBooleRule

protected double integrateBooleRule(PDF pdf,
                                    double min,
                                    double max)

Integrate the PDF over a finite interval using Boole's rule

Parameters:

pdf - probability distribution function to integrate

min - Minimum value of x (in linear units)

max - Maximum value of x (in linear units)

Returns:

Value of the integral of the pdf in the interval (min,max)

integrateBooleRule

protected double integrateBooleRule(PDF pdf,
                                    double min,
                                    double max,
                                    double pdfAtMin,
                                    double pdfAtMax)

Integrate the PDF over a finite interval using Boole's rule This implementation uses the value of the PDF at the boundaries that has already been calculated externally, thus avoiding unnecesary calls to PDF.getUnnormalizedProbabilityAt(double)

Parameters:

pdf - probability distribution function to integrate

min - Minimum value of x (in linear units)

max - Maximum value of x (in linear units)

pdfAtMin - Value of the PDF at min

pdfAtMax - Value of the PDF at max

Returns:

Value of the integral of the pdf in the interval (min,max)

integrateToInfinitum

protected double integrateToInfinitum(PDF pdf,
                                      boolean logAxis,
                                      double minR,
                                      double step)

getPercentiles

public double[] getPercentiles(double[][] cdf,
                               double... percentiles)

getDeltaThreshold

public double getDeltaThreshold()

setDeltaThreshold

public void setDeltaThreshold(double deltaThreshold)

getIterationLimit

public int getIterationLimit()

setIterationLimit

public void setIterationLimit(int iterationLimit)