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
% h( ^' L; b- [4 Q, F t=2.0
5 @+ X2 X+ A: e a0=1
p=0.1
q=0.9
1 x: b8 j8 h  a, ?9 \: `8 v/ x a1=0
a2=100000
r=dot_product(dir,gradt)
: P& k3 o. L7 k- ^1 p5 A. X1 y6 N f1=f(x,A,b)
2 Z1 B  t, B3 p! Y10 f2=f(x+a0*dir,A,b)
* ~8 y. r. k; A( }# y" w. M5 t if(f2&gt;f1+p*r*a0)then
4 T7 C' V. x" _+ i) i- v. K     a2=a0
  a0=a1+0.5*(a0-a1)/(1+(f1-f2)/(a0-a1)*r)
  goto 10
else
$ I6 `' }) J& A; J# k. O: F$ \+ o/ @     gradt1=matmul(A,x+a0*dir)+b
  r1=dot_product(gradt1,dir)
% _% Q; {8 X! ^/ d- |+ d1 u  if(r1&gt;=q*r)then
       wp_result=a0
  else
       a1=a0
6 I4 N, Y2 ?6 p. O    f1=f2
3 b4 d7 q* l0 z    r=r1
3 K' d+ C* F% z6 T/ h3 P    a0=a0+(a0-a1)*r1/(r-r1)
    goto 10
2 k! H7 \$ p# R, w4 u  ?' w+ d  endif
endif
3 w$ @2 }1 j7 C- N" M; c, c# v end function
</P>
- |, L7 |8 g' W% B: N6 @9 K<>本程序由Fortran 90语言编写，在Visual Fortran 5上编译通过！本程序由沙沙提供！</P>

ammon

真强，佩服，向您学习