anchor.diff:        stop
       endif
   
  -    TOM(1) = t2-t1
  -    TOM(2) = t3-t2
  -    TOM(3) = t4-t3
  -
       !* Solve
       !* Backsubstitute
       do i=1,hecMAT%NP*hecMESH%n_dof
  
positive.code:       !* Symbolic factorization
        call MATINI(NEQns,NTTbr,JCOl,IROw,LEN_colno,NSTop,LEN_dsln,IPErm,INVp,PARent,NCH,XLNzr,COLno,IALoc,RALoc,STAge,ir)
        hecMAT%IARRAY(98) = 0
        write (6,*) "symbolic fct done"
      endif
      call PTIME(t2)
      t3 = t2
  
negative.code:     !------
      integer(kind=kint):: i98
      integer(kind=kint):: i97
      integer(kind=kint):: timelog
      integer(kind=kint):: iterlog
      integer(kind=kint):: ir
      integer(kind=kint):: i
      real(kind=kreal):: t1
  