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

nf2ff: restore mesh lines after mirroring 

This fixes a bug in nf2ff. When using mirror planes, fields and mesh
lines are mirrored in place. Mesh lines must be restored after the
calculation of a single frequency point. Otherwise, the far field at
every even frequency point is incorrect.
pull/32/head
Detlef Conradin 2018-03-03 13:36:15 +01:00
parent 65ca6bfc44
commit 2799f19ae3
1 changed files with 12 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
nf2ff/nf2ff_calc.cpp
View File

@ -313,6 +313,10 @@ bool nf2ff_calc::AddPlane(float **lines, unsigned int* numLines, complex<float>*
if ((m_MirrorType[n]!=MIRROR_OFF) && (m_MirrorType[nP]==MIRROR_OFF) && (m_MirrorType[nPP]==MIRROR_OFF)) if ((m_MirrorType[n]!=MIRROR_OFF) && (m_MirrorType[nP]==MIRROR_OFF) && (m_MirrorType[nPP]==MIRROR_OFF))
{ {
this->AddMirrorPlane(n, lines, numLines, E_field, H_field, MeshType); this->AddMirrorPlane(n, lines, numLines, E_field, H_field, MeshType);
for (unsigned int i=0;i<numLines[n];++i)
lines[n][i] = 2.0*m_MirrorPos[n] - lines[n][i];
break; break;
} }
//check if two planes are on //check if two planes are on
@ -321,6 +325,10 @@ bool nf2ff_calc::AddPlane(float **lines, unsigned int* numLines, complex<float>*
this->AddMirrorPlane(nP, lines, numLines, E_field, H_field, MeshType); this->AddMirrorPlane(nP, lines, numLines, E_field, H_field, MeshType);
this->AddMirrorPlane(nPP, lines, numLines, E_field, H_field, MeshType); this->AddMirrorPlane(nPP, lines, numLines, E_field, H_field, MeshType);
this->AddMirrorPlane(nP, lines, numLines, E_field, H_field, MeshType); this->AddMirrorPlane(nP, lines, numLines, E_field, H_field, MeshType);
for (unsigned int i=0;i<numLines[nPP];++i)
lines[nPP][i] = 2.0*m_MirrorPos[nPP] - lines[nPP][i];
break; break;
} }
} }
@ -334,6 +342,10 @@ bool nf2ff_calc::AddPlane(float **lines, unsigned int* numLines, complex<float>*
this->AddMirrorPlane(0, lines, numLines, E_field, H_field, MeshType); this->AddMirrorPlane(0, lines, numLines, E_field, H_field, MeshType);
this->AddMirrorPlane(1, lines, numLines, E_field, H_field, MeshType); this->AddMirrorPlane(1, lines, numLines, E_field, H_field, MeshType);
this->AddMirrorPlane(0, lines, numLines, E_field, H_field, MeshType); this->AddMirrorPlane(0, lines, numLines, E_field, H_field, MeshType);
for (unsigned int i=0;i<numLines[2];++i)
lines[2][i] = 2.0*m_MirrorPos[2] - lines[2][i];
} }
//cleanup E- & H-Fields //cleanup E- & H-Fields