diff options
| -rw-r--r-- | lossdistrib.c | 38 |
1 files changed, 19 insertions, 19 deletions
diff --git a/lossdistrib.c b/lossdistrib.c index 3f8a1dbe..fef0016d 100644 --- a/lossdistrib.c +++ b/lossdistrib.c @@ -17,7 +17,7 @@ void lossdistrib(double *p, int *np, double *w, double *S, int *N, double *q) { double *qtemp;
lu = 1./(*N-1);
- qtemp = calloc( *N, sizeof(double));
+ qtemp = Calloc(*N, double);
q[0] = 1;
M = 1;
for(i=0; i<(*np); i++){
@@ -45,7 +45,7 @@ void lossdistrib(double *p, int *np, double *w, double *S, int *N, double *q) { }
q[*N-1] += 1-sum;
}
- free(qtemp);
+ Free(qtemp);
}
void lossdistrib_truncated(double *p, int *np, double *w, double *S, int *N, int *T, double *q) {
@@ -63,8 +63,8 @@ void lossdistrib_truncated(double *p, int *np, double *w, double *S, int *N, int double *q1, *q2;
lu = 1./(*N-1);
- q1 = calloc( *T, sizeof(double));
- q2 = calloc( *T, sizeof(double));
+ q1 = Calloc( *T, double);
+ q2 = Calloc( *T, double);
q[0] = 1;
M = 1;
for(i=0; i<(*np); i++){
@@ -84,10 +84,10 @@ void lossdistrib_truncated(double *p, int *np, double *w, double *S, int *N, int for(j=0; j < MIN(M, *T-d2); j++){
q[d2+j] += q2[j];
};
- M+=d2;
+ M += d2;
}
- free(q1);
- free(q2);
+ Free(q1);
+ Free(q2);
}
void lossdistrib_joint( double *p, int *np, double *w, double *S, int *N, double *q) {
@@ -108,7 +108,7 @@ void lossdistrib_joint( double *p, int *np, double *w, double *S, int *N, double double *qtemp;
lu = 1./(*N-1);
- qtemp = calloc( (*N) * (*N), sizeof(double));
+ qtemp = Calloc( (*N) * (*N), double);
q[0] = 1;
Mx=1;
My=1;
@@ -155,7 +155,7 @@ void lossdistrib_joint( double *p, int *np, double *w, double *S, int *N, double }
q[MIN(*N, Mx)*MIN(My,*N)-1] += 1 - sum;
}
- free(qtemp);
+ Free(qtemp);
}
void recovdist(double *dp, double *pp, int *n, double *w, double *S, int *N, double *q) {
@@ -175,7 +175,7 @@ void recovdist(double *dp, double *pp, int *n, double *w, double *S, int *N, dou double *qtemp;
lu = 1./(*N - 1);
- qtemp = calloc( (*N), sizeof(double));
+ qtemp = Calloc( (*N), double);
q[0] = 1;
M=1;
for(i=0; i<(*n); i++){
@@ -209,7 +209,7 @@ void recovdist(double *dp, double *pp, int *n, double *w, double *S, int *N, dou }
q[*N-1] += 1-sum;
}
- free(qtemp);
+ Free(qtemp);
}
void lossdistrib_prepay_joint(double *dp, double *pp, int *ndp, double *w, double *S, int *N, double *q) {
@@ -221,7 +221,7 @@ void lossdistrib_prepay_joint(double *dp, double *pp, int *ndp, double *w, doubl int Mx, My;
lu = 1./(*N-1);
- qtemp = calloc( (*N) * (*N), sizeof(double));
+ qtemp = Calloc((*N) * (*N), double);
q[0] = 1;
Mx=1;
My=1;
@@ -272,7 +272,7 @@ void lossdistrib_prepay_joint(double *dp, double *pp, int *ndp, double *w, doubl }
q[MIN(*N, Mx)*MIN(My,*N)-1] += 1 - sum;
}
- free(qtemp);
+ Free(qtemp);
}
double shockprob(double p, double rho, double Z, int give_log){
@@ -311,9 +311,9 @@ void BClossdist(double *SurvProb, int *dim1, int *dim2, double *issuerweights, double g;
double *gshocked, *Rshocked, *Sshocked, *Lw, *Rw;
- gshocked = calloc((*dim1), sizeof(double));
- Rshocked = calloc((*dim1), sizeof(double));
- Sshocked = calloc((*dim1), sizeof(double));
+ gshocked = Calloc((*dim1), double);
+ Rshocked = Calloc((*dim1), double);
+ Sshocked = Calloc((*dim1), double);
Lw = malloc((*N) * sizeof(double));
Rw = malloc((*N) * sizeof(double));
@@ -334,9 +334,9 @@ void BClossdist(double *SurvProb, int *dim1, int *dim2, double *issuerweights, addandmultiply(Rw, w[i], R + t * (*N), *N);
}
}
- free(gshocked);
- free(Rshocked);
- free(Sshocked);
+ Free(gshocked);
+ Free(Rshocked);
+ Free(Sshocked);
free(Lw);
free(Rw);
}
|