wolf_powell线搜索方法

ilikenba

<> function wolf_powell(A,b,dir,gradt,x) result(wp_result)
4 }: m7 l, P" P; W! j% q real,dimension(,intent(in)::b,dir,gradt,x
real,dimension(:,,intent(in)::A
, U8 \- X6 g9 F' r9 ^ real::a1,a2,f1,f2,a0,p,q,t,r,wp_result,gradt1(size(gradt)) ,r1
t=2.0
6 G/ {  `- O- o: n a0=1
7 b2 H! ~: X  _+ f1 d9 C p=0.1
7 ?6 @) j. {0 M6 C9 R$ I1 l q=0.9
a1=0
) J/ W  B6 b7 j$ t; R a2=100000
. G/ V5 c* S" {5 h  |; u* R& E& n) N r=dot_product(dir,gradt)
f1=f(x,A,b)
( ^3 i$ }1 J2 S10 f2=f(x+a0*dir,A,b)
. _" e8 c5 C" ]! |% ^ if(f2&gt;f1+p*r*a0)then
     a2=a0
3 f/ }, E* J2 I5 M  a0=a1+0.5*(a0-a1)/(1+(f1-f2)/(a0-a1)*r)
% l3 o4 o5 h4 Z  goto 10
else
" V1 ]) T$ V- ?9 k! H# k) N     gradt1=matmul(A,x+a0*dir)+b
' Y+ r+ O+ V- k9 |* I& x  r1=dot_product(gradt1,dir)
: }* H7 H, j+ G  if(r1&gt;=q*r)then
       wp_result=a0
2 U+ C6 W( B; a" I1 H  else
" b; g; s: l9 S% v0 ^       a1=a0
    f1=f2
3 y' ~# p+ p/ t( r  s" N, O& @    r=r1
    a0=a0+(a0-a1)*r1/(r-r1)
& h, S6 \2 }) @0 [7 Y    goto 10
  endif
endif
, n- i6 b% n% h3 F  R; G" h- o. M: } end function
. n. Q& ?+ K( D</P>
<>本程序由Fortran 90语言编写，在Visual Fortran 5上编译通过！本程序由沙沙提供！</P>

ammon

真强，佩服，向您学习