Linear Scaling DFT Calculations with Numerical Atomic Orbitals