fdtd_1d_ezhy_pml_loss.h

2维时域有限差分模拟

#pragma once

class CFDTD_1D_EZHY_PML_LOSS
{
private:
	double pi, eps_0, mu_0;
	
	double *Fz, *Ez, *Hy;
	double *K_Fz_a, *K_Fz_b, *K_Hy_a, *K_Hy_b;
	
	double K_Ez_a, K_Ez_b;
	
	int nx; //Dimension of the simulation space
	double dt, dx;
	int n_PML;

	double eps_r, mu_r, sigma;

	double E0, omega, phi, t0, tw, alfa;
	int Source_Type;

	int jel_source; // 1 - Gauss, 2 - Sin, 3 - SinGauss

public:
	CFDTD_1D_EZHY_PML_LOSS(void);
	~CFDTD_1D_EZHY_PML_LOSS(void);

	bool Init(int n_x, int n_pml, double d_t, double d_x, double eps_r, double mu_r, 
		      double sigma);

	void Init_PML_Param(void);

	void Init_Sin_Source(double E_0, double om, double Phi);
	void Init_Gauss_Source(double E_0, double t_0, double t_w);
	void Init_GaussSin(double E_0, double om, double Phi, double t_0, double t_w);

	void calc_Ez(double time);
	void calc_Hy(void);

	void Get_Data(double *&ez, double *&hy);

	void Free_Mem(void);
};