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

new: node interpolated dump in processfields

Browse Source
This commit is contained in:
Thorsten Liebig 2010-05-11 20:38:58 +02:00
parent 279b348724
commit 5e5a278ac7
4 changed files with 156 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

31
FDTD/processfields.cpp
View File

@ -138,6 +138,35 @@ void ProcessFields::DefineStartStopCoord(double* dstart, double* dstop)
discLines[n][i] = lines.at(i);
}
}
else if (m_DumpMode==NODE_INTERPOLATE)
{
if (Op->SnapToMesh(dstart,start)==false) cerr << "ProcessFields::DefineStartStopCoord: Warning: Snapping problem, check start value!!" << endl;
if (Op->SnapToMesh(dstop,stop)==false) cerr << "ProcessFields::DefineStartStopCoord: Warning: Snapping problem, check stop value!!" << endl;
//create mesh
for (int n=0;n<3;++n)
{
if (start[n]>stop[n])
{
unsigned int help = start[n];
start[n]=stop[n];
stop[n]=help;
}
if (stop[n]==Op->GetNumberOfLines(n)-1)
--stop[n];
// cerr << "start " << start[n] << "stop " << stop[n];
lines.clear();
for (unsigned int i=start[n];i<=stop[n];i+=subSample[n])
{
lines.push_back(Op->GetDiscLine(n,i));//0.5*(Op->discLines[n][i+1] + Op->discLines[n][i]));
}
numLines[n] = lines.size();
delete[] discLines[n];
discLines[n] = new double[numLines[n]];
for (unsigned int i=0;i<numLines[n];++i)
discLines[n][i] = lines.at(i);
}
}
else if (m_DumpMode==CELL_INTERPOLATE)
{
if (Op->SnapToMesh(dstart,start,true)==false) cerr << "ProcessFields::DefineStartStopCoord: Warning: Snapping problem, check start value!!" << endl;
@ -158,7 +187,7 @@ void ProcessFields::DefineStartStopCoord(double* dstart, double* dstop)
for (unsigned int i=start[n];i<stop[n];i+=subSample[n])
{
lines.push_back(Op->GetDiscLine(n,i,true));//0.5*(Op->discLines[n][i+1] + Op->discLines[n][i]));
}
}
numDLines[n] = lines.size();
delete[] discDLines[n];
discDLines[n] = new double[numDLines[n]];

4
FDTD/processfields.h
View File

@ -46,7 +46,9 @@ public:
//! Define the Dump-Mode
void SetDumpMode(DumpMode mode) {m_DumpMode=mode;}
//! This methode will dump all fields in the center of a main cell (dual-node) using 4 E-field and 2 H-fields per direction. (default)
//! This methode will dump all fields on a main cell node using 2 E-field and 4 H-fields per direction.
void SetDumpMode2Node() {m_DumpMode=NODE_INTERPOLATE;}
//! This methode will dump all fields in the center of a main cell (dual-node) using 4 E-field and 2 H-fields per direction.
void SetDumpMode2Cell() {m_DumpMode=CELL_INTERPOLATE;}
//! Set dump type: 0 for E-fields, 1 for H-fields, 2 for D-fields, 3 for B-fields, 4 for J-fields, etc...

122
FDTD/processfields_td.cpp
View File

