Commit 576178fa authored by Marcella Wijngaarden's avatar Marcella Wijngaarden

Start of output reduction (less waveform prints)

parent 9b58d5c5
Pipeline #160547 failed with stages
in 1 minute and 9 seconds
......@@ -519,7 +519,7 @@ void print_projected_cbc_waveform(double **SN, double Tobs, double ttrig, double
fac = ceil(sqrt(1.0/4.0));
for (id = 0; id < net->Nifo; ++id)
{
sprintf(command, "wavesCBC/projected_wavewhite_%d_%d_%d_%d.dat", iter, (int)(Tobs), (int)ttrig, net->labels[id]);
sprintf(command, "waveforms/projected_wavewhite_%d_%d_%d_%d.dat", iter, (int)(Tobs), (int)ttrig, net->labels[id]);
out = fopen(command,"w");
for (i = 0; i < N; ++i)
{
......@@ -527,7 +527,7 @@ void print_projected_cbc_waveform(double **SN, double Tobs, double ttrig, double
}
fclose(out);
sprintf(command, "wavesCBC/projected_wavedatawhite_%d_%d_%d_%d.dat", iter, (int)(Tobs), (int)ttrig, net->labels[id]);
sprintf(command, "waveforms/projected_wavedatawhite_%d_%d_%d_%d.dat", iter, (int)(Tobs), (int)ttrig, net->labels[id]);
out = fopen(command,"w");
for (i = 0; i < N; ++i)
{
......@@ -3796,11 +3796,11 @@ void BayesCBC_MCMC(int M, int N, double Tobs, double ttrig, gsl_rng *r, int ista
}
}
if(ic==0 && rundata->verbose==1)
{
// Save waveforms to file
printwaveall(rundata->net, N, rundata->freq, rundata->pallx[icbc], rundata->SN, Tobs, ttrig, rundata->mxc[0], rundata);
}
// if(ic==0 && rundata->verbose==1)
// {
// // Save waveforms to file
// printwaveall(rundata->net, N, rundata->freq, rundata->pallx[icbc], rundata->SN, Tobs, ttrig, rundata->mxc[0], rundata);
// }
// Save MAP likelihood. TODO: Do we want to have this?
// for(i = 0; i < NP; i++) rundata->pallx[ic][i] = pmax[i];
......@@ -6190,6 +6190,7 @@ void skymcmc(struct Net *net, int MCX, int *mxc, FILE *chain, double **paramx, d
id1 = net->labels[id1];
id2 = net->labels[id2];
// get new sky locations with the same time delays
Ring(skyx[q], skyy[q], id1, id2, r, rundata->gmst);
qyx = 1.0;
qxy = 1.0;
......
......@@ -2630,8 +2630,8 @@ void EvolveBayesCBCParameters(struct Data *data, struct Model **model, struct ba
// Get CBC waveforms projected on the detectors
projectCBCWaveform(model_x->cbcamphase, N, NI, fmin, data->Tobs, model_x->extParams, model_x->cbctemplate, projection->dtimes, projection->Fplus, projection->Fcross);
// Save waveform temporarily for checking residuals + template
if (ic==0 && bayescbc->debug == 1) print_projected_cbc_waveform(bayescbc->SN, data->Tobs, data->trigtime, model_x->cbctemplate, bayescbc->D, N, bayescbc->mxc[0], bayescbc);
// Save 200 waveforms temporarily for checking residuals + template
if (ic==0 && bayescbc->debug == 1 && chain->mc%(chain->count/200)==0) print_projected_cbc_waveform(bayescbc->SN, data->Tobs, data->trigtime, model_x->cbctemplate, bayescbc->D, N, bayescbc->mxc[0], bayescbc);
// Recompute likelihoods of current chain
model_x->logLnorm = 0.0;
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment