anchor.diff:    !                30  after LU decomposition
     !                40  after solving
     !======================================================================!
  -  subroutine MATINI(Ir)
  -    implicit none
  -    !------
  +  subroutine MATINI(NEQns,NTTbr,JCOl,IROw,LEN_colno,NSTop,LEN_dsln,IPErm,INVp,PARent,NCH,XLNzr,COLno,IALoc,RALoc,STAge,Ir)
  +    implicit none
  +    !------
  +    integer, intent(in):: NEQns
  +    integer, intent(in):: NTTbr
  +    integer(kind=4), intent(in) :: jcol(:)
  +    integer(kind=4), intent(in) :: irow(:)
  +    integer, intent(out):: LEN_colno
  +    integer, intent(out):: NSTop
  +    integer, intent(out):: LEN_dsln
  +    integer, pointer :: IPErm(:)
  +    integer, pointer :: INVp(:)
  +    integer, pointer :: PARent(:)
  +    integer, pointer :: NCH(:)
  +    integer, pointer :: XLNzr(:)
  +    integer, pointer :: COLno(:)
  +    integer, intent(out):: IALoc
  +    integer, intent(out):: RALoc
  +    integer, intent(out):: STAge
       integer, intent(out):: Ir
       !------
  +    !*Work arrays
  +    integer(kind=4), allocatable :: zpiv(:)
  +    integer(kind=4), allocatable :: jcpt(:)
  +    integer(kind=4), allocatable :: jcolno(:)
  +    integer(kind=4), allocatable :: ia(:)
  +    integer(kind=4), allocatable :: ja(:)
  +    integer(kind=4), allocatable :: deg(:)
  +    integer(kind=4), allocatable :: marker(:)
  +    integer(kind=4), allocatable :: rchset(:)
  +    integer(kind=4), allocatable :: nbrhd(:)
  +    integer(kind=4), allocatable :: qsize(:)
  +    integer(kind=4), allocatable :: qlink(:)
  +    integer(kind=4) :: nofsub
  +    integer(kind=4), allocatable :: adjncy(:)
  +    integer(kind=4), allocatable :: btree(:)
  +    integer(kind=4), allocatable :: pordr(:)
  +    integer(kind=4), allocatable :: adjncp(:)
  +    integer(kind=4), allocatable :: xleaf(:)
  +    integer(kind=4), allocatable :: leaf(:)
       integer:: ir1
  +    integer:: irr
       integer:: iv1
       integer:: izz
       integer:: izz0
  
positive.code:   !======================================================================!
    !> @brief SUM
    !======================================================================!
    subroutine sum(Ic,Xlnzr,Colno,Zln,Diag,Nch,Par,Temp,Indx)
      implicit none
      !------
      integer, intent(in):: Ic
  
negative.code:         deallocate (QLInk)
          deallocate (ADJncy)
          deallocate (ZPIv)
          !
          !   build up xlnzr,colno  (this is the symbolic fct.)
          !
  