anchor.diff:      do i=1,hecMAT%NP*hecMESH%n_dof
         hecMAT%X(i) = hecMAT%B(i)
       end do
  -    call NUSOL0(hecMAT%X,NEQns,NDEg,NSTop,IPErm,XLNzr,COLno,DIAg,ZLN,DSLn,STAge,ir)
  +    call NUSOL0(hecMAT%X,FCT%NEQns,FCT%NDEg,FCT%NSTop,FCT%IPErm,FCT%XLNzr,FCT%COLno,FCT%DIAg,FCT%ZLN,FCT%DSLn,FCT%STAge,ir)
       call PTIME(t5)
       !* Errors 4
       if ( ir/=0 ) then
  
positive.code:   integer(kind=kint), parameter :: IDBg_ini = 0
    integer(kind=kint), parameter :: IDBg_sym = 0
    integer(kind=kint), parameter :: IDBg_num = 0
  
  contains
    !----------------------------------------------------------------------
  
negative.code:     type (HECMWST_MATRIX), intent(inout)::hecMAT
      integer(kind=kint), intent(in):: Ifmsg
      !------
      type (cholesky_factor), save :: FCT
      !------
      integer(kind=kint):: i98
      integer(kind=kint):: i97
  