wolf_powell线搜索方法

ilikenba

<> function wolf_powell(A,b,dir,gradt,x) result(wp_result)
real,dimension(,intent(in)::b,dir,gradt,x
real,dimension(:,,intent(in)::A
( O& Z8 F3 U+ H1 U, U2 O0 `8 z real::a1,a2,f1,f2,a0,p,q,t,r,wp_result,gradt1(size(gradt)) ,r1
. a0 o, B' t8 ]+ [* O8 Q t=2.0
6 E/ n1 a, m, @, S  S1 p a0=1
p=0.1
0 \1 l2 R# F0 p q=0.9
a1=0
a2=100000
& c3 m4 M; v3 e* H8 {& y; y r=dot_product(dir,gradt)
6 e! e" u4 ^4 ]- w f1=f(x,A,b)
9 x# E! _' Y3 [; a- N7 a10 f2=f(x+a0*dir,A,b)
. b7 L* ]; L6 B% x if(f2&gt;f1+p*r*a0)then
. @" k% M5 @, l/ S; j3 w( N3 k. b     a2=a0
, z" z+ H% m& ?2 D  a0=a1+0.5*(a0-a1)/(1+(f1-f2)/(a0-a1)*r)
0 @4 b8 P7 W5 f/ }2 z/ p' H  n- b+ g# e* a  goto 10
else
5 s9 N3 @+ U- G+ r+ r     gradt1=matmul(A,x+a0*dir)+b
, Y6 b2 `' ]8 G( U' |  A  r1=dot_product(gradt1,dir)
  if(r1&gt;=q*r)then
) g2 [. E  A8 @  U$ ]& z+ D       wp_result=a0
  else
       a1=a0
! z; G# X: m" k3 p; u* w5 {: X    f1=f2
+ V9 i' f) n) C5 r    r=r1
' q/ u( T: }, B    a0=a0+(a0-a1)*r1/(r-r1)
4 S% ]. Y6 N. z" D( }' Z. P    goto 10
  endif
. ?) h( u% _, C' b endif
end function
0 Z# B) b" C6 L& J- W6 w</P>
( Y+ M% `6 X; _( o0 C; S2 A0 w% n<>本程序由Fortran 90语言编写，在Visual Fortran 5上编译通过！本程序由沙沙提供！</P>

ammon

真强，佩服，向您学习