openEMS/FDTD/extensions/operator_ext_mur_abc.h

/*
*	Copyright (C) 2010 Thorsten Liebig (Thorsten.Liebig@gmx.de)
*
*	This program is free software: you can redistribute it and/or modify
*	it under the terms of the GNU General Public License as published by
*	the Free Software Foundation, either version 3 of the License, or
*	(at your option) any later version.
*
*	This program is distributed in the hope that it will be useful,
*	but WITHOUT ANY WARRANTY; without even the implied warranty of
*	MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
*	GNU General Public License for more details.
*
*	You should have received a copy of the GNU General Public License
*	along with this program.  If not, see <http://www.gnu.org/licenses/>.
*/

#ifndef OPERATOR_EXT_MUR_ABC_H
#define OPERATOR_EXT_MUR_ABC_H

#include "FDTD/operator.h"
#include "operator_extension.h"

class Operator_Ext_Mur_ABC : public Operator_Extension
{
	friend class Engine_Ext_Mur_ABC;
public:
	Operator_Ext_Mur_ABC(Operator* op);
	~Operator_Ext_Mur_ABC();

	virtual Operator_Extension* Clone(Operator* op);

	//! Define the direction of this ABC: \a ny=0,1,2 -> x,y,z and if at bottom_ny -> e.g. x=0 or x=end
	void SetDirection(int ny, bool top_ny);

	//! Set (override) the expected phase velocity of the incoming wave
	void SetPhaseVelocity(double c_phase) {m_v_phase=c_phase;};

	virtual bool BuildExtension();

	virtual Engine_Extension* CreateEngineExtention();

	virtual bool IsCylinderCoordsSave(bool closedAlpha, bool R0_included) const;
	virtual bool IsCylindricalMultiGridSave(bool child) const;
	virtual bool IsMPISave() const {return true;}

	virtual string GetExtensionName() const {return string("Mur ABC Extension");}

	virtual void ShowStat(ostream &ostr) const;

protected:
	Operator_Ext_Mur_ABC(Operator* op, Operator_Ext_Mur_ABC* op_ext);
	void Initialize();
	int m_ny;
	int m_nyP,m_nyPP;
	bool m_top;
	unsigned int m_LineNr;
	int m_LineNr_Shift;

	double m_v_phase;

	unsigned int m_numLines[2];

	FDTD_FLOAT** m_Mur_Coeff_nyP;
	FDTD_FLOAT** m_Mur_Coeff_nyPP;
};

#endif // OPERATOR_EXT_MUR_ABC_H