cpmlfdtd3d.c

FDTD

		}

	    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++) {