use existing functions for wavelet SNR

parent 2cbf058b
......@@ -174,7 +174,9 @@ void draw_signal_amplitude(double *params, double *Snf, gsl_rng *seed, double To
}
//SNR defined with Sn(f) but Snf array holdes <n_i^2>
params[3] = SNR/sqrt((params[2]/(2.0*RT2PI*params[1]))/(Snf[i]*2.0/Tobs));
//params[3] = SNR/sqrt((params[2]/(2.0*RT2PI*params[1]))/(Snf[i]*2.0/Tobs));
params[3] *= SNR/SineGaussianSNR(params, Snf, Tobs);
// FILE *temp=fopen("prior.dat","a");
// fprintf(temp,"%lg\n",params[3]);
......
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