wolf_powell线搜索方法

ilikenba

<> function wolf_powell(A,b,dir,gradt,x) result(wp_result)
3 b% v, r# w, Y& v real,dimension(,intent(in)::b,dir,gradt,x
real,dimension(:,,intent(in)::A
real::a1,a2,f1,f2,a0,p,q,t,r,wp_result,gradt1(size(gradt)) ,r1
t=2.0
a0=1
; A6 k1 N( z* [7 b6 B* } p=0.1
q=0.9
& B9 o/ o7 A/ v/ R& q2 f- c a1=0
a2=100000
+ W0 I- L2 F9 X: }7 ~. w. v4 K r=dot_product(dir,gradt)
% Z. _7 O2 c: `" Q, B- j) x; [1 [0 \ f1=f(x,A,b)
10 f2=f(x+a0*dir,A,b)
if(f2&gt;f1+p*r*a0)then
. l  m. j% _- |     a2=a0
+ F. {0 l2 Q- U  a0=a1+0.5*(a0-a1)/(1+(f1-f2)/(a0-a1)*r)
  goto 10
else
     gradt1=matmul(A,x+a0*dir)+b
  r1=dot_product(gradt1,dir)
  if(r1&gt;=q*r)then
( T9 ]) k9 Y+ s+ w  j       wp_result=a0
9 [7 {/ N1 Y+ j7 m0 x* l  else
       a1=a0
  A) R$ q* J  |* R2 ~& a    f1=f2
8 d; k, ~5 {( _# i    r=r1
. e3 s! o5 i  C3 l    a0=a0+(a0-a1)*r1/(r-r1)
    goto 10
/ f3 p, x, p$ U3 w6 v3 N" o  endif
endif
% l. S' ]; D0 Y+ {" o end function
9 `& N6 o7 l6 r7 H' H! D</P>
<>本程序由Fortran 90语言编写，在Visual Fortran 5上编译通过！本程序由沙沙提供！</P>

ammon

真强，佩服，向您学习