wolf_powell线搜索方法

ilikenba

<> function wolf_powell(A,b,dir,gradt,x) result(wp_result)
1 |9 K6 u: S4 G4 b9 H real,dimension(,intent(in)::b,dir,gradt,x
% |  \4 k$ X: }# M; Y5 v6 P real,dimension(:,,intent(in)::A
2 V  r2 I. h6 X! ?. Y: A real::a1,a2,f1,f2,a0,p,q,t,r,wp_result,gradt1(size(gradt)) ,r1
  P* [9 e$ ~$ Q* T4 s, p" C t=2.0
9 s8 r% {+ P  O3 F0 G a0=1
6 m+ p1 R' ]+ F  E p=0.1
2 i, A6 @; a6 M" \ q=0.9
4 V0 s- f' V, t9 A3 h, e a1=0
a2=100000
r=dot_product(dir,gradt)
% b8 H! A: {4 w7 K f1=f(x,A,b)
2 N7 C4 i+ a1 Y3 C5 r10 f2=f(x+a0*dir,A,b)
if(f2&gt;f1+p*r*a0)then
* z7 R; u( ~& ?& G- a1 _     a2=a0
  a0=a1+0.5*(a0-a1)/(1+(f1-f2)/(a0-a1)*r)
* y7 R8 Z/ F! V  goto 10
1 F! G* [3 N) C% M/ m. G else
3 |8 Y1 c; [6 V, t* K     gradt1=matmul(A,x+a0*dir)+b
  r1=dot_product(gradt1,dir)
  if(r1&gt;=q*r)then
       wp_result=a0
) _5 U* u2 _6 X/ b; u4 c+ Q3 s  else
1 u; x) M  R( B; |; J( l- n0 B       a1=a0
    f1=f2
3 h2 Q8 Z1 k% K9 v    r=r1
    a0=a0+(a0-a1)*r1/(r-r1)
- H% j/ s% h2 Y2 t5 N0 ?/ {% i    goto 10
  endif
  a& F+ ?& h$ E# M' H$ V endif
end function
; f  {0 F9 u8 S- P) x, ~( u</P>
<>本程序由Fortran 90语言编写，在Visual Fortran 5上编译通过！本程序由沙沙提供！</P>

ammon

真强，佩服，向您学习