QPMS
/
qpms
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

Fix sign in Drude-Lorentz model. 

In our convention, where outgoing wave correspond to
Hankel functions of FIRST kind, the time-harmonic waves carry
(contrary to the more usual electrotechnical convention)
the factor exp(-iωt) and lossy materials are represented
by Im ε > 0, i.e. Im n > 0.


Former-commit-id: 19108c43a31899c2f9196c1a0fc44c5ecae69b5a
This commit is contained in:
Marek Nečada 2019-11-08 09:28:13 +02:00
parent 28f4e7f3d2
commit cfda120a20
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
qpms/materials.c
View File

@ -158,7 +158,7 @@ complex double qpms_lorentzdrude_eps(complex double omega, const qpms_ldparams_t
complex double eps = 0;
for(size_t j = 0; j < p->n; ++j) {
const qpms_ldparams_triple_t d = p->data[j];
eps += d.f * SQ(p->omega_p) / (SQ(d.omega) - SQ(omega) + I*omega*d.gamma );
eps += d.f * SQ(p->omega_p) / (SQ(d.omega) - SQ(omega) - I*omega*d.gamma );
}
return eps;
}