Calculation of fermionic Green functions with Grassmann higher-order tensor renormalization group

We develop a calculational method for fermionic Green functions in the framework of Grassmann higher-order tensor renormalization group. The validity of the method is tested by applying it to the three-dimensional free Wilson fermion system. We compare the numerical results for the chiral condensate and two-point correlation functions with the exact ones obtained by analytical methods.