float sdf = volumesample(0,0,@P); //0번인풋의 0번째 볼륨(prim)에 대해서 @P를 불러옵니다. sdf의 경우 표면이 0, 표면밖이 0보다 큰값(양수), 표면 안쪽이 0보다 작은 값(음수)이 되어, 표면과 멀 수록 표면(0)과의 값의 차가 더욱 커집니다.
vector grad = volumegradient(0,0,@P); //0번인풋의 0번째 볼륨(prim)에 대하여 @P를 그라디언트(노말라이징)합니다. 이 경우 표면에서 바깥과 안쪽에서 각각 표면으로 향하는 최단거리방향의 벡터가 모두 동일하게 생성됩니다.
vector inforce = -sdf*grad*ch('to_surface_mult'); //float의 sdf의 역수값을, 노말라이징된 vector인 grad와 곱함으로써, 표면에서 부터의 거리가 멀 수록 보다 표면에서 멀어지는 큰 벡터가 되도록 하여, inforce라는 벡터로 설정합니다. (ch값으로 곱하는 값을 추가 생성)
@vel *= ch('curl_mult'); //기존에 만들어진 @vel (노말방향과 curlnoise(@P+@Time)을 crossproduct 된 벡터를 @vel로 앞서 사용했습니다) 에, 추가의 곱하는 채널값을 만들어둡니다
@vel = lerp(@vel, inforce, abs(sdf)); //@vel과 상단의 inforce vector 두 벡터값을 lerp하여 원하는 비율로 믹싱합니다. 믹싱의 비율은 상단의 sdf float의 절대값을 따릅니다 (lerp의 경우 최소최대 0-1까지만 영향받습니다)
//즉, 표면에서의 거리가 멀어질수록 1이되거나 1을 넘게되어 @vel은 없어지고 inforce의 비중이 커진채로 @vel에 새롭게 갱신됩니다. (=표면에서 멀 수록 표면밖으로 향하는 벡터가 됩니다)
//반대로 표면에서의 거리가 가까울수록 0과 가까워 지기때문에, 기존에 만들어둔 @vel의 비중이 커지게 됩니다. 다만 표면을 제외한 나머지 표면에 근접한 부분은 inforce의 비중이 상대적으로 상당히 작을 뿐 영향을 주긴 합니다. (=표면에 가까울수록 @P의 curlnoise와 표면의 노말이 crossproduct한 벡터가 됩니다)
@density = clamp(fit(abs(sdf),0,.05,1,0),0,.05); //밀도 @density는 표면에 가까울수록 1이되고 멀 수록(최대 0.05까지를) 0이 되도록 새롭게 매핑하지만, 이 값을 0과 0.05로 다시 제한(clamp)합니다
@density *= fit(noise(@P*ch('density_noise_freq')),0.2,1,0,ch('density')); //그런 뒤 밀도 @density 에는 @P의 noise값에 대하여 새롭게 매핑한 값을 다시 곱합니다.
//정리하자면 @vel 필드는 표면과의 거리에 따른 @P 기준의 curlnoise의 영향력이 각기 다르도록 조정하였고, 그 전에 생성될 @density 필드는 미리 노이즈를 준 채로 시뮬에 적용되도록 합니다.
//이하 단계는 pyro solver의 영향으로 이어지게 됩니다.
원문 링크
참고자료 링크
