wolf_powell线搜索方法

ilikenba

<> function wolf_powell(A,b,dir,gradt,x) result(wp_result)
real,dimension(,intent(in)::b,dir,gradt,x
% x: `4 W7 U8 U9 f$ |' f! } real,dimension(:,,intent(in)::A
; i. ~- A0 F6 ^8 N  }- D real::a1,a2,f1,f2,a0,p,q,t,r,wp_result,gradt1(size(gradt)) ,r1
t=2.0
/ H  l, Z7 U8 @1 u a0=1
p=0.1
q=0.9
a1=0
a2=100000
r=dot_product(dir,gradt)
# h  N8 H8 s9 _- \1 K/ ]- X- ` f1=f(x,A,b)
+ v/ h" M+ F. _7 T( Z8 [10 f2=f(x+a0*dir,A,b)
6 s1 k( \2 {! m2 R2 g) N. ^, C if(f2&gt;f1+p*r*a0)then
8 x: g  s: M3 p. C2 r" u     a2=a0
  a0=a1+0.5*(a0-a1)/(1+(f1-f2)/(a0-a1)*r)
9 E$ @, i" B" C# n2 g  goto 10
else
% A0 d2 c( ?9 V7 q7 I     gradt1=matmul(A,x+a0*dir)+b
  r1=dot_product(gradt1,dir)
  if(r1&gt;=q*r)then
       wp_result=a0
  else
* X5 a" O& r+ \, w0 u0 b0 S       a1=a0
    f1=f2
1 k3 n( ?; ?. p# q- s( |% G    r=r1
- k+ b8 p. S- K# }0 y( u9 n0 c    a0=a0+(a0-a1)*r1/(r-r1)
8 T$ \' W' Z1 j5 m4 M6 {" Y4 U, e    goto 10
# d! m' g/ t) m1 I( @$ Z  endif
  E$ P# C8 H8 x' a endif
* n6 U& m2 Y! G4 a3 S4 C end function
& a0 K. T; g7 a) h3 H# X</P>
4 v+ R% G  r9 w# }) m<>本程序由Fortran 90语言编写，在Visual Fortran 5上编译通过！本程序由沙沙提供！</P>

ammon

真强，佩服，向您学习