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:       PARAMETER (infinity = 1.7976931348623157d308)      
        DOUBLE PRECISION PI
        PARAMETER (PI=3.141592653589793238462643d0) 
        DOUBLE PRECISION twopi_N, log_detC
  
        EXTERNAL DPOTRF
  ! DPOTRF( UPLO, N, A, LDA, INFO ) Cholesky factorization
        EXTERNAL DPOTRS
  ! DPOTRS( UPLO, N, NRHS, A, LDA, B, LDB, INFO ) Solves triangular system
        EXTERNAL DAXPY
  ! DAXPY(N,DA,DX,INCX,DY,INCY) Adding vectors
        EXTERNAL DCOPY
  ! DCOPY(N,DX,INCX,DY,INCY) copies x to y
        EXTERNAL DDOT
        
  !     C <- cholesky(C)      
        call DPOTRF( 'L', n, C, n, info )
  
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
  