Terrestrial laser scanning (TLS) can be used as a millimeter-level measurement tool for forest inventories. However, the stem diameter retrieval accuracy in sample plot scanning is not yet convincing. The errors in each step of stem diameter retrieval algorithms must be evaluated. In this study, six numerical calculation methods for the numerical calculation step, i.e., cylinder fitting (CYF), circle fitting (CF), convex hull line fitting (CLF), the proposed caliper simulation method (CSM), closure B-spline curve fitting (SP) and closure Bézier curve fitting with global convexity (SPC), were applied to stem diameter retrieval, and the similarities and differences were evaluated. The ovality, completeness and roughness were used to evaluate the stem slice point cloud quality. A total of 165 stem slice point clouds at breast height collected from three Larix kaempferi plots were used. Compared with the field-measured stem diameters at breast height (DBHs), the root mean square errors (RMSEs) of the CYF, CF, CLF, CSM, SP and SPC methods were 0.30 cm, 0.30 cm, 0.51 cm, 0.51 cm, 0.56 cm and 0.54 cm, respectively. Compared with the SPC method results, the RMSE of the CSM results was 0.05 cm. The results illustrated that the CYF and CF methods performed the same, as did the CLF and CSM methods. Most DBHs retrieved by the CYF and CF methods were smaller than the field-measured DBHs, and most DBHs retrieved by the CLF, CSM, SP and SPC methods were larger than the field-measured DBHs. This study demonstrated that the CYF and CF methods perform the best and are the most robust, and the measurements by a diameter tape and a caliper are similar enough for forestry inventories. Evaluating and preprocessing stem slice point clouds is a potential way to improve stem diameter retrieval accuracy.