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
' y0 Y3 P  N; m- Q5 m: R+ G* q real::a1,a2,f1,f2,a0,p,q,t,r,wp_result,gradt1(size(gradt)) ,r1
t=2.0
a0=1
p=0.1
, R% B/ D6 E* t. `0 ]" G7 i q=0.9
a1=0
; P" i/ J0 \  N* y( H a2=100000
r=dot_product(dir,gradt)
f1=f(x,A,b)
: j8 V3 H4 G, \3 m5 h  C& B$ t10 f2=f(x+a0*dir,A,b)
if(f2&gt;f1+p*r*a0)then
$ l: v% Z3 q- A+ v# m     a2=a0
  a0=a1+0.5*(a0-a1)/(1+(f1-f2)/(a0-a1)*r)
2 |+ C1 U5 w5 I5 X; W+ p  |* C  goto 10
else
     gradt1=matmul(A,x+a0*dir)+b
" t  `. }: g0 p8 [2 z  r# z  r1=dot_product(gradt1,dir)
  if(r1&gt;=q*r)then
$ l5 ?! f7 V& }  p5 U       wp_result=a0
1 R3 S. A. ^5 C  s. u( s) o  else
8 E: a5 e0 A; S) h8 i2 K$ B6 y- p       a1=a0
$ I; O6 Y$ n: U& D    f1=f2
4 N! s+ C* E. x' B  p+ T; ]. B    r=r1
    a0=a0+(a0-a1)*r1/(r-r1)
    goto 10
  endif
3 Y5 ^$ _- f* M. ^6 M% p& f: a endif
0 Q; N0 U3 K, i9 I* Q5 X end function
0 T$ J( X0 j4 B) j</P>
<>本程序由Fortran 90语言编写，在Visual Fortran 5上编译通过！本程序由沙沙提供！</P>

ammon

真强，佩服，向您学习