PyMieScatt > SU = (SR + SL) / 2 > unpolarized = (perpendicular + parallel) / 2

動作環境

GeForce GTX 1070 (8GB)
ASRock Z170M Pro4S [Intel Z170chipset]
Ubuntu 16.04 LTS desktop amd64
TensorFlow v1.2.1
cuDNN v5.1 for Linux
CUDA v8.0
Python 3.5.2
IPython 6.0.0 -- An enhanced Interactive Python.
gcc (Ubuntu 5.4.0-6ubuntu1~16.04.4) 5.4.0 20160609
GNU bash, version 4.3.48(1)-release (x86_64-pc-linux-gnu)
scipy v0.19.1
geopandas v0.3.0
MATLAB R2017b (Home Edition)
ADDA v.1.3b6
gnustep-gui-runtime v0.24.0-3.1

Reference

• Bohren C.F. and D. R. Huffman, "Absorption and Scattering of Light by Small Particles", John Wiley, New York, NY (1983)
• PDF: MATLAB Functions for Mie Scattering and Absorption, Research Report No. 2002-08, June 2002 by Christian Mätzler

SR, SL, SU

http://pymiescatt.readthedocs.io/en/latest/forward.html#ScatteringFunction
https://github.com/bsumlin/PyMieScatt/blob/master/PyMieScatt/Mie.py#L306

    SL[j] = (np.sum(np.conjugate(S1)*S1)).real
    SR[j] = (np.sum(np.conjugate(S2)*S2)).real
    SU[j] = (SR[j]+SL[j])/2

行列?

行列は@mpywさんとのやり取りでRuby実装を教えていただいたことがある。
https://qiita.com/mpyw/items/321c5485236615ca4664#comment-adb7fba3e87b85b4c18e
ideoneでは以下
https://ideone.com/XS2Mzc

L(左三角形）とR(右三角形）という定義はあっただろうか?

... then computes parallel, perpendicular, and unpolarized intensities by

SR(\theta) = | S_1 |^2

SL(\theta) = | S_2 |^2

SU(\theta) = \frac{1}{2}(SR + SL)

L(Parallel)とR(perpendicular)を足したのはU(unpolarized: 無偏光)になる。

L(左円偏光)とR(右円偏光)を足した場合は直線偏光になるだっただろうか。忘れた。

left-polarized (parallel) light

Returns
SL :numpy.ndarray
An array of the scattered intensity of left-polarized (parallel) light. Same size as the theta array.

SLは左偏光(parallel)という記載。

left-polarizedとparallel component of the electric field vectorは別の定義と認識していたが。。。

left-polarizedは電場の振動方向が時刻により変化する。
parallel componentはある軸を基準に、電場の振動方向はその軸にparallelのまま。
?