Skip to content
Snippets Groups Projects

fix fisher jump

1 file
+ 1
2
Compare changes
  • Side-by-side
  • Inline
+ 1
2
@@ -657,10 +657,9 @@ void intrinsic_fisher_proposal(UNUSED struct Model *model, double **range, doubl
// [0] t0 [1] f0 [2] Q [3] Amp [4] phi ([5] beta chirplet)
//perturb current location by N[0,scale*sigma]
double c = gaussian_draw(seed);
for(k=0; k<NW; k++)
{
paramsy[k] = paramsx[k] + scale*dparamsx[k]*c;
paramsy[k] = paramsx[k] + scale*dparamsx[k]*gaussian_draw(seed);
}
if(checkrange(paramsy, range, NW))
Loading