12.07.2015 Views

Metrics of curves in shape optimization and analysis - Andrea Carlo ...

Metrics of curves in shape optimization and analysis - Andrea Carlo ...

Metrics of curves in shape optimization and analysis - Andrea Carlo ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Def<strong>in</strong>ition 10.7 (Convolution) A arc-parameterized convolutional kernel Kalong the curve c <strong>of</strong> length L is a L-periodic function K : lR → lR. Given avector field f : S 1 → lR n <strong>and</strong> a kernel K, we def<strong>in</strong>e the convolution by arcparameter 17 formally as∫(f ⋆ K)(s) := K(s − ŝ)f(ŝ) dŝ . (10.6)By def<strong>in</strong><strong>in</strong>g the run-length function l : lR → lRl(τ) :=c∫ τ0|c ′ (x)| dxwe can rewrite the above eqn. (10.6) explicitly <strong>in</strong> θ parameter as(f ⋆ K)(θ) :=∫ 2π0K ( l(θ) − l(τ) ) f(τ)|c ′ (τ)| dτ . (10.7)Recall the def<strong>in</strong>ition 3.37 <strong>of</strong> gradient ∇E by means <strong>of</strong> the identity〈∇E, h〉 c = DE(c; h) ∀h ∈ T c M .Let f = ∇ H 0E, g = ∇ H 1E be the gradients w.r.to the <strong>in</strong>ner products H 0 <strong>and</strong>H 1 ; by the def<strong>in</strong>ition <strong>of</strong> gradient, we obta<strong>in</strong> that〈f, h〉 H0 ,c = 〈g, h〉 H1 ,c∀h ∈ T c Mthat is 18 ∫∫h · f ds = h · g + λL 2 (D s h · D s g) dsccby <strong>in</strong>tegrat<strong>in</strong>g by parts this becomes∫h · (f − g + λL 2 Dsg) 2 ds = 0 ∀hthen we conclude thatc∀h∇ H 0E = ∇ H 1E − λL 2 D 2 s(∇ H 1E) . (10.8)With similar computation, for ˜H 1 we obta<strong>in</strong> that∫∇ H 0E = ∇˜H1E ds − λL 2 D s(∇˜H1E) 2 . (10.9)cGiven ∇ H 0E, both equations can be solved for ∇ H 1E <strong>and</strong> ∇˜H1E, by means<strong>of</strong> suitable convolution kernels ˜K λ , K λ , that is, we have the formulas∇ H 1E = ∇ H 0E ⋆ K λ , ∇˜H1E = ∇ H 0E ⋆ ˜K λ ;17 Note this def<strong>in</strong>ition is different from the eqn. (13) <strong>in</strong> [55].18 We use the def<strong>in</strong>ition (10.4) <strong>of</strong> H 0 .64

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

Saved successfully!

Ooh no, something went wrong!