anchor.diff:        write (Ifmsg,*) 'error in nusol0. irr = ', ir
         stop
       endif
  +
  +    if ( timelog > 1 ) then
  +      write(*,*) 'sym fct time : ',t2-t1
  +      write(*,*) 'nuform time  : ',t3-t2
  +      write(*,*) 'num fct time : ',t4-t3
  +      write(*,*) 'solve time   : ',t5-t4
  +    elseif ( timelog > 0 ) then
  +      write(*,*) 'setup time : ',t4-t1
  +      write(*,*) 'solve time : ',t5-t4
  +    endif
  +
       call HECMW_MAT_DUMP_SOLUTION(hecMAT)
     end subroutine HECMW_SOLVE_DIRECT
   
  
positive.code:       call NUFCT0(NEQns,NDEg,NSTop,PARent,NCH,XLNzr,COLno,DIAg,ZLN,DSLn,STAge,ir)
        hecMAT%IARRAY(97) = 0
  
        !*Memory Details
        write (*,*) '*-----------------------------------*'
        write (*,*) '|   Direct  Solver  Memory  Usage   |'
        write (*,*) '*-----------------------------------*'
        write (*,*) 'INTEGER memory: ', real(IALoc*4)/real(1048576), 'MB'
        write (*,*) 'REAL*8  memory: ', real(RALoc*8)/real(1048576), 'MB'
        write (*,*) 'TOTAL   memory: ', real((RALoc*2+IALoc)*4)/real(1048576), 'MB'
        write (*,*) '*-----------------------------------*'
      endif
      call PTIME(t4)
  
  
negative.code:       stop
      endif
  
      !* Solve
      !* Backsubstitute
      do i=1,hecMAT%NP*hecMESH%n_dof
  