From 436db5b3aa46db1c6b245ea57ccf8acc97e03356 Mon Sep 17 00:00:00 2001 From: Katerina Chatziioannou Date: Tue, 26 May 2020 14:31:56 -0400 Subject: [PATCH] waves->wavesCBC and make it at runtime --- src/BayesCBC.c | 7 ++++--- src/BayesWave.c | 1 + 2 files changed, 5 insertions(+), 3 deletions(-) diff --git a/src/BayesCBC.c b/src/BayesCBC.c index ffccec0..102a672 100644 --- a/src/BayesCBC.c +++ b/src/BayesCBC.c @@ -849,6 +849,7 @@ void printwaveall(struct Net *net, int N, RealVector *freq, double *paramx, doub fs = fbegin(paramx); if(fs < fmin) fs = fmin; + printf("In printwaveall: "); printf("%d %f %f\n", iter, exp(paramx[0])/MSUN_SI, fs); @@ -916,7 +917,7 @@ void printwaveall(struct Net *net, int N, RealVector *freq, double *paramx, doub for (id = 0; id < net->Nifo; ++id) { - sprintf(command, "waves/wave_%d_%d_%d_%d.dat", iter, (int)(Tobs), (int)ttrig, net->labels[id]); + sprintf(command, "wavesCBC/wave_%d_%d_%d_%d.dat", iter, (int)(Tobs), (int)ttrig, net->labels[id]); out = fopen(command,"w"); for (i = 0; i < N; ++i) { @@ -927,7 +928,7 @@ void printwaveall(struct Net *net, int N, RealVector *freq, double *paramx, doub for (id = 0; id < net->Nifo; ++id) { - sprintf(command, "waves/wavef_%d_%d_%d_%d.dat", iter, (int)(Tobs), (int)ttrig, net->labels[id]); + sprintf(command, "wavesCBC/wavef_%d_%d_%d_%d.dat", iter, (int)(Tobs), (int)ttrig, net->labels[id]); out = fopen(command,"w"); for (i = 1; i < N/2; ++i) { @@ -949,7 +950,7 @@ void printwaveall(struct Net *net, int N, RealVector *freq, double *paramx, doub printf(" amplitude factor: %f\n\n",fac); for (id = 0; id < net->Nifo; ++id) { - sprintf(command, "waves/wavewhite_%d_%d_%d_%d.dat", iter, (int)(Tobs), (int)ttrig, net->labels[id]); + sprintf(command, "wavesCBC/wavewhite_%d_%d_%d_%d.dat", iter, (int)(Tobs), (int)ttrig, net->labels[id]); out = fopen(command,"w"); for (i = 0; i < N; ++i) { diff --git a/src/BayesWave.c b/src/BayesWave.c index 7f63947..63b5219 100644 --- a/src/BayesWave.c +++ b/src/BayesWave.c @@ -148,6 +148,7 @@ int main(int argc, char *argv[]) mode_t process_mask = umask(0); mkdir("chains", S_IRWXU | S_IRWXG | S_IROTH | S_IXOTH); mkdir("waveforms", S_IRWXU | S_IRWXG | S_IROTH | S_IXOTH); + mkdir("wavesCBC", S_IRWXU | S_IRWXG | S_IROTH | S_IXOTH); if(data->checkpointFlag) mkdir("checkpoint",S_IRWXU | S_IRWXG | S_IROTH | S_IXOTH); sprintf(filename,"%ssnr",data->runName); mkdir(filename, S_IRWXU | S_IRWXG | S_IROTH | S_IXOTH); -- GitLab