From 0f4a78cdbf66b3081941c63eed70cd6c42403a8d Mon Sep 17 00:00:00 2001
From: Thorsten Liebig <Thorsten.Liebig@gmx.de>
Date: Mon, 28 Jun 2010 19:45:19 +0200
Subject: [PATCH] new ProcessingIntegral as abstract base class to voltage and
 current processing

---
 FDTD/processcurrent.cpp  | 22 ++++----------------
 FDTD/processcurrent.h    | 11 +++-------
 FDTD/processing.cpp      | 10 +++++++++
 FDTD/processintegral.cpp | 44 ++++++++++++++++++++++++++++++++++++++++
 FDTD/processintegral.h   | 40 ++++++++++++++++++++++++++++++++++++
 FDTD/processvoltage.cpp  | 21 +++----------------
 FDTD/processvoltage.h    | 11 ++--------
 openEMS.pro              |  2 ++
 8 files changed, 108 insertions(+), 53 deletions(-)
 create mode 100644 FDTD/processintegral.cpp
 create mode 100644 FDTD/processintegral.h

diff --git a/FDTD/processcurrent.cpp b/FDTD/processcurrent.cpp
index f7517f2..49be042 100644
--- a/FDTD/processcurrent.cpp
+++ b/FDTD/processcurrent.cpp
@@ -20,13 +20,12 @@
 #include <iomanip>
 #include <complex.h>
 
-ProcessCurrent::ProcessCurrent(Operator* op, Engine* eng) : Processing(op, eng)
+ProcessCurrent::ProcessCurrent(Operator* op, Engine* eng) : ProcessIntegral(op, eng)
 {
 }
 
 ProcessCurrent::~ProcessCurrent()
 {
-	ProcessCurrent::FlushData();
 }
 
 void ProcessCurrent::DefineStartStopCoord(double* dstart, double* dstop)
@@ -46,16 +45,6 @@ void ProcessCurrent::DefineStartStopCoord(double* dstart, double* dstop)
 	}
 }
 
