diff options
| -rw-r--r-- | src/lossdistrib.c | 5 |
1 files changed, 1 insertions, 4 deletions
diff --git a/src/lossdistrib.c b/src/lossdistrib.c index d8257bd..bc39451 100644 --- a/src/lossdistrib.c +++ b/src/lossdistrib.c @@ -462,7 +462,6 @@ void lossdistrib_joint_Z(const double *dp, const double *pp, const int *ndp, #pragma omp parallel for for(int i = 0; i < *nZ; i++){ double* dpshocked = malloc(sizeof(double) * (*ndp)); - double* Sshocked = malloc(sizeof(double) * (*ndp)); double* ppshocked = NULL; if(pp) { ppshocked = malloc(sizeof(double) * (*ndp)); @@ -472,13 +471,11 @@ void lossdistrib_joint_Z(const double *dp, const double *pp, const int *ndp, if(pp) { ppshocked[j] = shockprob(pp[j], rho[j], -Z[i], 0); } - Sshocked[j] = shockseverity(S[j], Z[i], rho[j], dp[j]); } - lossdistrib_joint(dpshocked, ppshocked, ndp, w, Sshocked, + lossdistrib_joint(dpshocked, ppshocked, ndp, w, S + (*ndp) * i, N, defaultflag, qmat + N2 * i); free(dpshocked); - free(Sshocked); if(pp) { free(ppshocked); } |