150 REAL8 logZnoise = 0.;
151 struct timeval time1, time2;
152 gettimeofday( &time1, NULL );
169 if ( !testgausslike ) {
173 REAL8 h0val = 0., h0sigma = 2.5e-24, h0mean = 0.;
174 runState.
data = NULL;
177 runState.
data->likeli_counter = 0;
178 runState.
data->templa_counter = 0;
179 runState.
data->next = NULL;
205 if ( testgausslike ) {
215 if ( !testgausslike ) {
220 if ( !testgausslike ) {
225 if ( !testgausslike ) {
233 if ( !testgausslike ) {
238 if ( !testgausslike ) {
242 if ( !testgausslike ) {
247 if ( !testgausslike ) {
278 if ( !testgausslike ) {
283 if ( testgausslike ) {
289 gettimeofday( &time2, NULL );
290 REAL8 tottime = (
REAL8 )( ( time2.tv_sec + time2.tv_usec * 1.e-6 ) - ( time1.tv_sec + time1.tv_usec * 1.e-6 ) );
#define __func__
log an I/O error, i.e.
void LALInferenceAddVariable(LALInferenceVariables *vars, const char *name, const void *value, LALInferenceVariableType type, LALInferenceParamVaryType vary)
void * LALInferenceGetVariable(const LALInferenceVariables *vars, const char *name)
ProcessParamsTable * LALInferenceGetProcParamVal(ProcessParamsTable *procparams, const char *name)
ProcessParamsTable * LALInferenceParseCommandLine(int argc, char **argv)
int LALInferenceCheckVariable(LALInferenceVariables *vars, const char *name)
INT4 LALInferenceNestedSamplingOneStep(LALInferenceRunState *runState)
void * XLALMalloc(size_t n)
void * XLALCalloc(size_t m, size_t n)
XLALErrorHandlerType * XLALSetErrorHandler(XLALErrorHandlerType *newHandler)
void XLALExitErrorHandler(const char *func, const char *file, int line, int errnum)
void initialise_proposal(LALInferenceRunState *runState)
Initialise the MCMC proposal distribution for sampling new points.
void initialise_prior(LALInferenceRunState *runState)
Sets up the parameters to be searched over and their prior ranges.
void sum_data(LALInferenceRunState *runState)
Calculates the sum of the square of the data and model terms.
void initialise_threads(LALInferenceRunState *state, INT4 nthreads)
void initialise_algorithm(LALInferenceRunState *runState)
Initialises the nested sampling algorithm control.
void setup_live_points_array_wrapper(LALInferenceRunState *runState)
A wrapper around LALInferenceSetupLivePointsArray.
void nested_sampling_algorithm_wrapper(LALInferenceRunState *runState)
A wrapper around LALInferenceNestedSamplingAlgorithm.
Header file for the initialisation functions for the parameter estimation code for known pulsar searc...
void inject_signal(LALInferenceRunState *runState)
Inject a simulated signal into the data.
void get_loudest_snr(LALInferenceRunState *runState)
Get the signal-to-noise ratio of the maximum likelihood signal.
Header file for the signal injection functions for the parameter estimation code for known pulsar sea...
REAL8 priorFunction(LALInferenceRunState *runState, LALInferenceVariables *params, UNUSED LALInferenceModel *model)
The prior function.
REAL8 noise_only_likelihood(LALInferenceRunState *runState)
Calculate the natural logarithm of the evidence that the data consists of only Gaussian noise The fun...
REAL8 pulsar_log_likelihood(LALInferenceVariables *vars, LALInferenceIFOData *data, LALInferenceModel *get_model)
The log likelihood function.
Header file for the likelihood and prior functions used in parameter estimation code for known pulsar...
void get_pulsar_model(LALInferenceModel *model)
Defines the pulsar model/template to use.
Header file for the signal models functions used in parameter estimation code for known pulsar search...
void setup_from_par_file(LALInferenceRunState *runState)
Reads in the parameters of the pulsar being searched for.
void read_pulsar_data(LALInferenceRunState *runState)
Reads in the input gravitational wave data files, or creates fake data sets.
Header file for the data reading functions for the parameter estimation code for known pulsar searche...
void generate_interpolant(LALInferenceRunState *runState)
Generate an orthonormal basis set of waveforms from a training set.
Header file for the reduced order quadrature generation used in parameter estimation code for known p...
void test_gaussian_output(LALInferenceRunState *runState)
Output the analytic evidence for the test Gaussian likelihood and a 95% upper limit.
REAL8 test_gaussian_log_likelihood(LALInferenceVariables *vars, LALInferenceIFOData *data, LALInferenceModel *get_model)
Test the sampler using a Gaussian likelihood.
void compare_likelihoods(LALInferenceRunState *rs)
Read in an ascii text file of nested samples and compare the log likelihoods.
void outputPriorSamples(LALInferenceRunState *runState)
Output a number of prior samples based on the initial live points.
void gridOutput(LALInferenceRunState *runState)
A test function to calculate a 1D posterior on a grid.
Header file for functions used in testing the parameter estimation code for known pulsar searches usi...
Header file for the helper functions for the parameter estimation code for known pulsar searches usin...
LALStringVector * corlist
INT4 main(INT4 argc, CHAR *argv[])
Header file for the parameter estimation code for known pulsar searches using the nested sampling alg...
struct tagLALInferenceIFOModel * next
LALInferenceVariables * params
REAL8 * ifo_loglikelihoods
LALInferenceIFOModel * ifo
LALInferenceTemplateFunction templt
ProcessParamsTable * commandLine
LALInferenceVariables * algorithmParams
struct tagLALInferenceIFOData * data
LALInferenceThreadState * threads
LALInferenceAlgorithm algorithm
LALInferenceEvolveOneStepFunction evolve
LALInferencePriorFunction prior
LALInferenceLikelihoodFunction likelihood
LALInferenceVariables * currentParams
LALInferenceModel * model