From e4ba0bb6aa38f352cf33be92c66f9c4dd4dec576 Mon Sep 17 00:00:00 2001
From: Guillaume Horel <guillaume.horel@gmail.com>
Date: Fri, 5 Apr 2019 13:53:29 -0400
Subject: abort early

---
 src/lossdistrib.c | 2 ++
 1 file changed, 2 insertions(+)

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) -
-- 
cgit v1.2.3-70-g09d2