From 664496f99c15f94975ab90627ded8b7648c663ee Mon Sep 17 00:00:00 2001 From: Guillaume Horel Date: Fri, 12 Oct 2018 13:02:57 -0400 Subject: shortcuts --- src/lossdistrib.c | 12 ++++++++---- 1 file changed, 8 insertions(+), 4 deletions(-) (limited to 'src') diff --git a/src/lossdistrib.c b/src/lossdistrib.c index 47867b9..e472a44 100644 --- a/src/lossdistrib.c +++ b/src/lossdistrib.c @@ -393,9 +393,11 @@ void recovdist(const double *dp, const double *pp, const int *n, const double *w } double shockprob(double p, double rho, double Z, int give_log){ - if(rho==1){ + if(rho == 1.) { return((double)(Z<=qnorm(p, 0, 1, 1, 0))); - }else{ + } else if (rho == 0.) { + return p; + } else { return( pnorm( (qnorm(p, 0, 1, 1, 0) - sqrt(rho) * Z)/sqrt(1 - rho), 0, 1, 1, give_log)); } } @@ -419,9 +421,11 @@ void shockprobvec2(const double p, const double rho, const double* Z, const int } double shockseverity(double S, double Z, double rho, double p){ - if(p==0){ + if(p == 0.) { return 0; - }else{ + } else if (rho == 0.) { + return S; + } else { return( exp(shockprob(S * p, rho, Z, 1) - shockprob(p, rho, Z, 1)) ); } } -- cgit v1.2.3-70-g09d2