  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.
