wolf_powell线搜索方法

ilikenba

<> function wolf_powell(A,b,dir,gradt,x) result(wp_result)
7 [! }6 E0 f1 C% R, G real,dimension(,intent(in)::b,dir,gradt,x
1 M; N. S2 m8 @ real,dimension(:,,intent(in)::A
' b1 H! [% @. Z; g2 ?6 A* b real::a1,a2,f1,f2,a0,p,q,t,r,wp_result,gradt1(size(gradt)) ,r1
8 m" @8 E  f" E% H+ [: l9 O t=2.0
  z; e; ?  V$ P' d/ [' ?1 N! L a0=1
p=0.1
q=0.9
a1=0
( \& V: |# z. c5 H$ r2 p. w- z a2=100000
4 @+ M0 C% m  T. { r=dot_product(dir,gradt)
3 L3 |. ]( R1 n f1=f(x,A,b)
10 f2=f(x+a0*dir,A,b)
' S" c" l. [  G6 a0 a, @  n if(f2&gt;f1+p*r*a0)then
* C7 ~7 g/ w( G, v     a2=a0
  a0=a1+0.5*(a0-a1)/(1+(f1-f2)/(a0-a1)*r)
6 T5 N  {7 B7 f  s) Z  goto 10
else
     gradt1=matmul(A,x+a0*dir)+b
  r1=dot_product(gradt1,dir)
  if(r1&gt;=q*r)then
7 O$ }2 `9 {4 X" I# a       wp_result=a0
1 L3 J/ \2 F1 K0 y% L+ O2 W  else
  ^( t$ ]$ `% {; M" h$ b       a1=a0
7 X' U( L( h! F% u, \& p+ f    f1=f2
    r=r1
    a0=a0+(a0-a1)*r1/(r-r1)
& F' _) t9 }2 v; K    goto 10
, H# U* y% D1 |5 e% {  endif
endif
end function
& X7 F. ^2 ?$ |) L0 h</P>
1 O. ?5 I7 R! W<>本程序由Fortran 90语言编写，在Visual Fortran 5上编译通过！本程序由沙沙提供！</P>

ammon

真强，佩服，向您学习