?? cpmlfdtd3d.c
字號:
}
psi_Hyz_1 = (float ***)malloc((Imax-1) * sizeof(float **));
for(i = 0; i < Imax-1; i++) {
psi_Hyz_1[i] = (float **)malloc(Jmax * sizeof(float *));
for(j = 0; j < Jmax; j++) {
psi_Hyz_1[i][j] = (float *)malloc((nzPML_1-1) * sizeof(float));
for(k = 0; k < nzPML_1-1; k++) {
psi_Hyz_1[i][j][k] = 0.0;
}
}
}
psi_Hyz_2 = (float ***)malloc((Imax-1) * sizeof(float **));
for(i = 0; i < Imax-1; i++) {
psi_Hyz_2[i] = (float **)malloc(Jmax * sizeof(float *));
for(j = 0; j < Jmax; j++) {
psi_Hyz_2[i][j] = (float *)malloc((nzPML_2-1) * sizeof(float));
for(k = 0; k < nzPML_2-1; k++) {
psi_Hyz_2[i][j][k] = 0.0;
}
}
}
psi_Exz_1 = (float ***)malloc((Imax-1) * sizeof(float **));
for(i = 0; i < Imax-1; i++) {
psi_Exz_1[i] = (float **)malloc(Jmax * sizeof(float *));
for(j = 0; j < Jmax; j++) {
psi_Exz_1[i][j] = (float *)malloc(nzPML_1 * sizeof(float));
for(k = 0; k < nzPML_1; k++) {
psi_Exz_1[i][j][k] = 0.0;
}
}
}
psi_Exz_2 = (float ***)malloc((Imax-1) * sizeof(float **));
for(i = 0; i < Imax-1; i++) {
psi_Exz_2[i] = (float **)malloc(Jmax * sizeof(float *));
for(j = 0; j < Jmax; j++) {
psi_Exz_2[i][j] = (float *)malloc(nzPML_2 * sizeof(float));
for(k = 0; k < nzPML_2; k++) {
psi_Exz_2[i][j][k] = 0.0;
}
}
}
psi_Eyz_1 = (float ***)malloc((Imax-1) * sizeof(float **));
for(i = 0; i < Imax; i++) {
psi_Eyz_1[i] = (float **)malloc((Jmax-1) * sizeof(float *));
for(j = 0; j < Jmax-1; j++) {
psi_Eyz_1[i][j] = (float *)malloc(nzPML_1 * sizeof(float));
for(k = 0; k < nzPML_1; k++) {
psi_Eyz_1[i][j][k] = 0.0;
}
}
}
psi_Eyz_2 = (float ***)malloc((Imax-1) * sizeof(float **));
for(i = 0; i < Imax; i++) {
psi_Eyz_2[i] = (float **)malloc((Jmax-1) * sizeof(float *));
for(j = 0; j < Jmax-1; j++) {
psi_Eyz_2[i][j] = (float *)malloc(nzPML_2 * sizeof(float));
for(k = 0; k < nzPML_2; k++) {
psi_Eyz_2[i][j][k] = 0.0;
}
}
}
psi_Hzx_1 = (float ***)malloc((nxPML_1-1) * sizeof(float **));
for(i = 0; i < nxPML_1-1; i++) {
psi_Hzx_1[i] = (float **)malloc((Jmax-1) * sizeof(float *));
for(j = 0; j < Jmax-1; j++) {
psi_Hzx_1[i][j] = (float *)malloc((Kmax-1) * sizeof(float));
for(k = 0; k < Kmax-1; k++) {
psi_Hzx_1[i][j][k] = 0.0;
}
}
}
psi_Hzx_2 = (float ***)malloc((nxPML_2-1) * sizeof(float **));
for(i = 0; i < nxPML_2-1; i++) {
psi_Hzx_2[i] = (float **)malloc((Jmax-1) * sizeof(float *));
for(j = 0; j < Jmax-1; j++) {
psi_Hzx_2[i][j] = (float *)malloc((Kmax-1) * sizeof(float));
for(k = 0; k < Kmax-1; k++) {
psi_Hzx_2[i][j][k] = 0.0;
}
}
}
psi_Eyx_1 = (float ***)malloc(nxPML_1 * sizeof(float **));
for(i = 0; i < nxPML_1; i++) {
psi_Eyx_1[i] = (float **)malloc((Jmax-1) * sizeof(float *));
for(j = 0; j < Jmax-1; j++) {
psi_Eyx_1[i][j] = (float *)malloc((Kmax-1) * sizeof(float));
for(k = 0; k < Kmax-1; k++) {
psi_Eyx_1[i][j][k] = 0.0;
}
}
}
psi_Eyx_2 = (float ***)malloc(nxPML_2 * sizeof(float **));
for(i = 0; i < nxPML_2; i++) {
psi_Eyx_2[i] = (float **)malloc((Jmax-1) * sizeof(float *));
for(j = 0; j < Jmax-1; j++) {
psi_Eyx_2[i][j] = (float *)malloc((Kmax-1) * sizeof(float));
for(k = 0; k < Kmax-1; k++) {
psi_Eyx_2[i][j][k] = 0.0;
}
}
}
psi_Hzy_1 = (float ***)malloc((Imax-1) * sizeof(float **));
for(i = 0; i < Imax-1; i++) {
psi_Hzy_1[i] = (float **)malloc((nyPML_1-1) * sizeof(float *));
for(j = 0; j < nyPML_1-1; j++) {
psi_Hzy_1[i][j] = (float *)malloc((Kmax-1) * sizeof(float));
for(k = 0; k < Kmax-1; k++) {
psi_Hzy_1[i][j][k] = 0.0;
}
}
}
psi_Hzy_2 = (float ***)malloc((Imax-1) * sizeof(float **));
for(i = 0; i < Imax-1; i++) {
psi_Hzy_2[i] = (float **)malloc((nyPML_2-1) * sizeof(float *));
for(j = 0; j < nyPML_2-1; j++) {
psi_Hzy_2[i][j] = (float *)malloc((Kmax-1) * sizeof(float));
for(k = 0; k < Kmax-1; k++) {
psi_Hzy_2[i][j][k] = 0.0;
}
}
}
psi_Exy_1 = (float ***)malloc((Imax-1) * sizeof(float **));
for(i = 0; i < Imax-1; i++) {
psi_Exy_1[i] = (float **)malloc(nyPML_1 * sizeof(float *));
for(j = 0; j < nyPML_1; j++) {
psi_Exy_1[i][j] = (float *)malloc((Kmax-1) * sizeof(float));
for(k = 0; k < Kmax-1; k++) {
psi_Exy_1[i][j][k] = 0.0;
}
}
}
psi_Exy_2 = (float ***)malloc((Imax-1) * sizeof(float **));
for(i = 0; i < Imax-1; i++) {
psi_Exy_2[i] = (float **)malloc(nyPML_2 * sizeof(float *));
for(j = 0; j < nyPML_2; j++) {
psi_Exy_2[i][j] = (float *)malloc((Kmax-1) * sizeof(float));
for(k = 0; k < Kmax-1; k++) {
psi_Exy_2[i][j][k] = 0.0;
}
}
}
den_ex = (float *)malloc((Imax-1) * sizeof(float));
for(i = 0; i < Imax-1; i++) {
den_ex[i] = 0.0;
}
den_hx = (float *)malloc((Imax-1) * sizeof(float));
for(i = 0; i < Imax-1; i++) {
den_hx[i] = 0.0;
}
den_ey = (float *)malloc((Jmax-1) * sizeof(float));
for(i = 0; i < Jmax-1; i++) {
den_ey[i] = 0.0;
}
den_hy = (float *)malloc((Jmax-1) * sizeof(float));
for(i = 0; i < Jmax-1; i++) {
den_hy[i] = 0.0;
}
den_ez = (float *)malloc((Kmax-1) * sizeof(float));
for(i = 0; i < Kmax-1; i++) {
den_ez[i] = 0.0;
}
den_hz = (float *)malloc((Kmax-1) * sizeof(float));
for(i = 0; i < Kmax-1; i++) {
den_hz[i] = 0.0;
}
be_x_1 = (float *)malloc((nxPML_1) * sizeof(float));
for(i = 0; i < nxPML_1; i++) {
be_x_1[i] = 0.0;
}
ce_x_1 = (float *)malloc((nxPML_1) * sizeof(float));
for(i = 0; i < nxPML_1; i++) {
ce_x_1[i] = 0.0;
}
alphae_x_PML_1 = (float *)malloc((nxPML_1) * sizeof(float));
for(i = 0; i < nxPML_1; i++) {
alphae_x_PML_1[i] = 0.0;
}
sige_x_PML_1 = (float *)malloc((nxPML_1) * sizeof(float));
for(i = 0; i < nxPML_1; i++) {
sige_x_PML_1[i] = 0.0;
}
kappae_x_PML_1 = (float *)malloc((nxPML_1) * sizeof(float));
for(i = 0; i < nxPML_1; i++) {
kappae_x_PML_1[i] = 0.0;
}
bh_x_1 = (float *)malloc((nxPML_1-1) * sizeof(float));
for(i = 0; i < nxPML_1-1; i++) {
bh_x_1[i] = 0.0;
}
ch_x_1 = (float *)malloc((nxPML_1-1) * sizeof(float));
for(i = 0; i < nxPML_1-1; i++) {
ch_x_1[i] = 0.0;
}
alphah_x_PML_1 = (float *)malloc((nxPML_1-1) * sizeof(float));
for(i = 0; i < nxPML_1-1; i++) {
alphah_x_PML_1[i] = 0.0;
}
sigh_x_PML_1 = (float *)malloc((nxPML_1-1) * sizeof(float));
for(i = 0; i < nxPML_1-1; i++) {
sigh_x_PML_1[i] = 0.0;
}
kappah_x_PML_1 = (float *)malloc((nxPML_1-1) * sizeof(float));
for(i = 0; i < nxPML_1-1; i++) {
kappah_x_PML_1[i] = 0.0;
}
be_x_2 = (float *)malloc((nxPML_2) * sizeof(float));
for(i = 0; i < nxPML_2; i++) {
be_x_2[i] = 0.0;
}
ce_x_2 = (float *)malloc((nxPML_2) * sizeof(float));
for(i = 0; i < nxPML_2; i++) {
ce_x_2[i] = 0.0;
}
alphae_x_PML_2 = (float *)malloc((nxPML_2) * sizeof(float));
for(i = 0; i < nxPML_2; i++) {
alphae_x_PML_2[i] = 0.0;
}
sige_x_PML_2 = (float *)malloc((nxPML_2) * sizeof(float));
for(i = 0; i < nxPML_2; i++) {
sige_x_PML_2[i] = 0.0;
}
kappae_x_PML_2 = (float *)malloc((nxPML_2) * sizeof(float));
for(i = 0; i < nxPML_2; i++) {
kappae_x_PML_2[i] = 0.0;
}
bh_x_2 = (float *)malloc((nxPML_2-1) * sizeof(float));
for(i = 0; i < nxPML_2-1; i++) {
bh_x_2[i] = 0.0;
}
ch_x_2 = (float *)malloc((nxPML_2-1) * sizeof(float));
for(i = 0; i < nxPML_2-1; i++) {
ch_x_2[i] = 0.0;
}
alphah_x_PML_2 = (float *)malloc((nxPML_2-1) * sizeof(float));
for(i = 0; i < nxPML_2-1; i++) {
alphah_x_PML_2[i] = 0.0;
}
sigh_x_PML_2 = (float *)malloc((nxPML_2-1) * sizeof(float));
for(i = 0; i < nxPML_2-1; i++) {
sigh_x_PML_2[i] = 0.0;
}
kappah_x_PML_2 = (float *)malloc((nxPML_2-1) * sizeof(float));
for(i = 0; i < nxPML_1-1; i++) {
kappah_x_PML_2[i] = 0.0;
}
be_y_1 = (float *)malloc((nyPML_1) * sizeof(float));
for(i = 0; i < nyPML_1; i++) {
be_y_1[i] = 0.0;
}
ce_y_1 = (float *)malloc((nyPML_1) * sizeof(float));
for(i = 0; i < nyPML_1; i++) {
ce_y_1[i] = 0.0;
}
alphae_y_PML_1 = (float *)malloc((nyPML_1) * sizeof(float));
for(i = 0; i < nyPML_1; i++) {
alphae_y_PML_1[i] = 0.0;
}
sige_y_PML_1 = (float *)malloc((nyPML_1) * sizeof(float));
for(i = 0; i < nyPML_1; i++) {
sige_y_PML_1[i] = 0.0;
}
kappae_y_PML_1 = (float *)malloc((nyPML_1) * sizeof(float));
for(i = 0; i < nyPML_1; i++) {
kappae_y_PML_1[i] = 0.0;
}
bh_y_1 = (float *)malloc((nyPML_1-1) * sizeof(float));
for(i = 0; i < nyPML_1-1; i++) {
bh_y_1[i] = 0.0;
}
ch_y_1 = (float *)malloc((nyPML_1-1) * sizeof(float));
for(i = 0; i < nyPML_1-1; i++) {
ch_y_1[i] = 0.0;
}
alphah_y_PML_1 = (float *)malloc((nyPML_1-1) * sizeof(float));
for(i = 0; i < nyPML_1-1; i++) {
alphah_y_PML_1[i] = 0.0;
}
sigh_y_PML_1 = (float *)malloc((nyPML_1-1) * sizeof(float));
for(i = 0; i < nyPML_1-1; i++) {
sigh_y_PML_1[i] = 0.0;
}
kappah_y_PML_1 = (float *)malloc((nyPML_1-1) * sizeof(float));
for(i = 0; i < nyPML_1-1; i++) {
kappah_y_PML_1[i] = 0.0;
}
be_y_2 = (float *)malloc((nyPML_2) * sizeof(float));
for(i = 0; i < nyPML_2; i++) {
be_y_2[i] = 0.0;
}
ce_y_2 = (float *)malloc((nyPML_2) * sizeof(float));
for(i = 0; i < nyPML_2; i++) {
ce_y_2[i] = 0.0;
}
alphae_y_PML_2 = (float *)malloc((nyPML_2) * sizeof(float));
for(i = 0; i < nyPML_2; i++) {
alphae_y_PML_2[i] = 0.0;
}
sige_y_PML_2 = (float *)malloc((nyPML_2) * sizeof(float));
for(i = 0; i < nyPML_2; i++) {
sige_y_PML_2[i] = 0.0;
}
kappae_y_PML_2 = (float *)malloc((nyPML_2) * sizeof(float));
for(i = 0; i < nyPML_2; i++) {
kappae_y_PML_2[i] = 0.0;
}
bh_y_2 = (float *)malloc((nyPML_2-1) * sizeof(float));
for(i = 0; i < nyPML_2-1; i++) {
bh_y_2[i] = 0.0;
}
ch_y_2 = (float *)malloc((nyPML_2-1) * sizeof(float));
for(i = 0; i < nyPML_2-1; i++) {
ch_y_2[i] = 0.0;
}
alphah_y_PML_2 = (float *)malloc((nyPML_2-1) * sizeof(float));
for(i = 0; i < nyPML_2-1; i++) {
alphah_y_PML_2[i] = 0.0;
}
sigh_y_PML_2 = (float *)malloc((nyPML_2-1) * sizeof(float));
for(i = 0; i < nyPML_2-1; i++) {
sigh_y_PML_2[i] = 0.0;
}
kappah_y_PML_2 = (float *)malloc((nyPML_2-1) * sizeof(float));
for(i = 0; i < nyPML_1-1; i++) {
kappah_y_PML_2[i] = 0.0;
}
be_z_1 = (float *)malloc((nzPML_1) * sizeof(float));
for(i = 0; i < nzPML_1; i++) {
be_z_1[i] = 0.0;
}
ce_z_1 = (float *)malloc((nzPML_1) * sizeof(float));
for(i = 0; i < nzPML_1; i++) {
ce_z_1[i] = 0.0;
}
alphae_z_PML_1 = (float *)malloc((nzPML_1) * sizeof(float));
for(i = 0; i < nzPML_1; i++) {
alphae_z_PML_1[i] = 0.0;
}
sige_z_PML_1 = (float *)malloc((nzPML_1) * sizeof(float));
for(i = 0; i < nzPML_1; i++) {
sige_z_PML_1[i] = 0.0;
}
kappae_z_PML_1 = (float *)malloc((nzPML_1) * sizeof(float));
for(i = 0; i < nzPML_1; i++) {
kappae_z_PML_1[i] = 0.0;
}
bh_z_1 = (float *)malloc((nzPML_1-1) * sizeof(float));
for(i = 0; i < nzPML_1-1; i++) {
?? 快捷鍵說明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -