Closed 10110111 closed 2 years ago
Thank you for reporting! We will try to manage this by the next release. Unfortunately, Mie theory for near-field is not developed to the same state of the art as far-field computations...
It seems that changes from #42 also fixed this.
$ ./fieldnlay-dp -l 1 113 1.3346 0 -p 242.1 242.1 1 -187.1 -187.1 1 -0.001 0.001 9
Near-field early convergence at nmax = 148
Near-field early convergence at nmax = 148
Near-field early convergence at nmax = 148
Near-field early convergence at nmax = 148
Near-field early convergence at nmax = 148
Near-field early convergence at nmax = 148
Near-field early convergence at nmax = 148
Near-field early convergence at nmax = 148
Near-field early convergence at nmax = 148
Number of multipoles used in Mie series nmax=183
X, Y, Z, Ex.r, Ex.i, Ey.r, Ey.i, Ez.r, Ez.i, Hx.r, Hx.i, Hy.r, Hy.i, Hz.r, Hz.i
242.1000000, -187.1000000, -0.0010000, +9.90348e-01, -1.83520e-02, -9.49330e-03, -2.51073e-02, -6.98874e-03, -1.24879e-03, +2.31173e-05, +2.27662e-05, +2.66144e-03, -1.54915e-05, -3.36443e-05, -7.72585e-05
242.1000000, -187.1000000, -0.0007500, +9.90350e-01, -1.81027e-02, -9.49134e-03, -2.51085e-02, -6.98838e-03, -1.24841e-03, +2.31152e-05, +2.27663e-05, +2.66144e-03, -1.48295e-05, -3.36381e-05, -7.72618e-05
242.1000000, -187.1000000, -0.0005000, +9.90352e-01, -1.78533e-02, -9.48937e-03, -2.51096e-02, -6.98802e-03, -1.24802e-03, +2.31132e-05, +2.27664e-05, +2.66144e-03, -1.41674e-05, -3.36320e-05, -7.72650e-05
242.1000000, -187.1000000, -0.0002500, +9.90353e-01, -1.76039e-02, -9.48740e-03, -2.51108e-02, -6.98766e-03, -1.24764e-03, +2.31112e-05, +2.27665e-05, +2.66144e-03, -1.35053e-05, -3.36258e-05, -7.72683e-05
242.1000000, -187.1000000, 0.0000000, +9.90355e-01, -1.73546e-02, -9.48544e-03, -2.51120e-02, -6.98730e-03, -1.24725e-03, +2.31091e-05, +2.27666e-05, +2.66144e-03, -1.28432e-05, -3.36197e-05, -7.72716e-05
242.1000000, -187.1000000, 0.0002500, +9.90356e-01, -1.71052e-02, -9.48347e-03, -2.51132e-02, -6.98694e-03, -1.24687e-03, +2.31071e-05, +2.27667e-05, +2.66144e-03, -1.21811e-05, -3.36135e-05, -7.72749e-05
242.1000000, -187.1000000, 0.0005000, +9.90357e-01, -1.68558e-02, -9.48150e-03, -2.51143e-02, -6.98658e-03, -1.24649e-03, +2.31051e-05, +2.27669e-05, +2.66144e-03, -1.15190e-05, -3.36074e-05, -7.72781e-05
242.1000000, -187.1000000, 0.0007500, +9.90359e-01, -1.66065e-02, -9.47953e-03, -2.51155e-02, -6.98621e-03, -1.24610e-03, +2.31031e-05, +2.27670e-05, +2.66144e-03, -1.08569e-05, -3.36012e-05, -7.72814e-05
242.1000000, -187.1000000, 0.0010000, +9.90360e-01, -1.63571e-02, -9.47757e-03, -2.51167e-02, -6.98585e-03, -1.24572e-03, +2.31010e-05, +2.27671e-05, +2.66144e-03, -1.01948e-05, -3.35951e-05, -7.72847e-05
Add a related test here https://github.com/ovidiopr/scattnlay/blob/quasar/tests/test_near_field.cc#L7-L42 , it passes OK now in quasar branch
Consider the following command:
EM field in the range of z∈[-0.001,0.001] should be very smooth, but the calculated values exhibit a jump at z=0:
Here's a plot of these values (51 points instead of the 9 above) to make the jump more obvious: