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
; c3 W& q- K7 A! y1 R6 q9 Q9 N: y+ h t=2.0
a0=1
p=0.1
q=0.9
a1=0
. r( u' P! H4 j* x- u a2=100000
r=dot_product(dir,gradt)
4 D. W; ^% P8 T* C' D6 [& M f1=f(x,A,b)
  I; A0 g$ k5 T' ]) y# A. j10 f2=f(x+a0*dir,A,b)
if(f2&gt;f1+p*r*a0)then
     a2=a0
4 \  [0 l0 J8 C; [% h# ^; B  a0=a1+0.5*(a0-a1)/(1+(f1-f2)/(a0-a1)*r)
* v' D: w0 x3 P  e  goto 10
; f9 ]. T; C) D7 f; o else
     gradt1=matmul(A,x+a0*dir)+b
  ?" U6 L- ~" I, l: c$ [  r1=dot_product(gradt1,dir)
  if(r1&gt;=q*r)then
! y+ V0 _; ?( J! M       wp_result=a0
  else
* C! q' @2 @. S/ I5 K( n       a1=a0
    f1=f2
5 ^; S- x" R* j* K1 d9 x    r=r1
3 b0 g- Q& d% g& l% R    a0=a0+(a0-a1)*r1/(r-r1)
& g6 W6 [$ D3 B7 z* y1 |  m    goto 10
8 ^9 x  |0 a1 K7 M( I. M' l4 J# M  endif
' ?2 `5 D, H6 U endif
end function
/ ?/ j7 g/ X+ {; N7 A/ u</P>
<>本程序由Fortran 90语言编写，在Visual Fortran 5上编译通过！本程序由沙沙提供！</P>

ammon

真强，佩服，向您学习