powerbudget.h

矩量法仿真电磁辐射和散射的源代码（c++）

/*
	Copyright (C) 2004-2005  Timothy C.A. Molteno
	
	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 2 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, write to the Free Software
	Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
*/
#ifndef __PowerBudget__
#define __PowerBudget__

#include <vector>
#include "BaseInput.h"

/*
			---------- POWER BUDGET ---------
			INPUT POWER   =  3.3203E-05 Watts
			RADIATED POWER=  3.3203E-05 Watts
			STRUCTURE LOSS=  0.0000E+00 Watts
			NETWORK LOSS  =  0.0000E+00 Watts
			EFFICIENCY    =  100.00 Percent
*/
class PowerBudget : public BaseInput
{
public:	
	vector<double> input_power, radiated_power, structure_loss, network_loss, efficiency;
	long n_items;
	
	double get_power_line()
	{
		string line = readline();
		if (line.length() < 2)
			line = readline();
			
		// delete up to '=' character (+1 includes the = character)
		line.erase(0, line.find("=",0) + 1);
		stringstream ss(line);
		double ret = read_sci(ss);
		return ret;
	}
	
	PowerBudget(std::string& filename)
		: BaseInput(filename)
	{
		n_items = 0;
		string searchString("POWER BUDGET");
		while (m_stream.good())
		{
		string line = readline();

		if (line.find(searchString,0) != string::npos)
		{
			// check for a blank line here.
			
			{
				double p = get_power_line();
				input_power.push_back(p);
			}
			{
				double p = get_power_line();
				radiated_power.push_back(p);
			}
			{
				double p = get_power_line();
				structure_loss.push_back(p);
			}
			{
				double p = get_power_line();
				network_loss.push_back(p);
			}
			{
				double p = get_power_line();
				efficiency.push_back(p);
			}
			n_items++;
		}	
		}
	}

	bool equalto(const PowerBudget& pb)
	{
		if (difference(pb) > 1e-4)
			return false;

		return true;
	}

	double difference(const PowerBudget& pb)
	{
		double ret = 0;

		if (n_items != pb.n_items)
			return 1.0;

		try
		{
		for (long i=0; i<n_items; i++)
		{
			ret += diff(pb.input_power[i],input_power[i]);
			ret += diff(pb.radiated_power[i],radiated_power[i]);
			ret += diff(pb.structure_loss[i],structure_loss[i]);
			ret += diff(pb.network_loss[i],network_loss[i]);
			ret += diff(pb.efficiency[i],efficiency[i]);
		}
		}
		catch(string message)
		{
			cout << "diff : " << message << endl;
		}
		return ret;
	};
};

#endif /* __PowerBudget__ */