@ -29,6 +29,126 @@ ProcessFieldsTD::~ProcessFieldsTD()
{
}
void ProcessFieldsTD::DumpNodeInterpol(string filename)
{
if (m_DumpType==H_FIELD_DUMP)
{
//create array
FDTD_FLOAT**** H_T = Create_N_3DArray(numLines);
unsigned int pos[3] = {start[0],start[1],start[2]};
unsigned int OpPos[3];
double delta;
// cerr << "processing e-fields... " << endl;
for (pos[0]=0;pos[0]<numLines[0];++pos[0])
{
OpPos[0]=start[0]+pos[0]*subSample[0];
for (pos[1]=0;pos[1]<numLines[1];++pos[1])
{
OpPos[1]=start[1]+pos[1]*subSample[1];
for (pos[2]=0;pos[2]<numLines[2];++pos[2])
{
OpPos[2]=start[2]+pos[2]*subSample[2];
//in x
delta = Op->GetMeshDelta(0,OpPos,true);
if (delta)
{
H_T[0][pos[0]][pos[1]][pos[2]] = Eng->GetCurr(0,OpPos[0],OpPos[1],OpPos[2]) + Eng->GetCurr(0,OpPos[0],OpPos[1]+1,OpPos[2]) + Eng->GetCurr(0,OpPos[0],OpPos[1],OpPos[2]+1) + Eng->GetCurr(0,OpPos[0],OpPos[1]+1,OpPos[2]+1);
H_T[0][pos[0]][pos[1]][pos[2]] /= (4*delta);
}
//in y
delta = Op->GetMeshDelta(1,OpPos,true);
if (delta)
{
H_T[1][pos[0]][pos[1]][pos[2]] = Eng->GetCurr(1,OpPos[0],OpPos[1],OpPos[2]) + Eng->GetCurr(1,OpPos[0]+1,OpPos[1],OpPos[2]) + Eng->GetCurr(1,OpPos[0],OpPos[1],OpPos[2]+1) + Eng->GetCurr(1,OpPos[0]+1,OpPos[1],OpPos[2]+1);
H_T[1][pos[0]][pos[1]][pos[2]] /= (4*delta);//*Op->gridDelta);
}
//in z
delta = Op->GetMeshDelta(2,OpPos,true);
if (delta)
{
H_T[2][pos[0]][pos[1]][pos[2]] = Eng->GetCurr(2,OpPos[0],OpPos[1],OpPos[2]) + Eng->GetCurr(2,OpPos[0],OpPos[1]+1,OpPos[2]) + Eng->GetCurr(2,OpPos[0]+1,OpPos[1],OpPos[2]) + Eng->GetCurr(2,OpPos[0]+1,OpPos[1]+1,OpPos[2]);
H_T[2][pos[0]][pos[1]][pos[2]] /= (4*delta);
}
}
}
}
if (m_fileType==VTK_FILETYPE)
{
ofstream file(filename.c_str());
if (file.is_open()==false) { cerr << "ProcessFieldsTD::Process: can't open file '" << filename << "' for writing... abort! " << endl;};
DumpVectorArray2VTK(file,string("H-Field"),H_T,discLines,numLines,m_precision);
file.close();
}
else if (m_fileType==HDF5_FILETYPE)
{
stringstream ss;
ss << std::setw( pad_length ) << std::setfill( '0' ) << Eng->GetNumberOfTimesteps();
DumpVectorArray2HDF5(filename.c_str(),string( ss.str() ),H_T,numLines,Eng->GetNumberOfTimesteps()*Op->GetTimestep());
}
else
cerr << "ProcessFieldsTD::DumpNodeInterpol: unknown File-Type" << endl;
Delete_N_3DArray(H_T,numLines);
H_T = NULL;
}
if (m_DumpType==E_FIELD_DUMP)
{
//create array
FDTD_FLOAT**** E_T = Create_N_3DArray(numLines);
unsigned int pos[3] = {start[0],start[1],start[2]};
unsigned int OpPos[3];
unsigned int OpPosUp[3];
double delta, deltaUp, deltaRel;
// cerr << "processing h-fields... " << endl;
for (pos[0]=0;pos[0]<numLines[0];++pos[0])
{
OpPos[0]=start[0]+pos[0]*subSample[0];
OpPosUp[0]=start[0]+pos[0]*subSample[0];
for (pos[1]=0;pos[1]<numLines[1];++pos[1])
{
OpPos[1]=start[1]+pos[1]*subSample[1];
OpPosUp[1]=start[1]+pos[1]*subSample[1];
for (pos[2]=0;pos[2]<numLines[2];++pos[2])
{
OpPos[2]=start[2]+pos[2]*subSample[2];
OpPosUp[2]=start[2]+pos[2]*subSample[2];
for (int n=0;n<3;++n)
{
delta = Op->GetMeshDelta(n,OpPos);
++OpPosUp[n];
deltaUp = Op->GetMeshDelta(n,OpPos);
deltaRel = delta / (delta+deltaUp);
if (delta*deltaUp)
{
E_T[n][pos[0]][pos[1]][pos[2]] = Eng->GetVolt(n,OpPos)*(1-deltaRel)/delta + Eng->GetVolt(n,OpPosUp)/deltaUp*deltaRel;
}
--OpPosUp[n];
}
}
}
}
if (m_fileType==VTK_FILETYPE)
{
ofstream file(filename.c_str());
if (file.is_open()==false) { cerr << "ProcessFieldsTD::Process: can't open file '" << filename << "' for writing... abort! " << endl;};
DumpVectorArray2VTK(file,string("E-Field"),E_T,discLines,numLines,m_precision);
file.close();
}
else if (m_fileType==HDF5_FILETYPE)
{
stringstream ss;
ss << std::setw( pad_length ) << std::setfill( '0' ) << Eng->GetNumberOfTimesteps();
DumpVectorArray2HDF5(filename.c_str(),string( ss.str() ),E_T,numLines,(0.5+Eng->GetNumberOfTimesteps())*Op->GetTimestep());
}
else
cerr << "ProcessFieldsTD::DumpCellInterpol: unknown File-Type" << endl;
Delete_N_3DArray(E_T,numLines);
E_T = NULL;
}
}
void ProcessFieldsTD::DumpCellInterpol(string filename)
{
if (m_DumpType==E_FIELD_DUMP)
@ -259,6 +379,8 @@ int ProcessFieldsTD::Process()
ss << ".vtk";
if (m_DumpMode==NO_INTERPOLATION)
DumpNoInterpol(ss.str());
if (m_DumpMode==NODE_INTERPOLATE)
DumpNodeInterpol(ss.str());
if (m_DumpMode==CELL_INTERPOLATE)
DumpCellInterpol(ss.str());
}

1
FDTD/processfields_td.h
View File

@ -35,6 +35,7 @@ protected:
int pad_length;
void DumpNoInterpol(string filename);
void DumpNodeInterpol(string filename);
void DumpCellInterpol(string filename);
};