HomerReid / scuff-em

A comprehensive and full-featured computational physics suite for boundary-element analysis of electromagnetic scattering, fluctuation-induced phenomena (Casimir forces and radiative heat transfer), nanophotonics, RF device engineering, electrostatics, and more. Includes a core library with C++ and python APIs as well as many command-line applications.
http://www.homerreid.com/scuff-em
GNU General Public License v2.0
126 stars 50 forks source link

Nan output for E-field of an infinitely PEC beam #243

Open SaihangQie opened 1 year ago

SaihangQie commented 1 year ago

I am trying to compute the scattered field of an infinitely PEC beam. I use the command-line tool scuff-scatter to complete the calculation. The gmsh geometry file for unit-cell geometry is from /scuff-em/examples/SiliconBeams. Lamnbda = 1 micron. The incident wave is planewave. I set pwDirection 0 0 -1, pwPolarization 1 0 0, EPFile observePoints.

However, the result of scttered filed is all nan. I am not sure why this is the case. I will appreciate any suggestions on how to resolve this error. Thank you.

//Source https://github.com/SaihangQie/Issue //Command scuff-scatter --geometry RoundedBeamUnitCell.scuffgeo --Lambda 1 --pwDirection 0 0 -1 --pwPolarization 1 0 0 --EPFile ObservePoints.txt