PyMieScatt > if ndp[-1]>np.max(ndp)/100 or ndp[0]>np.max(ndp)/100: > not compact?

動作環境
GeForce GTX 1070 (8GB)
ASRock Z170M Pro4S [Intel Z170chipset]
Ubuntu 16.04.4 LTS desktop amd64
TensorFlow v1.7.0
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

ndpチェック

https://github.com/bsumlin/PyMieScatt/blob/master/PyMieScatt/Mie.py#L470

    ndp = numberOfParticles*ithPart(1,dp,geoMean,geoStdDev)
  if ndp[-1]>np.max(ndp)/100 or ndp[0]>np.max(ndp)/100:
    warnings.warn("Warning: distribution may not be compact on the specified interval. Consider using a higher upper bound.")

ndpというのは数密度分布だろうか。
http://dust.ess.uci.edu/facts/psd/psd.pdf
のp2 Eq(1) The total number concentration of particles N0...の$n_n(r)$
(ref. statistics > mean size and variance > 緯度経度ごとの値から緯度のmeanとvarianceを求める)

ndpの最小インデックスと最大インデックスに位置する要素値がndpの最大値/100を超えていると警告を出している。

100というマジックナンバーの意味が不明である。

数密度分布のオーダーが2(=log10(100))を超えた場合に、compactでないという定義にしているようだ。

数密度分布のオーダーが2(=log10(100))未満の時に、compactでないという定義にしているようだ。

オーダーを広げる方が「compact」という定義がよくわからない。

Sign up for free and join this conversation.
Sign Up
If you already have a Qiita account log in.