set S="33 34 35"

echo "-1 1   -1 1   -1 1   -1 1" | unu 2op x - 1.6 \
 | unu reshape -s 2 2 2 \
 | unu resample -s $S -k tent -c node \
 | unu 2op pow - 2 -o xx.nrrd

echo "-1 1   -1 1   -1 1   -1 1" | unu 2op x - 1.6 \
 | unu reshape -s 2 2 2 \
 | unu swap -a 0 1 \
 | unu resample -s $S -k tent -c node \
 | unu 2op pow - 2 -o yy.nrrd

echo "-1 1   -1 1   -1 1   -1 1" | unu 2op x - 1.6 \
 | unu reshape -s 2 2 2 \
 | unu swap -a 0 2 \
 | unu resample -s $S -k tent -c node \
 | unu 2op pow - 2 -o zz.nrrd

unu 2op spow xx.nrrd 3 -o xxxx.nrrd
unu 2op spow yy.nrrd 3 -o yyyy.nrrd
unu 2op pow zz.nrrd 2 -o zzzz.nrrd

unu 2op + xx.nrrd yy.nrrd \
 | unu 2op - - xxxx.nrrd \
 | unu 2op - - yyyy.nrrd \
 | unu 2op x - 0.6 \
 | unu 2op - - zz.nrrd \
 | unu axinfo -a 0 -mm -1 1 \
 | unu axinfo -a 1 -mm -1 1 \
 | unu axinfo -a 2 -mm -1 1 \
 | unu 3op max_sm -0.3 -0.1 - \
 | unu 2op - - 0.17 \
 | ../../bin/dnorm -i - -o - \
 | unu pad -min -2 -2 0 -max M+2 M+2 M \
 | unu crop -min 0 0 2 -max M M M-2 -o ring.nrrd

unu resample -i ring.nrrd -s x1 x1 x1 -k bspln3ai -o ring-bspln3.nrrd

rm -f {xx,yy,zz}.nrrd {xxxx,yyyy,zzzz}.nrrd