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