anchor.diff:      logical, intent(in) :: is_sym
   
       if (INITIALIZED) then
  -       hecMAT%Iarray(98) = 1
  -       hecMAT%Iarray(97) = 1
          INITIALIZED = .false.
       endif
  +
  +    hecMAT%Iarray(98) = 1
  +    hecMAT%Iarray(97) = 1
   
       if (is_sym) then
         SymType = 1
  
positive.code:     ! solve
      call hecmw_solve_33(hecMESHtmp, hecTKT)
      if (DEBUG > 0) write(0,*) myrank, 'DEBUG: solver finished', hecmw_wtime()-t1
  
      ! calc u_s
      call hecmw_localmat_mulvec(BT_all, hecTKT%X, hecMAT%X) !!!<== maybe, BT_all should be BTmat ???
  
negative.code: 
      call hecmw_solve_33(hecMESH, hecMATLag)
  
      hecMAT%Iarray=hecMATLag%Iarray
      hecMAT%Rarray=hecMATLag%Rarray
  
      if (hecMESH%n_neighbor_pe > 0) then
        do i = 1, hecMESH%import_index(hecMESH%n_neighbor_pe)
          hecMESH%import_item(i) = hecMESH%import_item(i) - nb_lag
  