anchor.diff:    subroutine GNCLNO(Parent,Xleaf,Leaf,Xlnzr,Colno,Neqns,Nstop,Lncol,Ir)
       implicit none
       !------
  -    integer, intent(in):: Neqns
  -    integer, intent(in):: Nstop
  -    integer, intent(in):: Parent(*)
  -    integer, intent(in):: Xleaf(*)
  -    integer, intent(in):: Leaf(*)
  -    integer, intent(out):: Ir
  -    integer, intent(out):: Lncol
  -    integer, intent(out):: Xlnzr(*)
  -    integer, intent(out):: Colno(*)
  -    !------
  -    integer:: i
  -    integer:: j
  -    integer:: k
  -    integer:: ke
  -    integer:: ks
  -    integer:: l
  -    integer:: nc
  -    integer:: nxleaf
  +    integer(kind=kint), intent(in):: Neqns
  +    integer(kind=kint), intent(in):: Nstop
  +    integer(kind=kint), intent(in):: Parent(*)
  +    integer(kind=kint), intent(in):: Xleaf(*)
  +    integer(kind=kint), intent(in):: Leaf(*)
  +    integer(kind=kint), intent(out):: Ir
  +    integer(kind=kint), intent(out):: Lncol
  +    integer(kind=kint), intent(out):: Xlnzr(*)
  +    integer(kind=kint), intent(out):: Colno(*)
  +    !------
  +    integer(kind=kint):: i
  +    integer(kind=kint):: j
  +    integer(kind=kint):: k
  +    integer(kind=kint):: ke
  +    integer(kind=kint):: ks
  +    integer(kind=kint):: l
  +    integer(kind=kint):: nc
  +    integer(kind=kint):: nxleaf
   
       nc = 0
       Ir = 0
  
positive.code:   subroutine INV33(Zln,Zz,Diag)
      implicit none
      !------
      real(kind=8), intent(in):: Diag(6)
      real(kind=8), intent(in):: Zz(9)
      real(kind=8), intent(out):: Zln(9)
      !------
      Zln(4) = Zz(4) - Zz(1)*Diag(2)
      Zln(7) = Zz(7) - Zz(1)*Diag(4) - Zln(4)*Diag(5)
  
negative.code:   subroutine D3SDOT(Wi,A,B,N)
      implicit none
      !------
      integer(kind=kint), intent(in):: N
      real(kind=kreal), intent(in):: A(3,*)
      real(kind=kreal), intent(in):: B(9,*)
      real(kind=kreal), intent(inout):: Wi(3)
      !------
      integer(kind=kint):: jj
  
      do jj = 1, N
        Wi(1) = Wi(1) - A(1,jj)*B(1,jj) - A(2,jj)*B(4,jj) - A(3,jj)*B(7,jj)
  