wolf_powell线搜索方法

ilikenba

<> function wolf_powell(A,b,dir,gradt,x) result(wp_result)
real,dimension(,intent(in)::b,dir,gradt,x
) n( E" o$ n# m7 X3 D4 ` real,dimension(:,,intent(in)::A
  j6 k+ H" [' r7 J  n real::a1,a2,f1,f2,a0,p,q,t,r,wp_result,gradt1(size(gradt)) ,r1
t=2.0
a0=1
" v9 U" l3 a" c* [3 U1 x9 l3 K p=0.1
q=0.9
+ W6 |6 y$ X9 S' b/ u4 U" n6 k a1=0
# U' C/ I# O" n. s6 O a2=100000
. V; I. R! a/ A$ u r=dot_product(dir,gradt)
f1=f(x,A,b)
10 f2=f(x+a0*dir,A,b)
, N  Y, j+ N: t if(f2&gt;f1+p*r*a0)then
! X# F) s6 g' S9 A) d     a2=a0
7 T" V2 {" g# p& `  a0=a1+0.5*(a0-a1)/(1+(f1-f2)/(a0-a1)*r)
: d; w" \4 \$ C( o  goto 10
6 k) U. S0 x" ]& u+ G4 P; _ else
     gradt1=matmul(A,x+a0*dir)+b
  r1=dot_product(gradt1,dir)
  if(r1&gt;=q*r)then
       wp_result=a0
0 }3 n+ D' s$ J/ Q1 B* v; m  |  _9 p  else
# a- {0 e- d- V! B) I       a1=a0
    f1=f2
" m, Q5 h( z$ h$ K7 B& N8 _    r=r1
    a0=a0+(a0-a1)*r1/(r-r1)
9 m6 q. q9 O2 w: F" h5 C    goto 10
  endif
endif
end function
9 J: e+ ?. Q- V8 T( r. _</P>
<>本程序由Fortran 90语言编写，在Visual Fortran 5上编译通过！本程序由沙沙提供！</P>

ammon

真强，佩服，向您学习