anchor.diff:  
   c        first Cholesky factor (1,1) block of wn to get LL'
   c                          with L' stored in the upper triangle of wn.
  -      call dpotrf('U',col,wn,m2,info)
  +      call dpofa(wn,m2,col,info)
         if (info .ne. 0) then
   	 info = -1
   	 return
  
positive.code: c
  c     Subprograms called:
  c
  c       CHANGED to use dpotrs from LAPACK 
  c       Linpack ... dtrsl.
  c
  c
  
negative.code: 
  c     Cholesky factorization of (2,2) block of wn.
  
        call dpofa(wn(col+1,col+1),m2,col,info)
        if (info .ne. 0) then
  	 info = -2
  	 return
  