anchor.diff:          i = hecMAT%ITEML(k)
           kk = kk + 1
           call VLCPY(val,hecMAT%AL(ndof2*(k-1)+1:ndof2*k),ndof)
  -        call STAIJ1(0,j,i,val,FCT%NEQns,FCT%NDEg,FCT%LEN_colno,FCT%NSTop,FCT%LEN_dsln,FCT%INVp,FCT%XLNzr,FCT%COLno,&
  -             FCT%DIAg,FCT%ZLN,FCT%DSLn,FCT%RALoc,FCT%STAge,Ir)
  +        call STAIJ1(0,j,i,val,FCT,Ir)
         enddo
       enddo
   
  
positive.code:       !*Diagonal
        kk = kk + 1
        call VLCPY(val,hecMAT%D(ndof2*(j-1)+1:ndof2*j),ndof)
        call STAIJ1(0,j,j,val,FCT%NEQns,FCT%NDEg,FCT%LEN_colno,FCT%NSTop,FCT%LEN_dsln,FCT%INVp,FCT%XLNzr,FCT%COLno,&
             FCT%DIAg,FCT%ZLN,FCT%DSLn,FCT%RALoc,FCT%STAge,Ir)
  
        do i = 1, ndof
          if ( val((i-1)*ndof+i)<=0 ) write (idbg,*) 'j,j,val:', j, i, val((i-1)*ndof+i)
  
negative.code:   !          aij      value
    !      (o)
    !          iv       communication array
    !  sets DIAg, ZLN, DSLn
    !  updates RALoc, STAge
    !======================================================================!
    subroutine STAIJ1(Isw,I,J,Aij,FCT,Ir)
      implicit none
      !------
      integer(kind=kint), intent(in):: I
      integer(kind=kint), intent(in):: Isw
      integer(kind=kint), intent(in):: J
      real(kind=kreal), intent(in):: Aij(:)
      type (cholesky_factor), intent(inout):: FCT
      integer(kind=kint), intent(out):: Ir
      !------
      integer(kind=kint):: ndeg2
  