anchor.diff:          if ( izz0/=izz ) then
             call BRINGU(ZPIv,IPErm,INVp,MARker,izz,NEQns,IRR)
           else
  -          lwk4 = last - neqns1
  -          lwk5 = lwk4 - neqns1
  -          last = lwk5
             call ROTATE(IA,JA,INVp,IPErm,MARker,BTRee,izz,NEQns,NBRhd,QSIze,IRR)
           endif
         endif
  
positive.code:         !
          !   generate skelton graph
          !
          lleaf = last - NTTbr
          lxleaf = lleaf - neqns1
          ladp = lxleaf - neqns1
          last = ladp
          !rmem
          allocate (ADJncp(NEQns+1),stat=IERror)
          if ( IERror/=0 ) stop "ALLOCATION ERROR, adjncp: SUB. matini.f"
          allocate (XLEaf(NEQns+1),stat=IERror)
  
negative.code:       !
        !   build up the binary tree
        !
        allocate (BTRee(2*(NEQns+1)),stat=IERror)
        if ( IERror/=0 ) stop "ALLOCATION ERROR, btree: SUB. matini"
        call GENBTQ(INVp,MARker,BTRee,ZPIv,izz,NEQns)
  