64,654
社区成员
发帖
与我相关
我的任务
分享
#include <stdio.h>
#include <malloc.h>
double ****PsiDDiffj,
****PsiDDiffi,
****sDDiff,
r_h = 10;
int ****NorDiff,
numv = 100;
Generate_PsiDuopolyDiff() {
int aa, bb, cc, dd;
double vv_ja, vv_jb, vv_ia, vv_ib, deltav = r_h / numv, ll;
PsiDDiffj = (double ****)malloc((numv + 2) * sizeof(double***));
PsiDDiffi = (double ****)malloc((numv + 2) * sizeof(double***));
sDDiff = (double ****)malloc((numv + 2) * sizeof(double***));
NorDiff = (int ****)malloc((numv + 2) * sizeof(int***));
for (aa = 0; aa <= numv; aa++) {
PsiDDiffj[aa] = (double ***)malloc((numv + 2) * sizeof(double**));
PsiDDiffi[aa] = (double ***)malloc((numv + 2) * sizeof(double**));
sDDiff[aa] = (double ***)malloc((numv + 2) * sizeof(double**));
NorDiff[aa] = (int ***)malloc((numv + 2) * sizeof(int**));
for (bb = 0; bb <= numv; bb++) {
PsiDDiffj[aa][bb] = (double **)malloc((numv + 2) * sizeof(double*));
PsiDDiffi[aa][bb] = (double **)malloc((numv + 2) * sizeof(double*));
sDDiff[aa][bb] = (double **)malloc((numv + 2) * sizeof(double*));
NorDiff[aa][bb] = (int **)malloc((numv + 2) * sizeof(int*));
for (cc = 0; cc <= numv; cc++) {
PsiDDiffj[aa][bb][cc] = (double *)malloc((numv + 2) * sizeof(double )); //where the error happened
PsiDDiffi[aa][bb][cc] = (double *)malloc((numv + 2) * sizeof(double ));
sDDiff[aa][bb][cc] = (double *)malloc((numv + 2) * sizeof(double ));
NorDiff[aa][bb][cc] = (int *)malloc((numv + 2) * sizeof(int ));
}
}
}
}
void
main() {
Generate_PsiDuopolyDiff();
}