anchor.diff:  c Created 29/01/2007 DH.
   
   cf2py integer dimension(n),intent(in) :: x
  -cf2py double precision dimension(np),intent(in) :: p
  +cf2py real dimension(np),intent(in) :: p
   cf2py integer intent(hide),depend(x) :: n=len(x)
   cf2py integer intent(hide),depend(p,n),check(np==1 || np==n) :: np=len(p)
  -cf2py double precision intent(out) :: like      
  +cf2py real intent(out) :: like      
   
         IMPLICIT NONE
         INTEGER n,np,i
         INTEGER x(n)
  -      DOUBLE PRECISION p(np), p_tmp
  -      DOUBLE PRECISION like
  +      REAL p(np), p_tmp
  +      REAL like
   
         p_tmp = p(1)
         like = 0.0
  
positive.code: C
  C     COMPUTE THE LIKELIHOOD OF THE GENERALIZED EXTREME VALUE DISTRIBUTION.
  C
  Cf2py double precision dimension(n), intent(in):: x
  Cf2py double precision dimension(nxi), intent(in):: xi
  Cf2py double precision dimension(nmu), intent(in):: mu
  Cf2py double precision dimension(nsigma), intent(in):: sigma
  Cf2py integer intent(hide), depend(x) :: n=len(x)
  Cf2py integer intent(hide), depend(xi,n),check(nxi==1||nxi==n) :: nxi=len(xi)
  Cf2py integer intent(hide), depend(mu,n),check(nmu==1||nmu==n) :: nmu=len(mu)
  Cf2py integer intent(hide), depend(sigma,n),check(nsigma==1||nsigma==n) :: nsigma=len(sigma)
  Cf2py double precision intent(out):: like
  
        INTEGER n, nmu, nxi, nsigma, i
        DOUBLE PRECISION x(n), xi(nxi), mu(nmu), sigma(nsigma), like
        DOUBLE PRECISION Z(N), EX(N), PEX(N)
        DOUBLE PRECISION XI_tmp, SIGMA_tmp
          
        CALL standardize(x,mu,sigma,n,nmu,nsigma,z)
  
  
negative.code: 
  cf2py integer dimension(nx,k),intent(in) :: x
  cf2py integer dimension(nn), intent(in) :: n
  cf2py real dimension(np,k),intent(in) :: p
  cf2py integer intent(hide),depend(x) :: nx=shape(x,0)
  cf2py integer intent(hide),depend(n) :: nn=shape(n,0)
  cf2py integer intent(hide),depend(p) :: np=shape(p,0)
  cf2py integer intent(hide),depend(x,p),check(k==shape(p,1)) :: k=shape(x,1)
  cf2py real intent(out) :: like      
  
        REAL like,sump
        REAL p(np,k), p_tmp
        INTEGER i,j,ll,n(nn),sumx, n_tmp
        INTEGER x(nx,k)
  
  