10.03.2015 Views

v2009.01.01 - Convex Optimization

v2009.01.01 - Convex Optimization

v2009.01.01 - Convex Optimization

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

D.1. DIRECTIONAL DERIVATIVE, TAYLOR SERIES 623<br />

Rotating our perspective, we get several views of the second-order gradient:<br />

⎡<br />

∇ 2 g(X) = ⎢<br />

⎣<br />

∇ 2 g 11 (X) ∇ 2 g 12 (X) · · · ∇ 2 g 1N (X)<br />

∇ 2 g 21 (X) ∇ 2 g 22 (X) · · · ∇ 2 g 2N (X)<br />

. .<br />

∇ 2 g M1 (X) ∇ 2 g M2 (X) · · ·<br />

.<br />

∇ 2 g MN (X)<br />

⎤<br />

⎥<br />

⎦ ∈ RM×N×K×L×K×L (1695)<br />

⎡<br />

∇ 2 g(X) T 1<br />

=<br />

⎢<br />

⎣<br />

∇ ∂g(X)<br />

∂X 11<br />

∇ ∂g(X)<br />

∂X 21<br />

.<br />

∇ ∂g(X)<br />

∂X K1<br />

∇ ∂g(X)<br />

∂X 12<br />

· · · ∇ ∂g(X)<br />

∂X 1L<br />

∇ ∂g(X)<br />

∂X 22<br />

.<br />

· · · ∇ ∂g(X)<br />

∂X 2L<br />

.<br />

∇ ∂g(X)<br />

∂X K2<br />

· · · ∇ ∂g(X)<br />

∂X KL<br />

⎤<br />

⎥<br />

⎦ ∈ RK×L×M×N×K×L (1696)<br />

⎡<br />

∇ 2 g(X) T 2<br />

=<br />

⎢<br />

⎣<br />

∂∇g(X)<br />

∂X 11<br />

∂∇g(X)<br />

∂X 21<br />

.<br />

∂∇g(X)<br />

∂X K1<br />

∂∇g(X)<br />

∂X 12<br />

· · ·<br />

∂∇g(X)<br />

∂X 22<br />

.<br />

· · ·<br />

∂∇g(X)<br />

∂X K2<br />

· · ·<br />

∂∇g(X)<br />

∂X 1L<br />

∂∇g(X)<br />

∂X 2L<br />

.<br />

∂∇g(X)<br />

∂X KL<br />

⎤<br />

⎥<br />

⎦ ∈ RK×L×K×L×M×N (1697)<br />

Assuming the limits exist, we may state the partial derivative of the mn th<br />

entry of g with respect to the kl th and ij th entries of X ;<br />

∂ 2 g mn(X)<br />

∂X kl ∂X ij<br />

=<br />

g mn(X+∆t e<br />

lim<br />

k e T l +∆τ e i eT j )−gmn(X+∆t e k eT l )−(gmn(X+∆τ e i eT)−gmn(X))<br />

j<br />

∆τ,∆t→0<br />

∆τ ∆t<br />

(1698)<br />

Differentiating (1678) and then scaling by Y ij<br />

∂ 2 g mn(X)<br />

∂X kl ∂X ij<br />

Y kl Y ij = lim<br />

∆t→0<br />

∂g mn(X+∆t Y kl e k e T l )−∂gmn(X)<br />

∂X ij<br />

Y<br />

∆t<br />

ij<br />

(1699)<br />

g mn(X+∆t Y kl e<br />

= lim<br />

k e T l +∆τ Y ij e i e T j )−gmn(X+∆t Y kl e k e T l )−(gmn(X+∆τ Y ij e i e T)−gmn(X))<br />

j<br />

∆τ,∆t→0<br />

∆τ ∆t

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!