zcy
/
openEMS
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

write H-field excitation into file "ht" 

If we have a separate signal for H-field excitation, we need to dump it.
The question is, do we really need a separate Signal?
pull/1/head
Sebastian Held 2010-07-16 10:33:54 +02:00
parent f7edb2f1fd
commit 5e904b112c
2 changed files with 12 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
FDTD/engine.cpp
View File

@ -54,6 +54,7 @@ void Engine::Init()
curr = Create_N_3DArray<FDTD_FLOAT>(numLines);
file_et.open( "et" );
file_ht.open( "ht" );
InitExtensions();
}
@ -80,6 +81,7 @@ void Engine::Reset()
curr=NULL;
file_et.close();
file_ht.close();
for (size_t n=0;n<m_Eng_exts.size();++n)
delete m_Eng_exts.at(n);
@ -132,7 +134,7 @@ void Engine::ApplyVoltageExcite()
GetVolt(Op->Exc->E_dir[n],Op->Exc->E_index[0][n],Op->Exc->E_index[1][n],Op->Exc->E_index[2][n]) += Op->Exc->E_amp[n]*Op->Exc->Signal_volt[exc_pos];
}
// write the first excitation into the file "et1"
// write the first excitation into the file "et"
if (numTS < Op->Exc->Length)
file_et << numTS * Op->GetTimestep() << "\t" << Op->Exc->Signal_volt[numTS] << "\n"; // do not use std::endl here, because it will do an implicit flush
else
@ -177,7 +179,15 @@ void Engine::ApplyCurrentExcite()
exc_pos *= (exc_pos>0 && exc_pos<=(int)Op->Exc->Length);
// if (n==0) cerr << numTS << " => " << Op->ExciteSignal[exc_pos] << endl;
GetCurr(Op->Exc->Curr_dir[n],Op->Exc->Curr_index[0][n],Op->Exc->Curr_index[1][n],Op->Exc->Curr_index[2][n]) += Op->Exc->Curr_amp[n]*Op->Exc->Signal_curr[exc_pos];
if (numTS == 0)
cout << "Engine::ApplyCurrentExcite(): n=" << n << " dir=" << Op->Exc->Curr_dir[n] << " coords=" << Op->Exc->Curr_index[0][n] << "," << Op->Exc->Curr_index[1][n] << "," << Op->Exc->Curr_index[2][n] << " amp=" << Op->Exc->Curr_amp[n] << endl;
}
// write the first excitation into the file "ht"
if (numTS < Op->Exc->Length)
file_ht << (numTS+0.5) * Op->GetTimestep() << "\t" << Op->Exc->Curr_amp[0]*Op->Exc->Signal_curr[numTS] << "\n"; // do not use std::endl here, because it will do an implicit flush
else
file_ht << (numTS+0.5) * Op->GetTimestep() << "\t0" << "\n"; // do not use std::endl here, because it will do an implicit flush
}
bool Engine::IterateTS(unsigned int iterTS)

2
FDTD/engine.h
View File

@ -77,7 +77,7 @@ protected:
virtual void InitExtensions();
vector<Engine_Extension*> m_Eng_exts;
ofstream file_et; //excite signal dump file
ofstream file_et, file_ht; //excite signal dump file
friend class NS_Engine_Multithread::thread; // evil hack to access numTS from multithreading context
};