电磁场完美匹配层PML资料汇总

// FDTD_1D_HzEy.cpp: implementation of the FDTD_1D_HzEy class.
//
//////////////////////////////////////////////////////////////////////

#include “stdafx.h“
#include “fdtd_2D_TE_PML_a.h“
#include “FDTD_1D_HzEy.h“
#include “Math.h“

#ifdef _DEBUG
#undef THIS_FILE
static char THIS_FILE[]=__FILE__;
#define new DEBUG_NEW
#endif

//////////////////////////////////////////////////////////////////////
// Construction/Destruction
//////////////////////////////////////////////////////////////////////
FDTD_1D_HzEy::FDTD_1D_HzEy（）
{
	Ey_1D = NULL; Hz_1D = NULL;
	K_E1_a = NULL; K_E1_b = NULL; K_E2_a = NULL; K_E2_b = NULL;
	pi = 3.1415926535897932384626433832795;
	eps_0 = 8.854e-12; // [F/m]
	mu_0 = 4*pi*1e-7; // [H/m]
}

FDTD_1D_HzEy::~FDTD_1D_HzEy（）
{
	if（Hz_1D）
	    free（Hz_1D）; 	
	if（Ey_1D）
	    free（Ey_1D）; 
	if（K_E1_a）
		free（K_E1_a）;
	if（K_E1_b）
		free（K_E1_b）;
	if（K_E2_a）
		free（K_E2_a）;
	if（K_E2_b）
		free（K_E2_b）;
}

///////////////////////////////////////////////////////////////////////////////////////
//Init Main Parameters
///////////////////////////////////////////////////////////////////////////////////////
BOOL FDTD_1D_HzEy::Init_Main_Param_1D（int n_x int n_pml double Eps_r double Mu_r double d_t 
									  double d_x）
{
	//dimension of the computational space
    nx = n_x;
	eps_r = Eps_r;
	mu_r = Mu_r;
	dt = d_t;
	dx = d_x;
	n_PML = n_pml;
	
	Hz_1D = （double *） calloc（nxsizeof（double））;
	if（!Hz_1D）
	{
		return FALSE;
	}
	
	Ey_1D = （double *） calloc（nx-1sizeof（double））;
	if（!Ey_1D）
	{
		free（Hz_1D）;
		return FALSE;
	}
	
	K_E1_a = （double *） calloc（nxsizeof（double））;
	if（!K_E1_a）
	{
		free（Hz_1D）;
		free（Ey_1D）;
		return FALSE;
	}
	
	K_E1_b = （double *） calloc（nxsizeof（double））;
	if（!K_E1_b）
	{
		free（Hz_1D）;
		free（Ey_1D）;
		free（K_E1_a）;
		return FALSE;
	}
	
	K_E2_a = （double *） calloc（nxsizeof（double））;
	if（!K_E2_a）
	{
		free（Hz_1D）;
		free（Ey_1D）;
		free（K_E1_a）;
		free（K_E1_b）;
		return FALSE;
	}
	
	K_E2_b = （double *） calloc（nxsizeof（double））;
	if（!K_E2_b）
	{
		free（Hz_1D）;
		free（Ey_1D）;
		free（K_E1_a）;
		free（K_E1_b）;
		free（K_E2_a）;
		return FALSE;
	}
	
	return TRUE;

}

//////////////////////////////////////////////////////////////////////
//Initialize the parameters of a Gaussian pulse
//////////////////////////////////////////////////////////////////////
void FDTD_1D_HzEy::Init_Gauss_1D（double H_0 double t_0 double t_w）
{
	source_type = 1;
	H0 = H_0;
	t0 = t_0;
	tw = t_w;
}

//////////////////////////////////////////////////////////////////////
//Initialize the parameters of a Sinusoidal plane wave
//////////////////////////////////////////////////////////////////////
void FDTD_1D_HzEy::Init_Sinus_1D（double H_0 double om double Phi）
{
	source_type = 2;
	H0 = H_0;
	omega = om;
	phi = Phi;
}


//////////////////////////////////////////////////////////////////////
//Initialize the param