wolf_powell线搜索方法

ilikenba

<> function wolf_powell(A,b,dir,gradt,x) result(wp_result)
real,dimension(,intent(in)::b,dir,gradt,x
' q/ ]0 z  {* w  A2 m& K real,dimension(:,,intent(in)::A
1 Z1 J* z5 |+ p0 A6 q, \ real::a1,a2,f1,f2,a0,p,q,t,r,wp_result,gradt1(size(gradt)) ,r1
$ B  R3 N  k. ?* L* { t=2.0
8 H" c5 j' D" U a0=1
p=0.1
, R+ R0 }4 G& E$ [ q=0.9
a1=0
/ X+ q  k, C# c8 { a2=100000
r=dot_product(dir,gradt)
! J) ]5 P  U0 o1 L& W8 c f1=f(x,A,b)
) l: G/ K3 f6 c4 W1 l10 f2=f(x+a0*dir,A,b)
) D3 P7 M& A) v  e if(f2&gt;f1+p*r*a0)then
8 L# d1 N$ P; ?     a2=a0
9 j  N. M: ]" w; G" C) M# n  a0=a1+0.5*(a0-a1)/(1+(f1-f2)/(a0-a1)*r)
  goto 10
7 y3 H7 g) I! ]" `1 q- |8 N9 Y* R else
" C, D' d6 B0 I. E4 o9 o- G2 E1 G     gradt1=matmul(A,x+a0*dir)+b
  r1=dot_product(gradt1,dir)
; G  W. v( ?4 [0 v! |. X  if(r1&gt;=q*r)then
4 s9 q1 y; _/ L& M/ t3 `6 r       wp_result=a0
  else
2 I! R) m- h7 t; r9 Z) B( |7 b0 P0 q       a1=a0
# \. o% ^, ^0 t+ g3 O    f1=f2
" M+ p- z4 r9 c2 }4 d9 }    r=r1
    a0=a0+(a0-a1)*r1/(r-r1)
0 E: k8 r! Y  ~% K( i# _    goto 10
  endif
endif
end function
</P>
<>本程序由Fortran 90语言编写，在Visual Fortran 5上编译通过！本程序由沙沙提供！</P>

ammon

真强，佩服，向您学习