systemic
Class Kernel

java.lang.Object
  extended by systemic.Kernel
All Implemented Interfaces:
java.io.Serializable

public class Kernel
extends java.lang.Object
implements java.io.Serializable

See Also:
Serialized Form

Field Summary
static int OESIZE
          Number of orbital elements per planet.
 java.lang.Object tagData
           
 
Constructor Summary
Kernel()
           
 
Method Summary
 int addPlanet(double[] plels)
           
 void addRVDataSet(double[][] rvSet, java.lang.String name)
          Adds a new RV dataset
 void addRVDataSet(java.lang.String fileName)
          Adds a new RV dataset from ASCII table fileName.
 boolean calculate()
          Calculates Chi^2 and residuals for the current fit, given the loaded data.
 void clearRVDataSets()
          Removes all RV datasets.
 Kernel clone()
          Creates a new, independent copy of the object.
 KernelList errorEst(int method, int steps, boolean[][] icFit, boolean[] voFit, java.lang.Object[] opts)
          Launches one of the error estimation algorithms.
static double[] ftest(Kernel k1, Kernel k2, int what)
          Computes the result of f-test test for two models, represented by kernel objects k1 and k2 (see Numerical Recipes).
 double[][] getAltElements()
          Returns an matrix containing extra orbital elements of the system (K, n, a, Ptr, Tperi) for each planet.
 double getChiSqr()
          Returns the combined (RV + transits) χ^2 (call after calculate()).
 double[] getChiSqrAll()
          Returns an array of χ^2 statistics.
 boolean getCloseEncounter()
          Returns true if a close encounter happened during the last integration.
 double[][] getCompiledData()
          Returns a table containing all the RV datasets compiled with their respective velocity offsets and sorted by JD.
 java.lang.String[] getDataNames()
          Returns a list of all the loaded RV datasets names.
 java.util.ArrayList<double[][]> getDataSets()
          Returns a list of all the loaded RV datasets.
 java.util.ArrayList<double[][]> getDataSetsCopy()
          Returns a copy of all the loaded RV datasets.
 double[] getElementColumn(int i)
          Returns an array containing the i-th orbital element for all planets.
 double[][] getElements()
          Returns an matrix containing the orbital elements of the system (P, M, MA, e, LOP, i, node, Rp) for each planet.
 double[] getElementsA()
          Returns an array containing the orbital elements of the system (P, M, MA, e, LOP, i, node, Rp, ...) for each planet.
