Skip to content

Commit

Permalink
rm unused acceptance counter
Browse files Browse the repository at this point in the history
  • Loading branch information
fabian-s committed Sep 3, 2024
1 parent 2474d6f commit 8e6901c
Showing 1 changed file with 9 additions and 9 deletions.
18 changes: 9 additions & 9 deletions src/sampler.c
View file
Original file line number Diff line number Diff line change
Expand Up @@ -96,7 +96,7 @@ SEXP sampler(
double *p1 =R_Calloc(*pPen, double);

double infV = 100000, oneV = 1.0, zeroV = 0.0, minusOneV =-1.0;
double *one=&oneV, *zero=&zeroV, *minusOne=&minusOneV, *inf=&infV, acceptance=0;
double *one=&oneV, *zero=&zeroV, *minusOne=&minusOneV, *inf=&infV; // , acceptance=0;
double invSigma2 = 1 / *sigma2, sqrtInvSigma2 = R_pow(invSigma2, 0.5);
double *penAlphaSq, *alphaLong, *varAlpha, *priorMeanAlpha, *modeAlpha, *offsetAlpha;;
penAlphaSq = R_Calloc(*pPen, double);
Expand Down Expand Up @@ -396,14 +396,14 @@ SEXP sampler(
PutRNGstate();

if(*verbose) Rprintf(".");
if(*family > 0) {
acceptance = 0.0;
for(j=0; j<*blocksAlpha; j++) acceptance += acceptAlpha[j];
acceptance = 0.0;
if(qKsiNoUpdate < *q){
for(j=0; j<*blocksKsi; j++) acceptance += acceptKsi[j];
}
}
// if(*family > 0) {
// acceptance = 0.0;
// for(j=0; j<*blocksAlpha; j++) acceptance += acceptAlpha[j];
// acceptance = 0.0;
// if(qKsiNoUpdate < *q){
// for(j=0; j<*blocksKsi; j++) acceptance += acceptKsi[j];
// }
// }

R_Free(etaOffset); R_Free(XKsiUpdate); R_Free(XAlpha); R_Free(resid); R_Free(eta);
R_Free(offsetKsi); R_Free(modeKsi); R_Free(priorMeanKsi); R_Free(ksiUpdate);
Expand Down

0 comments on commit 8e6901c

Please sign in to comment.