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

Custom python EpsMuGenerators seem to work 

Former-commit-id: b3a38dde88b019182fcee4e62028b93fe0bf7aa5
This commit is contained in:
Marek Nečada 2019-08-11 19:11:43 +03:00
parent 4572322553
commit d00e46b908
1 changed files with 4 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
qpms/cymaterials.pyx
View File

@ -3,7 +3,7 @@
import numpy as np import numpy as np
import cmath import cmath
from .qpms_cdefs cimport qpms_permittivity_interpolator_from_yml, qpms_permittivity_interpolator_free, qpms_permittivity_interpolator_omega_min, qpms_permittivity_interpolator_omega_max, gsl_interp_type, qpms_permittivity_interpolator_t, gsl_interp_cspline, qpms_permittivity_interpolator_eps_at_omega, qpms_epsmu_const_g, qpms_permittivity_interpolator_epsmu_g, qpms_epsmu_const_g, qpms_lorentzdrude_epsmu_g, qpms_ldparams_triple_t, qpms_lorentzdrude_eps from .qpms_cdefs cimport qpms_permittivity_interpolator_from_yml, qpms_permittivity_interpolator_free, qpms_permittivity_interpolator_omega_min, qpms_permittivity_interpolator_omega_max, gsl_interp_type, qpms_permittivity_interpolator_t, gsl_interp_cspline, qpms_permittivity_interpolator_eps_at_omega, qpms_epsmu_const_g, qpms_permittivity_interpolator_epsmu_g, qpms_epsmu_const_g, qpms_lorentzdrude_epsmu_g, qpms_ldparams_triple_t, qpms_lorentzdrude_eps, cdouble
from .cycommon cimport make_c_string from .cycommon cimport make_c_string
cimport cython cimport cython
import enum import enum
@ -76,7 +76,7 @@ lorentz_drude = {
} }
cdef qpms_epsmu_t python_epsmu_generator(cdouble omega, const void *params): cdef qpms_epsmu_t python_epsmu_generator(cdouble omega, const void *params):
object fun = <object> params cdef object fun = <object> params
cdef qpms_epsmu_t em cdef qpms_epsmu_t em
em.eps, em.mu = fun(omega) em.eps, em.mu = fun(omega)
return em return em
@ -97,7 +97,7 @@ cdef class EpsMuGenerator:
self.g.params = (<MaterialInterpolator?>self.holder).rawpointer() self.g.params = (<MaterialInterpolator?>self.holder).rawpointer()
elif isinstance(what, EpsMuGenerator): # Copy constructor elif isinstance(what, EpsMuGenerator): # Copy constructor
self.holder = what.holder self.holder = what.holder
self.g = what.g self.g = (<EpsMuGenerator?>what).g
elif callable(what): elif callable(what):
warnings.warn("Custom python (eps,mu) generators are an experimental feature") warnings.warn("Custom python (eps,mu) generators are an experimental feature")
self.holder = what self.holder = what
@ -126,7 +126,6 @@ cdef class EpsMuGenerator:
if(omega < i[0] or omega > i[1]): if(omega < i[0] or omega > i[1]):
raise ValueError("Input frequency %g is outside the interpolator domain (%g, %g)." raise ValueError("Input frequency %g is outside the interpolator domain (%g, %g)."
% (omega, i[0], i[1])) % (omega, i[0], i[1]))
with nogil:
em = self.g.function(omega, self.g.params) em = self.g.function(omega, self.g.params)
return EpsMu(em.eps, em.mu) return EpsMu(em.eps, em.mu)