Implementation of the meta-GGA exchange-correlation functional in numerical atomic orbital basis: With systematic testing on SCAN, rSCAN, and r2SCAN functionals.