<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Sat, Feb 7, 2015 at 2:39 AM, Lulu Liu <span dir="ltr"><<a href="mailto:lulu.liu@kaust.edu.sa" target="_blank">lulu.liu@kaust.edu.sa</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">X is the current approximate solution, Y is the Newton direction<div><br clear="all"><div>Step 1. the unknown X=[phi,omega,S], the third component S in [0,1]. Hence, we firstly compute a multiplier alpha such that</div><div><br></div><div>X1= X-alpha*Y, and the component S of X1 is still in [0,1].</div><div>In step1, we need access the component S to calculate multiplier.<br></div><div><br></div><div>Step 2. Let the damping factor damp=1;</div><div><br></div><div>for n=1:N</div><div><br></div><div>    lambda=alpha*damp;</div><div>    X1=X-lambda*Y;</div><div><br></div><div>     Solve J(X1)Y1 = F(X1);</div><div><br></div><div>     if ||Y1||<||Y||</div><div>           break;</div><div>      else</div><div>            damp = damp/2;</div><div>       end</div><div><span>endfor</span></div><div><br></div><div>the norm || || is weighted norm, and we need access each component and compute some weights based on each component.</div></div></div></blockquote><div><br></div><div>This is just a line search, so first it should be implemented as a SNESLineSearch. Second, you should really</div><div>try the new nleqerr line search that Patrick Farrell just implemented first, since it may already solve your problem.</div><div>Third, if that weight norm is important, I think you just use that norm in nleqerr.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>How to implement the damped Newton method above? Is it a easier way to do it? </div><div><br></div><div>Thanks!</div><div><br></div><br>
</div></div>

<br>
<div><hr></div><font face="Arial" size="1">This message and its contents, including attachments are intended solely for the original recipient. If you are not the intended recipient or have received this message in error, please notify me immediately and delete this message from your computer system. Any unauthorized use or distribution is prohibited. Please consider the environment before printing this email.</font></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div></div>