A new classical tearing mode (TM) stability simulation workflow has been developed that solves the resistive inner-layer equations in a plasma slab to yield a linear, quasi-toroidal TM growth rate and mode rotation frequency . This workflow combines two-fluid and drift MHD effects in a slab approximation of the resistive inner layer (SLAYER) with an effective tearing stability index as . SLAYER is used to calculate the inner-layer , the STRIDE code is used to calculate a toroidal that includes shaping effects, and the toroidal incorporates effects of thermal conduction on Glasser stabilization. This workflow is rapid and numerically robust across reactor-relevant plasma conditions, and yields growth rates that closely align with analytic predictions in well-documented linear growth rate regimes. Using synthetic equilibria, TM stability was calculated across scans of plasma , inverse aspect-ratio, and toroidal current profile gradient. These scans effectively benchmarked this STRIDE+SLAYER workflow against existing models and showed reliable stability predictions in shaped H-mode-like plasmas. This capability to quickly and robustly predict classical tearing stability in tokamaks will facilitate the mapping of TM stable operational regimes and design of safe discharge trajectories in future devices.