We study the computational aspects of two-temperature heat conduction in case of linear transport coefficients related to electron and lattice subsystems. Several popular finite-difference schemes were compared with commercial finite-element solver revealing the significant dependence of maximum temperature on the spatial-step error.