wolf_powell线搜索方法

ilikenba

<> function wolf_powell(A,b,dir,gradt,x) result(wp_result)
real,dimension(,intent(in)::b,dir,gradt,x
, R$ l  [0 D4 z$ b& b$ i; z3 g real,dimension(:,,intent(in)::A
  L' J  O' _  @# V. l* [ real::a1,a2,f1,f2,a0,p,q,t,r,wp_result,gradt1(size(gradt)) ,r1
  i( M, E2 s" C8 x# a$ ]2 Z+ r t=2.0
. w- r1 b1 O+ D; m% m. r6 `# K a0=1
1 N" E6 c3 h% X p=0.1
q=0.9
a1=0
a2=100000
0 M. M& P5 @& [! X* C( G8 _' T, _ r=dot_product(dir,gradt)
: ~! L7 d( Q% A f1=f(x,A,b)
10 f2=f(x+a0*dir,A,b)
. k: D# Q5 h- a8 \6 ~$ c# v& U: n if(f2&gt;f1+p*r*a0)then
     a2=a0
7 U4 L3 g* ^" m  a0=a1+0.5*(a0-a1)/(1+(f1-f2)/(a0-a1)*r)
  goto 10
else
     gradt1=matmul(A,x+a0*dir)+b
+ l+ V2 J! P, a# e% x3 s  r1=dot_product(gradt1,dir)
  if(r1&gt;=q*r)then
       wp_result=a0
  else
       a1=a0
    f1=f2
+ z2 h0 _: m3 m* O' y; d    r=r1
" z, ~+ U9 \; X9 e3 W) n    a0=a0+(a0-a1)*r1/(r-r1)
7 q! O& c/ J4 t; W' i1 a. }    goto 10
* a7 K( n3 u- z9 O6 X5 ^2 X  endif
endif
end function
. K" A7 w& [! q1 ^) `/ n</P>
& E3 W* m* s0 j" c<>本程序由Fortran 90语言编写，在Visual Fortran 5上编译通过！本程序由沙沙提供！</P>

ammon

真强，佩服，向您学习