Skip to content
Snippets Groups Projects

Fix PSD prior error

Merged Katerina Chatziioannou requested to merge katerina.chatziioannou/bayeswave:master into master
1 file
+ 5
2
Compare changes
  • Side-by-side
  • Inline
+ 5
2
@@ -2752,6 +2752,7 @@ void BayesLineBurnin(struct BayesLineParams *bayesline, double *timeData, double
double *Sl = malloc(sizeof(double)*(N+1));
int imin = data->fmin*Tobs;
for(i=0; i<N/2; i++)
{
@@ -2759,10 +2760,12 @@ void BayesLineBurnin(struct BayesLineParams *bayesline, double *timeData, double
logSn[i] = 1.0;
Sl[i] = 0.0;
if(i<N/2-imin){
bayesline->Sbase[i] = 1.0;
bayesline->Sline[i] = 0.0;
}
}
//interpolate cubic spline points
CubicSplineGSL(spline->n,spline->points,spline->data,N/2-imin,f+imin,logSn+imin);
@@ -2801,7 +2804,7 @@ void BayesLineBurnin(struct BayesLineParams *bayesline, double *timeData, double
bayesline->priors->lower[i] = bayesline->Sbase[i]/100.;
bayesline->priors->upper[i] = bayesline->Snf[i]*100.;
}
free(f);
free(logSn);
free(Sl);
Loading