wolf_powell线搜索方法

ilikenba

<> function wolf_powell(A,b,dir,gradt,x) result(wp_result)
& o: B5 f( V+ g& @/ I8 q, K( | real,dimension(,intent(in)::b,dir,gradt,x
real,dimension(:,,intent(in)::A
3 n, f0 M- _0 Q4 {) m4 \& [ real::a1,a2,f1,f2,a0,p,q,t,r,wp_result,gradt1(size(gradt)) ,r1
0 z9 M& t: ]) V6 v4 ] t=2.0
a0=1
p=0.1
2 e+ }% L9 ~, J0 a- y q=0.9
a1=0
a2=100000
8 n$ D5 x5 Z, O0 O" [7 F& e6 o9 s r=dot_product(dir,gradt)
' O" c0 [7 r0 h. }# @) j f1=f(x,A,b)
10 f2=f(x+a0*dir,A,b)
9 E# Y. n1 N2 S: |5 E" ]9 P! q if(f2&gt;f1+p*r*a0)then
     a2=a0
  a0=a1+0.5*(a0-a1)/(1+(f1-f2)/(a0-a1)*r)
  goto 10
2 s8 Q* T( X  D else
4 n* t: H2 G# i/ U     gradt1=matmul(A,x+a0*dir)+b
  r1=dot_product(gradt1,dir)
; h6 e6 d. J0 T8 O4 ^' p# Z3 ^  if(r1&gt;=q*r)then
6 F3 P2 G, I' x: X8 w/ m2 _       wp_result=a0
) |% g4 k2 p/ _0 }. u- ]  else
' x1 k. y2 t4 t5 d       a1=a0
$ E0 N4 g) f' z1 \. b    f1=f2
    r=r1
    a0=a0+(a0-a1)*r1/(r-r1)
% @+ I/ T) T% a4 n& L0 r    goto 10
  endif
endif
, ?# d+ |  N6 p( y, g( U7 O end function
, Q- u4 b( ^, s; c3 c</P>
& ]% v8 h0 y# v- u2 m1 g% w<>本程序由Fortran 90语言编写，在Visual Fortran 5上编译通过！本程序由沙沙提供！</P>

ammon

真强，佩服，向您学习