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
real::a1,a2,f1,f2,a0,p,q,t,r,wp_result,gradt1(size(gradt)) ,r1
( u" E8 p7 K- \5 x t=2.0
( H. M& g; Y; q9 E* E a0=1
p=0.1
. X! Q: h. U* k6 Z" b& q$ J q=0.9
$ F# K' S, D, A) e4 Q. A% j a1=0
a2=100000
3 ]5 U0 [( X7 p% `) S r=dot_product(dir,gradt)
# |) s" d( q$ B( g: C+ w; d# F: R f1=f(x,A,b)
8 ?& S# I& c1 q1 ^5 @# s- d10 f2=f(x+a0*dir,A,b)
/ Q3 T3 W, i2 j& V if(f2&gt;f1+p*r*a0)then
: @/ ^) w) R- n8 W* p0 T, v     a2=a0
$ d/ T! k, y5 q0 Y1 X, x7 z  a0=a1+0.5*(a0-a1)/(1+(f1-f2)/(a0-a1)*r)
  goto 10
else
% V0 m0 ]4 t& g9 I! T     gradt1=matmul(A,x+a0*dir)+b
  r1=dot_product(gradt1,dir)
  if(r1&gt;=q*r)then
! q% A+ [8 Q8 s       wp_result=a0
  else
. l6 |8 }9 ]+ g  b4 d  E( F$ A$ t       a1=a0
% o- U1 \, O' c4 P* L- @& ]7 @    f1=f2
    r=r1
! \9 \; }2 D  n& S    a0=a0+(a0-a1)*r1/(r-r1)
) r9 n& }% B+ M" g. Z  d/ m    goto 10
  endif
$ T% ~0 ^2 t5 N3 u endif
/ T& Y8 R: p7 e! k7 r: \; s0 {' m end function
$ a. E( z1 i0 E: b6 X</P>
4 G/ s, G! d& e6 W5 ?# L<>本程序由Fortran 90语言编写，在Visual Fortran 5上编译通过！本程序由沙沙提供！</P>

ammon

真强，佩服，向您学习