anchor.diff:  
         call make_new_b(hecMESH, hecMAT, BTtmat, iwS, wSL, &
              fstrMAT%num_lagrange, hecTKT%B)
  +
  +      ! use CG when the matrix is symmetric
  +      if (SymType == 1) call hecmw_mat_set_method(hecTKT, 1)
   
         ! solve
         call hecmw_solve_33(hecMESH, hecTKT)
  
positive.code: 
      if (fstrMAT%num_lagrange == 0) then
  
        call hecmw_solve_33(hecMESH, hecMAT)
  
      else
  
  
negative.code:     type (hecmwST_local_mesh), intent(in) :: hecMESH
      type (hecmwST_matrix    ), intent(inout) :: hecMAT
      type (fstrST_matrix_contact_lagrange), intent(inout) :: fstrMAT !< type fstrST_matrix_contact_lagrange
      integer :: ndof, method_org
      integer, allocatable :: iw2(:), iwS(:)
      real(kind=kreal), allocatable :: wSL(:), wSU(:)
      type(T_matrix) :: Tmat, Ttmat
  