Atom style variable의 연산

LIGGGHTS에서 variable은 equal style과 atom style로 구분된다.

예를들어 atom의 radius를 compute한 결과는 atom style (열벡터)로 얻어진다.

compute                    cal_radius all property/atom update_on_run_end yes radius

variable                   ak2 atom c_cal_radius
variable                   ak3 atom (c_cal_radius)^2

이렇게 계산된 atom variable은 custom dump style로 출력 가능하다.

dump                       dump_imsi2         all custom 100 ak1-*.txt v_ak2 v_ak3

출력 결과를 확인해 보면,

ITEM: ATOMS v_ak2 v_ak3 
0.00025 6.25e-08 
0.0005 2.5e-07 
0.0003 9e-08 
0.00055 3.025e-07 
0.00035 1.225e-07 
0.0006 3.6e-07 
0.0004 1.6e-07 
0.00065 4.225e-07 
0.00045 2.025e-07 
0.0007 4.9e-07

atom variable의 제곱은 각각의 값을 제곱하여 다시 atom style로 저장됨을 알 수 있다.

 

어디에 응용가능한가? atom의 부피를 계산하여, compute reduce를 사용하면 총 부피를 계산할 수 있다.