anchor.diff:      DOUBLE PRECISION, INTENT(OUT):: Dsln(Ndeg,Ndeg,*)
       DOUBLE PRECISION, INTENT(OUT):: Temp(Ndeg,Ndeg,*)
       DOUBLE PRECISION, INTENT(OUT):: Zln(Ndeg,Ndeg,*)
  -
  +    !------
       INTEGER:: ic
       INTEGER:: j
       INTEGER:: j1
  
positive.code: 
      joc = 0
      DO ic = Nstop, Neqns
         DO i = 1, Nstop
            Temp(i) = 0.0D0
         ENDDO
         ks = Xlnzr(ic)
         ke = Xlnzr(ic+1) - 1
         DO k = ks, ke
  
negative.code:        ENDDO
      ENDDO
    END SUBROUTINE SETIJ
  
    !======================================================================!
    !> @brief MATINI initializes storage for sparse matrix solver.
  