I want to calculate the directional derivative of $f(x,y,z) = e^{xyz}, v = (1,1,-1)$ without recourse to the Jacobian matrix (in which case this is a very easy exercise). Here's what I've got so far: $$\frac{\partial}{\partial_v}f=\lim_{h \to 0} \frac{e^{(x+h)(y+h)(z-h)}-e^{xyz}}{h}=e^{xyz} \cdot \lim_{h \to 0} \frac{e^{h(-xy+xz+yz)+h^2(-x-y+z)-h^3}-1}{h}.$$ The answer is now intuitively clear, i.e. the higher order $h$ terms go to zero faster than the linear $h$ term, thus I'm left with$$ e^{xyz} \cdot \lim_{h \to 0} \frac{e^{h(-xy+xz+yz)}-1}{h} = e^{xyz} (-xy+xz+yz). $$ But two things to be made precise:
$1.$ why is $\lim_{h \to 0} \frac{e^{hk}-1}{h} =k $
$2.$ why exactly can I let the higher order $h$ terms in the exponent go to zero?
By definition the directional derivative you are after is the quantity $${d\over dh}e^{(x+h)(y+h)(z-h)}\biggr|_{h=0}\ .$$ Now $${d\over dh}\bigl((x+h)(y+h)(z-h)\bigr)=(y+h)(z-h)+(x+h)(z-h)-(x+h)(y+h)\ ,$$ and letting $h=0$ here gives $yz+xz-xy$ on the RHS. It follows that $${d\over dh}e^{(x+h)(y+h)(z-h)}\biggr|_{h=0}=(yz+xz-xy)e^{xyz}\ .$$ Concerning your last questions: