anchor.diff:    !======================================================================!
     !> @brief GENQMD
     !======================================================================!
  -  subroutine GENQMD(Neqns,Xadj,Adj0,Perm,Invp,Deg,Marker,Rchset,Nbrhd,Qsize,Qlink,Nofsub,Adjncy)
  +  subroutine GENQMD(Neqns,Nttbr,Xadj,Adj0,Perm,Invp)
       implicit none
       !------
       integer(kind=kint), intent(in):: Neqns
  +    integer(kind=kint), intent(in):: Nttbr
       integer(kind=kint), intent(in):: Adj0(:)
       integer(kind=kint), intent(in):: Xadj(:)
  -    integer(kind=kint), intent(out):: Nofsub
  -    integer(kind=kint), intent(out):: Adjncy(:)
       integer(kind=kint), intent(out):: Perm(:)
       integer(kind=kint), intent(out):: Invp(:)
  -    integer(kind=kint), intent(out):: Deg(:)
  -    integer(kind=kint), intent(out):: Marker(:)
  -    integer(kind=kint), intent(out):: Rchset(:)
  -    integer(kind=kint), intent(out):: Nbrhd(:)
  -    integer(kind=kint), intent(out):: Qsize(:)
  -    integer(kind=kint), intent(out):: Qlink(:)
  -    !------
  +    !------
  +    integer(kind=kint), allocatable:: Deg(:)
  +    integer(kind=kint), allocatable:: Marker(:)
  +    integer(kind=kint), allocatable:: Rchset(:)
  +    integer(kind=kint), allocatable:: Nbrhd(:)
  +    integer(kind=kint), allocatable:: Qsize(:)
  +    integer(kind=kint), allocatable:: Qlink(:)
  +    integer(kind=kint), allocatable:: Adjncy(:)
       integer(kind=kint):: inode
       integer(kind=kint):: ip
       integer(kind=kint):: irch
  
positive.code:     integer(kind=kint), intent(in):: Invp(:)
      integer(kind=kint), intent(in):: Iperm(:)
      integer(kind=kint), intent(out):: Parent(:)
      integer(kind=kint), intent(out):: Ancstr(:)
      !------
      integer(kind=kint):: i
      integer(kind=kint):: ip
      integer(kind=kint):: it
      integer(kind=kint):: k
      integer(kind=kint):: l
  
      do i = 1, Neqns
        Parent(i) = 0
  
negative.code:     integer(kind=kint), intent(in):: Nch(:)
      integer(kind=kint), intent(in):: Invp(:)
      integer(kind=kint), intent(in):: Iperm(:)
      integer(kind=kint), intent(out):: Xleaf(:)
      integer(kind=kint), intent(out):: Leaf(:)
      !------
      integer(kind=kint), allocatable:: Adjncp(:)
      integer(kind=kint):: Lnleaf
      integer(kind=kint):: i
      integer(kind=kint):: ik
      integer(kind=kint):: ip
  