Skip to content
Snippets Groups Projects

Changes to BayesWaveCleanFrame (and bayeswave_pipe) to work with scitokens

Merged Changes to BayesWaveCleanFrame (and bayeswave_pipe) to work with scitokens
sophie.hourihane/bayeswave:scitokens-for-igwn-pool into master
Merged Sophie Hourihane requested to merge sophie.hourihane/bayeswave:scitokens-for-igwn-pool into master
Viewing commit 9fbd2db0
Prev Next
Show latest version
1 file
+ 41
− 42
Compare changes
  • Side-by-side
  • Inline
src/BayesWaveCleanFrame.c
+ 41
− 42
@@ -108,30 +108,29 @@ void printProgress (double percentage)
/* ============================ MAIN PROGRAM ============================ */
int main(int argc, char *argv[])
{
fprintf(stdout,"\n");
fprintf(stdout,"#############################################################\n");
fprintf(stdout," BayesWave Glitch Subtraction \n");
fprintf(stdout," * \n");
fprintf(stdout," * \n");
fprintf(stdout," * \n");
fprintf(stdout," * \n");
fprintf(stdout," ** \n");
fprintf(stdout," * * \n");
fprintf(stdout,"**************************** * ***************************\n");
fprintf(stdout," * ** \n");
fprintf(stdout," ** \n");
fprintf(stdout," * \n");
fprintf(stdout," * \n");
fprintf(stdout," * \n");
fprintf(stdout," * \n");
fprintf(stdout," \n");
fprintf(stdout,"#############################################################\n");
fprintf(stdout,"\n");
int main(int argc, char *argv[]) {
fprintf(stdout, "\n");
fprintf(stdout, "#############################################################\n");
fprintf(stdout, " BayesWave Glitch Subtraction \n");
fprintf(stdout, " * \n");
fprintf(stdout, " * \n");
fprintf(stdout, " * \n");
fprintf(stdout, " * \n");
fprintf(stdout, " ** \n");
fprintf(stdout, " * * \n");
fprintf(stdout, "**************************** * ***************************\n");
fprintf(stdout, " * ** \n");
fprintf(stdout, " ** \n");
fprintf(stdout, " * \n");
fprintf(stdout, " * \n");
fprintf(stdout, " * \n");
fprintf(stdout, " * \n");
fprintf(stdout, " \n");
fprintf(stdout, "#############################################################\n");
fprintf(stdout, "\n");
print_version(stdout);
fprintf(stdout,"\n");
fprintf(stdout,"\n");
fprintf(stdout, "\n");
fprintf(stdout, "\n");
/* Variable declaration */
int i, j, k, NGWF;
@@ -186,44 +185,44 @@ int main(int argc, char *argv[])
outdir = data->outdir;
char cleandir[MAXSTRINGSIZE];
sprintf(cleandir,"%s/clean_frame", data->outdir);
mkdir(cleandir,S_IRWXU | S_IRWXG | S_IROTH | S_IXOTH);
sprintf(cleandir, "%s/clean_frame", data->outdir);
mkdir(cleandir, S_IRWXU | S_IRWXG | S_IROTH | S_IXOTH);
char outframeType[MAXSTRINGSIZE];
char outframeChannel[MAXSTRINGSIZE];
char outframeGlitchChannel[MAXSTRINGSIZE];
snprintf(outframeType, MAXSTRINGSIZE, "%s_%s", data->fr_type, version);
snprintf(outframeType, MAXSTRINGSIZE, "%s_%s", data->fr_type, version);
snprintf(outframeChannel, MAXSTRINGSIZE, "%s_%s", data->fr_chanl, version);
snprintf(outframeGlitchChannel, MAXSTRINGSIZE, "%s_glitch", data->fr_chanl);
/* set the channel name */
strncpy(timeData->name, data->fr_chanl, LALNameLength);
strncpy(timeRes->name, outframeChannel, LALNameLength);
strncpy(timeData->name, data->fr_chanl, LALNameLength);
strncpy(timeRes->name, outframeChannel, LALNameLength);
strncpy(timeGlitch->name, outframeGlitchChannel, LALNameLength);
fprintf(stdout,"Creating glitch-subtracted frame:\n");
fprintf(stdout," Type: %s\n",outframeType);
fprintf(stdout," Channel: %s\n",outframeChannel);
fprintf(stdout,"\n");
fprintf(stdout, "Creating glitch-subtracted frame:\n");
fprintf(stdout, " Type: %s\n", outframeType);
fprintf(stdout, " Channel: %s\n", outframeChannel);
fprintf(stdout, "\n");
// shortcuts to original frame size/cadence
dT = timeData->deltaT;
NGWF = timeData->data->length;
dT = timeData->deltaT;
NGWF = timeData->data->length;
// Length of the up-sampled glitch data
NBW = (int)(data->bw_seglen * frame_sample_rate);
NBW = (int) (data->bw_seglen * frame_sample_rate);
double *glitch_model = malloc(sizeof(double)*NBW);
double *glitch_model_padded = malloc(sizeof(double)*NGWF);
double **glitch_model_posterior = malloc(sizeof(double *)*NBW);
// stores time series of glitch model
double *glitch_model = malloc(sizeof(double) * NBW);
double *glitch_model_padded = malloc(sizeof(double) * NGWF);
double **glitch_model_posterior = malloc(sizeof(double *) * NBW);
FILE *paramsfile = fopen(data->bw_model,"r");
FILE *paramsfile = fopen(data->bw_model, "r");
int glitch_model_size;
double *glitch_model_params = malloc(NP*sizeof(double));
double *glitch_model_params = malloc(NP * sizeof(double));
//count lines in file
int N_chain_samples = 0;
while(!feof(paramsfile))