Implementation of bond weighting method for the Grassmann tensor renormalization group

We demonstrate the efficiency of the bond weighting method for the Grassmann tensor renormalization group (TRG). Benchmarking with the two-dimensional Gross-Neveu model with the Wilson fermion at finite density, we show that the bond weighting method improves the accuracy of the original Grassmann TRG. We also provide a sample code of the bond-weighted TRG that can be applied to the two-dimensional models including fermions on a square lattice.