routines.inc

这是一个实数编码的遗传算法

interface
   RECURSIVE FUNCTION Differ(func_name,as,param_no,          &
                             eps_of_param_no)
      real(8), dimension(:,:)                ::as
      real(8), dimension(size(as,dim=1))     ::differ
      real(8)  eps_of_param_no
      integer  param_no
      external func_name
   end function Differ
end interface

interface
   FUNCTION Differentiate(func_name,as,param_no,             &
                          eps_of_param_no)
      real(8), dimension(:,:)                ::as
      real(8), dimension(size(as,dim=1))     ::differentiate
      real(8)  eps_of_param_no
      integer  param_no
      external func_name
   end function Differentiate
end interface

interface
   RECURSIVE FUNCTION Integ(func_name,as,param_no,lower,     &
                         upper,half_no_partitions)
      integer  half_no_partitions,param_no
      real(8), dimension(:,:)            ::as
      real(8), dimension(size(as,dim=1)) ::integ,lower,upper
      external func_name
   END function Integ
end interface

interface
   RECURSIVE FUNCTION Integ_over_partition(func_name,as,     &
                               param_no,partition)
      integer, dimension(:)             ::param_no
      real(8), dimension(:,:)           ::as
      real(8), dimension(:,:)           ::partition
      real(8), dimension(size(as,dim=1))::integ_over_partition
      external func_name
   END function Integ_over_partition
end interface

interface
   FUNCTION Integrate(func_name,as,param_no,lower,upper,     &
                   half_no_partitions)
      integer  half_no_partitions,param_no
      real(8), dimension(:,:)            ::as
      real(8), dimension(size(as,dim=1)) ::integrate
      real(8), dimension(size(as,dim=1)) ::lower,upper
      external func_name
   END function Integrate
end interface

interface
   FUNCTION Integrate_over_partition(func_name,as,           &
                                     param_no,partition)
      integer, dimension(:)               ::param_no
      real(8), dimension(:,:)             ::as
      real(8), dimension(:,:)             ::partition
      real(8), dimension(size(as,dim=1))  ::                 &
                                 integrate_over_partition
      external func_name
   END function Integrate_over_partition
end interface

interface
   SUBROUTINE MINIMIZE(spare_params,a,maxnumas,func_name,    &
                       HESSINreset,itmax,schlength,          &
                       dfunctype,eps_vector,num_derivs,      &
                       conv1,conv2,conv3,FRET,nerr)
      real(8), dimension(:)   ::a,eps_vector
      real(8), intent(in)     ::schlength
      integer, intent(in)     ::HESSINreset,itmax,num_derivs
      integer, dimension(:)      ::spare_params
      real(8)  conv1,conv2,conv3,FRET
      integer nerr,maxnumas,dfunctype
      external func_name
   END SUBROUTINE MINIMIZE
end interface

interface
   SUBROUTINE N01RNG(seed,random)
      real(8), dimension(:)   ::random
      integer, dimension(:)   ::seed
   END SUBROUTINE N01RNG
end interface