anchor.diff:  ! DAXPY(N,DA,DX,INCX,DY,INCY) Adding vectors
         EXTERNAL DCOPY
   ! DCOPY(N,DX,INCX,DY,INCY) copies x to y
  -      EXTERNAL DDOT
   
         twopi_N = 0.5D0 * N * dlog(2.0D0*PI)
   
  
positive.code: ! DAXPY(N,DA,DX,INCX,DY,INCY) Adding vectors
        EXTERNAL DCOPY
  ! DCOPY(N,DX,INCX,DY,INCY) copies x to y
        EXTERNAL DDOT
        
  !     x <- (x-mu)      
        call DAXPY(n, -1.0D0, mu, 1, x, 1)
  
negative.code: cf2py integer intent(hide) :: info
  
        DOUBLE PRECISION sig(n,n), x(n), mu(n), like
        INTEGER n, info, i
        DOUBLE PRECISION twopi_N, log_detC, gd
        DOUBLE PRECISION infinity
        PARAMETER (infinity = 1.7976931348623157d308)      
        DOUBLE PRECISION PI
        PARAMETER (PI=3.141592653589793238462643d0) 
  
  
        EXTERNAL DPOTRS
  ! DPOTRS( UPLO, N, NRHS, A, LDA, B, LDB, INFO ) Solves triangular system
  