-void ProcessCurrent::InitProcess()
-{
-	m_filename = m_Name;
-	OpenFile(m_filename);
-	FD_currents.clear();
-	for (size_t n=0;n<m_FD_Samples.size();++n)
-		FD_currents.push_back(0);
-}
-
-
 int ProcessCurrent::Process()
 {
 	if (Enabled==false) return -1;
@@ -157,11 +146,12 @@ int ProcessCurrent::Process()
 	//	cerr << "ts: " << Eng->numTS << " i: " << current << endl;
 	current*=m_weight;
 
+
 	if (ProcessInterval)
 	{
 		if (Eng->GetNumberOfTimesteps()%ProcessInterval==0)
 		{
-			v_current.push_back(current);
+			TD_Values.push_back(current);
 			//current is sampled half a timestep later then the voltages
 			file  << setprecision(m_precision) << (0.5 + (double)Eng->GetNumberOfTimesteps())*Op->GetTimestep() << "\t" << current << endl;
 		}
@@ -174,7 +164,7 @@ int ProcessCurrent::Process()
 			double T = ((double)Eng->GetNumberOfTimesteps() + 0.5) * Op->GetTimestep();
 			for (size_t n=0;n<m_FD_Samples.size();++n)
 			{
-				FD_currents.at(n) += current * cexp( -2.0 * 1.0i * M_PI * m_FD_Samples.at(n) * T );
+				FD_Values.at(n) += current * cexp( -2.0 * 1.0i * M_PI * m_FD_Samples.at(n) * T );
 			}
 			++m_FD_SampleCount;
 			if (m_Flush)
@@ -191,7 +181,3 @@ void ProcessCurrent::DumpBox2File( string vtkfilenameprefix, bool /*dualMesh*/ )
 	Processing::DumpBox2File( vtkfilenameprefix, true );
 }
 
-void ProcessCurrent::FlushData()
-{
-	Dump_FD_Data(FD_currents,1.0/(double)m_FD_SampleCount,m_filename + "_FD");
-}
diff --git a/FDTD/processcurrent.h b/FDTD/processcurrent.h
index 2c30cec..b3bf1c5 100644
--- a/FDTD/processcurrent.h
+++ b/FDTD/processcurrent.h
@@ -18,9 +18,9 @@
 #ifndef PROCESSCURRENT_H
 #define PROCESSCURRENT_H
 
-#include "processing.h"
+#include "processintegral.h"
 
-class ProcessCurrent : public Processing
+class ProcessCurrent : public ProcessIntegral
 {
 public:
 	ProcessCurrent(Operator* op, Engine* eng);
@@ -28,16 +28,11 @@ public:
 
 	virtual void DefineStartStopCoord(double* dstart, double* dstop);
 
-	virtual void InitProcess();
-	virtual void FlushData();
 	virtual int Process();
 
 	virtual void DumpBox2File( string vtkfilenameprefix, bool dualMesh = false ) const; //!< dump geometry to file
-protected:
-	vector<FDTD_FLOAT> v_current;
 
-	vector<_Complex double> FD_currents;
-	void WriteFDCurrents();
+protected:
 };
 
 #endif // PROCESSCURRENT_H
diff --git a/FDTD/processing.cpp b/FDTD/processing.cpp
index 1856670..ee7bbb5 100644
--- a/FDTD/processing.cpp
+++ b/FDTD/processing.cpp
@@ -275,11 +275,21 @@ void Processing::DumpBox2File( string vtkfilenameprefix, bool dualMesh ) const
 
 void Processing::Dump_FD_Data(vector<_Complex double> value, double factor, string filename)
 {
+	if (value.size()==0) 
+		return;
+	if (value.size()!=m_FD_Samples.size()) 
+	{
+		cerr << "Processing::Dump_FD_Data: Error: Complex value and frequency vector have different size! This should never happend!!!" << endl;
+		return;
+	}
 	ofstream file;
 	file.open( filename.c_str() );
 	if (!file.is_open())
 		cerr << "Can't open file: " << filename << endl;
 
+	time_t rawTime;
+	time(&rawTime);
+	file << "%dump by openEMS @" << ctime(&rawTime) << "%frequency\treal\timag\n";
 	for (size_t n=0;n<value.size();++n)
 	{
 		file << m_FD_Samples.at(n) << "\t" << 2.0 * creal(value.at(n))*factor << "\t" << 2.0 * cimag(value.at(n))*factor << "\n";
diff --git a/FDTD/processintegral.cpp b/FDTD/processintegral.cpp
new file mode 100644
index 0000000..8980711
--- /dev/null
+++ b/FDTD/processintegral.cpp
@@ -0,0 +1,44 @@
+/*
+*	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/>.
+*/
+
+#include "processintegral.h"
+
+ProcessIntegral::ProcessIntegral(Operator* op, Engine* eng)  : Processing(op, eng)
+{
+}
+
+ProcessIntegral::~ProcessIntegral()
+{
+	ProcessIntegral::FlushData();
+}
+
+
+void ProcessIntegral::InitProcess()
+{
+	m_filename = m_Name;
+	OpenFile(m_filename);
+	FD_Values.clear();
+	for (size_t n=0;n<m_FD_Samples.size();++n)
+		FD_Values.push_back(0);
+}
+
+void ProcessIntegral::FlushData()
+{
+	if (m_FD_Samples.size())
+		Dump_FD_Data(FD_Values,1.0/(double)m_FD_SampleCount,m_filename + "_FD");
+}
+
diff --git a/FDTD/processintegral.h b/FDTD/processintegral.h
new file mode 100644
index 0000000..b73f35b
--- /dev/null
+++ b/FDTD/processintegral.h
@@ -0,0 +1,40 @@
+/*
+*	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 PROCESSINTEGRAL_H
+#define PROCESSINTEGRAL_H
+
+#include "processing.h"
+
+//! Abstract base class for integral parameter processing
+class ProcessIntegral : public Processing
+{
+public:
+	virtual ~ProcessIntegral();
+
+	virtual void InitProcess();
+	virtual void FlushData();
+
+protected:
+	ProcessIntegral(Operator* op, Engine* eng);
+
+	vector<FDTD_FLOAT> TD_Values;
+
+	vector<_Complex double> FD_Values;
+};
+
+#endif // PROCESSINTEGRAL_H
diff --git a/FDTD/processvoltage.cpp b/FDTD/processvoltage.cpp
index 20e1805..014a3c0 100644
--- a/FDTD/processvoltage.cpp
+++ b/FDTD/processvoltage.cpp
@@ -19,22 +19,12 @@
 #include <complex.h>
 #include <iomanip>
 
-ProcessVoltage::ProcessVoltage(Operator* op, Engine* eng) : Processing(op, eng)
+ProcessVoltage::ProcessVoltage(Operator* op, Engine* eng) : ProcessIntegral(op, eng)
 {
 }
 
 ProcessVoltage::~ProcessVoltage()
 {
-	ProcessVoltage::FlushData();
-}
-
-void ProcessVoltage::InitProcess()
-{
-	m_filename = m_Name;
-	OpenFile(m_filename);
-	FD_voltages.clear();
-	for (size_t n=0;n<m_FD_Samples.size();++n)
-		FD_voltages.push_back(0);
 }
 
 int ProcessVoltage::Process()
@@ -49,7 +39,7 @@ int ProcessVoltage::Process()
 	{
 		if (Eng->GetNumberOfTimesteps()%ProcessInterval==0)
 		{
-			voltages.push_back(voltage);
+			TD_Values.push_back(voltage);
 			file << setprecision(m_precision) << (double)Eng->GetNumberOfTimesteps()*Op->GetTimestep() << "\t" << voltage << endl;
 		}
 	}
@@ -61,7 +51,7 @@ int ProcessVoltage::Process()
 			double T = (double)Eng->GetNumberOfTimesteps() * Op->GetTimestep();
 			for (size_t n=0;n<m_FD_Samples.size();++n)
 			{
-				FD_voltages.at(n) += voltage * cexp( -2.0 * 1.0i * M_PI * m_FD_Samples.at(n) * T );
+				FD_Values.at(n) += voltage * cexp( -2.0 * 1.0i * M_PI * m_FD_Samples.at(n) * T );
 			}
 			++m_FD_SampleCount;
 			if (m_Flush)
@@ -72,8 +62,3 @@ int ProcessVoltage::Process()
 
 	return GetNextInterval();
 }
-
-void ProcessVoltage::FlushData()
-{
-	Dump_FD_Data(FD_voltages,1.0/(double)m_FD_SampleCount,m_filename + "_FD");
-}
diff --git a/FDTD/processvoltage.h b/FDTD/processvoltage.h
index b6cb344..3c3dc03 100644
--- a/FDTD/processvoltage.h
+++ b/FDTD/processvoltage.h
@@ -18,25 +18,18 @@
 #ifndef PROCESSVOLTAGE_H
 #define PROCESSVOLTAGE_H
 
-#include "processing.h"
+#include "processintegral.h"
 
 //! Process voltage along a line from start to stop coordinates. ATM integration along the axis e.g.: in x, then y then z direction (Future: diagonal integration)
-class ProcessVoltage : public Processing
+class ProcessVoltage : public ProcessIntegral
 {
 public:
 	ProcessVoltage(Operator* op, Engine* eng);
 	virtual ~ProcessVoltage();
 
-	virtual void InitProcess();
-	virtual void FlushData();
-
 	virtual int Process();
 
 protected:
-	vector<FDTD_FLOAT> voltages;
-
-	vector<_Complex double> FD_voltages;
-	void WriteFDVoltages();
 };
 
 #endif // PROCESSVOLTAGE_H
diff --git a/openEMS.pro b/openEMS.pro
index 3df9bad..e4a8e48 100644
--- a/openEMS.pro
+++ b/openEMS.pro
@@ -38,6 +38,7 @@ SOURCES += main.cpp \
     tools/array_ops.cpp \
     FDTD/processvoltage.cpp \
     FDTD/processing.cpp \
+    FDTD/processintegral.cpp \
     FDTD/processfields.cpp \
     FDTD/processfields_td.cpp \
     FDTD/processcurrent.cpp \
@@ -66,6 +67,7 @@ HEADERS += tools/ErrorMsg.h \
     tools/array_ops.h \
     FDTD/processvoltage.h \
     FDTD/processing.h \
+    FDTD/processintegral.h \
     FDTD/processfields.h \
     FDTD/processfields_td.h \
     FDTD/processcurrent.h \