anchor.diff:    subroutine QQSORT(Iw,Ik)
       implicit none
       !------
  -    integer, intent(in):: Ik
  -    integer, intent(inout):: Iw(*)
  -    !------
  -    integer:: itemp
  -    integer:: l
  -    integer:: m
  +    integer(kind=kint), intent(in):: Ik
  +    integer(kind=kint), intent(inout):: Iw(*)
  +    !------
  +    integer(kind=kint):: itemp
  +    integer(kind=kint):: l
  +    integer(kind=kint):: m
   
       if ( Ik<=1 ) return
       do l = 1, Ik - 1
  
positive.code:   subroutine NUFCTX(Xlnzr,Colno,Dsln,Zln,Diag,Indx,Temp,Neqns,Parent,Nch,Nstop,Ndeg,Ndegl,Ir)
      implicit none
      !------
      integer, intent(in):: Ndeg
      integer, intent(in):: Ndegl
      integer, intent(in):: Neqns
      integer, intent(in):: Nstop
      integer, intent(in):: Xlnzr(*)
      integer, intent(in):: Colno(*)
      integer, intent(in):: Parent(*)
      integer, intent(out):: Indx(*)
      integer, intent(inout):: Nch(*)
      real(kind=8), intent(inout):: Zln(Ndeg*Ndeg,*)
      real(kind=8), intent(inout):: Diag(Ndegl,*)
      real(kind=8), intent(out):: Temp(Ndeg*Ndeg,*)
      real(kind=8), intent(inout):: Dsln(Ndeg*Ndeg,*)
      !------
      integer:: ic
      integer:: Ir
      integer:: l
      real(kind=8):: t1
      real(kind=8):: t2
      real(kind=8):: t3
      real(kind=8):: t4
      real(kind=8):: t5
      real(kind=8):: tt
      real(kind=8):: zz(100)
      real(kind=8):: t(100)
  
      !
      ! phase I
  
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)
  