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
This article is related to ADDA (light scattering simulator based on the discrete dipole approximation).
関連: Superposition T-matrix codes > bisphere.f > 2つの球の光散乱計算コード (analytical orientation averaging)
ADDA側でbisphereの計算を実行しようとしているが、実行引数について不明な点が出ている。
@ manual.pdfのpage 15
bisphere - two identical spheres with outer diameter d and center-to-center distance R cc
(along the z-axis). It describes both separate and sintered spheres.
実行時ヘルプを見てみると
$ ./adda -shape bisphere
ERROR: Illegal number of arguments (0) to '-shape bisphere' option (1 expected)
Usage: -shape bisphere <R_cc/d>
Type 'adda -h shape bisphere' for details
<R_cc/d>に指定する値により、２つの球の位置情報が変化する。
接触状態にするには0.5だろうか。
位置情報の可視化
Jupyter | Matplotlib > 凹凸の激しい形状の可視化 > v0.4, v0.5: ADDA InitField-Yファイル読込み > Sphere | Prism (N=3) | Chevyshev particles
で実装したコードを用いて接触状態かどうかの確認をする。
InitField-Yの生成
ADDAの実行時引数「-store_int_field」を指定することで、IntField-Yを作る。
$ ./adda -shape bisphere 0.5 -store_int_field
...
$ ls run22099_bisphere_g16_m1.5/
CrossSec-Y IntField-Y log mueller
上記のリンク先のコード: showChebyshev_170910.ipynbを実行する。
描画までに30秒程度かかる。
結果
$./adda -shape bisphere 0.5 -store_int_field
$./adda -shape bisphere 1.0 -store_int_field
2つの球が接触状態になるbisphereパラメータ
上記の画像より、以下のパラメータにすることで２つの球が接触状態になることが分かった。
$./adda -shape bisphere 1.0
Possible correction to the manual.pdf page 15
bisphere - two identical spheres with outer diameter d and center-to-center distance Rcc
(along the z-axis). It describes both separate and sintered spheres.
shall better be:
"bisphere - two identical spheres with composing sphere diameter d and center-to-center distance Rcc
(along the z-axis). It describes both separate and sintered spheres."