garrekstemo / transfermatrix.jl Goto Github PK

Right now the transfer matrix method is used. While the current research in transfer matrices tries to mitigate its shortfalls (singularities and numerical instabilities), scattering matrices avoid these from the start. It would probably be more efficient and definitely more broadly useful to the photonics community to use scattering matrices instead of transfer matrices. In any case, once the transfer matrix is calculated, the scattering matrix elements are still calculated in the end to get transmission and reflection. (see this discussion)

Proposal: keep the current transfer matrix code, but write a scattering matrix implementation and have both available.

Calculate the refractive index from a dispersion formula using RefractiveIndex.jl when loading data.

Looking through your code, it looks like in your 4x4 transfer matrix formulation has a opposite ordering to Yeh' s book p246 - 247. (Yeh Pochi. 2005. Optical Waves in Layered Media. Hoboken NJ: Wiley.) In Yeh's book, he ordered the s-pol in the upper half of the vector and p-pol in the lower, while your implementation seems to reverse such ordering and hence, your s, p subscripts are swapped compared with Eq. (9.7-23) - (9.7-30) in Yeh's book. I guess it may be helpful to include that in the documentation to avoid confusion? Maybe something like Eq. (9.7-22) in Yeh's book but with your ordering scheme would suffice as a clarification? Thank you.

The database at refractiveindex.info has been updated, changing the table structure when downloading refractive index data. This has caused the read_refractive() function to break.

Either warn the user to reformat the file, fix this function, or use RefractiveIndex.jl.

@JuliaRegistrator register

StaticArrays.jl would be more efficient than native Julia arrays.

Make this the default instead of requiring the user to type 1.

@JuliaRegistrator register

• Use interpolate
• Error if out of range

Argument information is missing in the docs for the Layer and Structure types.

@JuliaRegistrator register

Use PrecompileTools.jl

For easy transport of complicated Structures

This issue is used to trigger TagBot; feel free to unsubscribe.

If you haven't already, you should update your TagBot.yml to include issue comment triggers.
Please see this post on Discourse for instructions and more details.

If you'd like for me to do this for you, comment TagBot fix on this issue.
I'll open a PR within a few hours, please be patient!