#!/bin/csh -f
# RCSid $Id: pmdblur.csh,v 3.1 2005/01/18 03:59:41 greg Exp $
#
# Generate views for motion and depth blurring on picture
#
if ($#argv != 5) then
	echo "Usage: $0 speed aperture nsamp v0file v1file"
	exit 1
endif
set s = "$1"
set a = "$2"
set n = "$3"
set vc = "$4"
set vn = "$5"
if (`ev "if($s-.01,0,1)"`) then
	pdfblur $a $n $vc
	exit
endif
if (`ev "if($a,0,1)"`) then
	pmblur $s $n $vc $vn
	exit
endif
cnt $n | rcalc -e `vwright C < $vc` -e `vwright N < $vn` \
-e "t=$s/$n"'*($1+rand($1))' \
-e "r=$a/2"'*sqrt(rand(182+7*$1));theta=2*PI*rand(-10-$1)' \
-e 'rcost=r*cos(theta);rsint=r*sin(theta)' \
-e 'opx= (1-t)*(Cpx+rcost*Chx+rsint*Cvx) + t*(Npx+rcost*Nhx+rsint*Nvx)' \
-e 'opy= (1-t)*(Cpy+rcost*Chy+rsint*Cvy) + t*(Npy+rcost*Nhy+rsint*Nvy)' \
-e 'opz= (1-t)*(Cpz+rcost*Chz+rsint*Cvz) + t*(Npz+rcost*Nhz+rsint*Nvz)' \
-e 'odx= (1-t)*Cdx*Cd + t*Ndx*Nd' \
-e 'ody= (1-t)*Cdy*Cd + t*Ndy*Nd' \
-e 'odz= (1-t)*Cdz*Cd + t*Ndz*Nd' \
-e 'oux=(1-t)*Cux+t*Nux;ouy=(1-t)*Cuy+t*Nuy;ouz=(1-t)*Cuz+t*Nuz' \
-e 'oh=(1-t)*Ch+t*Nh;ov=(1-t)*Cv+t*Nv' \
-e 'os= (1-t)*(Cs-rcost/(Cd*Chn)) + t*(Ns-rcost/(Nd*Nhn))' \
-e 'ol= (1-t)*(Cl-rsint/(Cd*Cvn)) + t*(Nl-rsint/(Nd*Nvn))' \
-o 'VIEW= -vp ${opx} ${opy} ${opz} -vd ${odx} ${ody} ${odz} -vu ${oux} ${ouy} ${ouz} -vh ${oh} -vv ${ov} -vs ${os} -vl ${ol}'