Merge branch 'FisherProp' into 'master'

Fisher proposals

See merge request lscsoft/bayeswave!194

src/BayesWaveProposal.c

@@ -490,7 +490,8 @@ void fisher_matrix_proposal(struct FisherMatrix *fisher, double *params, gsl_rng
   for (j=0; j<N; j++) jump[j] += Amps[i]*fisher->evector[j][i]*sqD;
 
   //jump from current position
-  for(i=0; i<N; i++) y[i] = params[i] + gaussian_draw(seed)*jump[i];
+  double c = gaussian_draw(seed);
+  for(i=0; i<N; i++) y[i] = params[i] + c*jump[i];
 }
 
 void intrinsic_fisher_update(double *params, double *dparams, double *Snf, double Tobs, int NW, int TFQFLAG)
@@ -656,9 +657,10 @@ 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]*gaussian_draw(seed);
+    paramsy[k] = paramsx[k] + scale*dparamsx[k]*c;
   }
 
   if(checkrange(paramsy, range, NW))