18 lines
445 B
Text
18 lines
445 B
Text
# collect samples for pdf estimation from ic1 case
|
|
# sampling region (xs, xe) hard-coded for ic1 case
|
|
|
|
[u, v, w, y]
|
|
|
|
# field_name (arg_name=arg (, pair)* ) = expression
|
|
|
|
c (export=true, xs=176, xe=400) = 1.0 - y
|
|
|
|
fsd (export=true, xs=176, xe=400) = sqrt (sqr(ddx(c)) + sqr(ddy(c)) + sqr(ddz(c)))
|
|
|
|
nx = - ddx(c) / fsd
|
|
ny = - ddy(c) / fsd
|
|
nz = - ddz(c) / fsd
|
|
|
|
divn (export=true, xs=176, xe=400) = ddx(nx) + ddy(ny) + ddz(nz)
|
|
|
|
avg { c, fsd, divn }
|