static Kernel getGUIState()
          If the systemic GUI is open, returns a new independent Kernel object representing the state currently loaded in the GUI (e.g.
 double getInitialEpoch()
           
 int getIntMethod()
          Returns the integration method in use
 double getJitter()
          Returns the required jitter for the current model (call after calculate()).
 double getLambda()
          Returns the weight factor λ.
 double getMstar()
          Returns stellar mass (grams)
 java.lang.String[] getNames()
           
 int getNOffsets()
          Returns the number of RV datasets currently loaded.
 java.lang.String[] getNotes()
           
 int getNPlanets()
          Returns the number of planets in the current system.
 transitsBag getObsTransits()
          Returns a transitsBag collection containing all loaded transits.
 transitsBag getObsTransitsCopy()
          Returns a transitsBag collection containing all loaded transits.
 double[][] getPowerSpectrum(boolean useWeights)
           
static double[][] getPowerSpectrum(double[][] data, boolean useWeights)
          Returns the power spectrum of the given double array.
 double[][] getPredRV()
           
 transitsBag getPredTransits()
           
 double[][] getPredTrPos()
           
 java.lang.String getPrimaryTele()
          Returns the name of the primary telescope used as global stellar offset.
 double getRMS()
          Returns the radial velocity RMS for the current model (call after calculate()).
 double[][] getRVResiduals()
           
static double[][] getSamplingSpectrum(double[][] data)
          Returns the power spectrum of the window
 StellarSystem getSys()
          Undocumented
 SystemData getSystemData()
          Returns information about the stellar host.
 double[] getVoffs()
          Returns an array containing the velocity offsets.
 double[][] getXYZ()
          Returns the current cartesian coordinates.
 double[][] getXYZ(boolean cgs)
           
 double[][] integrate(double startTime, double endTime, double dt, boolean continued)
          Integrates a system forward in time.
 double[][] integrate(double startTime, double endTime, double dt, boolean continued, java.lang.Object[] opts)
           
 double[][] integrateAdv(double startTime, double endTime, double dt, double orbSample, boolean continued, java.lang.Object[] options)
          UNDOCUMENTED
 transitsBag integrateTransits(double startTime, double endTime, double dt, boolean continued)
          Returns the predicted transits between startTime and endTime in a transitsBag object.
 boolean isCoplanar()
           
 boolean isFitWithTransits()
          Whether the chi^2 budget in calculate() and the minimization routines will include transits.
 boolean isPredictTransits()
          Whether calculate() will predict transits (getPredTransits() will return the predicted transits within the window of the observed transits).
 boolean lm(int[][] icFit, int[] voFit)
          Deprecated. Use minimize(int, boolean[][], boolean[], java.lang.Object[]) instead
 java.lang.Object[][] loadCustomIntegrators()
          Loads all the integrators specified by integrators.ini
 boolean loadData(java.lang.String fileName)
          Loads observations (.sys files and associated .vels and .tds files) from fileName.
 boolean loadFit(java.lang.String fileNm)
          Opens a fit.
 boolean loadFit(java.lang.String fileNm, java.lang.String tab)
          UNDOCUMENTED
 java.lang.Object minimize(int method, boolean[][] icFit, boolean[] voFit, java.lang.Object[] pars)
          Attempts to locate a best-fit model using the current orbital elements as the starting point.
 void printElements()
          Prints all orbital elements to the standard output.
 void printElements(java.io.Writer w)
          Prints all orbital elements to the specified Writer object.
 void printState()
          Prints a summary of the loaded data and orbital elements.
protected  void readBinary(java.io.FileInputStream fos)
           
protected  void readBinary(java.io.FileInputStream fos, java.io.ObjectInputStream xenc)
           
protected  void readXML(java.io.FileInputStream fos)
           
protected  void readXML(java.io.FileInputStream fos, java.beans.XMLDecoder xenc)
           
protected  void reduce()
           
 void restore(java.lang.String fn, int type)
          Restores a Kernel object from a previously saved snapshot (created using save())
 MCMC sa(int steps, double T, double alpha, boolean[][] icFit, boolean[] voFit, double[] bm, double[] bmv, double[] cmin, double[] cmax, java.lang.Object[] pars)
          Deprecated. Use minimize(int, boolean[][], boolean[], java.lang.Object[]) instead
 void save(java.lang.String fn, int type)
          Saves the entire state of the Kernel object into a file.
protected  void saveBinary(java.io.FileOutputStream fos)
           
protected  void saveBinary(java.io.FileOutputStream fos, java.io.ObjectOutputStream xenc)
           
protected  void saveXML(java.io.FileOutputStream fos)
           
protected  void saveXML(java.io.FileOutputStream fos, java.beans.XMLEncoder xenc)
           
 void setCartesian(double[][] xyzEls)
           
 void setCoplanar(boolean c)
           
 void setElements(double[][] elements)
          Sets the orbital elements of the system for each planet (see getElements for description).
 void setFitWithTransits(boolean fitWithTransits)
          Sets whether the chi^2 budget in calculate() and the minimization routines will include transits.
 void setInitialEpoch(double epoch)
          Sets the epoch for the orbital elements.
 void setIntMethod(int integrationMethod)
          Sets the integration method used by all the methods in this class.
 void setLambda(double lambda)
          Sets the weight factor for χ^2_TR: χ^2 = 1/NDATA (χ^2_RV + λχ^2_TR).
 void setMstar(double starMass)
          Sets stellar mass (grams).
 void setNames(java.lang.String[] planetNames)
           
 void setNotes(java.lang.String[] planetNotes)
           
 void setObsTransits(java.lang.String fileName, int flag)
          Loads a transitsBag collection of transits from a file.
 void setObsTransits(transitsBag tr)
          Loads a transitsBag collection of transits.
 void setPredictTransits(boolean predictTransits)
          Activates transit prediction when calculate() is called.
 void setPrimaryTele(java.lang.String primaryTele)
          Sets the name of the primary telescope used as global stellar offset.
 void setSystemData(SystemData data)
           
 void setTransitingPlanet(int tp)
          Sets the transiting planet.
 Bootstrap setupBootstrap(int steps)
          Deprecated. (use errorEst(int, int, boolean[][], boolean[], java.lang.Object[]) instead)
 void setVoffs(double[] vOffsets)
          Sets the velocity offsets for each RV dataset.
 void sortElements()
           
 void sortElementsByPnum()
           
 
Methods inherited from class java.lang.Object
equals, finalize, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait
 

Field Detail

tagData

public java.lang.Object tagData

OESIZE

public static int OESIZE
Number of orbital elements per planet.

Constructor Detail

Kernel

public Kernel()
Method Detail

isPredictTransits

public boolean isPredictTransits()
Whether calculate() will predict transits (getPredTransits() will return the predicted transits within the window of the observed transits).

Returns:

getGUIState

public static Kernel getGUIState()
If the systemic GUI is open, returns a new independent Kernel object representing the state currently loaded in the GUI (e.g. orbital elements, loaded datasets, velocity offsets, transits, etc.).

Returns:
A Kernel object representing the GUI state, or null if the GUI is not open.

save

public void save(java.lang.String fn,
                 int type)
          throws java.io.FileNotFoundException,
                 java.io.IOException
Saves the entire state of the Kernel object into a file.

Parameters:
fn - Name of the file
type - 0 -> binary format, 1 -> XML
Throws:
java.io.FileNotFoundException
java.io.IOException

saveBinary

protected void saveBinary(java.io.FileOutputStream fos)
                   throws java.io.IOException
Throws:
java.io.IOException

reduce

protected void reduce()

saveBinary

protected void saveBinary(java.io.FileOutputStream fos,
                          java.io.ObjectOutputStream xenc)
                   throws java.io.IOException
Throws:
java.io.IOException

saveXML

protected void saveXML(java.io.FileOutputStream fos)

saveXML

protected void saveXML(java.io.FileOutputStream fos,
                       java.beans.XMLEncoder xenc)

readXML

protected void readXML(java.io.FileInputStream fos)

readXML

protected void readXML(java.io.FileInputStream fos,
                       java.beans.XMLDecoder xenc)

readBinary

protected void readBinary(java.io.FileInputStream fos)
                   throws java.io.IOException,
                          java.lang.ClassNotFoundException
Throws:
java.io.IOException
java.lang.ClassNotFoundException

readBinary

protected void readBinary(java.io.FileInputStream fos,
                          java.io.ObjectInputStream xenc)
                   throws java.io.IOException,
                          java.lang.ClassNotFoundException
Throws:
java.io.IOException
java.lang.ClassNotFoundException

restore

public void restore(java.lang.String fn,
                    int type)
             throws java.io.FileNotFoundException,
                    java.io.IOException,
                    java.lang.ClassNotFoundException
Restores a Kernel object from a previously saved snapshot (created using save())

Parameters:
fn - File name
type - type of the snapshot (0 = binary, 1 = XML)
Throws:
java.io.FileNotFoundException
java.io.IOException
java.lang.ClassNotFoundException

setPredictTransits

public void setPredictTransits(boolean predictTransits)
Activates transit prediction when calculate() is called. Transits can be retrieved with getPredTransits().

Parameters:
predictTransits -

getPrimaryTele

public java.lang.String getPrimaryTele()
Returns the name of the primary telescope used as global stellar offset.

Returns:
Name of the primary telescope.

setPrimaryTele

public void setPrimaryTele(java.lang.String primaryTele)
Sets the name of the primary telescope used as global stellar offset.

Parameters:
primaryTele - String indicating the name of the primary telescope (must be equal to one of the strings returned by getDataNames())

ftest

public static double[] ftest(Kernel k1,
                             Kernel k2,
                             int what)
Computes the result of f-test test for two models, represented by kernel objects k1 and k2 (see Numerical Recipes). When the returned probability is close to 1.0, the two models are statistically equivalent.

Parameters:
k1 - First model
k2 - Second model
what -
  • what = 0 computes the statistic for the RV residuals;
  • what = 1 computes the statistic for the transit residuals
Returns:
an array containing {F, prob}.

getInitialEpoch

public double getInitialEpoch()
Returns:
returns the current epoch.

getIntMethod

public int getIntMethod()
Returns the integration method in use

Returns:
integration method id

setInitialEpoch

public void setInitialEpoch(double epoch)
Sets the epoch for the orbital elements. This is automatically set to the JD of the first RV observation when you load a .fit (loadFit), a .sys (loadData) file or a new RV dataset (addRVDataSet).

Parameters:
epoch - in JD.

isFitWithTransits

public boolean isFitWithTransits()
Whether the chi^2 budget in calculate() and the minimization routines will include transits.

Returns:

setFitWithTransits

public void setFitWithTransits(boolean fitWithTransits)
Sets whether the chi^2 budget in calculate() and the minimization routines will include transits.

Parameters:
fitWithTransits -

getElements

public double[][] getElements()
Returns an matrix containing the orbital elements of the system (P, M, MA, e, LOP, i, node, Rp) for each planet. Size of the matrix: npl x OESIZE.

Returns:
matrix containing the orbital elements of the system (P, M, MA, e, LOP, i, node, Rp) for each planet.

getAltElements

public double[][] getAltElements()
Returns an matrix containing extra orbital elements of the system (K, n, a, Ptr, Tperi) for each planet.

Returns:

getElementsA

public double[] getElementsA()
Returns an array containing the orbital elements of the system (P, M, MA, e, LOP, i, node, Rp, ...) for each planet. Size of the array: 1 x (npl x OESIZE).

Returns:
matrix containing the orbital elements of the system (P, M, MA, e, LOP, i, node, Rp, ...) for each planet.

getElementColumn

public double[] getElementColumn(int i)
Returns an array containing the i-th orbital element for all planets.

Parameters:
i - column number (see getElements)
Returns:
array containing the i-th orbital element for all planets.

getVoffs

public double[] getVoffs()
Returns an array containing the velocity offsets.

Returns:
array containing the velocity offsets.

setElements

public void setElements(double[][] elements)
Sets the orbital elements of the system for each planet (see getElements for description). Notice that each row must have at least OESIZE columns (e.g. the matrix must be sized npl x OESIZE).

Parameters:
elements - the matrix

setCartesian

public void setCartesian(double[][] xyzEls)

sortElements

public void sortElements()

sortElementsByPnum

public void sortElementsByPnum()

setNames

public void setNames(java.lang.String[] planetNames)

getNames

public java.lang.String[] getNames()

setNotes

public void setNotes(java.lang.String[] planetNotes)

getNotes

public java.lang.String[] getNotes()

setVoffs

public void setVoffs(double[] vOffsets)
Sets the velocity offsets for each RV dataset.

Parameters:
vOffsets -

setMstar

public void setMstar(double starMass)
Sets stellar mass (grams).

Parameters:
starMass - (grams)

getMstar

public double getMstar()
Returns stellar mass (grams)

Returns:
stellar mass

setIntMethod

public void setIntMethod(int integrationMethod)
Sets the integration method used by all the methods in this class. Possible values: 0 = Keplerian, 4 = Embedded Runge-Kutta, 5 = Hermite, 6 = Fast Hermite, 7 = Bulirsch-Stoer. Values above 10 indicate external integrators.

Parameters:
integrationMethod -

getNOffsets

public int getNOffsets()
Returns the number of RV datasets currently loaded.

Returns:
number of RV datasets

getNPlanets

public int getNPlanets()
Returns the number of planets in the current system.

Returns:
number of planets

getSystemData

public SystemData getSystemData()
Returns information about the stellar host.

Returns:
data object, containing data.Mstar, data.Rstar, data.nickNm, etc.

setSystemData

public void setSystemData(SystemData data)

loadData

public boolean loadData(java.lang.String fileName)
Loads observations (.sys files and associated .vels and .tds files) from fileName.

Parameters:
fileName - the path of the .sys file
Returns:
true on success.

getPowerSpectrum

public static double[][] getPowerSpectrum(double[][] data,
                                          boolean useWeights)
Returns the power spectrum of the given double array.

Parameters:
data - (time, measurement, error)
useWeights - include RV errors in the calculation
Returns:
matrix of rows (period, power)

getSamplingSpectrum

public static double[][] getSamplingSpectrum(double[][] data)
Returns the power spectrum of the window

Parameters:
data -
Returns:
power spectrum

getPowerSpectrum

public double[][] getPowerSpectrum(boolean useWeights)
Parameters:
useWeights -
Returns:

getSys

public StellarSystem getSys()
Undocumented

Returns:

integrate

public double[][] integrate(double startTime,
                            double endTime,
                            double dt,
                            boolean continued)
Integrates a system forward in time. If continued is false, the starting orbital elements will be those specified with setElements() (or loaded from a fit). If continued is true, it will use the orbital elements reached after the first call to integrate.

Parameters:
startTime -
endTime -
dt - (only for RK -- for other integrators stepsize is determined automatically)
continued -
Returns:
returns the orbital elements at endTime

integrate

public double[][] integrate(double startTime,
                            double endTime,
                            double dt,
                            boolean continued,
                            java.lang.Object[] opts)

integrateAdv

public double[][] integrateAdv(double startTime,
                               double endTime,
                               double dt,
                               double orbSample,
                               boolean continued,
                               java.lang.Object[] options)
UNDOCUMENTED

Parameters:
startTime -
endTime -
dt -
orbSample -
continued -
options -
Returns:

setTransitingPlanet

public void setTransitingPlanet(int tp)
Sets the transiting planet.

Parameters:
tp -

integrateTransits

public transitsBag integrateTransits(double startTime,
                                     double endTime,
                                     double dt,
                                     boolean continued)
Returns the predicted transits between startTime and endTime in a transitsBag object.

Parameters:
startTime -
endTime -
dt -
continued -
Returns:
a transitsBag (collection of transit times)

getXYZ

public double[][] getXYZ()
Returns the current cartesian coordinates.

Returns:

getXYZ

public double[][] getXYZ(boolean cgs)

loadFit

public boolean loadFit(java.lang.String fileNm)
Opens a fit.

Parameters:
fileNm -
Returns:
true on success.

loadFit

public boolean loadFit(java.lang.String fileNm,
                       java.lang.String tab)
UNDOCUMENTED

Parameters:
fileNm -
tab -
Returns:

printElements

public void printElements()
Prints all orbital elements to the standard output.


printElements

public void printElements(java.io.Writer w)
                   throws java.io.IOException
Prints all orbital elements to the specified Writer object.

Throws:
java.io.IOException

printState

public void printState()
Prints a summary of the loaded data and orbital elements.


addPlanet

public final int addPlanet(double[] plels)

addRVDataSet

public void addRVDataSet(double[][] rvSet,
                         java.lang.String name)
Adds a new RV dataset

Parameters:
rvSet - table of RVs; each row contains time, RV and sigma.
name -

addRVDataSet

public final void addRVDataSet(java.lang.String fileName)
Adds a new RV dataset from ASCII table fileName. Each row in the table represents a single RV measurement in the format
JDDate of measurement [JD]
RVRadial velocity measurement [m/s]
sigmaSigma of measurement [m/s]

Parameters:
fileName - ASCII table containing a RV dataset

clearRVDataSets

public final void clearRVDataSets()
Removes all RV datasets.


setupBootstrap

public final Bootstrap setupBootstrap(int steps)
Deprecated. (use errorEst(int, int, boolean[][], boolean[], java.lang.Object[]) instead)

Parameters:
steps -
Returns:

errorEst

public final KernelList errorEst(int method,
                                 int steps,
                                 boolean[][] icFit,
                                 boolean[] voFit,
                                 java.lang.Object[] opts)
Launches one of the error estimation algorithms.

Parameters:
method - 0 = bootstrap, 1 = MCMC
steps - number of models to be computed
icFit - boolean matrix of size npl x OESIZE; elements set to false correspond to orbital elements to keep fixed.
voFit - boolean array of size ndatasets, as above.
opts - Array of property Name/Value pairs:
BOOTSTRAP
BURN(int) number of models to use for burn-in
TRIALS(int) number of trial fits to compute for each bootstrap realization; to be used in conjunction with BURN
Returns:
a KernelList object containing an ensamble of fits. See documentation for KernelList for details

getCompiledData

public final double[][] getCompiledData()
Returns a table containing all the RV datasets compiled with their respective velocity offsets and sorted by JD.

Returns:

calculate

public final boolean calculate()
Calculates Chi^2 and residuals for the current fit, given the loaded data.

Returns:
true on success, false on error.

getCloseEncounter

public boolean getCloseEncounter()
Returns true if a close encounter happened during the last integration.

Returns:

setCoplanar

public void setCoplanar(boolean c)

isCoplanar

public boolean isCoplanar()

minimize

public java.lang.Object minimize(int method,
                                 boolean[][] icFit,
                                 boolean[] voFit,
                                 java.lang.Object[] pars)
Attempts to locate a best-fit model using the current orbital elements as the starting point. Upon completion, the new orbital elements will be set to those of the best-fitting model.

Parameters:
method - (0 = Levenberg-Marquardt, 1 = Simulated annealing, 2 = AMOEBA)
icFit - boolean matrix of size npl x OESIZE; elements set to true correspond to orbital elements to minimize on.
voFit - boolean array of size ndatasets, as above.
pars - Array of property Name/Value pairs:
SIMULATED ANNEALING
"STEPS"(integer) Number of iterations
"T"(double) initial temperature
"ALPHA"(double) slope of temperature scheduler function
"BM"(double[]) flat array of length npl x OESIZE; each element corresponds to the βμ (characteristic step size) for the orbital element.
"BMV"(double[]) flat array of length ndata, as above.
"CMIN"/"CMAX"(double[]) flat array of length npl x OESIZE; gives min/max constraint for the corresponding orbital element. A value of Double.NaN implies no constraint.
"LOG"(boolean[]) flat array of length npl x OESIZE; each element set to true implies the step will be logarithmic for the corresponding orbital element
AMOEBA
"MAXITER"(double) maximum number of iterations

Example: kernel.minimize(1, icFit, voFit, new Object[] { "STEPS", 1000, "ALPHA", 1., "BM", bm });

Returns:
reserved for future use.

sa

public MCMC sa(int steps,
               double T,
               double alpha,
               boolean[][] icFit,
               boolean[] voFit,
               double[] bm,
               double[] bmv,
               double[] cmin,
               double[] cmax,
               java.lang.Object[] pars)
Deprecated. Use minimize(int, boolean[][], boolean[], java.lang.Object[]) instead

Parameters:
steps -
T -
alpha -
icFit -
voFit -
bm -
bmv -
cmin -
cmax -
pars -
Returns:

lm

public boolean lm(int[][] icFit,
                  int[] voFit)
Deprecated. Use minimize(int, boolean[][], boolean[], java.lang.Object[]) instead

Parameters:
icFit -
voFit -

getObsTransits

public transitsBag getObsTransits()
Returns a transitsBag collection containing all loaded transits.

Returns:
Returns a transitsBag collection containing all loaded transits.

getObsTransitsCopy

public transitsBag getObsTransitsCopy()
Returns a transitsBag collection containing all loaded transits.

Returns:
Returns a transitsBag collection containing all loaded transits.

getDataSets

public java.util.ArrayList<double[][]> getDataSets()
Returns a list of all the loaded RV datasets.

Returns:
List of RV datasets (each row in JD, RV, σ format).

getDataNames

public java.lang.String[] getDataNames()
Returns a list of all the loaded RV datasets names.

Returns:
Array of RV names.

getDataSetsCopy

public java.util.ArrayList<double[][]> getDataSetsCopy()
Returns a copy of all the loaded RV datasets.

Returns:
List of RV datasets (each row in JD, RV, σ format).

setObsTransits

public void setObsTransits(transitsBag tr)
Loads a transitsBag collection of transits.

Parameters:
tr - transitsBag collection.

setObsTransits

public void setObsTransits(java.lang.String fileName,
                           int flag)
Loads a transitsBag collection of transits from a file.

Parameters:
fileName - table of central transits (in the format JD sigma flag, where flag = 1 for primary transit or 2 for secondary eclipse).
flag - reserved for future use

getPredRV

public double[][] getPredRV()
Returns:
Returns an array of predicted RVs at each observed RV time (call after calculate()).

getPredTrPos

public double[][] getPredTrPos()
Returns:
Returns a table of projected separations at the observed central times (call after calculate()).

getPredTransits

public transitsBag getPredTransits()
Returns:
Returns a transitsBag collection of predicted transits within the observational window (call after calculate()).

getRVResiduals

public double[][] getRVResiduals()
Returns:
Returns an array of predicted RV residuals at each observed RV time (call after calculate()).

getChiSqr

public final double getChiSqr()
Returns the combined (RV + transits) χ^2 (call after calculate()).

Returns:
Returns total χ^2

getRMS

public final double getRMS()
Returns the radial velocity RMS for the current model (call after calculate()).

Returns:
RV rms [m/s]

getJitter

public final double getJitter()
Returns the required jitter for the current model (call after calculate()).

Returns:
RV jitter [m/s]

getChiSqrAll

public final double[] getChiSqrAll()
Returns an array of χ^2 statistics.

Returns:
Returns an array containing: combined chi^2, unnormalized combined chi^2, chi^2_RV, chi^2_transits.

setLambda

public void setLambda(double lambda)
Sets the weight factor for χ^2_TR: χ^2 = 1/NDATA (χ^2_RV + λχ^2_TR).

Parameters:
lambda -

getLambda

public double getLambda()
Returns the weight factor λ.

Returns:
Returns lambda (see setLambda).

clone

public Kernel clone()
Creates a new, independent copy of the object. The entire state of the kernel (e.g. orbital elements, voffs, datasets, integration method, etc.) is copied into the new object.

Overrides:
clone in class java.lang.Object
Returns:
A copy of the kernel state.

loadCustomIntegrators

public java.lang.Object[][] loadCustomIntegrators()
Loads all the integrators specified by integrators.ini

Returns:
list of integrators (see docs)