anchor.diff:  !C | calc. QSI and ETA |
   !C +-------------------+
   !C===
  -      call hecmw_InnerProduct_R_nocomm(hecMESH, NDOF, WW(:,Y ), WW(:,Y ), CG(1))
  -      call hecmw_InnerProduct_R_nocomm(hecMESH, NDOF, WW(:,TT), WW(:,T ), CG(2))
  -      call hecmw_InnerProduct_R_nocomm(hecMESH, NDOF, WW(:,Y ), WW(:,T ), CG(3))
  -      call hecmw_InnerProduct_R_nocomm(hecMESH, NDOF, WW(:,TT), WW(:,Y ), CG(4))
  -      call hecmw_InnerProduct_R_nocomm(hecMESH, NDOF, WW(:,TT), WW(:,TT), CG(5))
  -      S_TIME= HECMW_WTIME()
  -      call hecmw_allreduce_R(hecMESH, CG, 5, HECMW_SUM)
  -      E_TIME= HECMW_WTIME()
  -      Tcomm =  Tcomm + E_TIME - S_TIME
  -
  -      if (iter.eq.1) then
  -        EQ(1)= CG(2)/CG(5)
  -        EQ(2)= 0.d0
  -       else
  -        EQ(1)= (CG(1)*CG(2)-CG(3)*CG(4)) / (CG(5)*CG(1)-CG(4)*CG(4))
  -        EQ(2)= (CG(5)*CG(3)-CG(4)*CG(2)) / (CG(5)*CG(1)-CG(4)*CG(4))
  -      endif
  -
  -      QSI= EQ(1)
  -      ETA= EQ(2)
  +      !call pol_coef(iter, WW, T, TT, Y, QSI, ETA)
  +      !call pol_coef_vanilla(iter, WW, T, TT, Y, QSI, ETA)
  +      call pol_coef_vanilla2(iter, WW, T, TT, Y, QSI, ETA)
   !C===
   
   !C
  
positive.code: !C***
  !
        module hecmw_solver_GPBiCG_33
        contains
  !C
  !C*** hecmw_solve_GPBiCG_33
  
negative.code: 
        real(kind=kreal), dimension(:,:),  allocatable       ::  WW
  
        real(kind=kreal), dimension(2) :: RR
  
        integer(kind=kint ) :: MAXIT
  