anchor.diff:           p(i2) = v(i2) + sum
     20  continue  
   c     Solve the triangular system
  -c      call dtrsl(wt,m,col,p(col+1),11,info)
  -      call dpotrs('L',col,1,wt,m,p(col+1),col,info)
  +      call dtrsl(wt,m,col,p(col+1),11,info)
         if (info .ne. 0) return
    
   c     	solve D^(1/2)p1=v1.
  
positive.code: c   Modified for SciPy by removing dependency on linpack
  c================    L-BFGS-B (version 2.1)   ==========================
   
        subroutine setulb(n, m, x, l, u, nbd, f, g, factr, pgtol, wa, iwa,
  
negative.code: c
  c     Subprograms called:
  c
  c       Linpack ... dcopy, dpofa, dtrsl.
  c
  c
  