From 0ad23ce6ee6118c3c795f88918169bd2faf571cc Mon Sep 17 00:00:00 2001 From: Thierry Gosselin Date: Wed, 9 Sep 2026 14:26:45 +1000 Subject: [PATCH] fix SAF invalid-likelihood exit status --- abcSaf.cpp | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/abcSaf.cpp b/abcSaf.cpp index 5c3ba3f..d6432c8 100644 --- a/abcSaf.cpp +++ b/abcSaf.cpp @@ -713,7 +713,7 @@ void abcSaf::algoJointPost(double **post, if(liks[i] < 0.0) { fprintf(stderr, "Negative genotype probability in %s will exit (site=%d,ind=%d,lik=%f)\n", __FUNCTION__, s, i, liks[i]); - exit(0); + exit(EXIT_FAILURE); } } memcpy(hj, liks, 3*sizeof(double)); @@ -1775,8 +1775,8 @@ void abcSaf::algoGeno(int refId,double **liks,char *major,char *minor,int nsites shouldBeOne = exp(res[0]-log(denominator))+exp(res[1]-log(denominator))+exp(res[2]-log(denominator)); //fprintf(stderr,"this should be one: %f\n",shouldBeOne); if((fabs(shouldBeOne-1)>0.000001)){ - fprintf(stderr,"this should be one: %f at position:(%s,%d)\n",shouldBeOne,header->target_name[refId],posi[it+1]); - exit(0); + fprintf(stderr,"this should be one: %f at position:(%s,%d)\n",shouldBeOne,header->target_name[refId],posi[it+1]); + exit(EXIT_FAILURE); }