“…The recent literature relevant to the construction and analysis of discretizations for (1.1) and closely related problems is very diverse. For instance, numerical methods specifically designed to capture boundary layers in singularly perturbed generalized Burgers-Huxley equations have been studied in [18], different types of finite differences have been used in [23,19,25,28], spectral, B-spline and Chebyshev wavelet collocation methods have been advanced in [1,15,31,7], numerical solutions obtained with the so-called adomain decomposition were analyzed in [14], homotopy perturbation techniques were used in [20], Strang splittings were proposed in [8], meshless radial basis functions were studied in [17], generalized finite differences and finite volume schemes have been analyzed in [9,33] for the restriction of (1.1) to the diffusive Nagumo (or bistable) model, and a finite element method satisfying a discrete maximum principle was introduced in [12] (the latter reference is closer to the present study). Although there is a growing interest in developing numerical techniques for the generalized Burgers-Huxley equation, it appears that the aspects of error analysis for finite element discretizations have not been yet thoroughly addressed.…”