anchor.diff:    !C*** hecmw_mpc_mat_finalize
     !C***
     !C
  -  subroutine hecmw_mpc_mat_finalize(hecMESH, hecMAT, hecMATmpc)
  +  subroutine hecmw_mpc_mat_finalize(hecMESH, hecMAT, hecMESHmpc, hecMATmpc)
       implicit none
       type (hecmwST_local_mesh), intent(in) :: hecMESH
       type (hecmwST_matrix), intent(in) :: hecMAT
  +    type (hecmwST_local_mesh), pointer :: hecMESHmpc
       type (hecmwST_matrix), pointer :: hecMATmpc
       integer(kind=kint) :: totalmpc, MPC_METHOD
   
  
positive.code:     call hecmw_allreduce_I1 (hecMESH, totalmpc, hecmw_sum)
  
      if (totalmpc == 0) then
        hecMATmpc => hecMAT
        return
      endif
  
negative.code:       !if (hecMESH%my_rank.eq.0) write(0,*) "MPC Method: MPC-CG"
      case (3)  ! elimination
        !if (hecMESH%my_rank.eq.0) write(0,*) "MPC Method: Elimination"
        call hecmw_trimatmul_TtKT_mpc(hecMESHmpc, hecMAT, hecMATmpc)
      end select
  
    end subroutine hecmw_mpc_mat_ass
  