diff options
Diffstat (limited to 'src/lossdistrib.c')
| -rw-r--r-- | src/lossdistrib.c | 2 |
1 files changed, 2 insertions, 0 deletions
diff --git a/src/lossdistrib.c b/src/lossdistrib.c index ce4b662..a9dc82c 100644 --- a/src/lossdistrib.c +++ b/src/lossdistrib.c @@ -3,6 +3,7 @@ #include <string.h> #include <omp.h> #include "lossdistrib.h" +#include <assert.h> #ifdef USE_BLAS #include <cblas.h> @@ -493,6 +494,7 @@ void stochasticrecov(double* R, double* Rtilde, double* Z, double* w, int* nZ, }else{ ptemp = (1 - *R) / (1 - *Rtilde) * *porig; fitprob(Z, w, nZ, rho, &ptemp, &ptilde); + assert(ptilde > *pmod); #pragma omp parallel for for(int i = 0; i < *nZ; i++){ q[i] = fabs(1 - (1 - *Rtilde) * exp( shockprob(ptilde, *rho, Z[i], 1) - |