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
9 {% b4 g5 F. Z real::a1,a2,f1,f2,a0,p,q,t,r,wp_result,gradt1(size(gradt)) ,r1
t=2.0
a0=1
p=0.1
1 i( k: W) K  v4 o" f9 ^5 M q=0.9
" s' q( ]  `, F7 T0 q  U a1=0
: ?. d- X. p# Z. z2 _1 J a2=100000
, c* a& u/ B; u. T. `; y r=dot_product(dir,gradt)
2 a3 L( V0 }, y! t f1=f(x,A,b)
10 f2=f(x+a0*dir,A,b)
  y) E7 |9 W% N7 _, e$ a& l2 t) r if(f2&gt;f1+p*r*a0)then
     a2=a0
* S6 C6 H  t" i  a0=a1+0.5*(a0-a1)/(1+(f1-f2)/(a0-a1)*r)
  goto 10
else
$ \8 v. u: e4 J     gradt1=matmul(A,x+a0*dir)+b
: T" F8 R9 y/ E  r1=dot_product(gradt1,dir)
  if(r1&gt;=q*r)then
& g8 K+ M0 H- W5 m8 f       wp_result=a0
2 \( d6 g) N- k$ a  R* D8 X  else
: i% r: u/ B2 L4 X) {" z+ H3 b; S; f       a1=a0
# v0 y" U, d% I+ B( N+ o" a) G7 i    f1=f2
  R3 q% e) N; i( k* ?% I; U1 O# e    r=r1
    a0=a0+(a0-a1)*r1/(r-r1)
# _1 ?7 f- S; s    goto 10
  endif
endif
. z" h+ X% o+ v) h end function
( ?" T) G6 p. C# M3 g</P>
<>本程序由Fortran 90语言编写，在Visual Fortran 5上编译通过！本程序由沙沙提供！</P>

ammon

真强，佩服，向您学习