wolf_powell线搜索方法

ilikenba

<> function wolf_powell(A,b,dir,gradt,x) result(wp_result)
# ]/ U/ t1 |  o# r$ J real,dimension(,intent(in)::b,dir,gradt,x
1 ?; g$ e7 c7 m/ V6 z& N 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
/ V" B. s8 J) Z8 z: Y- T0 A5 }( p p=0.1
  q  W" g" T( D/ `; | q=0.9
a1=0
/ l# i1 k: X0 @; x" d a2=100000
% ]' f) f1 P0 @) P# l- J r=dot_product(dir,gradt)
& [3 K2 r3 v! h" r f1=f(x,A,b)
10 f2=f(x+a0*dir,A,b)
if(f2&gt;f1+p*r*a0)then
. V1 j! n' k. p3 R! }     a2=a0
9 p. y3 [0 v; M1 n  a0=a1+0.5*(a0-a1)/(1+(f1-f2)/(a0-a1)*r)
  }2 ?  E. A" J: M! U3 K7 K! r/ u  goto 10
else
     gradt1=matmul(A,x+a0*dir)+b
0 A6 ]' P6 u5 |& ?/ a6 j  r1=dot_product(gradt1,dir)
' j2 \7 |" ^3 F$ W2 o' r  if(r1&gt;=q*r)then
       wp_result=a0
  else
' e) {- v- l7 F3 ~$ w$ M       a1=a0
    f1=f2
    r=r1
. _9 C8 F/ S2 R/ [6 }    a0=a0+(a0-a1)*r1/(r-r1)
    goto 10
5 e/ M4 v- Y8 b- [  endif
, @+ m* O* I7 u+ Z8 I% s endif
" o1 R! t, Z, ? end function
</P>
1 O) }) Y7 l& B<>本程序由Fortran 90语言编写，在Visual Fortran 5上编译通过！本程序由沙沙提供！</P>

ammon

真强，佩服，向您学习