wolf_powell线搜索方法

ilikenba

<> function wolf_powell(A,b,dir,gradt,x) result(wp_result)
- M- v! |- |1 ~+ W# [ real,dimension(,intent(in)::b,dir,gradt,x
real,dimension(:,,intent(in)::A
/ y3 i4 U: G! A( T3 s- y real::a1,a2,f1,f2,a0,p,q,t,r,wp_result,gradt1(size(gradt)) ,r1
" A0 r5 L+ U, j; F' u/ h2 S t=2.0
2 w) E+ v" ]0 Q a0=1
p=0.1
/ M# k' U) u  t' |: R! r; ~9 U3 a# R q=0.9
( j0 d% |3 w+ O) w* `% \9 f% P- [: a a1=0
1 r2 b% S( r: \0 r' g9 W a2=100000
; p, u  k9 C6 K$ _# V  [: F" n5 ]6 B r=dot_product(dir,gradt)
2 R5 c& e# T) U( K+ f* n6 g/ S& ^ f1=f(x,A,b)
8 V2 r1 d) V/ E7 W1 I) [6 B6 N6 b6 V6 ^10 f2=f(x+a0*dir,A,b)
if(f2&gt;f1+p*r*a0)then
     a2=a0
# c5 E: E! q/ o! P6 C0 Q3 J( x  a0=a1+0.5*(a0-a1)/(1+(f1-f2)/(a0-a1)*r)
9 d" I  q( Z" N2 N8 C  goto 10
else
( D& ]1 L: [- I- {     gradt1=matmul(A,x+a0*dir)+b
( d; G' U' r& {' a  r1=dot_product(gradt1,dir)
. v- Z* N: m. r) ^$ I  if(r1&gt;=q*r)then
       wp_result=a0
2 u7 G: T% r1 r2 ~2 u: l' `  else
) ^  w, q0 G+ [: K5 y       a1=a0
! y8 [7 e0 J8 W. G# u; _    f1=f2
9 |8 [0 H8 ]( M2 L0 d: `    r=r1
    a0=a0+(a0-a1)*r1/(r-r1)
    goto 10
+ R; J; y1 h# A6 P$ j, Y% Y% a  endif
endif
; Z% y( m8 X6 W9 W: j! i/ \ end function
</P>
<>本程序由Fortran 90语言编写，在Visual Fortran 5上编译通过！本程序由沙沙提供！</P>

ammon

真强，佩服，向您学习