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

bug fix in interpolated h-field dump

pull/1/head
Thorsten Liebig 2010-05-07 22:58:00 +02:00
parent 065c05eeae
commit 5b3ccd3f53
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
FDTD/processfields_td.cpp
View File

@ -112,7 +112,7 @@ void ProcessFieldsTD::DumpCellInterpol(string filename)
//in x //in x
// if (OpPos[0]==0) delta = Op->discLines[0][OpPos[0]+1] - Op->discLines[0][OpPos[0]]; // if (OpPos[0]==0) delta = Op->discLines[0][OpPos[0]+1] - Op->discLines[0][OpPos[0]];
// else delta = 0.5* (Op->discLines[0][OpPos[0]+1] - Op->discLines[0][OpPos[0]-1]); // else delta = 0.5* (Op->discLines[0][OpPos[0]+1] - Op->discLines[0][OpPos[0]-1]);
delta = Op->GetDiscLine(0,OpPos[0],true); delta = Op->GetMeshDelta(0,OpPos,true);
if (delta) 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]+1,OpPos[1],OpPos[2]); H_T[0][pos[0]][pos[1]][pos[2]] = Eng->GetCurr(0,OpPos[0],OpPos[1],OpPos[2]) + Eng->GetCurr(0,OpPos[0]+1,OpPos[1],OpPos[2]);
@ -120,7 +120,7 @@ void ProcessFieldsTD::DumpCellInterpol(string filename)
} }
//in y //in y
// delta = Op->discLines[1][OpPos[1]+1] - Op->discLines[1][OpPos[1]]; // delta = Op->discLines[1][OpPos[1]+1] - Op->discLines[1][OpPos[1]];
delta = Op->GetDiscLine(1,OpPos[1],true); delta = Op->GetMeshDelta(1,OpPos,true);
if (delta) 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],OpPos[1]+1,OpPos[2]); H_T[1][pos[0]][pos[1]][pos[2]] = Eng->GetCurr(1,OpPos[0],OpPos[1],OpPos[2]) + Eng->GetCurr(1,OpPos[0],OpPos[1]+1,OpPos[2]);
@ -128,7 +128,7 @@ void ProcessFieldsTD::DumpCellInterpol(string filename)
} }
//in z //in z
// delta = Op->discLines[2][OpPos[2]+1] - Op->discLines[2][OpPos[2]]; // delta = Op->discLines[2][OpPos[2]+1] - Op->discLines[2][OpPos[2]];
delta = Op->GetDiscLine(2,OpPos[2],true); delta = Op->GetMeshDelta(2,OpPos,true);
if (delta) 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],OpPos[2]+1); 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],OpPos[2]+1);