wolf_powell线搜索方法

ilikenba

<> function wolf_powell(A,b,dir,gradt,x) result(wp_result)
2 @8 D; |" ?2 P8 `: I. T real,dimension(,intent(in)::b,dir,gradt,x
real,dimension(:,,intent(in)::A
4 p7 i( }' Q1 \5 J: N; b" n4 s real::a1,a2,f1,f2,a0,p,q,t,r,wp_result,gradt1(size(gradt)) ,r1
t=2.0
) A4 r& p4 e: y' K3 G; K% o a0=1
5 ?, r  A  u3 n' J+ P p=0.1
q=0.9
2 e# r$ Y% B% n/ h& j) O2 P a1=0
a2=100000
r=dot_product(dir,gradt)
( j9 o3 e. y8 U: h* d, I3 \ f1=f(x,A,b)
10 f2=f(x+a0*dir,A,b)
if(f2&gt;f1+p*r*a0)then
* x) a. J( m- J0 B     a2=a0
  a0=a1+0.5*(a0-a1)/(1+(f1-f2)/(a0-a1)*r)
  goto 10
& r' r. B" V3 q$ X else
     gradt1=matmul(A,x+a0*dir)+b
7 A9 M0 O1 m2 Q! K  r1=dot_product(gradt1,dir)
! R" u* }) y: u  if(r1&gt;=q*r)then
/ M4 D1 H+ w+ L; Q8 ?7 c       wp_result=a0
  else
$ `1 E8 e' l5 c' D4 ]: }& _       a1=a0
7 ?5 @3 L* U1 Y1 J: ?; V    f1=f2
    r=r1
! a' I7 ^: F4 G. a- S7 k    a0=a0+(a0-a1)*r1/(r-r1)
    goto 10
  endif
endif
end function
& o" v2 R, E& ^2 A8 N3 d</P>
<>本程序由Fortran 90语言编写，在Visual Fortran 5上编译通过！本程序由沙沙提供！</P>

ammon

真强，佩服，向您学习